Studies on the degradation pathway of iron-sulfur centers during unfolding of a hyperstable ferredoxin: cluster dissociation,iron release and protein stability

The ferredoxin from the thermoacidophile Acidianus ambivalens is a representative of the archaeal family of di-cluster [3Fe-4S][4Fe-4S] ferredoxins. Previous studies have shown that these ferredoxins are intrinsically very stable and led to the suggestion that upon protein unfolding the iron-sulfur clusters degraded via linear three-iron sulfur center species, with 610 and 520 nm absorption bands, resembling those observed in purple aconitase. In this work, a kinetic and spectroscopic investigation on the alkaline chemical denaturation of the protein was performed in an attempt to elucidate the degradation pathway of the iron-sulfur centers in respect to protein unfolding events. For this purpose we investigated cluster dissociation, iron release and protein unfolding by complementary biophysical techniques. We found that shortly after initial protein unfolding, iron release proceeds monophasically at a rate comparable to that of cluster degradation, and that no typical EPR features of linear three-iron sulfur centers are observed. Further, it was observed that EDTA prevents formation of the transient bands and that sulfide significantly enhances its intensity and lifetime, even after protein unfolding. Altogether, our data suggest that iron sulfides, which are formed from the release of iron and sulfide resulting from cluster degradation during protein unfolding in alkaline conditions, are in fact responsible for the observed intermediate spectral species, thus disproving the hypothesis suggesting the presence of a linear three-iron center intermediate. Kinetic studies monitored by visible, fluorescence and UV second-derivative spectroscopies have elicited that upon initial perturbation of the tertiary structure the iron-sulfur centers start decomposing and that the presence of EDTA accelerates the process. Also, the presence of EDTA lowers the observed melting temperature in thermal ramp experiments and the midpoint denaturant concentration in equilibrium chemical unfolding experiments, further suggesting that the clusters also play a structural role in the maintenance of the conformation of the folded state.     

Molecular mechanism underlying the thermal stability and pH-induced unfolding of CHABII

The 37-residue alpha/beta protein CHABII was previously demonstrated to undergo a gradual pH-induced unfolding. It has been shown that even at pH 4.0 CHABII still retained a highly native-like secondary structure and tertiary topology although its tight side-chain packing was severely disrupted, typical of the molten globule state. Here, we have expressed and refolded the recombinant proteins of CHABII and its mutant [Phe21]-CHABII, and subsequently conducted extensive CD and NMR characterizations. The results indicated: (1) replacement of His21 by Phe in [Phe21]-CHABII eliminated the pH-induced unfolding from pH 6.5 to 4.0, indicating that His21 was responsible for the observed pH-induced unfolding of CHABII. Further examinations revealed that although the pH-induced unfolding of CHABII was also triggered by the protonation of the His residue as previously uncovered for apomyoglobin, their molecular mechanisms are different. (2) Monitoring the pH-induced unfolding by 1H-15N HSQC spectroscopy allowed us to visualize the gradual development of the CHABII molten globule. At pH 4.0, the HSQC spectrum of CHABII was poorly dispersed with dispersions of approximately 1 ppm over proton dimension and 10 ppm over 15N dimension, characteristic of severely or even "completely unfolded" proteins. One the other hand, unambiguous assignments of the NOESY spectra of CHABII led to the identification of the persistent medium and long-range NOEs at pH 4.0, which define a highly native-like secondary structure and tertiary packing. This implies that the degree of the native-like topology might be underestimated in the previous characterization of partially folded and even completely unfolded proteins. (3) Replacement of His21 by Phe with higher side-chain hydrophobicity only caused a minor structural rearrangement but considerably enhanced the packing interaction of the hydrophobic core, as evident from a dramatic increase in NOE contacts in [Phe21]-CHABII. The enhancement led to an increase of the thermal stability of [Phe21]-CHABII by approximately 17 deg. C.     

Thermal effects in stretching of Go-like models of titin and secondary structures

