Dry molten globule intermediates and the mechanism of protein unfolding

Determinants of Unfolding Cooperativity and Binding Are Decoupled in a DNA Binding Domain

The relative magnitudes of noncovalent stabilization energies or the coupling free energies in folded proteins are anisotropically distributed, uniquely influencing folding and functional behaviors. In this regard, the fructose repressor (FruR) DBD belonging to the LacR repressor family harbors a three-residue insertion─KQY─between the canonical second and third helices. This sequence insertion promotes a strong Tyr-Tyr stacking interaction that is not observed in related homologues. Combining experiments with simulations, we show that the Tyr-Tyr stacking contributes to a decoupled unfolding due to the localization of a large part of the stabilization energy in this specific structural region. This leads to melting temperatures from different probes spanning nearly 10 K, while concomitantly stabilizing a partially structured intermediate state. Disruption of the aromatic stacking interaction via an alanine mutation promotes a molten-globular state whose native ensemble is replete with non-native interactions while displaying enhanced thermodynamic fluctuations and minimal calorimetric cooperativity. Surprisingly, the molten-globular variant of FruR DBD binds to the operator site on DNA with an affinity similar to that of the wild-type but with altered secondary-structure characteristics in the bound state, underscoring the chaperone-like role of DNA through its large negative electrostatic potential. FruR DBD thus appears to be at the verge of disorder as expected of an entropically destabilizing three-residue insertion but is rescued by the aromatic stacking interaction that distinctly dictates the finer details of stability, cooperativity, and binding.

Dry molten globule conformational state of CYP11A1 (SCC) regulates the first step of steroidogenesis in the mitochondrial matrix

Multiple metabolic events occur in mitochondria. Mitochondrial protein translocation from the cytoplasm across compartments depends on the amino acid sequence within the precursor. At the mitochondria associated-ER membrane, misfolding of a mitochondrial targeted protein prior to import ablates metabolism. CYP11A1, cytochrome P450 cholesterol side chain cleavage enzyme (SCC), is imported from the cytoplasm to mitochondrial matrix catalyzing cholesterol to pregnenolone, an essential step for metabolic processes and mammalian survival. Multiple steps regulate the availability of an actively folded SCC; however, the mechanism is unknown. We identified that a dry molten globule state of SCC exists in the matrix by capturing intermediate protein folding steps dictated by its C-terminus. The intermediate dry molten globule state in the mitochondrial matrix of living cells is stable with a limited network of interaction and is inactive. The dry molten globule is activated with hydrogen ions availability, triggering cleavage of cholesterol sidechain, and initiating steroidogenesis.

Protein folding as a jamming transition

Proteins fold to a specific functional conformation with a densely packed hydrophobic core that controls their stability. We develop a geometric, yet all-atom model for proteins that explains the universal core packing fraction of ϕ c = 0.55 found in experimental measurements. We show that as the hydrophobic interactions increase relative to the temperature, a novel jamming transition occurs when the core packing fraction exceeds ϕ c . The model also recapitulates the global structure of proteins since it can accurately refold to native-like structures from partially unfolded states.

A pressure-jump EPR system to monitor millisecond conformational exchange rates of spin-labeled proteins

Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) using nitroxide spin labels is a well-established technology for mapping site-specific secondary and tertiary structure and for monitoring conformational changes in proteins of any degree of complexity, including membrane proteins, with high sensitivity. SDSL-EPR also provides information on protein dynamics in the time scale of ps-µs using continuous wave lineshape analysis and spin lattice relaxation time methods. However, the functionally important time domain of µs-ms, corresponding to large-scale protein motions, is inaccessible to those methods. To extend SDSL-EPR to the longer time domain, the perturbation method of pressure-jump relaxation is implemented. Here, we describe a complete high-pressure EPR system at Q-band for both static pressure and millisecond-timescale pressure-jump measurements on spin-labeled proteins. The instrument enables pressure jumps both up and down from any holding pressure, ranging from atmospheric pressure to the maximum pressure capacity of the system components (~3500 bar). To demonstrate the utility of the system, we characterize a local folding-unfolding equilibrium of T4 lysozyme. The results illustrate the ability of the system to measure thermodynamic and kinetic parameters of protein conformational exchange on the millisecond timescale.

The attenuating ability of deep eutectic solvents towards the carboxylated multiwalled carbon nanotubes induced denatured β-lactoglobulin structure

