Equilibrium unfolding of dimeric desulfoferrodoxin involves a monomeric intermediate: iron cofactors dissociate after polypeptide unfolding

A decision tree model for the prediction of homodimer folding mechanism

The formation of protein homodimer complexes for molecular catalysis and regulation is fascinating. The homodimer formation through 2S (2 state), 3SMI (3 state with monomer intermediate) and 3SDI (3 state with dimer intermediate) folding mechanism is known for 47 homodimer structures. Our dataset of forty-seven homodimers consists of twenty-eight 2S, twelve 3SMI and seven 3SDI. The dataset is characterized using monomer length, interface area and interface/total (I/T) residue ratio. It is found that 2S are often small in size with large I/T ratio and 3SDI are frequently large in size with small I/T ratio. Nonetheless, 3SMI have a mixture of these features. Hence, we used these parameters to develop a decision tree model. The decision tree model produced positive predictive values (PPV) of 72% for 2S, 58% for 3SMI and 57% for 3SDI in cross validation. Thus, the method finds application in assigning homodimers with folding mechanism.

Sequential dissociation of subunits from bovine heart cytochrome C oxidase by urea

The quaternary stability of purified, detergent-solubilized, cytochrome c oxidase (CcO) was probed using two chemical denaturants, urea and guanidinium chloride (GdmCl). Each chaotrope induces dissociation of five subunits in a concentration-dependent manner. These five subunits are not scattered over the surface of CcO but are clustered together in close contact at the dimer interface. Increasing the concentration of urea selectively dissociates subunits from CcO in the following order: VIa and VIb, followed by III and VIIa, and finally Vb. After incubation in urea for 10 min at room temperature, the sigmoidal dissociation transitions were centered at 3.7, 4.6, and 7.0 M urea, respectively. The secondary structure of CcO was only minimally perturbed, indicating that urea causes disruption of subunit interactions without urea-induced conformational changes. Incubation of CcO in urea for 120 min produced similar results but shifted the sigmoidal dissociation curves to lower urea concentrations. Incubation of CcO with increasing concentrations of GdmCl produces an analogous effect; however, the GdmCl-induced dissociation of subunits occurs at lower concentrations and with a narrower concentration range. Thermodynamic parameters for each subunit dissociation were evaluated from the sigmoidal dissociation data by assuming a single transition from bound to dissociated subunit. The free energy change accompanying urea-induced dissociation of each subunit ranged from 18.0 to 29.7 kJ/mol, which corresponds to 0.32-0.59 kJ/mol per 100 A(2) of newly exposed solvent-accessible surface area. These values are 30-50-fold smaller than previously reported for the unfolding of soluble or membrane proteins.

Reversible unfolding of dimeric phosphofructokinase-2 from Escherichia coli reveals a dominant role of inter-subunit contacts for stability

Escherichia coli phosphofructokinase-2 (Pfk-2) is a homodimer whose subunits consist of a large domain and an additional beta-sheet that provides the interfacial contacts between the subunits, creating a beta-barrel flattened-like structure with the adjacent subunit's beta-sheet. To determine how the structural organization of Pfk-2 determines its stability, the reversible unfolding of the enzyme was characterized under equilibrium conditions by enzymatic activity, circular dichroism, fluorescence and hydrodynamic measurements. Pfk-2 undergoes a cooperative unfolding/dissociation process with the accumulation of an expanded and unstructured monomeric intermediate with a marginal stability and a large solvent accessibility with respect to the native dimer.

Structural features for homodimer folding mechanism

