Snapshots of protein folding. A study on the multiple transition state pathway of cytochrome c(551) from Pseudomonas aeruginosa

Cytochrome c(551) (cyt c(551)) from Pseudomonas aeruginosa is a small protein (82 residues) that folds via a three-state pathway with the accumulation in the microsecond time-range of a compact collapsed intermediate. The presence of a single His residue, at position 16, permits the study of the refolding at pH 7.0 in the absence of miscoordination events. Here, we report on folding kinetics in the millisecond time-range as a function of urea under different pH conditions. Analysis of this process (over-and-above proline cis-trans isomerization) at pH 7.0, suggests the existence of a multiple transition state pathway in which we postulate three transition states. Taking advantage of site-directed mutagenesis we propose that the first "unfolded-like" transition state (t(1)) originates from the electrostatic properties of the collapsed state, while the second transition state (t(2)) involves the interaction between the N and C-terminal helices and is stabilized by the salt bridge between Lys10 and Glu70 ( approximately 1 kcal mol(-1)). Our results suggest that, contrary to other cytochromes c, the roll-over effect observed for cyt c(551) at low denaturant concentration can be interpreted in terms of a broad energy barrier without population of any intermediates. The third and more "native-like" transition state (M) can be associated with the breaking/formation of the Fe(3+)-Met61 bond. This strong interaction is stabilized by the hydrogen bond between Trp56 and heme propionate 17 (HP-17) as suggested by the increase in the unfolding rate at high denaturant concentration of the Trp56Phe site-directed mutant.     

Pseudomonas aeruginosa cytochrome C(551): probing the role of the hydrophobic patch in electron transfer

Cytochrome c(551) from Pseudomonas aeruginosa is a monomeric redox protein of 82 amino-acid residues, involved in dissimilative denitrification as the physiological electron donor of cd(1) nitrite reductase. The distribution of charged residues on the surface of c(551) is very anisotropic: one side is richer in acidic residues whereas the other shows a ring of positive side chains, mainly lysines, located at the border of an hydrophobic patch which surrounds the heme crevice. In order to map in cytochrome c(551) the surface involved in electron transfer, we have introduced specific mutations in three residues belonging to the hydrophobic patch, namely Val23-->Asp, Pro58-->Ala and Ile59-->Glu. The effect of these mutations was analyzed studying both the self-exchange rate and the electron-transfer activity towards P. aeruginosa cd(1) nitrite reductase, the physiological partner and P. aeruginosa azurin, a copper protein often used as a model redox partner in vitro. Our results show that introduction of a negative charge in the hydrophobic patch severely hampers both homonuclear and heteronuclear electron transfer.     

Nuclear-magnetic-resonance studies of Pseudomonas aeruginosa cytochrome c-551.

Nuclear magnetic resonance (NMR) spectroscopy was used to study Pseudomonas aeruginosa cytochrome c-551. Assignments of resonances to specific residues have been made. A low-resolution X-ray structure was used to aid assignments. A structural comparison was made between P. aeruginosa cytochrome c-551 and mammalian cytochrome c, based on comparisons of NMR data.     

Electron transfer between azurin from Alcaligenes faecalis and cytochrome c551 from Pseudomonas aeruginosa

The electron transfer equilibrium and kinetics between azurin from Alcaligenes faecalis and cytochrome c551 from Pseudomonas aeruginosa have been studied. The equilibrium constant K = ([Cyt(III)] . [Az(I)])/([Cyt(II)] . [Az(II))]) = 0.5 at 25 degrees C is about seven times smaller than that observed between the cytochrome c551 and the titrations confirmed a 43-mV difference between the mid-point potentials of +266 mV and +309 mV for the Alcaligenes and Pseudomonas azurins respectively. The kinetics of the reaction between Alcaligenes azurin and Pseudomonas cytochrome c551 were investigated by the temperature-jump chemical relaxation method. Only a single relaxation mode was observed throughout the range of concentrations and temperatures examined. Thus, the slow relaxation time observed in the reaction between P. aeruginosa azurin and cytochrome c551 is not observed with the Alcaligenes azurin. The simplest mechanism that can therefore be ascribed to the investigated system is: [formula: see text]. This scheme is similar to that proposed earlier for the reaction between P. aeruginosa azurin and cytochrome c551 but does not involve the conformational transition proposed for azurin. The specific rates for the electron transfer are still fast: 1.8 x 10(6) M-1 . s-1 and 3.0 x 10(6) M-1 . s-1 respectively at 25 degrees C.     

Circular-dichroic properties and secondary structure of Pseudomonas aeruginosa soluble cytochrome c oxidase.