The effect of temperature on mechanical unfolding of proteins is studied using a Go-like model with a realistic contact map and Lennard-Jones contact interactions. The behavior of the I27 domain of titin and its serial repeats is contrasted to that of simple secondary structures. In all cases, thermal fluctuations accelerate the unraveling process, decreasing the unfolding force nearly linearly at low temperatures. However, differences in bonding geometry lead to different sensitivity to temperature and different changes in the unfolding pattern. Due to its special native-state geometry, titin is much more thermally and elastically stable than the secondary structures. At low temperatures, serial repeats of titin show a parallel unfolding of all domains to an intermediate state, followed by serial unfolding of the domains. At high temperatures, all domains unfold simultaneously, and the unfolding distance decreases monotonically with the contact order, that is, the sequence distance between the amino acids that form the native contact.     

Hierarchy of structure loss in MD simulations of src SH3 domain unfolding.

To complement experimental studies of the src SH3 domain folding, we studied 30 independent, high-temperature, molecular dynamics simulations of src SH3 domain unfolding. These trajectories were observed to differ widely from each other. Thus, rather than analyzing individual trajectories, we sought to identify the recurrent features of the high-temperature unfolding process. The conformations from all simulations were combined and then divided into groups based on the number of native contacts. Average occupancies of each side-chain hydrophobic contact and hydrogen bond in the protein were then determined. In the symmetric funnel limit, the occupancies of all contacts should decrease in concert with the loss in total number of native contacts. If there is a lack of symmetry or hierarchy to the unfolding process, the occupancies of some contacts should decrease more slowly, and others more rapidly. Despite the heterogeneity of the individual trajectories, the ensemble averaging revealed an order to the unfolding process: contacts between the N and C-terminal strands are the first to disappear, whereas contacts within the distal beta-hairpin and a hydrogen-bonding network involving the distal loop beta-turn and the diverging turn persist well after the majority of the native contacts are lost. This hierarchy of events resembles but is somewhat less pronounced than that observed in our experimental studies of the folding of src SH3 domain.     

Two-dimensional infrared correlation spectroscopy study of sequential events in the heat-induced unfolding and aggregation process of myoglobin

Unfolding and aggregation are basic problems in protein science with serious biotechnological and medical implications. Probing the sequential events occurring during the unfolding and aggregation process and the relationship between unfolding and aggregation is of particular interest. In this study, two-dimensional infrared (2D IR) correlation spectroscopy was used to study the sequential events and starting temperature dependence of Myoglobin (Mb) thermal transitions. Though a two-state model could be obtained from traditional 1D IR spectra, subtle noncooperative conformational changes were observed at low temperatures. Formation of aggregation was observed at a temperature (50-58 degrees C) that protein was dominated by native structures and accompanied with unfolding of native helical structures when a traditional thermal denaturation condition was used. The time course NMR study of Mb incubated at 55 degrees C for 45 h confirmed that an irreversible aggregation process existed. Aggregation was also observed before fully unfolding of the Mb native structure when a relative high starting temperature was used. These findings demonstrated that 2D IR correlation spectroscopy is a powerful tool to study protein aggregation and the protein aggregation process observed depends on the different environmental conditions used.     

Structural heterogeneity of the various forms of apomyoglobin: implications for protein folding.

Temperature-induced denaturation transitions of different structural forms of apomyoglobin were studied monitoring intrinsic tryptophan fluorescence. It was found that the tryptophans are effectively screened from solvent both in native and acid forms throughout most of the temperature range tested. Thus, the tryptophans' surrounding do not show a considerable change in structure where major protein conformational transitions have been found in apomyoglobin using other techniques. At high temperatures and under strong destabilizing conditions, the tryptophans' fluorescence parameters show sigmoidal thermal denaturation. These results, combined with previous studies, show that the structure of this protein is heterogeneous, including native-like (tightly packed) and molten globule-like substructures that exhibit conformation (denaturation) transitions under different conditions of pH and temperature (and denaturants). The results suggest that the folding of this protein proceeds via two "nucleation" events whereby native-like contacts are formed. One of these events, which involves AGH "core" formation, appears to occur very early in the folding process, even before significant hydrophobic collapse in the rest of the protein molecule. From the current studies and other results, a rather detailed picture of the folding of myoglobin is presented, on the level of specific structures and their thermodynamical properties as well as formation kinetics.     

