Heme orientation affects holo-myoglobin folding and unfolding kinetics

Water stabilizes an alternate turn conformation in horse heart myoglobin

Comparison of myoglobin structures reveals that protein isolated from horse heart consistently adopts an alternate turn conformation in comparison to its homologues. Analysis of hundreds of high-resolution structures discounts crystallization conditions or the surrounding amino acid protein environment as explaining this difference, that is also not captured by the AlphaFold prediction. Rather, a water molecule is identified as stabilizing the conformation in the horse heart structure, which immediately reverts to the whale conformation in molecular dynamics simulations excluding that structural water.

Changes in the structure and digestibility of myoglobin treated with sodium chloride

Myoglobin is a protein not easily broken down by digestive enzymes due to its rigid structure. This study evaluated the structural characteristics of myoglobin under various sodium chloride treatments (0.4-0.8 mol/L for 5-10 h) and the impacts on its digestibility using spectroscopic and molecular dynamics simulation techniques. Myoglobin digestibility was 40% following pepsin digestion and 60% after being sequentially digested by pepsin and trypsin. The α-helix content of myoglobin did not change significantly following sodium chloride treatment but hydrophobic amino acids were exposed and the binding of phenylalanine targeted by some digestive enzymes became more stable, leading to the reduced digestibility.

Heme-bound SiaA from Streptococcus pyogenes: Effects of mutations and oxidation state on protein stability

The protein SiaA (HtsA) is part of a heme uptake pathway in Streptococcus pyogenes. In this report, we present the heme binding of the alanine mutants of the axial histidine (H229A) and methionine (M79A) ligands, as well as a lysine (K61A) and cysteine (C58A) located near the heme propionates (based on homology modeling) and a control mutant (C47A). pH titrations gave pKa values ranging from 9.0 to 9.5, close to the value of 9.7 for WT SiaA. Resonance Raman spectra of the mutants suggested that the ferric heme environment may be distinct from the wild-type; spectra of the ferrous states were similar. The midpoint reduction potential of the K61A mutant was determined by spectroelectrochemical titration to be 61±3mV vs. SHE, similar to the wild-type protein (68±3mV). The addition of guanidine hydrochloride showed two processes for protein denaturation, consistent with heme loss from protein forms differing by the orientation of the heme in the binding pocket (the half-life for the slower process ranged from less than half a day to two days). The ease of protein unfolding was related to the strength of interaction of the residues with the heme. We hypothesize that kinetically facile but only partial unfolding, followed by a very slow approach to the completely unfolded state, may be a fundamental attribute of heme trafficking proteins. Small motions to release/transfer the heme accompanied by resistance to extensive unfolding may preserve the three dimensional form of the protein for further uptake and release.

Effects of heme binding on myoglobin folding: Coarse grained molecular simulations

Many proteins contain cofactors, such as heme, ATP and metal ions. Binding of cofactors is not only essential for their biological functions, but also can reshape the intrinsic energy landscape of protein molecules and modulate the folding and stability. However, the molecular mechanism of cofactor coupled protein folding is not well understood. In this work, we study the cofactor coupled folding of myoglobin, which is a typical cofactor (heme) containing protein, by performing molecular dynamics simulations with a structure-based protein model developed based on the energy landscape theory. We showed that the heme binding increases the stability of the myoglobin. More importantly, the heme binding tends to increase the protein folding cooperativity, and switch the folding process from a “three-state” mechanism to a “two-state” mechanism. We also showed that the folding pathways of the myoglobin can be modulated by the heme binding. By performing comparative simulations, we revealed that the above effects of heme binding are resulted from the heme induced folding of F-helix, which is otherwise unstructured at apo state, and the heme mediated contacting interactions around the heme binding site. The simulation results are consistent with available experimental data, and provide insights into the molecular mechanism of the effects of cofactor binding on the protein folding and stability.

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.

Spectroscopic studies on the effect of high pressure treatment on the soluble protein fraction of porcine longissimus dorsi