The stabilization of proteins has been a major challenge for their practical utilization in industrial applications. Proteins can easily lose their native conformation in the presence of denaturants, which unfolds the protein structure. Since the introduction of deep eutectic solvents (DESs), there are numerous studies in which DESs act as promising co-solvents that are biocompatible with biomolecules. DESs have emerged as sustainable biocatalytic media and an alternative to conventional organic solvents and ionic liquids (ILs). However, the superiority of DESs over the deleterious influence of denaturants on proteins is often neglected. To address this, we present the counteracting ability of biocompatible DESs, namely, choline chloride-glycerol (DES-1) and choline chloride-urea (DES-2), against the structural changes induced in β-lactoglobulin (Blg) by carboxylated multiwalled carbon nanotubes (CA-MWCNTs). The work is substantiated with various spectroscopic and thermal studies. The spectroscopic results revealed that the fluorescence emission intensity enhances for the protein in DESs. Contrary to this, the emission intensity extremely quenches in the presence of CA-MWCNTs. However, in the mixture of DESs and CA-MWCNTs, there was a slight increase in the fluorescence intensity. Circular dichroism spectral studies reflect the reappearance of the native band that was lost in the presence of CA-MWCNTs, which is a good indicator of the counteraction ability of DESs. Further, thermal fluorescence studies showed that the protein exhibited extremely great thermal stability in both DESs as well as in the mixture of DES-CA-MWCNTs compared to the protein in buffer. This study is also supported by dynamic light scattering and zeta potential measurements; the results reveal that DESs were successfully able to maintain the protein structure. The addition of CA-MWCNTs results in complex formation with the protein, which is indicated by the increased hydrodynamic size of the protein. The presence of DESs in the mixture of CA-MWCNTs and DESs was quite successful in eliminating the negative impact of CA-MWCNTs on protein structural alteration. DES-1 proved to be superior to DES-2 over counteraction against CA-MWCNTs and maintained the native conformation of the protein. Overall, both DESs act as recoiling media for both native and unfolded (denatured by CA-MWCNTs) Blg structures. Both the DESs can be described as potential co-solvents for Blg with increased structural and thermal stability of the protein. To the best of our knowledge, this study for the first time has demonstrated the role of choline-based DESs in the mixture with CA-MWCNTs in the structural transition of Blg. The DESs in the mixture successfully enhance the stability of the protein by reducing the perturbation caused by CA-MWCNTs and then amplifying the advantages of the DESs present in the mixture. Overall, these results might find implications for understanding the role of DES-CA-MWCNT mixtures in protein folding/unfolding and pave a new direction for the development of eco-friendly protein-protective solvents.

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.

Pre-Molten, Wet, and Dry Molten Globules en Route to the Functional State of Proteins

Transitions between the unfolded and native states of the ordered globular proteins are accompanied by the accumulation of several intermediates, such as pre-molten globules, wet molten globules, and dry molten globules. Structurally equivalent conformations can serve as native functional states of intrinsically disordered proteins. This overview captures the characteristics and importance of these molten globules in both structured and intrinsically disordered proteins. It also discusses examples of engineered molten globules. The formation of these intermediates under conditions of macromolecular crowding and their interactions with nanomaterials are also reviewed.

Dry Molten Globule-Like Intermediates in Protein Folding, Function, and Disease

The performance of a protein depends on its correct folding to the final functional native form. Hence, understanding the process of protein folding has remained an important field of research for the scientific community for the past five decades. Two important intermediate states, namely, wet molten globule (WMG) and dry molten globule (DMG), have emerged as critical milestones during protein folding-unfolding reactions. While much has been discussed about WMGs as a common unfolding intermediate, the evidence for DMGs has remained elusive owing to their near-native features, which makes them difficult to probe using global structural probes. This Review puts together the available literature and new evidence on DMGs to give a broader perspective on the universality of DMGs and discuss their significance in protein folding, function, and disease.

The Molten Globule State of a Globular Protein in a Cell Is More or Less Frequent Case Rather than an Exception

Quite a long time ago, Oleg B. Ptitsyn put forward a hypothesis about the possible functional significance of the molten globule (MG) state for the functioning of proteins. MG is an intermediate between the unfolded and the native state of a protein. Its experimental detection and investigation in a cell are extremely difficult. In the last decades, intensive studies have demonstrated that the MG-like state of some globular proteins arises from either their modifications or interactions with protein partners or other cell components. This review summarizes such reports. In many cases, MG was evidenced to be functionally important. Thus, the MG state is quite common for functional cellular proteins. This supports Ptitsyn’s hypothesis that some globular proteins may switch between two active states, rigid (N) and soft (MG), to work in solution or interact with partners.

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

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

Reframing the Protein Folding Problem: Entropy as Organizer