The homodimers have essential role in catalysis and regulation. The homodimer folding mechanism through 2-state without stable intermediate (2S), 3-state with monomer intermediate (3SMI) and 3-state with dimer intermediate (3SDI) is fascinating. 23MI and 3SDI constitute 3-state (3S). Hence, it is important to differentiate 2S, 3SMI and 3SDI homodimers using structural features. We used the dataset of Li et al. [L. Li, K. Gunasekaran, J.G. Gan, C. Zhanhua, P. Shapshak, M.K. Sakharkar, P. Kangueane, Structural features differentiate the mechanisms between 2S and 3S folding of homodimers, Bioinformation 1 (2005) 42-49] consisting of twenty-five 2S, ten 3SMI and six 3SDI homodimer structures for the study. Interface to total (I/T) residues ratio is large for 2S than 3SMI and 3SDI. Interface to total residues ratio is similar for 3SMI (mean monomer length (ML)=208) and 3SDI (mean monomer length (ML)=404) despite difference in mean monomer size. Interface residues correlate with monomer size in 2S (Pearson's correlation coefficient (r); r(2)=0.41) and 3SMI (r(2)=0.52). This is not true for 3SDI with interface residues and monomer length (r(2)=0.17). Interface area (B/2) does not correlate with interface residues (r(2)<0.001) and monomer size (r(2)=0.023) in 2S. This is despite a relationship with interface residues and monomer size (r(2)=0.41) in 2S. However, this is not true for 3SMI (r(2)=0.61 with interface residues and r(2)=0.25 with monomer size). In 3SDI, a different relationship is seen (r(2)=0.28 with interface residues and r(2)=0.09 with size). The mean hydrophobicity factor (H(f)) is 3-fold less in 3S than 2S. H(f) does not correlate with interface area in 2S (r(2)=0.03) and 3SDI (r(2)=0.0). However, a weak causal relation is seen in 3SMI (r(2)=0.23). Hydrophilic amino acid residues (E, R, K, S and Q) are prominent in 2S than 3S. Charged negative amino acid residues (D, E) are more than positive amino acid residues (R, K, H) in 2S and charged positive amino acid residues (R, K, H) are more than negative amino acid residues (D, E) in 3S. These features help to distinguish 2S, 3SMI and 3SDI providing insights to homodimer folding and binding.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Kinetic Folding and Assembly Mechanisms Differ for Two Homologous Heptamers

Here we investigate the time-resolved folding and assembly mechanism of the heptameric co-chaperonin protein 10 (cpn10) in vitro. The structure of cpn10 is conserved throughout nature: seven beta-barrel subunits are non-covalently assembled through beta-strand pairings in an overall doughnut-like shape. Kinetic folding/assembly experiments of chemically denatured cpn10 from Homo sapiens (hmcpn10) and Aquifex aeolicus (Aacpn10) were monitored by far-UV circular dichroism and fluorescence. We find the processes to be complex, involving several kinetic steps, and to differ between the mesophilic and hyper-thermophilic proteins. The hmcpn10 molecules partition into two parallel pathways, one involving polypeptide folding before protein-protein assembly and another in which inter-protein interactions take place prior to folding. In contrast, the Aacpn10 molecules follow a single sequential path that includes initial monomer misfolding, relaxation to productive intermediates and, subsequently, final folding and heptamer assembly. An A. aeolicus variant lacking the unique C-terminal extension of Aacpn10 displays the same kinetic mechanism as Aacpn10, signifying that the tail is not responsible for the rapid misfolding step. This study demonstrates that molecular details can overrule similarity of native-state topology in defining apparent protein-biophysical properties.

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.

Structural features differentiate the mechanisms between 2S (2 state) and 3S (3 state) folding homodimers

The formation of homodimer complexes for interface stability, catalysis and regulation is intriguing. The mechanisms of homodimer complexations are even more interesting. Some homodimers form without intermediates (two-state (2S)) and others through the formation of stable intermediates (three-state (3S)). Here, we analyze 41 homodimer (25 2S and 16 3S) structures determined by X-ray crystallography to estimate structural differences between them. The analysis suggests that a combination of structural properties such as monomer length, subunit interface area, ratio of interface to interior hydrophobicity can predominately distinguish 2S and 3S homodimers. These findings are useful in the prediction of homodimer folding and binding mechanisms using structural data.

The importance of being dimeric

Why are there so many dimeric proteins and enzymes? While for heterodimers a functional explanation seems quite reasonable, the case of homodimers is more puzzling. The number of homodimers found in all living organisms is rapidly increasing. A thorough inspection of the structural data from the available literature and stability (measured from denaturation-renaturation experiments) allows one to suggest that homodimers can be divided into three main types according to their mass and the presence of a (relatively) stable monomeric intermediate in the folding-unfolding pathway. Among other explanations, we propose that an essential advantage for a protein being dimeric may be the proper and rapid assembly in the cellular milieu.

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.

An engineered two-iron superoxide reductase lacking the [Fe(SCys)4] site retains its catalytic properties in vitro and in vivo

Superoxide reductases (SORs) contain a characteristic square-pyramidal [Fe(NHis)(4)(SCys)] active site that catalyzes reduction of superoxide to hydrogen peroxide in several anaerobic bacteria and archaea. Some SORs, referred to as two-iron SORs (2Fe-SORs), also contain a lower-potential [Fe(SCys)(4)] site that is presumed to have an electron transfer function. However, the intra- and inter-subunit distances between [Fe(SCys)(4)] and [Fe(NHis)(4)(SCys)] iron centers within the 2Fe-SOR homodimer seem too long for efficient electron transfer between these sites. The possible role of the [Fe(SCys)(4)] site in 2Fe-SORs was addressed in this work by examination of an engineered Desulfovibrio vulgaris 2Fe-SOR variant, C13S, in which one ligand residue of the [Fe(SCys)(4)] site, cysteine 13, was changed to serine. This single amino acid residue change destroyed the native [Fe(SCys)(4)] site with complete loss of its iron, but left the [Fe(NHis)(4)(SCys)] site and the protein homodimer intact. The spectroscopic, redox and superoxide reactivity properties of the [Fe(NHis)(4)(SCys)] site in the C13S variant were nearly indistinguishable from those of the wild-type 2Fe-SOR. Aerobic growth complementation of a superoxide dismutase (SOD)-deficient Escherichia coli strain showed that the presence of the [Fe(NHis)(4)(SCys)] site in C13S 2Fe-SOR was apparently sufficient to catalyze reduction of the intracellular superoxide to nonlethal levels. As is the case for the wild-type protein, C13S 2Fe-SOR did not show any detectable SOD activity, i.e., destruction of the [Fe(SCys)(4)] site did not unmask latent SOD activity of the [Fe(NHis)(4)(SCys)] site. Possible alternative roles for the [Fe(SCys)(4)] site in 2Fe-SORs are considered.

Thermodynamic characterization of yeast triosephosphate isomerase refolding: insights into the interplay between function and stability as reasons for the oligomeric nature of the enzyme

The reasons underlying the oligomeric nature of some proteins such as triosephosphate isomerase (TIM) are unclear. It has been proposed that this enzyme is an oligomer, mainly because of its stability rather than for functional reasons. To address this issue, the reversible denaturation and renaturation of the homodimeric TIM from baker's yeast ( Saccharomyces cerevisiae ) induced by guanidinium chloride and urea have been characterized by spectroscopic, functional and hydrodynamic techniques. The unfolding and refolding of this enzyme are not coincident after 'conventional' equilibrium times. Unfolding experiments did not reach equilibrium, owing to a very slow dissociation and/or unfolding process. By contrast, equilibrium was reached in the refolding direction. The simplest equilibrium pathway compatible with the obtained data was found to be a three-state process involving an inactive and expanded monomer. The Gibbs energy changes for monomer folding (delta G (0)(fold) = -16.6+/-0.7 kJ x mol(-1)) and monomer association (delta G (0)(assoc) = -70.3+/-1.1 kJ x mol(-1)) were calculated from data obtained in the two denaturants. From an analysis of the present data and data from the literature on the stability of TIM from different species and for other beta/alpha barrels, and model simulations on the effect of stability in the catalytic activity of the enzyme, it is concluded that the low stability of the monomers is neither the only, nor the main, cause for the dimeric nature of TIM. There is interplay between function and stability.

Equilibrium unfolding of dimeric human prostatic acid phosphatase involves an inactive monomeric intermediate

Guanidine hydrochloride (GdnHCl)-induced unfolding of human prostatic acid phosphatase (hPAP), a homodimer of 50 kDa subunit molecular weight, was investigated with activity measurements, size exclusion HPLC, tryptophan fluorescence, 1-anilinonaphtalene-8-sulfonate (ANS) binding and reactivity with 2-(4'-maleimidoanilino)naphthalene-6-sulfonate (MIANS). Equilibrium analysis was performed to shed light on the role of dimerization in the folding and stability of the catalytically active oligomeric protein. Unfolding was reversible, as verified by activity measurements and tryptophan fluorescence. The noncoincidence of the unfolding curves obtained by different techniques suggests the occurrence of a multiphasic process. The reaction of hPAP inactivation is accompanied by dissociation of the dimer into two monomers. The midpoint of this transition is at 0.65 M GdnHCl with 4.24+/-0.12 kcalmol(-1) free energy change. Binding of ANS to the inactive phosphatase monomer, especially remarkable in the region from 0.8 to 1.25M GdnHCl, suggests that the hydrophobic probe indicates exposition of the intersubunit hydrophobic surface and a loosening of the monomer's tertiary structure. Strong fluorescence of thiol group derivatives, the products of their reaction with MIANS, appears in a limited range of GdnHCl concentrations (1.2-1.6M). This shows that in the relaxed structure of the intermediate, the reagent is allowed to penetrate into the hydrophobic environment of the partially hidden thiol groups. The equilibrium unfolding reaction of hPAP, as monitored by tryptophan fluorescence, does not depend on the protein concentration and displays a single transition curve with a midpoint at 1.7 M GdnHCl and value of DeltaG(unf)(H(2)O)=3.38+/-0.08 kcalmol(-1) per monomer, a result implying that this transition is related to the conformational change of the earlier dissociated and already inactive subunit of the protein.

Folding, stability, and secondary structure of a new dimeric cysteine proteinase inhibitor

Clitocypin, a new type of cysteine proteinase inhibitor from the mushroom Clitocybe nebularis, is a 34-kDa homodimer lacking disulphide bonds, reported to have unusual stability properties. Sequence similarity is limited solely to certain proteins from mushrooms. Infrared spectroscopy shows that clitocypin is a high beta-structure protein which was lost at high temperatures. The far UV circular dichroism spectrum is not that of classical beta-structure, but similar to those of a group of small beta-strand proteins, with a peak at 189nm and a trough at 202nm. An aromatic peak at 232nm and infrared bands at 1633 and 1515cm(-1) associated with the peptide backbone and the tyrosine microenvironment, respectively, were used to characterize the thermal unfolding. The reversible transition has a midpoint at 67 degrees C, with DeltaG=34kJ/mol and DeltaH=300kJ/mol, and is, unusually, independent of protein concentration. The kinetics of thermal unfolding and refolding are slow, with activation energies of 167 and 44kJ/mol, respectively. A model for folding and assembly is discussed.

Reversible unfolding of bovine odorant binding protein induced by guanidinium hydrochloride at neutral pH

An analysis of the unfolding and refolding curves at equilibrium of dimeric bovine odorant binding protein (bOBP) has been performed. Unfolding induced by guanidinium chloride (GdnHCl) is completely reversible as far as structure and ligand binding capacity are concerned. The transition curves, as obtained by fluorescence and ellipticity measurements, are very similar and have the same protein concentration-independent midpoint (C1/2 approximately 2.6 M). This result implies a sequential, rather than a concerted, unfolding mechanism, with the involvement of an intermediate. However, since it has not been detected, this intermediate must be present in small amounts or have the same optical properties of either native or denatured protein. The thermodynamic best fit parameters, obtained according to a simple two-state model, are: deltaG degrees un,w = 5.0 +/- 0.6 kcal mol(-1), m = 1.9 +/- 0.2 kcal mol(-1) M(-1) and C1/2 = 2.6 +/- 0.1 M. The presence of the ligand dihydromyrcenol has a stabilising effect against unfolding by GdnHCl, with an extrapolated deltaG degrees un,w of 22.2 +/- 0.9 kcal mol(-1), a cooperative index of 3.2 +/- 0.3 and a midpoint of 4.6 +/- 0.4 M. The refolding curves, recorded after 24 h from dilution of denaturant are not yet at equilibrium: they show an apparently lower midpoint (C1/2 = 2.2 M), but tend to overlap the unfolding curve after several days. In contrast to chromatographic unfolding data, which fail to reveal the presence of folded intermediates, chromatographic refolding data as a function of time clearly show a rapid formation of folded monomers, followed by a slower step leading to folded dimers. Therefore, according to this result, we believe that the preferential unfolding/refolding mechanism is one in which dimer dissociation occurs before unfolding rather than the reverse.

Impact of domain interchange on conformational stability and equilibrium folding of chimeric class micro glutathione transferases

Rat micro class glutathione transferases M1-1 and M2-2 are homodimers that share a 78% sequence identity but display differences in stability. M1-1 is more stable at the secondary and tertiary structural levels, whereas its quaternary structure is less stable. Each subunit in these proteins consists of two structurally distinct domains with intersubunit contacts occurring between domain 1 of one subunit and domain 2 of the other subunit. The chimeric subunit variants M(12), which has domain 1 of M1 and domain 2 of M2, and its complement M(21), were used to investigate the conformational stability of the chimeric homodimers M(12)-(12) and M(21)-(21) to determine the contribution of each domain toward stability. Exchanging entire domains between class micro GSTs is accommodated by the GST fold. Urea-induced equilibrium unfolding data indicate that whereas the class micro equilibrium unfolding mechanism (i.e., N(2) <--> 2I <--> 2U) is not altered, domain exchanges impact significantly on the conformational stability of the native dimers and monomeric folding intermediates. Data for the wild-type and chimeric proteins indicate that the order of stability for the native dimer (N(2)) is M2-2 > M(12)-(12) M1-1 approximately M(21)-(21), and that the order of stability of the monomeric intermediate (I) is M1 > M2 approximately M(12) > M(21). Interactions involving Arg 77, which is topologically conserved in GSTs, appear to play an important role in the stability of both the native dimeric and folding monomeric structures.

Role of cofactors in protein folding

Although cofactors are essential components of many proteins to attain biological activity, the role of cofactors in protein folding is not well understood. Biophysical characterization of four types of cofactor-binding proteins (with copper, flavin moiety, iron-sulfur cluster, and heme cofactors, respectively) provides the following insights. (1) The presence of the cofactor often stabilizes the native protein. (2) The cofactor has the ability to interact specifically with the unfolded polypeptide. (3) The presence of the cofactor is sometimes essential for the polypeptide to fold. (4) Coordination of the cofactor prior to polypeptide folding can dramatically accelerate formation of the functional protein.