The effect of high pressure (HP) treatment (400 MPa, 10 min) of porcine longissimus dorsi was investigated using reflectance spectroscopy and by UV-vis and circular dichroism spectroscopy for the soluble protein fraction. The soluble protein content was expectedly lowered significantly by HP treatment, whereas the solid state fraction of the meat responded to HP by exhibition of characteristic spectral changes in the visible reflectance data with a temporal evolution over the course of 2 days. However, the soluble protein fraction did not exhibit the same altered spectral characteristics in the visible region as seen in the solid state following HP, and there were no indications of altered folding of the proteins that remain in solution.

Role of heme in the unfolding and assembly of myoglobin

The unfolding of wild-type holomyoglobin in the ferric state (metMb) appears to be a simple two-state process, even though hemichrome spectra are often observed and apoMb denaturation involves an intermediate. To resolve these discrepancies, we measured GuHCl-induced, equilibrium unfolding of five sperm whale metMb variants, which were selected to examine the relative importance of apoglobin stability and hemin affinity. Combined analysis of CD, Trp fluorescence, and Soret absorbance titration curves for all the variants requires a six-state mechanism containing native (N), intermediate (I), and unfolded (U) states of apoMb and their hemin-bound counterparts, NH (holoMb), IH, and UH, respectively. The unfolding parameters for the apoMbs were obtained in independent experiments and then fixed in the analysis of the holoprotein data, where only the affinities of the apoglobin states for hemin were allowed to vary. This cofactor binding analysis applies generally to all globins and led to three specific conclusions. (1) The stability of holo-metMb is determined primarily by the high affinity (K(d) approximately 10(-13) M) of native apoMb (N) for hemin. (2) The partially unfolded intermediate with hemin bound (IH) has a hemichrome spectrum indicative of a bis-histidyl axial coordination and is seen clearly when the stability of the N state or its affinity for hemin is reduced. (3) Although the affinity of the intermediate for hemin (K(d) approximately 10(-11) M) is approximately 100-fold lower than that for the native state, free hemin can bind to it and promote the assembly of the holoprotein.

Probing the role of hydration in the unfolding transitions of carbonmonoxy myoglobin and apomyoglobin

We show that the equilibrium unfolding transition of horse carbonmonoxy myoglobin monitored by the stretching vibration of the CO ligand, a local environmental probe, is very sharp and, thus, quite different from those measured by global conformational reporters. In addition, the denatured protein exhibits an A(0)-like CO band. We hypothesize that this sharp transition reports penetration of water into the heme pocket of the protein. Parallel experiments on horse apomyoglobin, wherein an environment-sensitive fluorescent probe, nile red, was used, also reveals a similar putative hydration event. Given the importance of dehydration in protein folding and also the recent debate over the interpretation of probe-dependent unfolding transitions, these results have strong implications on the mechanism of protein folding.

Dynamics of ligand rebinding to unfolded MbCO by guanidine HCl

Femtosecond vibrational spectroscopy was used to probe a functionally important dynamics and residual structure of myoglobin unfolded by 4 M guanidine HCl. The spectra of the dissociated CO indicated that the residual structure of unfolded myoglobin (Mb) forms a few hydrophobic cavities that could accommodate the dissociated ligand. Geminate rebinding (GR) of CO to the unfolded Mb is three-orders-of-magnitude faster and more efficient than the native Mb but similar to a model heme in a viscous solvent, suggesting that the GR of CO to heme is accelerated by the longer retention of the dissociated ligand near the Fe atom by the poorly-structured protein matrix of the unfolded Mb or viscous solvent. The inefficient GR of CO in native Mb, while dissociated CO is trapped in the primary heme pocket located near the active binding site, indicates that the tertiary structure of the pocket in native Mb plays a functionally significant role.

Novel properties of P450s in Streptomyces coelicolor

In the process of investigating the function of the 18 P450s that make up the CYPome in the soil bacterium Streptomyces coelicolor, we have discovered new and novel biochemical features of some of these monooxygenases. This article which is dedicated to Ron Estabrook summarizes these discoveries with emphasis on their novelty. It is concluded from this work that there are many interesting features of P450 enzymes waiting to be discovered. Since there are more than 6500 P450 gene sequences known today, we can expect that many new biochemical aspects of P450s will be discovered as these individual monooxygenases are investigated.

