Effect of H helix destabilizing mutations on the kinetic and equilibrium folding of apomyoglobin

Mapping Protein-Protein Interactions at Birth: Single-Particle Cryo-EM Analysis of a Ribosome-Nascent Globin Complex

Interactions between ribosome-bound nascent chains (RNCs) and ribosomal components are critical to elucidate the mechanism of cotranslational protein folding. Nascent protein-ribosome contacts within the ribosomal exit tunnel were previously assessed mostly in the presence of C-terminal stalling sequences, yet little is known about contacts taking place in the absence of these strongly interacting motifs. Further, there is nearly no information about ribosomal proteins (r-proteins) interacting with nascent chains within the outer surface of the ribosome. Here, we combine chemical cross-linking, single-particle cryo-EM, and fluorescence anisotropy decays to determine the structural features of ribosome-bound apomyoglobin (apoMb). Within the ribosomal exit tunnel core, interactions are similar to those identified in previous reports. However, once the RNC enters the tunnel vestibule, it becomes more dynamic and interacts with ribosomal RNA (rRNA) and the L23 r-protein. Remarkably, on the outer surface of the ribosome, RNCs interact mainly with a highly conserved nonpolar patch of the L23 r-protein. RNCs also comprise a compact and dynamic N-terminal region lacking contact with the ribosome. In all, apoMb traverses the ribosome and interacts with it via its C-terminal region, while N-terminal residues sample conformational space and form a compact subdomain before the entire nascent protein sequence departs from the ribosome.

From Immunogenic Peptides to Intrinsically Disordered Proteins

It is hard to evaluate the role of individual mentors in the genesis of important ideas. In the case of our realization that proteins do not have to be stably folded to be functional, the influence of Richard Lerner and our collaborative work in the 1980s on the conformations of immunogenic peptides provided a base level of thinking about the nature of polypeptides in water solutions that led us to formulate and develop our ideas on the importance of intrinsic disorder in proteins. This review describes how the insights gained into the behavior of peptides led directly to the realization that proteins were not only capable of being functional while disordered, but also that disorder provided a distinct functional advantage in many important cellular processes.

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

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

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.

Dimer domain swapping versus monomer folding in apo-myoglobin studied by molecular simulations

Using a coarse-grained symmetrized Go model, we performed a series of folding simulations of two apo-myoglobin molecules restrained at a high density, addressing competition of formation of a domain-swapped dimer with folding to two monomer structures.

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.

Independent of their localization in protein the hydrophobic amino acid residues have no effect on the molten globule state of apomyoglobin and the disulfide bond on the surface of apomyoglobin stabilizes this intermediate state

At present it is unclear which interactions in proteins reveal the presence of intermediate states, their stability and formation rate. In this study, we have investigated the effect of substitutions of hydrophobic amino acid residues in the hydrophobic core of protein and on its surface on a molten globule type intermediate state of apomyoglobin. It has been found that independent of their localization in protein, substitutions of hydrophobic amino acid residues do not affect the stability of the molten globule state of apomyoglobin. It has been shown also that introduction of a disulfide bond on the protein surface can stabilize the molten globule state. However in the case of apomyoglobin, stabilization of the intermediate state leads to relative destabilization of the native state of apomyoglobin. The result obtained allows us not only to conclude which mutations can have an effect on the intermediate state of the molten globule type, but also explains why the introduction of a disulfide bond (which seems to “strengthen” the protein) can result in destabilization of the protein native state of apomyoglobin.

Sequence analysis on the information of folding initiation segments in ferredoxin-like fold proteins

BACKGROUND

While some studies have shown that the 3D protein structures are more conservative than their amino acid sequences, other experimental studies have shown that even if two proteins share the same topology, they may have different folding pathways. There are many studies investigating this issue with molecular dynamics or Go-like model simulations, however, one should be able to obtain the same information by analyzing the proteins' amino acid sequences, if the sequences contain all the information about the 3D structures. In this study, we use information about protein sequences to predict the location of their folding segments. We focus on proteins with a ferredoxin-like fold, which has a characteristic topology. Some of these proteins have different folding segments.

RESULTS