It has been a long-standing conviction that a protein's native fold is selected from a vast number of conformers by the optimal constellation of enthalpically favorable interactions. In marked contrast, this Perspective introduces a different mechanism, one that emphasizes conformational entropy as the principal organizer in protein folding while proposing that the conventional view is incomplete. This mechanism stems from the realization that hydrogen bond satisfaction is a thermodynamic necessity. In particular, a backbone hydrogen bond may add little to the stability of the native state, but a completely unsatisfied backbone hydrogen bond would be dramatically destabilizing, shifting the U(nfolded) ⇌ N(ative) equilibrium far to the left. If even a single backbone polar group is satisfied by solvent when unfolded but buried and unsatisfied when folded, that energy penalty alone, approximately +5 kcal/mol, would rival almost the entire free energy of protein stabilization, typically between -5 and -15 kcal/mol under physiological conditions. Consequently, upon folding, buried backbone polar groups must form hydrogen bonds, and they do so by assembling scaffolds of α-helices and/or strands of β-sheet, the only conformers in which, with rare exception, hydrogen bond donors and acceptors are exactly balanced. In addition, only a few thousand viable scaffold topologies are possible for a typical protein domain. This thermodynamic imperative winnows the folding population by culling conformers with unsatisfied hydrogen bonds, thereby reducing the entropy cost of folding. Importantly, conformational restrictions imposed by backbone···backbone hydrogen bonding in the scaffold are sequence-independent, enabling mutation─and thus evolution─without sacrificing the structure.

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.

Protein folding - seeing is deceiving

This Perspective is intended to raise questions about the conventional interpretation of protein folding. According to the conventional interpretation, developed over many decades, a protein population can visit a vast number of conformations under unfolding conditions, but a single dominant native population emerges under folding conditions. Accordingly, folding comes with a substantial loss of conformational entropy. How is this price paid? The conventional answer is that favorable interactions between and among the side chains can compensate for entropy loss, and moreover, these interactions are responsible for the structural particulars of the native conformation. Challenging this interpretation, the Perspective introduces a proposal that high energy (i.e., unfavorable) excluding interactions winnow the accessible population substantially under physical-chemical conditions that favor folding. Both steric clash and unsatisfied hydrogen bond donors and acceptors are classified as excluding interactions, so called because conformers with such disfavored interactions will be largely excluded from the thermodynamic population. Both excluding interactions and solvent factors that induce compactness are somewhat nonspecific, yet together they promote substantial chain organization. Moreover, proteins are built on a backbone scaffold consisting of α-helices and strands of β-sheet, where the number of hydrogen bond donors and acceptors is exactly balanced. These repetitive secondary structural elements are the only two conformers that can be both completely hydrogen-bond satisfied and extended indefinitely without encountering a steric clash. Consequently, the number of fundamental folds is limited to no more than ~10,000 for a protein domain. Once excluding interactions are taken into account, the issue of "frustration" is largely eliminated and the Levinthal paradox is resolved. Putting the "bottom line" at the top: it is likely that hydrogen-bond satisfaction represents a largely under-appreciated parameter in protein folding models.

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.

Impact of denaturing agents on surface properties of myoglobin solutions

The addition of denaturants strongly influences the surface properties of aqueous myoglobin solutions. The effect differs from the results for mixed solutions of the denaturants and other globular proteins, for example, bovine serum albumin (BSA), lysozyme and β-lactoglobulin (BLG), although the surface properties of the solutions of the pure proteins are similar. The kinetic dependencies of the dynamic surface elasticity of myoglobin solutions with guanidine hydrochloride (GuHCl) reveal at least two adsorption steps at denaturant concentrations higher than 1 M: a very fast increase of the dynamic surface elasticity to approximately 30 mN/m at the beginning of adsorption, and a slower growth to abnormally high values of 250-300 mN/m. At the same time, the surface elasticity of BSA/GuHCl, BLG/GuHCl and lysozyme/GuHCl solutions is a non-monotonic function of the surface age, and does not exceed 50 mN/m close to equilibrium. The high surface elasticity of myoglobin/GuHCl solutions may be associated with protein aggregation in the surface layer. The formation of aggregates is confirmed by ellipsometry and Brewster angle microscopy. The addition of ionic surfactants to protein solutions leads to the formation of myoglobin/surfactant complexes, and the kinetic dependencies of the dynamic surface elasticity display local maxima indicating multistep adsorption kinetics, unlike the corresponding results for solutions of other globular proteins mixed with ionic surfactants. Ellipsometry and infrared reflection-absorption spectroscopy allow tracing the adsorption of the complexes and their displacement from the interface at high surfactant concentrations.

FTIR Spectroscopy Study of the Secondary Structure Changes in Human Serum Albumin and Trypsin under Neutral Salts