Flow injection analysis with diode array absorbance detection and dynamic surface tension detection for studying denaturation and surface activity of globular proteins

In this article, a multidimensional dynamic surface tension detector (DSTD), in a parallel configuration with a UV-visible diode array absorbance detector, is presented in a novel flow injection analysis (FIA) application to study the effects of chemical denaturants urea, guanidinium hydrochloride (GdmHCl), and guanidinium thyocyanate (GdmSCN) on the surface activity of globular proteins at the liquid-air interface. The DSTD signal is obtained by measuring the changing pressure across the liquid-air interface of 4-mul drops repeatedly forming at the end of a capillary using FIA. The sensitivity and selectivity of the DSTD signal is related to the surface-active protein concentration in aqueous solution combined with the thermodynamics and kinetics of protein interaction at a liquid-air drop interface. Rapid on-line calibration and measurement of dynamic surface tension is applied, with the surface tension converted into surface pressure results. Continuous surface tension measurement throughout the entire drop growth is achieved, providing insight into kinetic behavior of protein interactive processes at the liquid-air drop interface. Specifically, chemical denaturation of 12 commercial globular proteins-chicken egg albumin, bovine serum albumin, human serum albumin, alpha-lactalbumin (alpha-Lac), myoglobin, cytochrome c, hemoglobin, carbonic anhydrase, alpha-chymotrypsinogen A, beta-lactoglobulin (beta-LG), lysozyme, and glyceraldehyde-3-phosphate-dehydrogenase-is studied in terms of surface pressure (i.e., surface activity) after treatment with increasing concentrations of urea, GdmHCl, and GdmSCN in the 0-8, 0-6, and 0-5 M ranges, respectively. For several of these proteins, the spectroscopic absorbance changes are monitored simultaneously to provide additional information prior to drop formation. Results show that surface pressure of proteins generally increases as the denaturant concentration increases and that effectiveness is GdmSCN > GdmHCl > urea. Protein unfolding curves obtained by plotting surface pressure as a function of denaturant concentration are presented and compared with respect to unfolding curves obtained by using UV absorbance and literature data. Kinetic information relative to the protein adsorption to the air-liquid interface of two proteins, alpha-Lac and beta-LG (chosen as representative proteins for comparison), denatured by the three denaturants is also studied and discussed.

Open reading frame yjbI of Bacillus subtilis codes for truncated hemoglobin

A hypothetical open reading frame from Bacillus subtilis genome, yjbI [NCBI genome database Accession No. ] having homology to many globin and globin-like proteins from different microbial genomes, was selectively amplified from the chromosomal DNA of B. subtilis strain DB104 based on genome sequence database of B. subtilis strain 168. The gene was cloned and over-expressed in Escherichia coli under the transcriptional control of tandem lambda P(L) and P(R) promoters, and the protein was purified to homogeneity. The single-chain monomeric hemoglobin-like protein is stable to the extent of 5.45 kcal/mol at 25 degrees C, binds carbon mono-oxide, and shows optical spectra characteristic of hemoproteins. The protein also exhibits peroxidase-like activity. This is the first report of a truncated bacterial globin endowed with peroxidase-like activity. The activity is enhanced in the presence of urea and guanidine hydrochloride, more so in the presence of the latter. Presumably, only a small portion of the protein is involved in peroxidase activity, which is exposed with increasing concentration of the denaturants.

Detection of coexisting protein conformations in capillary zone electrophoresis subsequent to transient contact with sodium dodecyl sulfate solutions

