Consequences of stabilizing the natively disordered f helix for the folding pathway of apomyoglobin

Non-Native Structures of Apomyoglobin and Apoleghemoglobin in Folding Intermediates Related to the Protein Misfolding

Protein folding is essential for a polypeptide chain to acquire its proper structure and function. Globins are a superfamily of ubiquitous heme-binding α-helical proteins whose function is principally to regulate oxygen homoeostasis. In this review, we explore the hierarchical helical formation in the globin proteins apomyoglobin and leghemoglobin, and we discuss the existence of non-native and misfolded structures occurring during the course of folding to its native state. This review summarizes the research aimed at characterizing and comparing the equilibrium and kinetic intermediates, as well as delineating the complete folding pathway at a molecular level, in order to answer the following questions: "What is the mechanism of misfolding via a folding intermediate? Does the non-native structure stabilize the contemporary intermediate structure? Does the non-native structure induce slower folding?" The role of the non-native structures in the folding intermediate related to misfolding is also discussed.

Complex Folding Landscape of Apomyoglobin at Acidic pH Revealed by Ultrafast Kinetic Analysis of Core Mutants

Under mildly acidic conditions (pH 4-4.5) apomyoglobin (apoMb) adopts a partially structured equilibrium state ( M-state) that structurally resembles a kinetic intermediate encountered at a late stage of folding to the native structure at neutral pH. We have previously reported that the M-state is formed rapidly (<1 ms) via a multistate process and thus offers a unique opportunity for exploring early stages of folding by both experimental and computational techniques. In order to gain structural insight into intermediates and barriers at the residue level, we studied the folding/unfolding kinetics of 12 apoMb mutants at pH 4.2 using fluorescence-detected ultrafast mixing techniques. Global analysis of the submillisecond folding/unfolding kinetics vs urea concentration for each variant, based on a sequential four-state mechanism ( U ⇔ I ⇔ L ⇔ M), allowed us to determine elementary rate constants and their dependence on urea concentration for most transitions. Comparison of the free energy diagrams constructed from the kinetic data of the mutants with that of wild-type apoMb yielded quantitative information on the effects of mutations on the free energy (ΔΔ G) of both intermediates and the first two kinetic barriers encountered during folding. Truncation of conserved aliphatic side chains on helices A, G, and H gives rise to a stepwise increase in ΔΔ G as the protein advances from U toward M, consistent with progressive stabilization of native-like contacts within the primary core of apoMb. Helix-helix contacts in the primary core contribute little to the first folding barrier ( U ⇔ I) and thus are not required for folding initiation but are critical for the stability of the late intermediate, L, and the M-state. Alanine substitution of hydrophobic residues at more peripheral helix-helix contact sites of the native structure, which are still absent or unstable in the M-state, shows both positive (destabilizing) and negative (stabilizing) ΔΔ G, indicating that non-native contacts are formed initially and weakened or lost as a result of subsequent structural rearrangement steps.

Tracing a protein's folding pathway over evolutionary time using ancestral sequence reconstruction and hydrogen exchange

The conformations populated during protein folding have been studied for decades; yet, their evolutionary importance remains largely unexplored. Ancestral sequence reconstruction allows access to proteins across evolutionary time, and new methods such as pulsed-labeling hydrogen exchange coupled with mass spectrometry allow determination of folding intermediate structures at near amino-acid resolution. Here, we combine these techniques to monitor the folding of the ribonuclease H family along the evolutionary lineages of T. thermophilus and E. coli RNase H. All homologs and ancestral proteins studied populate a similar folding intermediate despite being separated by billions of years of evolution. Even though this conformation is conserved, the pathway leading to it has diverged over evolutionary time, and rational mutations can alter this trajectory. Our results demonstrate that evolutionary processes can affect the energy landscape to preserve or alter specific features of a protein’s folding pathway.

Non-Native α-Helices in the Initial Folding Intermediate Facilitate the Ordered Assembly of the β-Barrel in β-Lactoglobulin