The effect of neutral salts on protein conformation was first analyzed by Hofmeister in 1888, however, even today this phenomenon is not completely understood. To clarify this effect, we studied changes in the secondary structure of two proteins: human serum albumin with predominantly α-helical structure and porcine pancreas β-trypsin with the typical β-structural arrangement in aqueous solutions of neutral salts (KSCN, KCl, (NH₄)₂SO₄). The changes in the secondary structure were studied at 23 °C and 80 °C by using the second derivative deconvolution method of the IR spectra. Our results demonstrated that the ability of the salts to stabilize/destabilize these two proteins correlates with the Hofmeister series of ions. At the same time, some exceptions were also observed. The destabilization of the native structures of both α-helical albumin and β-structural trypsin upon interaction with neutral salts leads to the formation of intermolecular β-sheets typical for amyloid fibrils or amorphous aggregates. Thus, our quantitative FTIR-spectroscopy analysis allowed us to further clarify the mechanisms and complexity of the neutral salt actions on protein structures which may lead to strategies preventing unwelcome misfolding of proteins.

Equilibrium partially folded states of B. licheniformis[Formula: see text]-lactamase

[Formula: see text]-Lactamases (penicillinases) facilitate bacterial resistance to antibiotics and are excellent theoretical and experimental models in protein structure, dynamics and evolution. Bacillus licheniformis exo-small penicillinase (ESP) is a Class A [Formula: see text]-lactamase with three tryptophan residues located one in each of its two domains and one in the interface between domains. The conformational landscape of three well-characterized ESP Trp[Formula: see text]Phe mutants was characterized in equilibrium unfolding experiments by measuring tryptophan fluorescence, far-UV CD, activity, hydrodynamic radius, and limited proteolysis. The Trp[Formula: see text]Phe substitutions had little impact on the native conformation, but changed the properties of the partially folded states populated at equilibrium. The results were interpreted in the framework of modern theories of protein folding.

CcdB at pH 4 Forms a Partially Unfolded State with a Dry Core

pH is an important factor that affects the protein structure, stability, and activity. Here, we probe the nature of the low-pH structural form of the homodimeric CcdB (controller of cell death B) protein. Characterization of CcdB protein at pH 4 and 300 K using circular dichroism spectroscopy, 8-anilino-1-naphthalene-sulphonate binding, and Trp solvation studies suggests that it forms a partially unfolded state with a dry core at equilibrium under these conditions. CcdB remains dimeric at pH 4 as shown by multiple techniques, such as size-exclusion chromatography coupled to multiangle light scattering, analytical ultracentrifugation, and electron paramagnetic resonance. Comparative analysis using two-dimensional ¹⁵N-¹H heteronuclear single-quantum coherence NMR spectra of CcdB at pH 4 and 7 suggests that the pH 4 and native state have similar but nonidentical structures. Hydrogen-exchange-mass-spectrometry studies demonstrate that the pH 4 state has substantial but anisotropic changes in local stability with core regions close to the dimer interface showing lower protection but some other regions showing higher protection relative to pH 7.

Slow Motion Protein Dance Visualized Using Red-Edge Excitation Shift of a Buried Fluorophore

It has been extremely challenging to detect protein structures with a dynamic core, such as dry molten globules, that remain in equilibrium with the tightly packed native (N) state and that are important for a myriad of entropy-driven protein functions. Here, we detect the higher entropy conformations of a human serum protein, using red-edge excitation shift experiments. We covalently introduced a fluorophore inside the protein core and observed that in a subset of native population, the side chains of the polar and buried residues have different spatial arrangements than the mean population and that they solvate the fluorophore on a timescale much slower than the nanosecond timescale of fluorescence. Our results provide direct evidence for the dense fluidity of protein core and show that alternate side-chain packing arrangements exist in the core that might be important for multiple binding functions of this protein.

Native mass spectrometry beyond ammonium acetate: effects of nonvolatile salts on protein stability and structure

Structures and stabilities of proteins investigated with native mass spectrometry can be affected by nonvolatile salts, including Tris buffer, in solution.

Mechanical Softening of a Small Ubiquitin-Related Modifier Protein Due to Temperature Induced Flexibility at the Core