Despite the simplicity of our methods, we are able to correctly determine the experimentally identified folding segments by predicting the location of the compact regions considered to play an important role in structural formation. We also apply our sequence analyses to some homologues of each protein and confirm that there are highly conserved folding segments despite the homologues' sequence diversity. These homologues have similar folding segments even though the homology of two proteins' sequences is not so high.

CONCLUSION

Our analyses have proven useful for investigating the common or different folding features of the proteins studied.

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.

Polymer uncrossing and knotting in protein folding, and their role in minimal folding pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold [Formula: see text], [Formula: see text], [Formula: see text], and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation "alignment". The consensus minimal pathway is constructed and shown schematically for representative cases of an [Formula: see text], [Formula: see text], and knotted protein. An overlap parameter is defined between pathways; we find that [Formula: see text] proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and [Formula: see text] proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.

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

The F helix region of sperm whale apomyoglobin is disordered, undergoing conformational fluctuations between a folded helical conformation and one or more locally unfolded states. To examine the effects of F helix stabilization on the folding pathway of apomyoglobin, we have introduced mutations to augment intrinsic helical structure in the F helix of the kinetic folding intermediate and to increase its propensity to fold early in the pathway, using predictions based on plots of the average area buried upon folding (AABUF) derived from the primary sequence. Two mutant proteins were prepared: a double mutant, P88K/S92K (F2), and a quadruple mutant, P88K/A90L/S92K/A94L (F4). Whereas the AABUF for F2 predicts that the F helix will not fold early in the pathway, the F helix in F4 shows a significantly increased AABUF and is therefore predicted to fold early. Protection of amide protons by formation of hydrogen-bonded helical structure during the early folding events has been analyzed by pH-pulse labeling. Consistent with the AABUF prediction, many of the F helix residues for F4 are significantly protected in the kinetic intermediate but are not protected in the F2 mutant. F4 folds via a kinetically trapped burst-phase intermediate that contains stabilized secondary structure in the A, B, F, G, and H helix regions. Rapid folding of the F helix stabilizes the central core of the misfolded intermediate and inhibits translocation of the H helix back to its native position, thereby decreasing the overall folding rate.

Protein folding at the exit tunnel

Over five decades of research have yielded a large body of information on how purified proteins attain their native state when refolded in the test tube, starting from a chemically or thermally denatured state. Nevertheless, we still know little about how proteins fold and unfold in their natural biological habitat: the living cell. Indeed, a variety of cellular components, including molecular chaperones, the ribosome, and crowding of the intracellular medium, modulate folding mechanisms in physiologically relevant environments. This review focuses on the current state of knowledge in protein folding in the cell with emphasis on the early stage of a protein's life, as the nascent polypeptide traverses and emerges from the ribosomal tunnel. Given the vectorial nature of ribosome-assisted translation, the transient degree of chain elongation becomes a relevant variable expected to affect nascent protein foldability, aggregation propensity and extent of interaction with chaperones and the ribosome.

pH-induced equilibrium unfolding of apomyoglobin: substitutions at conserved Trp14 and Met131 and non-conserved Val17 positions

A number of residues in globins family are well conserved but are not directly involved in the primary oxygen-carrying function of these proteins. A possible role for these conserved, non-functional residues has been suggested in promoting a rapid and correct folding process to the native tertiary structure. To test this hypothesis, we have studied pH-induced equilibrium unfolding of mutant apomyoglobins with substitutions of the conserved residues Trp14 and Met131, which are not involved in the function of myoglobin, by various amino acids. This allowed estimating their impact on the stability of various conformational states of the proteins and selecting conditions for a folding kinetics study. The results obtained from circular dichroism, tryptophan fluorescence, and differential scanning microcalorimetry for these mutant proteins were compared with those for the wild type protein and for a mutant with the non-conserved Val17 substituted by Ala. In the native folded state, all of the mutant apoproteins have a compact globular structure, but are destabilized in comparison to the wild type protein. The pH-induced denaturation of the mutant proteins occurs through the formation of a molten globule-like intermediate similar to that of the wild type protein. Thermodynamic parameters for all of the proteins were calculated using the three state model. Stability of equilibrium intermediates at pH ~4.0 was shown to be slightly affected by the mutations. Thus, all of the above substitutions influence the stability of the native state of these proteins. The cooperativity of conformational transitions and the exposed to solvent protein surface were also changed, but not for the substitution at Val17.