The c.d. spectra of Pseudomonas aeruginosa cytochrome c oxidase in the oxidized state and the reduced state are reported in the visible- and u.v. absorption regions. In the visible region the comparison between the spectra of reduced cytochrome c oxidase and ferrocytochrome c-551 allows the identification of the c.d. bands mainly due to the d1 haem chromophore in cytochrome c oxidase. In the near-u.v. region the assignment of some of the observed peaks to the haem groups and to the aromatic amino acid residues is proposed. A careful analysis of the data in the far-u.v. region leads to the determination of the relative amounts of alpha-helix and beta-sheet in the enzyme, giving for the first time a picture of its secondary structure. A significant difference in this respect between the reduced and the oxidized species is observed and discussed in the light of similar conclusions reported by other workers.     

The reaction of Pseudomonas aeruginosa cytochrome c-551 oxidase with oxygen.

The reaction of ascorbate-reduced Pseudomonas cytochrome oxidase with oxygen was studied by using stopped-flow techniques at pH 7.0 and 25 degrees C. The observed time courses were complex, the reaction consisting of three phases. Of these, only the fastest process, with a second-order rate constant of 3.3 X 10(4) M-1.S-1, was dependent on oxygen concentration. The two slower processes were first-order reactions with rates of 1.0 +/- 0.4s-1 and 0.1 +/- 0.03s-1. A kinetic titration experiment revealed that the enzyme had a relatively low affinity constant for oxygen, approx. 10(4)M-1. Kinetic difference spectra were determined for all three reaction phases, showing each to have different characteristics. The fast-phase difference spectrum showed that changes occurred at both the haem c and haem d1 components of the enzyme during this process. These changes were consistent with the haem c becoming oxidized, but with the haem d1 assuming a form that did not correspond to the normal oxidized state, a situation that was not restored even after the second kinetic phase, which reflected further changes in the haem d1 component. The results are discussed in terms of a kinetic scheme.     

Control of redox properties of cytochrome c by special electrostatic interactions

An assessment is made of the proposal: electrostatic interactions between the ferric ion of oxidised cytochrome c and its haem propionate sidechains assists in determining the value of the redox potential and plays an important role in the redox state conformation change. Differences between the properties of homologous cytochromes are proposed to be due to differences associated with the charge on their haem propionates.     

Sequential 1H NMR assignments of iron(II) cytochrome c551 from Pseudomonas aeruginosa

Sequence-specific 1H NMR resonance assignments for all but the C-terminal Lys 82 are reported for iron(II) cytochrome c551 from Pseudomonas aeruginosa at 25 degrees C and pH = 6.8. Spin systems were identified by using TOCSY and DQF-COSY spectra in 2H2O and 1H2O. Sequential assignments were made by using NOESY connectivities between adjacent amide, alpha, and beta protons. Resonances from several amino acids including His 16, Gly 24, Ile 48, and Met 61 experience strong ring-current shifts due to their placement near the heme. All heme protons, including the previously unassigned propionates, have been identified. Preliminary analysis of sequential and medium-range NOEs provides evidence for substantial amounts of helix in the solution structure. Long-range NOEs indicate that the folds in solution and crystal structures are similar. For one aromatic side chain (Tyr 27) that is close to the heme group we found a transition from hindered ring rotation at low temperature to rapid rotation at high temperature.     

Backbone dynamics and hydrogen exchange of Pseudomonas aeruginosa ferricytochrome c 551

A model-free analysis of Pseudomonas aeruginosa ferricytochrome c(551) dynamics based on (15)N R(1), (15)N R(2), and [(1)H]-(15)N heteronuclear nuclear Overhauser effect data is reported. The protein backbone is highly rigid (< S(2)>=0.924+/-0.005) and displays little variation in picosecond-nanosecond time scale dynamics over the structure. The loop structure containing the axial methionine ligand (loop 3) displays anomalous rigidity, which is attributed to its high proline content. Also reported are protection factors calculated from hydrogen-exchange rates. These data reveal that loop 3 residues, including the axial methionine, are protected from exchange as a result of long-range hydrogen-bonding interactions. These results are contrasted with data reported for Saccharomyces cerevisiae iso-1-ferricytochrome c, which displays higher overall flexibility (< S(2)>=0.80+/-0.07), greater variation of dynamics as a function of structure, and low protection factors for loop 3. This analysis reveals that heme proteins with similar functions and topologies may display diverse dynamical properties.     

Cloning and sequencing of the gene encoding cytochrome c-551 from Pseudomonas aeruginosa