Despite the growing interest in the thermal softening of proteins, the mechanistic details of it are far from understood. β-Grasp proteins have globular shape with compact structure and they are mechanically resilient. The β-clamp or mechanical clamp in them formed by the interactions between the terminal β-strands is generally associated with the protein mechanical resistance. Although previous studies showed that temperature can perturb the protein mechanical stability, the structural changes leading to the lowered mechanical resistance are not known. Here, we investigated the temperature dependent mechanical stability of small ubiquitin-related modifier 2 (SUMO2) using single-molecule force spectroscopy (SMFS) and the corresponding conformational changes using ensemble experiments. SMFS studies on the polyprotein of SUMO2 estimate a decrease in the spring constant of the protein from 4.50 to 1.35 N/m upon increasing the temperature from 5 to 45 °C. Interestingly, near-UV circular dichroism spectroscopy reveals a decrease in tertiary structure content while the overall secondary structure of the protein remains unchanged. Steady-state fluorescence and quenching studies on SUMO2 with a tryptophan mutation at the core (F60W) show that the nonpolar environment of the tryptophan is unchanged and the protein core is inaccessible to the bulk solvent, in the same temperature range. We attribute the thermal softening observed in atomic force microscopy (AFM) experiments to the reduction in tertiary structure of SUMO2. Our results provide evidence for the importance of the intramolecular interactions at the protein core along with the β-clamp or mechanical clamp in providing the mechanical resistance as well as in modulating the protein stiffness.

Aggregation-primed molten globule conformers of the p53 core domain provide potential tools for studying p53C aggregation in cancer

The functionality of the tumor suppressor p53 is altered in more than 50% of human cancers, and many individuals with cancer exhibit amyloid-like buildups of aggregated p53. An understanding of what triggers the pathogenic amyloid conversion of p53 is required for the further development of cancer therapies. Here, perturbation of the p53 core domain (p53C) with subdenaturing concentrations of guanidine hydrochloride and high hydrostatic pressure revealed native-like molten globule (MG) states, a subset of which were highly prone to amyloidogenic aggregation. We found that MG conformers of p53C, probably representing population-weighted averages of multiple states, have different volumetric properties, as determined by pressure perturbation and size-exclusion chromatography. We also found that they bind the fluorescent dye 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) and have a native-like tertiary structure that occludes the single Trp residue in p53. Fluorescence experiments revealed conformational changes of the single Trp and Tyr residues before p53 unfolding and the presence of MG conformers, some of which were highly prone to aggregation. p53C exhibited marginal unfolding cooperativity, which could be modulated from unfolding to aggregation pathways with chemical or physical forces. We conclude that trapping amyloid precursor states in solution is a promising approach for understanding p53 aggregation in cancer. Our findings support the use of single-Trp fluorescence as a probe for evaluating p53 stability, effects of mutations, and the efficacy of therapeutics designed to stabilize p53.

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.

pH-Dependent cooperativity and existence of a dry molten globule in the folding of a miniprotein BBL

Solution pH plays an important role in protein dynamics, stability, and folding; however, detailed mechanisms remain poorly understood. Here we use continuous constant pH molecular dynamics in explicit solvent with pH replica exchange to describe the pH profile of the folding cooperativity of a miniprotein BBL, which has drawn intense debate in the past. Our data reconciled the two opposing hypotheses (downhill vs. two-state) and uncovered a sparsely populated unfolding intermediate. As pH is lowered from 7 to 5, the folding barrier vanishes. As pH continues to decrease, the unfolding barrier lowers and denaturation is triggered by the protonation of Asp162, consistent with experimental evidence. Interestingly, unfolding proceeded via an intermediate, with intact secondary structure and a compact, unlocked hydrophobic core shielded from solvent, lending support to the recent hypothesis of a universal dry molten globule in protein folding. Our work demonstrates that constant pH molecular dynamics is a unique tool for testing this and other hypotheses to advance the knowledge in protein dynamics, stability, and folding.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Identification and Characterization of an Inside-Out Folding Intermediate of T4 Phage Sliding Clamp

Protein folding process involves formation of transiently occurring intermediates that are difficult to isolate and characterize. It is both necessary and interesting to characterize the structural conformations adopted by these intermediates, also called molten globules (MG), to understand protein folding. Here, we investigated the equilibrium (un)folding intermediate state of T4 phage gene product 45 (gp45, also known as DNA polymerase processivity factor or sliding clamp) obtained during chemical denaturation. We show that gp45 undergoes substantial conformational rearrangement during unfolding and forms an expanded dry-MG. By monitoring the fluorescence of tryptophans that were strategically introduced at various sites, we demonstrate that the urea-treated molecule has its surface residues flip inside the core, and closely placed residues move farther. We were also able to isolate and purify the MG form of gp45 in native condition (i.e., nondenaturing buffer, at physiological pH and temperature); characteristics of this purified molecule substantially match with urea-treated wild-type gp45. To the best of our knowledge, this is one of the few reports that demonstrate the isolation and purification of a protein folding intermediate in native condition. We believe that our work not only allows us to dissect the process of protein folding, but will also help in the designing of folding inhibitors against sliding clamps to treat a wide variety of diseases from bacterial infection to cancer, due to the vast presence of clamps in all the domains of life.

Macromolecular crowding induces molten globule state in the native myoglobin at physiological pH