Non-native conformations of proteins were generated by temporary contact with aqueous solutions of sodium dodecyl sulfate (SDS) and separated from the native state with capillary zone electrophoresis (CZE) in alkaline borate buffer deficient of SDS. Nine proteins at concentrations of 2.0 or 3.0 mg.L(-1) were compared in terms of their susceptibility to SDS. For superoxide dismutase and ferritin the tendency of unfolding was modest with < 25% of the protein being transformed to the non-native state at 10 mmol.L(-1) SDS. Highest susceptibility was observed for albumin, myoglobin (Mb), and hemoglobin with > 75% in the non-native state even at 2.0 mmol.L(-1) SDS. The influence of varying SDS concentrations on the conformational state of Mb was tested. Increasing the SDS concentration, circular dichroism revealed a reduction in alpha-helix, an increase in random coil, and an introduction of beta-sheet, which is absent in native structure. Modifications in the secondary structure were in agreement with distinct changes in the shape of the non-native Mb peak in CZE and make a gradual unfolding/refolding process with several coexisting molten globules instead of two-state transition of conformations most plausible for Mb. CZE was found to contribute to a further understanding of holo-Mb transformation towards a population of non-native conformations (i) by means of calculated peak area ratios of native to non-native states, which showed sigmoid transition, (ii) by detecting the release of the prosthetic heme group, and (iii) by changes in the effective electrophoretic mobility of the Mb-SDS peaks. Reconstituted holo-Mb forms differed in the Soret band around 410 nm, indicating diversity in the conformation of the heme pocket.

Synchrotron radiation circular dichroism spectroscopy applied to metmyoglobin and a 4-?-helix bundle carboprotein

The novel technique, synchrotron radiation-based circular dichroism (SR-CD), has been applied to the study of metmyoglobin and a carboprotein (carbohydrate-based peptide with protein tertiary structure) with 4-alpha-helix bundle structure, as well as a carbopeptide (carbohydrate-based peptide) with a truncated peptide sequence. The use of synchroton radiation (SR) enabled circular dichroism (CD) measurements in the vacuum ultraviolet (VUV) down to 168 nm in D(2)O and 160 nm in 2,2,2-trifluoroethanol (TFE). The band shape in the CD spectra in the low wavelength region was studied, comparing samples with two types of alpha-helical tertiary structure, namely the globin fold and the 4-alpha-helix bundle motif. No significant differences were found between the CD spectra of the alpha-helical samples (metmyoglobin and carboprotein) in D(2)O solution. The use of 2,2,2-TFE (TFE) as solvent clearly alters the VUV CD but the two samples have very similar CD spectra. The solvent-induced denaturing of metmyoglobin in TFE was observed using absorption and CD spectroscopy of the Soret band, with results indicating heme release. The VUV spectrum of TFE-denatured metmyoglobin exhibits dramatic differences in comparison with previous studies of the native enzyme in aqueous solution. The implications of this observation are discussed.

Comparison of the 1.85 A structure of CYP154A1 from Streptomyces coelicolor A3(2) with the closely related CYP154C1 and CYPs from antibiotic biosynthetic pathways

The genus Streptomyces produces two-thirds of microbially derived antibiotics. Polyketides form the largest and most diverse group of these natural products. Antibiotic diversity of polyketides is generated during their biosynthesis by several means, including postpolyketide modification performed by oxidoreductases, a broad group of enzymes including cytochrome P450 monooxygenases (CYPs). CYPs catalyze site-specific oxidation of macrolide antibiotic precursors significantly affecting antibiotic activity. Efficient manipulation of Streptomyces CYPs in generating new antibiotics will require identification and/or engineering of monooxygenases with activities toward a diverse array of chemical substrates. To begin to link structure to function of CYPs involved in secondary metabolic pathways of industrially important species, we determined the X-ray structure of Streptomyces coelicolor A3(2) CYP154A1 at 1.85 A and analyzed it in the context of the closely related CYP154C1 and more distant CYPs from polyketide synthase (EryF) and nonribosomal peptide synthetase (OxyB) biosynthetic pathways. In contrast to CYP154C1, CYP154A1 reveals an active site inaccessible from the molecular surface, and an absence of catalytic activities observed for CYP154C1. Systematic variations in the amino acid patterns and length of the surface HI loop correlate with degree of rotation of the F and G helices relative to the active site in CYP154A1-related CYPs, presumably regulating the degree of active site accessibility and its dimensions. Heme in CYP154A1 is in a 180 degrees flipped orientation compared with most other structurally determined CYPs.

Monomer topology defines folding speed of heptamer