Hierarchical folding mechanism of apomyoglobin revealed by ultra-fast H/D exchange coupled with 2D NMR

The earliest steps in the folding of proteins are complete on an extremely rapid time scale that is difficult to access experimentally. We have used rapid-mixing quench-flow methods to extend the time resolution of folding studies on apomyoglobin and elucidate the structural and dynamic features of members of the ensemble of intermediate states that are populated on a submillisecond time scale during this process. The picture that emerges is of a continuum of rapidly interconverting states. Even after only 0.4 ms of refolding time a compact state is formed that contains major parts of the A, G, and H helices, which are sufficiently well folded to protect amides from exchange. The B, C, and E helix regions fold more slowly and fluctuate rapidly between open and closed states as they search docking sites on this core; the secondary structure in these regions becomes stabilized as the refolding time is increased from 0.4 to 6 ms. No further stabilization occurs in the A, G, H core at 6 ms of folding time. These studies begin to time-resolve a progression of compact states between the fully unfolded and native folded states and confirm the presence an ensemble of intermediates that interconvert in a hierarchical sequence as the protein searches conformational space on its folding trajectory.

Thermodynamic and kinetic characterization of apoHmpH, a fast-folding bacterial globin

Despite the widespread presence of the globin fold in most living organisms, only eukaryotic globins have been employed as model proteins in folding/stability studies so far. This work introduces the first thermodynamic and kinetic characterization of a prokaryotic globin, that is, the apo form of the heme-binding domain of flavohemoglobin (apoHmpH) from Escherichia coli. This bacterial globin has a widely different sequence but nearly identical structure to its eukaryotic analogues. We show that apoHmpH is a well-folded monomeric protein with moderate stability at room temperature [apparent Delta G degrees (UN(w))=-3.1+/-0.3 kcal mol(-1); m(UN)=-1.7 kcal mol(-1) M(-1)] and predominant alpha-helical structure. Remarkably, apoHmpH is the fastest-folding globin known to date, as it refolds about 4- to 16-fold more rapidly than its eukaryotic analogues (e.g., sperm whale apomyoglobin and soybean apoleghemoglobin), populating a compact kinetic intermediate (beta(I)=0.9+/-0.2) with significant helical content. Additionally, the single Trp120 (located in the native H helix) becomes locked into a fully native-like environment within 6 ms, suggesting that this residue and its closest spatial neighbors complete their folding at ultrafast (submillisecond) speed. In summary, apoHmpH is a bacterial globin that shares the general folding scheme (i.e., a rapid burst phase followed by slower rate-determining phases) of its eukaryotic analogues but displays an overall faster folding and a kinetic intermediate with some fully native-like traits. This study supports the view that the general folding features of bacterial and eukaryotic globins are preserved through evolution while kinetic details differ.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Binding specificity of an alpha-helical protein sequence to a full-length Hsp70 chaperone and its minimal substrate-binding domain

Hsp70 chaperones are involved in the prevention of misfolding, and possibly the folding, of newly synthesized proteins. The members of this chaperone family are capable of interacting with polypeptide chains both co- and posttranslationally, but it is currently not clear how different structural domains of the chaperone affect binding specificity. We explored the interactions between the bacterial Hsp70, DnaK, and the sequence of a model all-alpha-helical globin (apoMb) by cellulose-bound peptide scanning. The binding specificity of the full-length chaperone was compared with that of its minimal substrate-binding domain, DnaK-beta. Six specific chaperone binding sites evenly distributed along the apoMb sequence were identified. Binding site locations are identical for the full-length chaperone and its substrate-binding domain, but relative affinities differ. The binding specificity of DnaK-beta is only slightly decreased relative to that of full-length DnaK. DnaK's binding motif is known to comprise hydrophobic regions flanked by positively charged residues. We found that the simple fractional mean buried area correlates well with Hsp70's binding site locations along the apoMb sequence. In order to further characterize the properties of the minimal binding host, the stability of DnaK-beta upon chemical denaturation by urea and protons was investigated. Urea unfolding titrations yielded an apparent folding DeltaG degrees of 3.1 +/- 0.9 kcal mol-1 and an m value of 1.7 +/- 0.4 kcal mol-1 M-1.