The roles of non-native α-helices frequently observed in the initial folding stage of β-sheet proteins have been examined for many years. We herein investigated the residue-level structures of several mutants of bovine β-lactoglobulin (βLG) in quenched-flow pH-pulse labeling experiments. βLG assumes a collapsed intermediate with a non-native α-helical structure (I₀) in the early stage of folding, although its native form is predominantly composed of β-structures. The protection profile in I₀ of pseudo-wild type (WT*) βLG was found to deviate from the pattern of the "average area buried upon folding" (AABUF). In particular, the level of protection at the region of strand A, at which non-native α-helices form in the I₀ state, was significantly low compared to AABUF. G17E, the mutant with an increased helical propensity, showed a similar protection pattern. In contrast, the protection pattern for I₀ of E44L, the mutant with an increased β-sheet propensity, was distinct from that of WT* and resembled the AABUF pattern. Transverse relaxation measurements demonstrated that the positions of the residual structures in the unfolded states of these mutants were consistent with those of the protected residues in the respective I₀ states. On the basis of the slower conversion of I₀ to the native state for E44L to that for WT*, non-native α-helices facilitate the ordered assembly of the β-barrel by preventing interactions that trap folding.

Exploring the Sequence-based Prediction of Folding Initiation Sites in Proteins

Protein folding is a complex process that can lead to disease when it fails. Especially poorly understood are the very early stages of protein folding, which are likely defined by intrinsic local interactions between amino acids close to each other in the protein sequence. We here present EFoldMine, a method that predicts, from the primary amino acid sequence of a protein, which amino acids are likely involved in early folding events. The method is based on early folding data from hydrogen deuterium exchange (HDX) data from NMR pulsed labelling experiments, and uses backbone and sidechain dynamics as well as secondary structure propensities as features. The EFoldMine predictions give insights into the folding process, as illustrated by a qualitative comparison with independent experimental observations. Furthermore, on a quantitative proteome scale, the predicted early folding residues tend to become the residues that interact the most in the folded structure, and they are often residues that display evolutionary covariation. The connection of the EFoldMine predictions with both folding pathway data and the folded protein structure suggests that the initial statistical behavior of the protein chain with respect to local structure formation has a lasting effect on its subsequent states.

How Does Your Protein Fold? Elucidating the Apomyoglobin Folding Pathway

Although each type of protein fold and in some cases individual proteins within a fold classification can have very different mechanisms of folding, the underlying biophysical and biochemical principles that operate to cause a linear polypeptide chain to fold into a globular structure must be the same. In an aqueous solution, the protein takes up the thermodynamically most stable structure, but the pathway along which the polypeptide proceeds in order to reach that structure is a function of the amino acid sequence, which must be the final determining factor, not only in shaping the final folded structure, but in dictating the folding pathway. A number of groups have focused on a single protein or group of proteins, to determine in detail the factors that influence the rate and mechanism of folding in a defined system, with the hope that hypothesis-driven experiments can elucidate the underlying principles governing the folding process. Our research group has focused on the folding of the globin family of proteins, and in particular on the monomeric protein apomyoglobin. Apomyoglobin (apoMb) folds relatively slowly (∼2 s) via an ensemble of obligatory intermediates that form rapidly after the initiation of folding. The folding pathway can be dissected using rapid-mixing techniques, which can probe processes in the millisecond time range. Stopped-flow measurements detected by circular dichroism (CD) or fluorescence spectroscopy give information on the rates of folding events. Quench-flow experiments utilize the differential rates of hydrogen-deuterium exchange of amide protons protected in parts of the structure that are folded early; protection of amides can be detected by mass spectrometry or proton nuclear magnetic resonance spectroscopy (NMR). In addition, apoMb forms an intermediate at equilibrium at pH ∼ 4, which is sufficiently stable for it to be structurally characterized by solution methods such as CD, fluorescence and NMR spectroscopies, and the conformational ensembles formed in the presence of denaturing agents and low pH can be characterized as models for the unfolded states of the protein. Newer NMR techniques such as measurement of residual dipolar couplings in the various partly folded states, and relaxation dispersion measurements to probe invisible states present at low concentrations, have contributed to providing a detailed picture of the apomyoglobin folding pathway. The research summarized in this Account was aimed at characterizing and comparing the equilibrium and kinetic intermediates both structurally and dynamically, as well as delineating the complete folding pathway at a residue-specific level, in order to answer the question: "What is it about the amino acid sequence that causes each molecule in the unfolded protein ensemble to start folding, and, once started, to proceed towards the formation of the correctly folded three-dimensional structure?"