Small monomeric proteins often fold in apparent two-state processes with folding speeds dictated by their native-state topology. Here we test, for the first time, the influence of monomer topology on the folding speed of an oligomeric protein: the heptameric cochaperonin protein 10 (cpn10), which in the native state has seven beta-barrel subunits noncovalently assembled through beta-strand pairing. Cpn10 is a particularly useful model because equilibrium-unfolding experiments have revealed that the denatured state in urea is that of a nonnative heptamer. Surprisingly, refolding of the nonnative cpn10 heptamer is a simple two-state kinetic process with a folding-rate constant in water (2.1 sec(-1); pH 7.0, 20 degrees C) that is in excellent agreement with the prediction based on the native-state topology of the cpn10 monomer. Thus, the monomers appear to fold as independent units, with a speed that correlates with topology, although the C and N termini are trapped in beta-strand pairing with neighboring subunits. In contrast, refolding of unfolded cpn10 monomers is dominated by a slow association step.

Characterization of denatured proteins by hydrophobic interaction chromatography: a preliminary study

In this preliminary study hydrophobic interaction chromatography (HIC) is proposed as a good tool in order to detect conformational changes induced by chemical denaturants in two globular proteins, cytochrome C (Cyt C) and myoglobin (MYO). Alterations in protein structure were manifested chromatographically by reproducible changes in peak heights, retention time, and appearance of multiple peaks. The HIC behavior of the two model proteins denatured by guanidinium thyocyanate (GdmSCN) was investigated, keeping constant various concentrations of urea in the mobile phase in a TSK-Gel Phenyl-5PW column (TosoBiosep). Suitable elution conditions provide evidence of the simultaneous presence of two denatured forms in the case of MYO, and sequential different denatured states of Cyt C.

Effect of gamma-irradiation on the molecular properties of myoglobin

To elucidate the effect of gamma-irradiation on the molecular properties of myoglobin, the secondary and tertiary structures, as well as the molecular weight size of the protein, were examined after irradiation at various irradiation doses. Gamma-irradiation of myoglobin solutions caused the disruption of the ordered structure of the protein molecules, as well as degradation, crosslinking, and aggregation of the polypeptide chains. A SDSPAGE study indicated that irradiation caused initial fragmentation of the proteins and subsequent aggregation, due to cross-linking of the protein molecules. The effect of irradiation on the protein was more significant at lower protein concentrations. Ascorbic acid protected against the degradation and aggregation of proteins by scavenging oxygen radicals that are produced by irradiation. A circular dichroism study showed that an increase of the irradiation decreased the alpha-helical content of myoglobin with a concurrent increase of the aperiodic structure content. Fluorescence spectroscopy indicated that irradiation increased the emission intensity that was excited at 280 nm.

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.

Toward Engineering the Stability and Hemin-Binding Properties of Microsomal Cytochromesb5into Rat Outer Mitochondrial Membrane Cytochromeb5: Examining the Influence of Residues 25 and 71†

As part of a larger effort to engineer the stability and hemin-binding properties of microsomal (Mc) cytochromes b(5) into rat liver outer mitochondrial membrane (OM) cytochrome (cyt) b(5), several mutants of rat OM cyt b(5) were prepared to study the effect of gradual and complete elimination of two extended hydrophobic networks, which are present in the structure of the mitochondrial protein and are absent in the structure of mammalian Mc cytochromes b(5). One of the hydrophobic networks, identified in a previous study [Altuve, A., Silchenko, S., Lee, K.-H., Kuczera, K., Terzyan, S., Zhang, X., Benson, D. R., and Rivera, M. (2001) Biochemistry 40, 9469-9483], encompasses the side chains of Ala-18, Ile-32, Leu-36, and Leu-47, whereas a second hydrophobic network, identified as part of this work, encompasses the side chains of Ile-25, Phe-58, Leu-71, and the heme. The X-ray structure of the A18S/I25L/I32L/L47R/L71S quintuple mutant of rat OM cyt b(5) demonstrates that both hydrophobic networks have been eliminated and that the corresponding structural elements of the Mc isoform have been introduced. The stability of the rat OM mutant proteins studied was found to decrease in the order wild type > I25L > A18S/I32L/L47R > L71S > A18S/I32L/L47R/L71S > 18S/I25L/I32L/L47R/L71S, indicating that the two hydrophobic networks do indeed contribute to the high stability of rat OM cyt b(5) relative to the bovine Mc isoform. Surprisingly, the quintuple mutant of rat OM cyt b(5) is less stable than bovine Mc cyt b(5), even though the former exhibits significantly slower rates of hemin release and hemin reorientation at pH 7.0. However, at pH 5.0 the bovine Mc and rat OM quintuple mutant proteins release hemin at comparable rates, suggesting that one or both of the His axial ligands in the rat OM protein are more resistant to protonation under physiological conditions. Results obtained from chemical denaturation experiments conducted with the apoproteins demonstrated that mutants containing L71S are significantly less stable than bovine Mc apocyt b(5), strongly suggesting that Leu-71 plays a pivotal role in the stabilization of rat OM apocyt b(5), presumably via hydrophobic interactions with Ile-25 and Phe-58. Because comparable interactions are absent in bovine Mc apocyt b(5), which contains Ser at position 71, it must resort to different interactions to stabilize its fold, thus highlighting yet another difference between rat OM and bovine Mc cyt b(5). During the course of these investigations we also discovered that rat OM cyt b(5) can be made to strongly favor hemin orientational isomer A (I32L) or isomer B (L71S) with a single point mutation and that release of hemin orientational isomers A and B can be kinetically resolved in certain rat OM mutants.