The role of hydrophobic interactions in initiation and propagation of protein folding

Globular proteins fold by minimizing the nonpolar surface that is exposed to water, while simultaneously providing hydrogen-bonding interactions for buried backbone groups, usually in the form of secondary structures such as alpha-helices, beta-sheets, and tight turns. A primary thermodynamic driving force for the formation of globular structure is thus the sequestration of nonpolar groups, but the correlation between the parts of proteins that are observed to fold first (termed folding initiation sites) and the "hydrophobicity" (as customarily defined) of the amino acids in these regions has been quite weak. It has previously been noted that many amino acid side chains contain considerable nonpolar sections, even if they also contain polar or charged groups. For example, a lysine side chain contains four methylenes, which may undergo hydrophobic interactions if the charged epsilon-NH(3)(+) group is salt-bridged or hydrogen-bonded. Folding initiation sites might therefore contain not only accepted "hydrophobic" amino acids, but also larger charged side chains. Recent experiments on the folding of mutant apomyoglobins provides corroboration for models based on the hypothesis that folding initiation sites arise from hydrophobic interactions. A near-perfect correlation was observed between the areas of the molecule that are present in the burst-phase kinetic intermediate and both the free energy of formation of hydrophobic initiation sites and the parameter "average area buried upon folding," which pinpoints large side chains, even those containing charged or polar portions. These results provide a putative mechanism for the control of protein-folding initiation and growth by polar/nonpolar sequence propensity alone.

Diffusion-collision of foldons elucidates the kinetic effects of point mutations and suggests control strategies of the folding process of helical proteins

In this article we use mutation studies as a benchmark for a minimal model of the folding process of helical proteins. The model ascribes a pivotal role to the collisional dynamics of a few crucial residues (foldons) and predicts the folding rates by exploiting information drawn from the protein sequence. We show that our model rationalizes the effects of point mutations on the kinetics of folding. The folding times of two proteins and their mutants are predicted. Stability and location of foldons have a critical role as the determinants of protein folding. This allows us to elucidate two main mechanisms for the kinetic effects of mutations. First, it turns out that the mutations eliciting the most notable effects alter protein stability through stabilization or destabilization of the foldons. Secondly, the folding rate is affected via a modification of the foldon topology by those mutations that lead to the birth or death of foldons. The few mispredicted folding rates of some mutants hint at the limits of the current version of the folding model proposed in the present article. The performance of our folding model declines in case the mutated residues are subject to strong long-range forces. That foldons are the critical targets of mutation studies has notable implications for design strategies and is of particular interest to address the issue of the kinetic regulation of single proteins in the general context of the overall dynamics of the interactome.

Residue-level NMR View of the Urea-driven Equilibrium Folding Transition of SUMO-1 (1-97): Native Preferences Do Not Increase Monotonously

SUMO-1 (1-97) is a crucial protein in the machinery of post-translational modifications. We observed by circular dichroism and fluorescence spectroscopy that urea-induced unfolding of this protein is a complex process with the possibility of occurrence of detectable intermediates along the way. The tertiary structure is completely lost around approximately 4.5 M urea with a transition mid-point at 2.53 M urea, while the secondary structure unfolding seems to show two transitions, with mid-points at 2.42 M and 5.69 M urea. We have elucidated by systematic urea titration, the equilibrium residue level structural and dynamics changes along the entire folding/unfolding transition by multidimensional NMR. With urea dilution, the protein is seen to progressively lose most of the broad beta-domain structural preferences present at 8 M urea, acquire some helical propensities at 5 M urea, and lose some of them again on further dilution of urea. Between 3 M and 2 M urea, the protein starts afresh to acquire native structural features. These observations are contrary to the conventional notion that proteins fold with monotonously increasing native-type preferences. For folding below approximately 3 M urea, the region around the alpha1 helix appears to be a potential folding initiation site. The folding seems to start with a collapse into native-like topologies, at least in parts, and is followed by formation of secondary and tertiary structure, perhaps by cooperative rearrangements. The motional characteristics of the protein show sequence-dependent variation as the concentration of urea is progressively reduced. At the sub-nanosecond level, the features are extremely unusual for denatured states, and only certain segments corresponding to the flexible regions in the native protein display these motions at the different concentrations of urea.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