Folding of apomyoglobin: Analysis of transient intermediate structure during refolding using quick hydrogen deuterium exchange and NMR

The structures of apomyoglobin folding intermediates have been widely analyzed using physical chemistry methods including fluorescence, circular dichroism, small angle X-ray scattering, NMR, mass spectrometry, and rapid mixing. So far, at least two intermediates (on sub-millisecond- and millisecond-scales) have been demonstrated for apomyoglobin folding. The combination of pH-pulse labeling and NMR is a useful tool for analyzing the kinetic intermediates at the atomic level. Its use has revealed that the latter-phase kinetic intermediate of apomyoglobin (6 ms) was composed of helices A, B, G and H, whereas the equilibrium intermediate, called the pH 4 molten-globule intermediate, was composed mainly of helices A, G and H. The improved strategy for the analysis of the kinetic intermediate was developed to include (1) the dimethyl sulfoxide method, (2) data processing with the various labeling times, and (3) a new in-house mixer. Particularly, the rapid mixing revealed that helices A and G were significantly more protected at the earlier stage (400 µs) of the intermediate (former-phase intermediate) than the other helices. Mutation studies, where each hydrophobic residue was replaced with an alanine in helices A, B, E, F, G and H, indicated that both non-native and native-like structures exist in the latter-phase folding intermediate. The N-terminal part of helix B is a weak point in the intermediate, and the docking of helix E residues to the core of the A, B, G and H helices was interrupted by a premature helix B, resulting in the accumulation of the intermediate composed of helices A, B, G and H. The prediction-based protein engineering produced important mutants: Helix F in a P88K/A90L/S92K/A94L mutant folded in the latter-phase intermediate, although helix F in the wild type does not fold even at the native state. Furthermore, in the L11G/W14G/A70L/G73W mutant, helix A did not fold but helix E did, which is similar to what was observed in the kinetic intermediate of apoleghemoglobin. Thus, this protein engineering resulted in a changed structure for the apomyoglobin folding intermediate.

Thermal unfolding of myoglobin in the Landau-Ginzburg-Wilson approach

The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature crystallographic Cα backbone structure of sperm whale myoglobin. The Glauber protocol is employed to simulate its response to an increase in ambient temperature. The myoglobin is found to unfold from its native state by a succession of α-helical intermediates, fully in line with the observed folding and unfolding patterns in denaturation experiments. In particular, a molten globule intermediate is identified with experimentally correct attributes. A detailed, experimentally testable contact map is constructed to characterize the specifics of the unfolding pathway, including the formation of long-range interactions. The results reveal how the unfolding process of a protein is driven by the interplay between, and a successive melting of, its modular secondary structure components.

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native-like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse-grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel.

Probing the non-native H helix translocation in apomyoglobin folding intermediates

Apomyoglobin folds via sequential helical intermediates that are formed by rapid collapse of the A, B, G, and H helix regions. An equilibrium molten globule with a similar structure is formed near pH 4. Previous studies suggested that the folding intermediates are kinetically trapped states in which folding is impeded by non-native packing of the G and H helices. Fluorescence spectra of mutant proteins in which cysteine residues were introduced at several positions in the G and H helices show differential quenching of W14 fluorescence, providing direct evidence of translocation of the H helix relative to helices A and G in both the kinetic and equilibrium intermediates. Förster resonance energy transfer measurements show that a 5-({2-[(acetyl)amino]ethyl}amino)naphthalene-1-sulfonic acid acceptor coupled to K140C (helix H) is closer to Trp14 (helix A) in the equilibrium molten globule than in the native state, by a distance that is consistent with sliding of the H helix in an N-terminal direction by approximately one helical turn. Formation of an S108C–L135C disulfide prevents H helix translocation in the equilibrium molten globule by locking the G and H helices into their native register. By enforcing nativelike packing of the A, G, and H helices, the disulfide resolves local energetic frustration and facilitates transient docking of the E helix region onto the hydrophobic core but has only a small effect on the refolding rate. The apomyoglobin folding landscape is highly rugged, with several energetic bottlenecks that frustrate folding; relief of any one of the major identified bottlenecks is insufficient to speed progression to the transition state.