Molecular dynamics simulations of a protein model in uniform and elongational flows

We present a molecular dynamics study of the conformational deformation of a minimalist beta-barrel protein model in two different types of hydrodynamic flows: uniform and elongational. We investigate the characteristics of protein stretching, paying special attention to the unfolding intermediate states and their relationship to the protein folding/unfolding problem. In the uniform flow simulations, one end of the modeled protein was tethered to a fixed point in space and the forced unfolding process was observed. The unfolding takes place via a few stages involving one or two intermediate states, depending on which end is tethered. The calculated force-extension curves show plateau regimes and hysteresis as the protein is stretched and refolded, in qualitative agreement with the experimental measurements. The physical behavior observed in our numerical simulations of the forced unfolding in an elongational flow is very different from that in uniform flow. The protein unfolds abruptly from the globular state to a fully stretched state without going through any observable intermediate states. From these observation, we stress that protein unfolding pathways under the influence of an external force are highly dependent on the mechanism of the exerted force.     

Surprisingly high correlation between early and late stages in non-two-state protein folding

We collected quantitative kinetic data on early and late stages of folding in non-two-state proteins from the literature, and studied the relationship between the kinetics of the two stages. There was a surprisingly high correlation between the rate constants of these stages. The correlation coefficient of the logarithmic rate constants was as high as 0.97, which could not be caused by chance. We also studied relationships of the logarithmic rate constants of the two stages with native three-dimensional structures represented by the residue-residue contact map. There were again surprisingly high correlations between the logarithmic rate constants and the number of non-local contact clusters obtained from the contact maps. Because the number of non-local contact clusters represents overall arrangement of substructures in a native protein, the results strongly suggested the importance of the arrangement of the substructures for the kinetics of both early and late stages of protein folding.     

Cell-cell contact promotes specific activity of pp60c-src protein kinase in human retinoblastoma cells.

Tyrosine-specific protein kinase activity of pp60c-src was examined in human Y79 retinoblastoma cells cultured in monolayers after clusters in suspension culture had been dissociated. The activity increased five- to six-fold between Days 1 and 7 in the monolayer cultures, with a concomitant increase in numbers of cellular contacts per cell. There was no effect of conditioned medium from high-density cultures in suspension on the activity of cultures with a low degree of contacts. The level of c-src protein in cell lysates was nearly constant irrespective of the degree of cellular contacts. These results suggest that the specific activity of pp60c-src is regulated by cell-cell contact.     

Topological frustration and the folding of interleukin-1 beta

The cytokine, interleukin-1beta (IL-1beta), adopts a beta-trefoil fold. It is known to be much slower folding than similarly sized proteins, despite having a low contact order. Proteins are sufficiently well designed that their folding is not dominated by local energetic traps. Therefore, protein models that encode only the folded structure and are energetically unfrustrated (Gō-type), can capture the essentials of the folding routes. We investigate the folding thermodynamics of IL-1beta using such a model and molecular dynamics (MD) simulations. We develop an enhanced sampling technique (a modified multicanonical method) to overcome the sampling problem caused by the slow folding. We find that IL-1beta has a broad and high free energy barrier. In addition, the protein fold causes intermediate unfolding and refolding of some native contacts within the protein along the folding trajectory. This "backtracking" occurs around the barrier region. Complex folds like the beta-trefoil fold and functional loops like the beta-bulge of IL-1beta can make some of the configuration space unavailable to the protein and cause topological frustration.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Prediction of optimal folding routes of proteins that satisfy the principle of lowest entropy loss: dynamic contact maps and optimal control