Removal of kinetic traps and enhanced protein folding by strategic substitution of amino acids in a model α-helical hairpin peptide

The presence of non-native kinetic traps in the free energy landscape of a protein may significantly lengthen the overall folding time so that the folding process becomes unreliable. We use a computational model alpha-helical hairpin peptide to calculate structural free energy landscapes and relate them to the kinetics of folding. We show how protein engineering through strategic changes in only a few amino acid residues along the primary sequence can greatly increase the speed and reliability of the folding process, as seen experimentally. These strategic substitutions also prevent the formation of long-lived misfolded configurations that can cause unwanted aggregations of peptides. These results support arguments that removal of kinetic traps, obligatory or nonobligatory, is crucial for fast folding.

In vitro expression and characterization of native apomyoglobin under low molecular crowding conditions

The labile nature of membranes and organelles poses serious challenges to in situ biomolecule characterization in intact cells. Cell-free in vitro systems provide an alternative promising medium for the expression and characterization of protein conformation and function in a biochemical context that bears several similarities to the cellular environment. In addition, cell-free transcription-translation has recently emerged as a convenient method for protein selective isotope labeling, providing significant advantages for detailed NMR analysis. We report the cell-free expression of the model protein apomyoglobin (apoMb) in an Escherichia coli cell-free system and the effect of polyethylene glycol (PEG) on the expression yields. In contrast with in vivo protein production under control of the strong T7 promoter, apoMb is expressed in vitro in 100% soluble form. In-gel tryptic digestion followed by mass spectrometry were performed to confirm the protein identity. In order to probe the conformation of the newly expressed protein and investigate the feasibility of in situ structural analysis, high resolution protein characterization was carried out by 2D NMR spectroscopy. In vitro apoMb expression in a PEG-free environment is a convenient method for the production of soluble native-like protein under conditions amenable to selective isotopic labeling. Yields can be easily scaled-up by dialysis-assisted cell-free expression.

Identification of native and non-native structure in kinetic folding intermediates of apomyoglobin

Site-directed mutagenesis has been used to probe the interactions that stabilize the equilibrium and burst phase kinetic intermediates formed by apomyoglobin. Nine bulky hydrophobic residues in the A, E, G and H helices were replaced by alanine, and the effects on protein stability and kinetic folding pathways were determined. Hydrogen exchange pulse-labeling experiments, with NMR detection, were performed for all mutants. All of the alanine substitutions resulted in changes in proton occupancy or an increased rate of hydrogen-deuterium exchange for amides in the immediate vicinity of the mutation. In addition, most mutations affected residues in distant parts of the amino acid sequence, providing insights into the topology of the burst phase intermediate and the interactions that stabilize its structure. Differences between the pH 4 equilibrium molten globule and the kinetic intermediate are evident: the E helix region plays no discernible role in the equilibrium intermediate, but contributes significantly to stabilization of the ensemble of compact intermediates formed during kinetic refolding. Mutations that interfere with docking of the E helix onto the preformed A/B/G/H helix core substantially decrease the folding rate, indicating that docking and folding of the E helix region occurs prior to formation of the apomyoglobin folding transition state. The results of the mutagenesis experiments are consistent with rapid formation of an ensemble of compact burst phase intermediates with an overall native-like topological arrangement of the A, B, E, G, and H helices. However, the experiments also point to disorder in docking of the E helix and to non-native contacts in the kinetic intermediate. In particular, there is evidence for translocation of the H helix by approximately one helical turn towards its N terminus to maximize hydrophobic interactions with helix G. Thus, the burst phase intermediate observed during kinetic refolding of apomyoglobin consists of an ensemble of compact, kinetically trapped states in which the helix docking appears to be topologically correct, but in which there are local non-native interactions that must be resolved before the protein can fold to the native structure.

Structural characterization of apomyoglobin self-associated species in aqueous buffer and urea solution