Here, we report the formation of molten globule state of the native myoglobin in crowded environment. We have used Soret absorption spectroscopy and far-UV circular dichroism to monitor changes in tertiary and secondary structures of myoglobin, respectively. Our results reveal that in the presence of ficoll 70, the secondary structure of myoglobin remains unchanged while tertiary structure is lost significantly. 1-anilinonaphthalene-8-sulfonate binding experiments showed that myoglobin in the presence of various concentrations of ficoll 70, has newly exposed hydrophobic surfaces. Dynamic light scattering measurements show that there is almost 1.5 times increase in the hydrodynamic volume of myoglobin in the crowded environment. These structural characteristics of myoglobin in the presence of 300mg/ml ficoll 70 resemble those of molten globule state. Isothermal titration calorimetric (ITC) measurements show that ficoll 70 binds to myoglobin, whereas it shows no interaction with apo form of the protein. ITC results indicate that the reason behind this unique behavior of ficoll 70 towards myoglobin may be interaction of ficoll 70 with the heme group of myoglobin, which was further confirmed by the docking studies. We hypothesize that the soft interactions between heme and ficoll 70 leads to the formation of molten globule in myoglobin.

Monitoring Hydrogen Exchange During Protein Folding by Fast Pressure Jump NMR Spectroscopy

A method is introduced that permits direct observation of the rates at which backbone amide hydrogens become protected from solvent exchange after rapidly dropping the hydrostatic pressure inside the NMR sample cell from denaturing (2.5 kbar) to native (1 bar) conditions. The method is demonstrated for a pressure-sensitized ubiquitin variant that contains two Val to Ala mutations. Increased protection against hydrogen exchange with solvent is monitored as a function of time during the folding process. Results for 53 backbone amides show narrow clustering with protection occurring with a time constant of ca. 85 ms, but slower protection is observed around a reverse turn near the C-terminus of the protein. Remarkably, the native NMR spectrum returns with this slower time constant of ca. 150 ms, indicating that the almost fully folded protein retains molten globule characteristics with severe NMR line broadening until the final hydrogen bonds are formed. Prior to crossing the transition state barrier, hydrogen exchange protection factors are close to unity, but with slightly elevated values in the β₁-β₂ hairpin, previously shown to be already lowly populated in the urea-denatured state.

Probing the Acid-Induced Packing Structure Changes of the Molten Globule Domains of a Protein near Equilibrium Unfolding

Using simultaneously scanning small-angle X-ray scattering (SAXS) and UV-vis absorption with integrated online size exclusion chromatography, supplemental with molecular dynamics simulations, we unveil the long-postulated global structure evolution of a model multidomain protein bovine serum albumin (BSA) during acid-induced unfolding. Our results differentiate three global packing structures of the three molten globule domains of BSA, forming three intermediates I₁, I₂, and E along the unfolding pathway. The I₁-I₂ transition, overlooked in all previous studies, involves mainly coordinated reorientations across interconnected molten globule subdomains, and the transition activates a critical pivot domain opening of the protein for entering into the E form, with an unexpectedly large unfolding free energy change of -9.5 kcal mol^-1, extracted based on the observed packing structural changes. The revealed local packing flexibility and rigidity of the molten globule domains in the E form elucidate how collective motions of the molten globule domains profoundly influence the folding-unfolding pathway of a multidomain protein.

A dry molten globule-like intermediate during the base-induced unfolding of a multidomain protein

An early intermediate during the base-induced unfolding of a multidomain protein resembles a dry molten globule state in which the structure is expanded without core hydration.

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.

The Effect of pH on Globular State of Lipase-3646; an Appropriate Model for Molten Globule Investigations

Secondary structure content of proteins in molten globule state is relatively constant while the quantity of tertiary structures clearly declines due to alterations in side-chain packing. In the present study, we analyze the MG state of lipase-3646 for the first time. We introduce lipase-3646 as an appropriate model for investigating the properties and behavior of a protein in MG state as well as folding pathway. Applying fluorescence spectroscopy we measured both intrinsic and extrinsic fluorescence of lipase-3646 in a pH range from 1.0 to 12.0. It was found that at pH 3.0 the protein acquires a MG state. Applying far-UV circular dichroism (CD), our analysis on the secondary structure of lipase-3646 revealed a slight change in the MG state intermediate (pH 3.0) compared to the native state (pH 8.5), which this amount of change is common for MG. Measurements in near-UV CD also showed a significant change in the enzyme conformation at pH 3.0 in comparison with the pH 8.5 wherein the protein acquires its native structure. Quenching the fluorescence by applying acrylamide, the amount 23 and 35 M(-1) were measured at pHs 8.5 and 3.0 respectively for stern-volmer constant (KSV). An increase in the enzyme molecular volume in the MG state was confirmed by gel filtration chromatography.