The cytochrome c-551 gene from Pseudomonas aeruginosa was cloned by using two oligonucleotide probes, which had been synthesized based on the known primary structure of the protein. The restriction map of the cloned DNA and sequence analysis showed that the cytochrome c-551 gene is located 50 bp downstream of the nitrite reductase gene, which has recently been cloned and sequenced. DNA sequence analysis also indicated that cytochrome c-551 is synthesized in vivo as a precursor having an amino-terminal signal sequence consisting of 22 amino acid residues.     

Redox-linked spin-state changes in the di-haem cytochrome c-551 peroxidase from Pseudomonas aeruginosa.

Magnetic-c.d., e.p.r. and optical-absorption spectra are reported for the half-reduced form of Pseudomonas aeruginosa cytochrome c-551 peroxidase, a di-haem protein, and its fluoride derivative. Comparison of this enzyme species with oxidized peroxidase shows the occurrence of spin-state changes at both haem sites. The high-potential haem changes its state from partially high-spin to low-spin upon reduction. This is linked to a structural alteration at the ferric low-potential haem group, causing it to change from low-spin to high-spin. Low-temperature spectra demonstrate photolysis of an endogenous ligand of the high-potential haem. In addition, an inactive form of enzyme is examined in which the structural change at the ferric low-potential haem does not occur on reduction of the high-potential haem.     