The biophysical characterization of nonfunctional protein aggregates at physiologically relevant temperatures is much needed to gain deeper insights into the kinetic and thermodynamic relationships between protein folding and misfolding. Dynamic and static laser light scattering have been employed for the detection and detailed characterization of apomyoglobin (apoMb) soluble aggregates populated at room temperature upon dissolving the purified protein in buffer at pH 6.0, both in the presence and absence of high concentrations of urea. Unlike the beta-sheet self-associated aggregates previously reported for this protein at high temperatures, the soluble aggregates detected here have either alpha-helical or random coil secondary structure, depending on solvent and solution conditions. Hydrodynamic diameters range from 80 to 130 nm, with semiflexible chain-like morphology. The combined use of low pH and high urea concentration leads to structural unfolding and complete elimination of the large aggregates. Even upon starting from this virtually monomeric unfolded state, however, protein refolding leads to the formation of severely self-associated species with native-like secondary structure. Under these conditions, kinetic apoMb refolding proceeds via two parallel routes: one leading to native monomer, and the other leading to a misfolded and heavily self-associated state bearing native-like secondary structure.

Three-state protein folding: experimental determination of free-energy profile

When considering protein folding with a transient intermediate, a difficulty arises as to determination of the rates of separate transitions. Here we overcome this problem, using the kinetic studies of the unfolding/refolding reactions of the three-state protein apomyoglobin as a model. Amplitudes of the protein refolding kinetic burst phase corresponding to the transition from the unfolded (U) to intermediate (I) state, that occurs prior to the native state (N) formation, allow us to estimate relative populations of the rapidly converting states at various final urea concentrations. On the basis of these proportions, a complicated experimental chevron plot has been deconvolved into the urea-dependent rates of the I<-->N and U<-->N transitions to give the dependence of free energies of the main transition state and of all three (N, I, and U) stable states on urea concentration.

Sequence determinants of a protein folding pathway

Local hydrophobic collapse of the polypeptide chain and transient long-range interactions in unfolded states of apomyoglobin appear to occur in regions of the amino acid sequence which, upon folding, bury an above-average area of hydrophobic surface. To explore the role of these interactions in protein folding, we prepared and characterized apomyoglobins with compensating point mutations designed to change the average buried surface area in local regions of the sequence, while conserving as much as possible the constitution of the hydrophobic core. The behavior of the mutants in quench-flow experiments to determine the folding pathway was exactly as predicted by the changes in the buried surface area parameter calculated from the amino acid sequence. In addition, spin label experiments with acid-unfolded mutant apomyoglobin showed that the transient long-range contacts that occur in the wild-type protein are abolished in the mutant, while new contacts are observed between areas that now have above-average buried surface area. We conclude that specific groupings of amino acid side-chains, which can be predicted from the sequence, are responsible for early hydrophobic interactions in the first phase of folding in apomyoglobin, and that these early interactions determine the subsequent course of the folding process.

Enhanced picture of protein-folding intermediates using organic solvents in H/D exchange and quench-flow experiments

Hydrogen/deuterium exchange followed by trapping of the labeled species in the aprotic solvent DMSO has been used to elucidate structure in both the burst-phase molten globule-folding intermediate of apomyoglobin and in an equilibrium intermediate that models the kinetic intermediate. Precise estimates can be made of exchange times in an interrupted exchange-out experiment at pH 4 followed by analysis in DMSO solution, giving extensive sequence-specific information about the structure of the equilibrium intermediate. In addition, the use of DMSO as a solvent for NMR measurements after quench-flow pH-pulse labeling experiments gives a greatly increased data set for the elucidation of the kinetic folding pathway. Interestingly, differences are observed in some regions of apomyoglobin between the equilibrium and kinetic intermediates. These differences are quantitative rather than qualitative; that is, the overall patterns of labeling and secondary structure formation remain similar between the two species. However, local differences are observed, which probably reflect the difference in the solution conditions for the equilibrium experiment (pH 4) vs. the kinetic experiment (pH 6) and the change in the status of the stabilizing hydrogen bond between the side chains of His-24 and His-119.

How Cytochrome c Folds, and Why: Submolecular Foldon Units and their Stepwise Sequential Stabilization