Evidence for Dry Molten Globule-Like Domains in the pH-Induced Equilibrium Folding Intermediate of a Multidomain Protein

The role of van der Waals (vdW) packing interactions compared to the hydrophobic effect in stabilizing the functional structure of proteins is poorly understood. Here we show, using fluorescence resonance energy transfer, dynamic fluorescence quenching, red-edge excitation shift, and near- and far-UV circular dichroism, that the pH-induced structural perturbation of a multidomain protein leads to the formation of a state in which two out of the three domains have characteristics of dry molten globules, that is, the domains are expanded compared to the native protein with disrupted packing interactions but have dry cores. We quantitatively estimate the energetic contribution of vdW interactions and show that they play an important role in the stability of the native state and cooperativity of its structural transition, in addition to the hydrophobic effect. Our results also indicate that during the pH-induced unfolding, side-chain unlocking and hydrophobic solvation occur in two distinct steps and not in a concerted manner, as commonly believed.

The Unfolding MD Simulations of Cyclophilin: Analyzed by Surface Contact Networks and Their Associated Metrics

Currently, considerable interest exists with regard to the dissociation of close packed aminoacids within proteins, in the course of unfolding, which could result in either wet or dry moltenglobules. The progressive disjuncture of residues constituting the hydrophobic core ofcyclophilin from L. donovani (LdCyp) has been studied during the thermal unfolding of the molecule, by molecular dynamics simulations. LdCyp has been represented as a surface contactnetwork (SCN) based on the surface complementarity (Sm) of interacting residues within themolecular interior. The application of Sm to side chain packing within proteins make it a very sensitive indicator of subtle perturbations in packing, in the thermal unfolding of the protein. Network based metrics have been defined to track the sequential changes in the disintegration ofthe SCN spanning the hydrophobic core of LdCyp and these metrics prove to be highly sensitive compared to traditional metrics in indicating the increased conformational (and dynamical) flexibility in the network. These metrics have been applied to suggest criteria distinguishing DMG, WMG and transition state ensembles and to identify key residues involved in crucial conformational/topological events during the unfolding process.

A hypothesis to reconcile the physical and chemical unfolding of proteins

High pressure (HP) or urea is commonly used to disturb folding species. Pressure favors the reversible unfolding of proteins by causing changes in the volumetric properties of the protein-solvent system. However, no mechanistic model has fully elucidated the effects of urea on structure unfolding, even though protein-urea interactions are considered to be crucial. Here, we provide NMR spectroscopy and 3D reconstructions from X-ray scattering to develop the "push-and-pull" hypothesis, which helps to explain the initial mechanism of chemical unfolding in light of the physical events triggered by HP. In studying MpNep2 from Moniliophthora perniciosa, we tracked two cooperative units using HP-NMR as MpNep2 moved uphill in the energy landscape; this process contrasts with the overall structural unfolding that occurs upon reaching a threshold concentration of urea. At subdenaturing concentrations of urea, we were able to trap a state in which urea is preferentially bound to the protein (as determined by NMR intensities and chemical shifts); this state is still folded and not additionally exposed to solvent [fluorescence and small-angle X-ray scattering (SAXS)]. This state has a higher susceptibility to pressure denaturation (lower p1/2 and larger ΔVu); thus, urea and HP share concomitant effects of urea binding and pulling and water-inducing pushing, respectively. These observations explain the differences between the molecular mechanisms that control the physical and chemical unfolding of proteins, thus opening up new possibilities for the study of protein folding and providing an interpretation of the nature of cooperativity in the folding and unfolding processes.

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.

All-atom molecular simulations have provided valuable insight into the workings of molecular machines and the folding and unfolding of proteins. However, commonly employed molecular dynamics simulations suffer from a limitation in accessible time scale, making it difficult to model large-scale unfolding events in a realistic amount of simulation time without employing unrealistically high temperatures. Here, we describe a rapid all-atom Monte Carlo simulation approach to simulate unfolding of the essential bacterial enzyme Dihydrofolate Reductase (DHFR) and all possible single point-mutants. We use these simulations to predict which mutants will be more thermodynamically stable (i.e., reside more often in the native folded state vs. the unfolded state) than the wild-type protein, and we confirm our predictions experimentally, creating several highly stable and catalytically active mutants. Thermally stable active engineered proteins can be used as a starting point in directed evolution experiments to evolve new functions on the background of this additional “reservoir of stability.” The stabilized enzyme may be able to accumulate a greater number of destabilizing yet functionally important mutations before unfolding, protease digestion, and aggregation abolish its activity.