An optimization model is introduced in which proteins try to evade high energy regions of the folding landscape, and prefer low entropy loss routes during folding. We make use of the framework of optimal control whose convenient solution provides practical and useful insight into the sequence of events during folding. We assume that the native state is available. As the protein folds, it makes different set of contacts at different folding steps. The dynamic contact map is constructed from these contacts. The topology of the dynamic contact map changes during the course of folding and this information is utilized in the dynamic optimization model. The solution is obtained using the optimal control theory. We show that the optimal solution can be cast into the form of a Gaussian Network that governs the optimal folding dynamics. Simulation results on three examples (CI2, Sso7d and Villin) show that folding starts by the formation of local clusters. Non-local clusters generally require the formation of several local clusters. Non-local clusters form cooperatively and not sequentially. We also observe that the optimal controller prefers "zipping" or small loop closure steps during folding. The folding routes predicted by the proposed method bear strong resemblance to the results in the literature.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.     

Folding rates and low-entropy-loss routes of two-state proteins

We develop a simple model for computing the rates and routes of folding of two-state proteins from the contact maps of their native structures. The model is based on the graph-theoretical concept of effective contact order (ECO). The model predicts that proteins fold by "zipping up" in a sequence of small-loop-closure events, depending on the native chain fold. Using a simple equation, with a few physical rate parameters, we obtain a good correlation with the folding rates of 24 two-state folding proteins. The model rationalizes data from Phi-value analysis that have been interpreted in terms of delocalized or polarized transition states. This model indicates how much of protein folding may take place in parallel, not along a single reaction coordinate or with a single transition state.     

Burst-phase expansion of native protein prior to global unfolding in SDS

Although numerous studies have been directed at understanding early folding events through the characterization of folding intermediates, there are few reports on the very late folding events, i.e. on the events taking place on the native side of the folding barrier and on alternative conformations of the folded state. To shed further light on these issues, we have characterized by protein engineering the structure of an expanded but native-like intermediate that accumulates transiently in the unfolding reaction of the small protein S6 in the presence of SDS. The results show that the SDS micelles attack the native protein in the dead-time of the denaturation experiment, causing an expansion of the hydrophobic core prior to the major unfolding transition. We distinguish two forms of the unfolding intermediate that are correlated with the micellar structure. With spherical micelles, the expansion is seen mainly as a weakening of the interactions which anchor the two alpha-helices to the core of the S6 structure. With cylindrical micelles, prevalent at higher SDS concentrations, the expansion is more global and produces a species which closely resembles the transition-state structure for unfolding in GdmCl. Despite the highly weakened core, the micelle-associated intermediate displays cooperative unfolding, indicating a significant structural plasticity of the species on the native side of the folding barrier in the presence of SDS.     

Thermal folding and mechanical unfolding pathways of protein secondary structures

Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force versus displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an alpha-helix and a beta-hairpin structure. In a tandem of two alpha-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous versus sequential unraveling of repeat units is presented.     

The importance of explicit chain representation in protein folding models: an examination of Ising-like models

A class of models that represents a protein chain as a sequence of "folded" and "unfolded" residues has recently been used to correlate rates and mechanisms of protein folding with the protein native structure. In order to better understand the conditions under which these "Ising-like" models apply, we compare results from this model to those obtained from an off-lattice model which uses the same potential function. We find that Ising-like models by construction impose folding via a highly sequential nucleation-condensation mechanism, which in turn leads to more rugged energy landscapes, fewer "pathways" to the native state, and in the specific case examined here, the cold shock protein A from Escherichia coli, a qualitative difference in the most likely order of events in folding.     

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high‐energy conformation susceptible to proteolysis (cleavable form) using native‐state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m‐value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F‐G and Met‐20 loops on one face of the central β‐sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near‐native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F‐G and Met‐20 loops, which is coupled with compaction of the rest of the protein.     

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.