Unfolding simulations of holomyoglobin from four mammals: identification of intermediates and β-sheet formation from partially unfolded states

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations.

Sub-millisecond chain collapse of the Escherichia coli globin ApoHmpH

Myoglobins are ubiquitous proteins that play a seminal role in oxygen storage, transport, and NO metabolism. The folding mechanism of apomyoglobins from different species has been studied to a fair extent over the last two decades. However, integrated investigations of the entire process, including both the early (sub-ms) and late (ms-s) folding stages, have been missing. Here, we study the folding kinetics of the single-Trp Escherichia coli globin apoHmpH via a combination of continuous-flow microfluidic and stopped-flow approaches. A rich series of molecular events emerges, spanning a very wide temporal range covering more than 7 orders of magnitude, from sub-microseconds to tens of seconds. Variations in fluorescence intensity and spectral shifts reveal that the protein region around Trp120 undergoes a fast collapse within the 8 μs mixing time and gradually reaches a native-like conformation with a half-life of 144 μs from refolding initiation. There are no further fluorescence changes beyond ca. 800 μs, and folding proceeds much more slowly, up to 20 s, with acquisition of the missing helicity (ca. 30%), long after consolidation of core compaction. The picture that emerges is a gradual acquisition of native structure on a free-energy landscape with few large barriers. Interestingly, the single tryptophan, which lies within the main folding core of globins, senses some local structural consolidation events after establishment of native-like core polarity (i.e., likely after core dedydration). In all, this work highlights how the main core of the globin fold is capable of becoming fully native efficiently, on the sub-millisecond time scale.

Direct ubiquitin independent recognition and degradation of a folded protein by the eukaryotic proteasomes-origin of intrinsic degradation signals

Eukaryotic 26S proteasomes are structurally organized to recognize, unfold and degrade globular proteins. However, all existing model substrates of the 26S proteasome in addition to ubiquitin or adaptor proteins require unstructured regions in the form of fusion tags for efficient degradation. We report for the first time that purified 26S proteasome can directly recognize and degrade apomyoglobin, a globular protein, in the absence of ubiquitin, extrinsic degradation tags or adaptor proteins. Despite a high affinity interaction, absence of a ligand and presence of only helices/loops that follow the degradation signal, apomyoglobin is degraded slowly by the proteasome. A short floppy F-helix exposed upon ligand removal and in conformational equilibrium with a disordered structure is mandatory for recognition and initiation of degradation. Holomyoglobin, in which the helix is buried, is neither recognized nor degraded. Exposure of the floppy F-helix seems to sensitize the proteasome and primes the substrate for degradation. Using peptide panning and competition experiments we speculate that initial encounters through the floppy helix and additional strong interactions with N-terminal helices anchors apomyoglobin to the proteasome. Stabilizing helical structure in the floppy F-helix slows down degradation. Destabilization of adjacent helices accelerates degradation. Unfolding seems to follow the mechanism of helix unraveling rather than global unfolding. Our findings while confirming the requirement for unstructured regions in degradation offers the following new insights: a) origin and identification of an intrinsic degradation signal in the substrate, b) identification of sequences in the native substrate that are likely to be responsible for direct interactions with the proteasome, and c) identification of critical rate limiting steps like exposure of the intrinsic degron and destabilization of an unfolding intermediate that are presumably catalyzed by the ATPases. Apomyoglobin emerges as a new model substrate to further explore the role of ATPases and protein structure in proteasomal degradation.