Structural characterization of MG and pre-MG states of proteins by MD simulations, NMR, and other techniques

Almost all proteins fold via a number of partially structured intermediates such as molten globule (MG) and pre-molten globule states. Understanding the structure of these intermediates at atomic level is often a challenge, as these states are observed under extreme conditions of pH, temperature, and chemical denaturants. Furthermore, several other processes such as chemical modification, site-directed mutagenesis (or point mutation), and cleavage of covalent bond of natural proteins often lead to MG like partially unfolded conformation. However, the dynamic nature of proteins in these states makes them unsuitable for most structure determination at atomic level. Intermediate states studied so far have been characterized mostly by circular dichroism, fluorescence, viscosity, dynamic light scattering measurements, dye binding, infrared techniques, molecular dynamics simulations, etc. There is a limited amount of structural data available on these intermediate states by nuclear magnetic resonance (NMR) and hence there is a need to characterize these states at the molecular level. In this review, we present characterization of equilibrium intermediates by biophysical techniques with special reference to NMR.

Evidence for close side-chain packing in an early protein folding intermediate previously assumed to be a molten globule

The molten globule, a conformational ensemble with significant secondary structure but only loosely packed tertiary structure, has been suggested to be a ubiquitous intermediate in protein folding. However, it is difficult to assess the tertiary packing of transiently populated species to evaluate this hypothesis. Escherichia coli RNase H is known to populate an intermediate before the rate-limiting barrier to folding that has long been thought to be a molten globule. We investigated this hypothesis by making mimics of the intermediate that are the ground-state conformation at equilibrium, using two approaches: a truncation to generate a fragment mimic of the intermediate, and selective destabilization of the native state using point mutations. Spectroscopic characterization and the response of the mimics to further mutation are consistent with studies on the transient kinetic intermediate, indicating that they model the early intermediate. Both mimics fold cooperatively and exhibit NMR spectra indicative of a closely packed conformation, in contrast to the hypothesis of molten tertiary packing. This result is important for understanding the nature of the subsequent rate-limiting barrier to folding and has implications for the assumption that many other proteins populate molten globule folding intermediates.

WW domain folding complexity revealed by infrared spectroscopy

Although the intrinsic tryptophan fluorescence of proteins offers a convenient probe of protein folding, interpretation of the fluorescence spectrum is often difficult because it is sensitive to both global and local changes. Infrared (IR) spectroscopy offers a complementary measure of structural changes involved in protein folding, because it probes changes in the secondary structure of the protein backbone. Here we demonstrate the advantages of using multiple probes, infrared and fluorescence spectroscopy, to study the folding of the FBP28 WW domain. Laser-induced temperature jumps coupled with fluorescence or infrared spectroscopy have been used to probe changes in the peptide backbone on the submillisecond time scale. The relaxation dynamics of the β-sheets and β-turn were measured independently by probing the corresponding IR bands assigned in the amide I region. Using these wavelength-dependent measurements, we observe three kinetics phases, with the fastest process corresponding to the relaxation kinetics of the turns. In contrast, fluorescence measurements of the wild-type WW domain and tryptophan mutants exhibit single-exponential kinetics with a lifetime that corresponds to the slowest phase observed by infrared spectroscopy. Mutant sequences provide evidence of an intermediate dry molten globule state. The slowest step in the folding of this WW domain is the tight packing of the side chains in the transition from the dry molten globule intermediate to the native structure. This study demonstrates that using multiple complementary probes enhances the interpretation of protein folding dynamics.

Unfolding of a small protein proceeds via dry and wet globules and a solvated transition state

Dissecting a protein unfolding process into individual steps can provide valuable information on the forces that maintain the integrity of the folded structure. Solvation of the protein core determines stability, but it is not clear when such solvation occurs during unfolding. In this study, far-UV circular dichroism measurements suggest a simplistic two-state view of the unfolding of barstar, but the use of multiple other probes brings out the complexity of the unfolding reaction. Near-UV circular dichroism measurements show that unfolding commences with the loosening of tertiary interactions in a native-like intermediate, N(∗). Fluorescence resonance energy transfer measurements show that N(∗) then expands rapidly but partially to form an early unfolding intermediate IE. Fluorescence spectral measurements indicate that both N(∗) and IE have retained native-like solvent accessibility of the core, suggesting that they are dry molten globules. Dynamic quenching measurements at the single tryptophan buried in the core suggest that the core becomes solvated only later in a late wet molten globule, IL, which precedes the unfolded form. Fluorescence anisotropy decay measurements show that tight packing around the core tryptophan is lost when IL forms. Of importance, the slowest step is unfolding of the wet molten globule and involves a solvated transition state.