Folding of horse cytochrome c in the reduced state

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10 degrees C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded DeltaG(H2O)=18.8(+/-1.45) kcal mol(-1), and C(m)=5.1(+/-0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of -164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U<-->N) of 80-90 % of unfolded molecules with a time constant in the range 190-550 micros estimated by linear extrapolation and model calculations. The remaining 10-20 % of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I<-->U<-->N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10-20 % of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pK(a) value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is approximately 300(+/-125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I<-->U<-->N scheme have yielded the microscopic rate constants.     

Oxidation by nitrite of azurin and cytochrome c-551 from Pseudomonas aeruginosa

The nitrite oxidizes reduced azurin and cytochrome c-551 from Pseudomonas aeruginosa. The effects of pH, ionic strength and concentrations of nitrite, EDTA and the protein on the oxidation were investigated. The results obtained indicate that nitrite interacts not only with the terminal electron carrier of the nitrite reducing chain (nitrite reductase, cytochrome cd1) but also with the intermediate electron carrier components of the chain (azurin and cytochrome c-551).     

The evolutionary stability of cytochrome c-551 in Pseudomonas aeruginosa and Pseudomonas fluorescens biotype C

Cytochrome c-551 was prepared from nine different strains of Pseudomonas aeruginosa and six of Pseudomonas fluorescens biotype C, and their amino acid sequences were compared with the sequences previously determined for the cytochromes of type strains of each species. The standard of sequence examination was such that all single amino acid substitutions, delections or insertions ought to have been detected. Balanced double changes in sites in the same part of the sequence might have escaped detection. The standard of some of the quantitative amino acid analyses was not as high as would be required for the investigation of completely unknown sequences. Eight of the Ps. aeruginosa sequences could not be distinguished from the type sequence, whereas the ninth had a single amino acid substitution. The sequences from Ps. fluorescens biotype C were more varied, differing in from zero to four substitutions from the type sequence, with the most diverse sequences differing in seven positions. The results for Ps. aeruginosa are interpreted as evidence that neutral mutations are not responsible for much molecular evolution. The superficially paradoxical differences in the results for the two species are discussed.     

Pseudomonas aeruginosa cytochrome c551 denaturation by five systematic urea derivatives that differ in the alkyl chain length

Reversible denaturation of Pseudomonas aeruginosa cytochrome c₅₅₁ (PAc₅₅₁) could be followed using five systematic urea derivatives that differ in the alkyl chain length, i.e. urea, N-methylurea (MU), N-ethylurea (EU), N-propylurea (PU), and N-butylurea (BU). The BU concentration was the lowest required for the PAc₅₅₁ denaturation, those of PU, EU, MU, and urea being gradually higher. Furthermore, the accessible surface area difference upon PAc₅₅₁ denaturation caused by BU was found to be the highest, those by PU, EU, MU, and urea being gradually lower. These findings indicate that urea derivatives with longer alkyl chains are stronger denaturants. In this study, as many as five systematic urea derivatives could be applied for the reversible denaturation of a single protein, PAc₅₅₁, for the first time, and the effects of the alkyl chain length on protein denaturation were systematically verified by means of thermodynamic parameters.     

Anaerobically induced expression of the nitrite reductase cytochrome c-551 operon from Pseudomonas aeruginosa.

The nitrite reductase gene (denA) and the cytochrome c-551 gene (denB) are located only 50 bp apart from each other in the Pseudomonas aeruginosa chromosome. We report evidence that these two genes are co-transcribed as an operon only under anaerobic (denitrifying) conditions. The nucleotide sequence of the promoter (regulatory) region of the operon is highly AT-rich and contains a sequence closely resembling the consensus FNR binding site in E. coli.     

Thermal aggregation of alpha-chymotrypsin: role of hydrophobic and electrostatic interactions

We have recently reported that electrostatic interactions may play a critical role in alcohol-induced aggregation of alpha-chymotrypsin (CT). In the present study, we have investigated the heat-induced aggregation of this protein. Thermal aggregation of CT obeyed a characteristic pattern, with a clear lag phase followed by a sharp rise in turbidity. Intrinsic and ANS fluorescence studies, together with fluorescence quenching by acrylamide, suggested that the hydrophobic patches are more exposed in the denatured conformation. Typical chaperone-like proteins, including alpha- and beta-caseins and alpha-crystalline could inhibit thermal aggregation of CT, and their inhibitory effect was nearly pH-independent (within the pH range of 7-9). This was partially counteracted by alpha-, beta- and especially gamma-cyclodextrins, suggesting that hydrophobic interactions may play a major role. Loss of thermal aggregation at extreme acidic and basic conditions, combined with changes in net charge/pH profile of aggregation upon chemical modification of lysine residues are taken to support concomitant involvement of electrostatic interactions.     

Stabilization of Pseudomonas aeruginosa cytochrome c(551) by systematic amino acid substitutions based on the structure of thermophilic Hydrogenobacter thermophilus cytochrome c(552)

A heterologous overexpression system for mesophilic Pseudomonas aeruginosa holocytochrome c(551) (PA c(551)) was established using Escherichia coli as a host organism. Amino acid residues were systematically substituted in three regions of PA c(551) with the corresponding residues from thermophilic Hydrogenobacter thermophilus cytochrome c(552) (HT c(552)), which has similar main chain folding to PA c(551), but is more stable to heat. Thermodynamic properties of PA c(551) with one of three single mutations (Phe-7 to Ala, Phe-34 to Tyr, or Val-78 to Ile) showed that these mutants had increased thermostability compared with that of the wild-type. Ala-7 and Ile-78 may contribute to the thermostability by tighter hydrophobic packing, which is indicated by the three dimensional structure comparison of PA c(551) with HT c(552). In the Phe-34 to Tyr mutant, the hydroxyl group of the Tyr residue and the guanidyl base of Arg-47 formed a hydrogen bond, which did not exist between the corresponding residues in HT c(552). We also found that stability of mutant proteins to denaturation by guanidine hydrochloride correlated with that against the thermal denaturation. These results and others described here suggest that significant stabilization of PA c(551) can be achieved through a few amino acid substitutions determined by molecular modeling with reference to the structure of HT c(552). The higher stability of HT c(552) may in part be attributed to some of these substitutions.     

Evolutionary change of the heme c electronic structure: ferricytochrome c-551 from Pseudomonas aeruginosa and horse heart ferricytochrome c.

Individual assignments of the ¹H n.m.r. lines of heme c in reduced and oxidized cytochrome c-551 from $\text{Pseudomonas aeruginosa}$ were obtained by nuclear Overhauser enhancement and saturation transfer experiments. Comparison with the corresponding data on horse heart cytochrome c showed that the locations of high spin density on the heme c periphery as well as the in-plane principal axes x and y of the electronic g-tensor are rotated by approximately 90° in ferricytochrome c-551 relative to horse ferricytochrome c. High spin density in ferricytochrome c-551 is thus localized on the pyrrole ring III. While this pyrrole ring is well shielded in the interior of mammalian-type cytochromes c, it is more easily accessible in cytochrome c-551. It is suggested that this evolutionary change of the heme c electronic structure would be compatible with the hypothesis that the electron transfer in both species is via solvent exposed peripheral ring carbon atoms.     

Solution structure of Fe(II) cytochrome c551 from Pseudomonas aeruginosa as determined by two-dimensional 1H NMR.

The solution structure of Fe(II) cytochrome c551 from Pseudomonas aeruginosa based on 2D 1H NMR data is reported. Two sets of structure calculations were completed with a combination of simulated annealing and distance geometry calculations: one set of 20 structures included the heme-peptide covalent linkages, and one set of 10 structures excluded them. The main-chain atoms were well constrained within the two structural ensembles (1.30 and 1.35 A average RMSD, respectively) except for two regions spanning residues 30-40 and 60-70. The results were essentially the same when global fold comparisons were made between the ensembles with an average RMSD of 1.33 A. In total, 556 constraints were used, including 479 NOEs, 53 volume constraints, and 24 other distances. This report represents the first solution structure determination of a heme protein by 2D 1H NMR and should provide a basis for the application of these techniques to other proteins containing large prosthetic groups or cofactors.