Unfolding of the loggerhead sea turtle (Caretta caretta) myoglobin: A (1)H-NMR and electronic absorbance study

The effect of urea concentration on the backbone solution structure of the cyanide derivative of ferric Caretta caretta myoglobin (at pH 5.4) is reported. By addition of urea, sequential and long-range nuclear Overhauser effects (NOEs) are gradually lost. By using the residual NOE constraints to build the molecular model, a picture of the unfolding pathway was obtained. When the urea concentration is raised to 2.2 M, helices A and B appear largely disordered; helices C, D, and F loose structural constraints at 3.0 M urea. At urea concentration >6 M, the protein appears to be fully unfolded, including the GH hairpin and helix E stabilizing the prosthetic group. Reversible and cooperative denaturation isotherms obtained by following NOE peaks are considerably different from those obtained by monitoring electronic absorption changes. The reversible and cooperative urea-dependent folding-unfolding process of C. caretta myoglobin follows the minimum three-state mechanism N long left and right arrow X long left and right arrow D, where X represents a disordered globin structure (occurring at approximately 4 M urea) that still binds the heme.

High stability of a ferredoxin from the hyperthermophilic archaeon A. ambivalens: involvement of electrostatic interactions and cofactors

The ferredoxin from the thermophilic archaeon Acidianus ambivalens is a small monomeric seven-iron protein with a thermal midpoint (T(m)) of 122 degrees C (pH 7). To gain insight into the basis of its thermostability, we have characterized unfolding reactions induced chemically and thermally at various pHs. Thermal unfolding of this ferredoxin, in the presence of various guanidine hydrochloride (GuHCl) concentrations, yields a linear correlation between unfolding enthalpies (DeltaH[T(m)]) and T(m) from which an upper limit for the heat capacity of unfolding (DeltaC(P)) was determined to be 3.15 +/- 0.1 kJ/(mole * K). Only by the use of the stronger denaturant guanidine thiocyanate (GuSCN) is unfolding of A. ambivalens ferredoxin at pH 7 (20 degrees C) observed ([GuSCN](1/2) = 3.1 M; DeltaG(U)[H(2)O] = 79 +/- 8 kJ/mole). The protein is, however, less stable at low pH: At pH 2.5, T(m) is 64 +/- 1 degrees C, and GuHCl-induced unfolding shows a midpoint at 2.3 M (DeltaG(U)[H(2)O] = 20 +/- 1 kJ/mole). These results support that electrostatic interactions contribute significantly to the stability. Analysis of the three-dimensional molecular model of the protein shows that there are several possible ion pairs on the surface. In addition, ferredoxin incorporates two iron-sulfur clusters and a zinc ion that all coordinate deprotonated side chains. The zinc remains bound in the unfolded state whereas the iron-sulfur clusters transiently form linear three-iron species (in pH range 2.5 to 10), which are associated with the unfolded polypeptide, before their complete degradation.