Native state hydrogen exchange experiments have shown that the cytochrome c (Cyt c) protein consists of five cooperative folding-unfolding units, called foldons. These are named, in the order of increasing unfolding free energy, the nested-Yellow, Red, Yellow, Green, and Blue foldons. Previous results suggest that these units unfold in a stepwise sequential way so that each higher energy partially unfolded form includes all of the previously unfolded lower free energy units. If this is so, then selectively destabilizing any given foldon should equally destabilize each subsequent unfolding step above it in the unfolding ladder but leave the lower ones before it unaffected. To perform this test, we introduced the mutation Glu62Gly, which deletes a salt link in the Yellow unit and destabilizes the protein by 0.8 kcal/mol. Native state hydrogen exchange and other experiments show that the stability of the Yellow unit and the states above it in the free energy ladder are destabilized by about the same amount while the lower lying states are unaffected. These results help to confirm the sequential stepwise nature of the Cyt c unfolding pathway and therefore a similar refolding pathway. The steps in the pathway are dictated by the concerted folding-unfolding property of the individual unit foldons; the order of steps is determined by the sequential stabilization of progressively added foldons in the native context. Much related information for Cyt c strongly conforms with this mechanism. Its generality is supported by available information for other proteins.

Unfolding of globular proteins: monte carlo dynamics of a realistic reduced model

Reduced lattice models of proteins and Monte Carlo dynamics were used to simulate the initial stages of the unfolding of several proteins of various structural types, and the results were compared to experiment. The models semiquantitatively reproduce the approximate order of events of unfolding as well as subtle mutation effects and effects resulting from differences in sequences of similar folds. The short-time mobility of particular residues, observed in simulations, correlates with the crystallographic temperature factor. The main factor controlling unfolding is the native state topology, with sequence playing a less important role. The correlation with various experiments, especially for sequence-specific effects, strongly suggests that properly designed reduced models of proteins can be used for qualitative studies (or prediction) of protein unfolding pathways.

Native structural propensity in cellular retinoic acid-binding protein I 64-88: the role of locally encoded structure in the folding of a beta-barrel protein

A central question in protein folding is the relative importance of locally encoded structure and cooperative interactions among residues distant in sequence. We have been exploring this question in a predominantly beta-sheet protein, since beta-structure formation clearly relies on both local and global sequence information. We present evidence that a 24-residue peptide corresponding to two linked hairpins of cellular retinoic acid-binding protein I (CRABP I) adopts significant native structure in aqueous solution. Prior work from our laboratory showed that the two turns contained in this fragment (turns III and IV) had the highest tendency of any of the eight turns in this anti-parallel beta-barrel to fold into native turns. In addition, the primary sequence of these two turns is well conserved throughout the structural family to which CRABP I belongs, and residues in the turns and their associated hairpins participate in a network of conserved long-range interactions. We propose that the strong local-sequence biases within the chain segment comprising turns III and IV favor longer-range interactions that are crucial to the folding and native-state stability of CRABP I, and may play a similar role in related intracellular lipid-binding proteins (iLBPs).

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.

Similarity of force-induced unfolding of apomyoglobin to its chemical-induced unfolding: an atomistic molecular dynamics simulation approach

We have compared force-induced unfolding with traditional unfolding methods using apomyoglobin as a model protein. Using molecular dynamics simulation, we have investigated the structural stability as a function of the degree of mechanical perturbation. Both anisotropic perturbation by stretching two terminal atoms and isotropic perturbation by increasing the radius of gyration of the protein show the same key event of force-induced unfolding. Our primary results show that the native structure of apomyoglobin becomes destabilized against the mechanical perturbation as soon as the interhelical packing between the G and H helices is broken, suggesting that our simulation results share a common feature with the experimental observation that the interhelical contact is more important for the folding of apomyoglobin than the stability of individual helices. This finding is further confirmed by simulating both helix destabilizing and interhelical packing destabilizing mutants.

Chain length dependence of apomyoglobin folding: structural evolution from misfolded sheets to native helices

Very little is known about how protein structure evolves during the polypeptide chain elongation that accompanies cotranslational protein folding. This in vitro model study is aimed at probing how conformational space evolves for purified N-terminal polypeptides of increasing length. These peptides are derived from the sequence of an all-alpha-helical single domain protein, Sperm whale apomyoglobin (apoMb). Even at short chain lengths, ordered structure is found. The nature of this structure is strongly chain length dependent. At relatively short lengths, a predominantly non-native beta-sheet conformation is present, and self-associated amyloid-like species are generated. As chain length increases, alpha-helix progressively takes over, and it replaces the beta-strand. The observed trends correlate with the specific fraction of solvent-accessible nonpolar surface area present at different chain lengths. The C-terminal portion of the chain plays an important role by promoting a large and cooperative overall increase in helical content and by consolidating the monomeric association state of the full-length protein. Thus, a native-like energy landscape develops late during apoMb chain elongation. This effect may provide an important driving force for chain expulsion from the ribosome and promote nearly-posttranslational folding of single domain proteins in the cell. Nature has been able to overcome the above intrinsic misfolding trends by modulating the composition of the intracellular environment. An imbalance or improper functioning by the above modulating factors during translation may play a role in misfolding-driven intracellular disorders.

Analysis of the differences in the folding kinetics of structurally homologous proteins based on predictions of the gross features of residue contacts

It is a general notion that proteins with very similar three-dimensional structures would show very similar folding kinetics. However, recent studies reveal that the folding kinetic properties of some proteins contradict this thought (i.e., the members in a same protein family fold through different pathways). For example, it has been reported that some beta-proteins in the intracellular lipid-binding protein family fold through quite different pathways (Burns et al., Proteins 1998;33:107-118). Similar differences in folding kinetics are also observed in the members of the globin family (Nishimura et al., Nat Struct Biol 2000;7:679-686). In our study, we examine the possibility of predicting qualitative differences in folding kinetics of the intracellular lipid-binding proteins and two globin proteins (i.e., myoglobin and leghemoglobin). The problem is tackled by means of a contact map based on the average distance statistics between residues, the Average Distance Map (ADM), as constructed from sequence. The ADMs for the three proteins show overall similarity, but some local differences among maps are also observed. Our results demonstrate that some properties of the protein folding kinetics are consistent with local differences in the ADMs. We also discuss the general possibility of predicting folding kinetics from sequence information.

Chain length is the main determinant of the folding rate for proteins with three-state folding kinetics

We demonstrate that chain length is the main determinant of the folding rate for proteins with the three-state folding kinetics. The logarithm of their folding rate in water (k(f)) strongly anticorrelates with their chain length L (the correlation coefficient being -0.80). At the same time, the chain length has no correlation with the folding rate for two-state folding proteins (the correlation coefficient is -0.07). Another significant difference of these two groups of proteins is a strong anticorrelation between the folding rate and Baker's "relative contact order" for the two-state folders and the complete absence of such correlation for the three-state folders.

Structural characterization of non-native states of sperm whale myoglobin in aqueous ethanol or 2,2,2-trifluoroethanol media

The effects of aqueous ethanol or 2,2,2-trifluoroethanol media on the structure of sperm whale myoglobin have been investigated by absorption, CD, and NMR spectra. The structural properties of myoglobin such as heme environments, helix contents, protein folding, and interactions between heme and the protein moiety have been sharply manifested in these spectra. The characterization demonstrated that alcohol-induced conformational change of myoglobin depends on the nature of alcohol and its concentration. It was shown for the first time that, upon the alcohol-induced denaturation of myoglobin, heme is released from partially denatured protein of which helix contents is altered by only about 20% relative to that of native state. Myoglobin has shown to unfold and refold reversibly by controlling the alcohol concentration. Novel methods for the preparation of apomyoglobin and in situ reconstitution of apomyoglobin with heme, based on the alcohol-induced denaturation of the protein, were presented.

A circularly permuted myoglobin possesses a folded structure and ligand binding similar to those of the wild-type protein but with a reduced thermodynamic stability

A circular permutein of sperm whale myoglobin in which the G helix is C-terminal, the H helix is N-terminal, and 16 amino acids link the H helix to the A helix has been expressed in Escherichia coli. The permutein sequence begins with Gly121 (using the numbering scheme for the wild-type protein) and terminates with Pro120. The ligand binding function of the permutein was assayed using stopped-flow methods and shown to be essentially identical to that of the wild-type protein. In addition, one- and two-dimensional NMR studies of the cyanomet isoform of the permutein show a nativelike structure with a heme binding pocket very similar to that of the wild-type myoglobin. Although the structure and function of the permutein resemble those of the wild-type myoglobin, the permutein is less stable to chemical denaturation by 5.2 kcal/mol.