EPR studies of iso-1-cytochrome c: effect of temperature on two-component spectra of spin label attached to cysteine at positions 102 and 47

Wild-type iso-1-cytochrome c from Saccharomyces cerevisiae containing naturally occurring cysteine at position 102 and mutated protein S47C (derived from the protein in which C102 had been replaced by threonine) were labeled with cysteine-specific methanethiosulfonate spin label. Continuous wave (CW) electron paramagnetic resonance (EPR) was used to examine the effect of temperature on the behavior of the spin label in the oxidized and reduced forms of wild-type cytochrome c and in the oxidized form of the mutated protein. The computer simulations revealed that the CW EPR spectrum for each form of cytochrome c consists of at least two components [a fast (F) and a slow (S) component], which differ in the values of the rotational correlation times tauRparallel (longitudinal rotational correlation time) and tauRperpendicular (transverse rotational correlation time) and that the relative contributions of the F and S components of the spectra change with temperature. In addition, the values of the rotational correlation times (tauRparallel and tauRperpendicular) for the F component appear to change much more dramatically with the temperature than the respective values for the S component. A large difference between the behavior of the oxidized and reduced wild-type spin-labeled cytochromes c indicates that the temperature-induced unfolding of the protein in the region around C102 progresses more rapidly when cytochrome c is in the oxidized form.     

EPR-detected folding kinetics of externally located cysteine-directed spin-labeled mutants of iso-1-cytochrome c.

We report the application of our newly developed dielectric resonator-based flow and stopped-flow kinetic EPR systematically to probe protein folding in yeast iso-1-cytochrome c at cysteine-directed spin-labeled locations. The locations studied have not been previously directly probed by other techniques, and we observe them on a time scale stretching from 50 micros to seconds. On the basis of crystal structure and homology information, the following mutation-tolerant, externally located cysteine labeling sites were chosen (in helices, T8C, E66C, and N92C; in loops, E21C, V28C, H39C, D50C, and K79C), and labeling at these sites was not destabilizing. Dilution of denaturant was used to induce folding and thereby to cause a change in the spin label EPR signal as folding altered the motion of the spin label. Under folding conditions, including the presence of imidazole to eliminate kinetic trapping due to heme misligation, a phase of folding on the 20-30 ms time scale was found. This phase occurred not only at the T8C and N92C labeling sites in the N- and C-terminal helices, where such a phase has been associated with folding in these helices, but overall at labeling sites throughout the protein. In the absence of imidazole the 20-30 ms phase disappeared, and another phase having the time scale of 1 s appeared throughout the protein. There was evidence under all conditions for a burst phase on a scale of less than several milliseconds which occurred at labeling positions V28C, H39C, D50C, E66C, and K79C in the middle of the protein sequence. At spin-labeled D50C rapid-mix flow EPR indicated a very short approximately 50 micros phase possibly associated with the prefolding or compaction of the loop to which D50 belongs. Spin labels have been criticized as perturbing the phenomena which they measure, but our spin labeling strategy has reported common kinetic themes and not perturbed, disconnected kinetic events.     

Mapping RNA-protein interactions in ribonuclease P from Escherichia coli using electron paramagnetic resonance spectroscopy

Ribonuclease P (RNase P) is a catalytic ribonucleoprotein (RNP) essential for tRNA biosynthesis. In Escherichia coli, this RNP complex is composed of a catalytic RNA subunit, M1 RNA, and a protein cofactor, C5 protein. Using the sulfhydryl-specific reagent (1-oxyl-2,2,5, 5-tetramethyl-Delta3-pyrroline-3-methyl)methanethiosulfonate (MTSL), we have introduced a nitroxide spin label individually at six genetically engineered cysteine residues (i.e., positions 16, 21, 44, 54, 66, and 106) and the native cysteine residue (i.e., position 113) in C5 protein. The spin label covalently attached to any protein is sensitive to structural changes in its microenvironment. Therefore, we expected that if the spin label introduced at a particular position in C5 protein was present at the RNA-protein interface, the electron paramagnetic resonance (EPR) spectrum of the spin label would be altered upon binding of the spin-labeled C5 protein to M1 RNA. The EPR spectra observed with the various MTSL-modified mutant derivatives of C5 protein indicate that the spin label attached to the protein at positions 16, 44, 54, 66, and 113 is immobilized to varying degrees upon addition of M1 RNA but not in the presence of a catalytically inactive, deletion derivative of M1 RNA. In contrast, the spin label attached to position 21 displays an increased mobility upon binding to M1 RNA. The results from this EPR spectroscopy-based approach together with those from earlier studies identify residues in C5 protein which are proximal to M1 RNA in the RNase P holoenzyme complex.     

Spin-labeling studies of the conformational changes in the vicinity of D36, D38, T46, and E161 of bacteriorhodopsin during the photocycle.

Electron paramagnetic resonance (EPR) spectroscopy of site-directed spin-labeled bacteriorhodopsin mutants is used to study structural changes during the photocycle. After exchange of the native amino acids D36 and D38 in the A-B loop, E161 in the E-F loop, and T46 in the putative proton channel by cysteines, these positions were modified by a methanethiosulfonate spin label. Time-resolved EPR spectroscopy reveals spectral changes during the photocycle for the mutants with spin labels attached to C36, C161, and C46. A comparison of the transient spectral amplitudes with simulated EPR difference spectra shows that the detected signals are due to changes in the spin label mobility and not to possible polarity changes in the vicinity of the attached spin label. The kinetic analysis of the EPR and the visible data with a global fitting procedure exhibits a structural rearrangement near position 161 in the E-F loop in the M state. The environmental changes at positions 36 and 46, however, occur during the M-to-N transition. All structural changes reverse with the recovery of the BR ground state. No structural changes are detected with a spin label attached to C38.     

Influence of the disulfide bond configuration on the dynamics of the spin label attached to cytochrome c

A series of multi-nanosecond molecular dynamics (MD) simulations of wild-type cytochrome c and its spin-labeled variants with the methanethiosulfonate moiety attached at position C102 were performed (1) to elucidate the effect of the spin probe presence on the protein structure and (2) to describe the structure and dynamics of the spin-label moiety. Comparisons with the reference crystal structure of cytochrome c (PDB entry: 1YCC) indicate that the protein secondary structure is well preserved during simulations of the wild-type cytochrome c but slightly changed in simulations of the cytochrome c labeled at position C102. At the time scale covered in our simulations, the spin label exhibits highly dynamical behavior. The number of observed distinct conformations of the spin label moiety is between 3 and 13. The spin probe was found to form short-lived hydrogen bonds with the protein. Temporary hydrophobic interactions between the probe and the protein were also detected. The MD simulations directly show that the disulfide bond in the tether linking a spin probe with a protein strongly influence the behavior of the nitroxide group. The conformational flexibility and interaction with the protein are different for each of the two low energy conformations of the disulfide bond.     

The Role of the Fast Motion of the Spin Label in the Interpretation of EPR Spectra for Spin-Labeled Macromolecules

Abstract

The spin label method was used to observe the nature of the fast motions of side chains in protein monocrystals. The EPR spectra of spin-labeled lysozyme monocrystals (with different orientations of the tetragonal protein crystal in relation to the direction of the magnetic field) were interpreted using the method of molecular dynamics (MD). Within the proposed simple model, MD calculations of the spin label motion trajectories are performed in a reasonable real time. The model regards the protein molecule as frozen as a whole and the spin labeled amino acid residue as unfrozen. To calculate the trajectories in vacuum, a model of spin-labeled lysozyme was assembled, and the parameters of the force fields were specified for atoms of the protein molecule, including the spin label. The calculations show that the protein environment sterically limits the area of the possible angular reorientations for the NO reporter group of the nitroxide (within the spin label), and this, in turn, affects the shape of the EPR spectrum. However, it turned out that the spread in the positions of the reporter group in the angle space strictly adheres to the Gaussian distribution. Using the coordinates of the spin label atoms obtained by the MD method within a selected time range and considering the distribution of the spin label states over the ensemble of spin-labeled macro- molecules in a crystal, the EPR spectra of spin-labeled lysozyme monocrystals were simulated. The resultant theoretical EPR spectra appeared to be similar to experimental ones.     

Membrane location of apocytochrome c and cytochrome c determined from lipid-protein spin exchange interactions by continuous wave saturation electron spin resonance.

Apocytochrome c derived from horse heart cytochrome c was spin-labeled on the cysteine residue at position 14 or 17 in the N-terminal region of the primary sequence, and cytochrome c from yeast was spin-labeled on the single cysteine residue at sequence position 102 in the C-terminal region. The spin-labeled apocytochrome c and cytochrome c were bound to fluid bilayers composed of different negatively charged phospholipids that also contained phospholipid probes that were spin-labeled either in the headgroup or at different positions in the sn-2 acyl chain. The location of the spin-labeled cysteine residues on the lipid-bound proteins was determined relative to the spin-label positions in the different spin-labeled phospholipids by the influence of spin-spin interactions on the microwave saturation properties of the spin-label electron spin resonance spectra. The enhanced spin relaxation observed in the doubly labeled systems arises from Heisenberg spin exchange, which is determined by the accessibility of the spin-label group on the protein to that on the lipid. It is found that the labeled cysteine groups in horse heart apocytochrome c are located closest to the 14-C atom of the lipid acyl chain when the protein is bound to dimyristoyl- or dioleoyl-phosphatidylglycerol, and to that of the 5-C atom when the protein is bound to a dimyristoylphosphatidylglycerol/dimyristoylphosphatidylcholine (15:85 mol/mol mixture. On binding to dioleoylphosphatidylglycerol, the labeled cysteine residue in yeast cytochrome c is located closest to the phospholipid headgroups but possibly between the polar group region and the 5-C atom of the acyl chains. These data determine the extent to which the different regions of the proteins are able to penetrate negatively charged phospholipid bilayers.     

Rational design of a more stable yeast iso-1-cytochrome c.

Yeast iso-1-cytochrome c is one of the least stable mitochondrial cytochromes c. We have used a coordinated approach, combining the known functional and structural properties of cytochromes c, to engineer mutations into yeast iso-1-cytochrome c with the goal of selectively increasing the stability of the protein. The two redox forms of the native protein and six different mutant forms of yeast iso-1-cytochrome c were analyzed by differential scanning calorimetry (DSC). The relative stability, expressed as the difference in the Gibb's free energy of denaturation at a given temperature between the native and mutant forms (DeltaDeltaG(Tref)), was determined for each of the proteins. In both oxidation states, the mutant proteins C102T, T69E/C102T, T96A/C102T, and T69E/T96A/C102T were more stable than the wild-type protein, respectively. The increased stability of the mutant proteins is proposed to be due to the removal of a rare surface cysteine and the stabilization of two distorted alpha-helices.     

The role of the fast motion of the spin label in the interpretation of EPR spectra for spin-labeled macromolecules

The spin label method was used to observe the nature of the fast motions of side chains in protein monocrystals. The EPR spectra of spin-labeled lysozyme monocrystals (with different orientations of the tetragonal protein crystal in relation to the direction of the magnetic field) were interpreted using the method of molecular dynamics (MD). Within the proposed simple model, MD calculations of the spin label motion trajectories are performed in a reasonable real time. The model regards the protein molecule as frozen as a whole and the spin-labeled amino acid residue as unfrozen. To calculate the trajectories in vacuum, a model of spin-labeled lysozyme was assembled, and the parameters of the force fields were specified for atoms of the protein molecule, including the spin label. The calculations show that the protein environment sterically limits the area of the possible angular reorientations for the NO reporter group of the nitroxide (within the spin label), and this, in turn, affects the shape of the EPR spectrum. However, it turned out that the spread in the positions of the reporter group in the angle space strictly adheres to the Gaussian distribution. Using the coordinates of the spin label atoms obtained by the MD method within a selected time range and considering the distribution of the spin label states over the ensemble of spin-labeled macromolecules in a crystal, the EPR spectra of spin-labeled lysozyme monocrystals were simulated. The resultant theoretical EPR spectra appeared to be similar to experimental ones.     

Role of rapid movement of spin labels in interpreting EPR spectra for spin-labelled macromolecules

The method of spin labeling was used to monitor quick movements of side residues in protein monocrystals. The EPR spectra of monocrystals of spin-labeled lysozyme at different orientations of the tetrahonal crystal relative to the direction of the magnetic field were interpreted using the molecular dynamics method. A simple model was proposed, which enables one to calculate the trajectory of movements of the spin label by the molecular dynamic method over a relatively short period of time. The entire "frozen" protein molecule and a "defrozen" spin-labeled amino acid residue were considered in the framework of the model. To calculate the trajectories in vacuum, a model of spin-labeled lysozyme was constructed, and the parameters of force potentials for the atoms of the protein molecule and the spin label were specified. It follows from the calculations that the protein environment sterically hinders the range of eventual angular reorientations of the reporter NO-group of nitroxyl incorporated into the spin label, thereby affecting the shape of the EPR spectrum. However, the scatter in the positions of the reporter group in the angular space turned out to correspond to the Gauss distribution. Using the atomic coordinates of the spin label, obtained in a chosen time interval by the method of molecular dynamics, and taking into account the distribution of the states of the spin label in the ensemble of spin-labeled macromolecules in the crystal, we simulated the EPR spectra of monocrystals of spin-labeled lysozyme. The theoretical EPR spectra coincide well with the experimental.     

Spin label studies on the interactions of bovine adrenodoxin with NADPH-adrenodoxin reductase and with cytochrome P-450scc.

Adrenodoxin of bovine adrenocortical mitochondria was spin-labeled with two different spin-labeling reagents, N-(2,2,5,5-tetramethyl-3-carbonylpyrroline-1-oxyl)imidazole (I) and N-(1-oxyl-2,2,6,6-tetramethyl-4-piperidinyl)maleimide (II), without major loss of its activity for electron transport from NADPH to cytochrome c. The EPR spectrum of adrenodoxin spin-labeled with either of the reagents showed a pattern typical of a moderately immobilized spin label. When adrenodoxin was treated with (I), approximately two amino acid residues per molecule were spin-labeled, whereas a single residue was labeled by (II). While assition of NADPH to adrenodoxin spin-labeled with (I) did not diminish the EPR signal intensity, addition of the reductant to the labeled adrenodoxin in the presence of adrenodoxin reductase caused slow reduction of the spin label, the rate of which was dependent on the aerobicity. Addition of adrenodoxin reductase to adrenodoxin spin-labeled with (I) or (II) resulted in the appearance of a more immobilized component in the EPR spectrum. The ratio of the more immobilized component to the less immobilized component was saturated at a molar ratio of one to one. Addition of cytochrome P-450scc to adrenodoxin labeled with (I) had similar effects on the EPR spectrum.     

Dynamic spin label study of the barstar-barnase complex

The dynamic spin label method was used to study protein-protein interactions in the model complex of the enzyme barnase (Bn) with its inhibitor barstar. The C40A mutant of barstar (Bs) containing a single cysteine residue was modified with two different spin labels varying in length and structure of a flexible linker. Each spin label was selectively bound to the Cys82 residue, located near the Bn-Bs contact site. The formation of the stable protein complex between Bn and spin labeled Bs was accompanied by a substantial restriction of spin label mobility, indicated by remarkable changes in the registered EPR spectra. Order parameter, S, as an estimate of rapid reorientation of spin label relative to protein molecule, was sharply increasing approaching 1. However, the rotational correlation time tau for spin-labeled Bs and its complex with Bn in solution corresponded precisely to their molecular weights. These data indicate that both Bs and its complex with Bn are rigid protein entities. Spin labels attached to Bs in close proximity to an interface of interaction with Bn, regardless of its structure, undergo significant restriction of mobility by the environment of the contact site of the two proteins. The results show that this approach can be used to investigate fusion proteins containing Bn or Bs.     

Structure determination and analysis of yeast iso-2-cytochrome c and a composite mutant protein.

As part of a study of protein folding and stability, the three-dimensional structures of yeast iso-2-cytochrome c and a composite protein (B-2036) composed of primary sequences of both iso-1 and iso-2-cytochromes c have been solved to 1.9 A and 1.95 A resolutions, respectively, using X-ray diffraction techniques. The sequences of iso-1 and iso-2-cytochrome c share approximately 84% identity and the B-2036 composite protein has residues 15 to 63 from iso-2-cytochrome c with the rest being derived form the iso-1 protein. Comparison of these structures reveals that amino acid substitutions result in alterations in the details of intramolecular interactions. Specifically, the substitution Leu98Met results in the filling of an internal cavity present in iso-1-cytochrome c. Further substitutions of Val20Ile and Cys102Ala alter the packing of secondary structure elements in the iso-2 protein. Blending the isozymic amino acid sequences in this latter area results in the expansion of the volume of an internal cavity in the B-2036 structure to relieve a steric clash between Ile20 and Cys102. Modification of hydrogen bonding and protein packing without disrupting the protein fold is illustrated by the His26Asn and Asn63Ser substitutions between iso-1 and iso-2-cytochromes c. Alternatively, a change in main-chain fold is observed at Gly37 apparently due to a remote amino acid substitution. Further structural changes occur at Phe82 and the amino terminus where a four residue extension is present in yeast iso-2-cytochrome c. An additional comparison with all other eukaryotic cytochrome c structures determined to date is presented, along with an analysis of conserved water molecules. Also determined are the midpoint reduction potentials of iso-2 and B-2036 cytochromes c using direct electrochemistry. The values obtained are 286 and 288 mV, respectively, indicating that the amino acid substitutions present have had only a small impact on the heme reduction potential in comparison to iso-1-cytochrome c, which has a reduction potential of 290 mV.     

Effects of ouabain on the rotational dynamics of renal Na,K-ATPase studied by saturation-transfer EPR.

The interaction of a nitroxide spin-labeled derivative of ouabain with sheep kidney Na,K-ATPase and the motional behavior of the ouabain spin label-Na,K-ATPase complex have been studied by means of electron paramagnetic resonance (EPR) and saturation-transfer EPR (ST-EPR). Spin-labeled ouabain binds with high affinity to the Na,K-ATPase with concurrent inhibition of ATPase activity. Enzyme preparations retain 0.61 +/- 0.1 mol of bound ouabain spin label per mole of ATP-dependent phosphorylation sites, even after repeated centrifugation and resuspension of the purified ATPase-containing membrane fragments. The conventional EPR spectrum of the ouabain spin label bound to the ATPase consists almost entirely (greater than 99%) of a broad resonance at 0 degrees C, characteristic of a tightly bound spin label which is strongly immobilized by the protein backbone. Saturation-transfer EPR measurements of the spin-labeled ATPase preparations yield effective correlation times for the bound labels significantly longer than 100 microseconds at 0 degrees C. Since the conventional EPR measurements of the ouabain spin-labeled Na,K-ATPase indicated the label was strongly immobilized, these rotational correlation times most likely represent the motion of the protein itself rather than the independent motion of mobile spin probes relative to a slower moving protein. Additional ST-EPR measurements of ouabain spin-labeled Na,K-ATPase (a) cross-linked with glutaraldehyde and (b) crystallized in two-dimensional arrays indicated that the observed rotational correlation times predominantly represented the motion of large Na,K-ATPase-containing membrane fragments, as opposed to the motion of individual monomeric or dimeric polypeptides within the membrane fragment. The results suggest that the binding of spin-labeled ouabain to the ATPase induces the protein to form large aggregates, implying that cardiac glycoside induced enzyme aggregation may play a role in the mechanism of action of the cardiac glycosides in inhibiting the Na,K-ATPase.     

Variable velocity liquid flow EPR applied to submillisecond protein folding

We have developed a variable velocity, rapid-mix, continuous-flow method for observing and delineating kinetics by dielectric resonator-based electron paramagnetic resonance (EPR). The technology opens a new facet for kinetic study of radicals in liquid at submillisecond time resolution. The EPR system (after Sienkiewicz, A., K. Qu, and C. P. Scholes. 1994. Rev. Sci. Instrum. 65:68-74) accommodated a miniature quartz capillary mixer with an approximately 0.5 microliter delivery volume to the midpoint of the EPR-active zone. The flow velocity was varied in a preprogrammed manner, giving a minimum delivery time of approximately 150 microseconds. The mixing was efficient, and we constructed kinetics in the 0.15-2. 1-ms time range by plotting the continuous wave EPR signal taken during flow versus the reciprocal of flow velocity. We followed the refolding kinetics of iso-1-cytochrome c spin-labeled at Cysteine 102. At 20 degrees C, upon dilution of guanidinium hydrochloride denaturant, a fast phase of refolding was resolved with an exponential time constant of 0.12 ms, which was consistent with the "burst" phase observed by optically detected flow techniques. At 7 degrees C the kinetic refolding time of this phase increased to 0.5 ms.     

Spin label method reveals barnase-barstar interaction: a temperature and viscosity dependence approach

Spin labeling was used to study the protein-protein interaction between the enzyme barnase (Bn) and its inhibitor barstar (Bs). A mutant of barstar (C40A), which contains only one cysteine, C82, located near the Bn-Bs contact region, was selectively modified by two spin labels having different lengths and structures of the flexible tether. The formation of a strong protein complex resulted in significant restriction of spin label mobility at the C82 residue of barstar, as indicated by notable changes in the recorded EPR spectra. The dependence of the separation between broad outer peaks of the EPR spectra on viscosity at constant temperature was used to evaluate the order parameter S and the rotational correlation time tau (a temperature-viscosity dependence approach). The order parameter S, which characterizes fast reorientation of a spin label relative to the protein molecule, sharply increases and approaches unity when Bs binds to Bn. In addition, formation of a Bs-Bs complex was observed; it is also accompanied by restriction of spin label mobility. At the same time, the rotational correlation times tau of spin-labeled Bs, its complex with Bn, and the Bs dimer in solution agree well with their molecular masses. This indicates that barstar, its complex with barnase, and barstar dimer are rigid protein entities. The described approach is suitable for studying any spin-labeled macromolecular complexes.     

Spin-Labeled Nucleotide Mobility in the Boundary of the EcoRI Endonuclease Binding Site

Abstract

A complex consisting of the EcoRI endonuclease site-specifically bound to spin-labeled DNA 26mers was prepared to provide a model system for studying possible conformational changes resulting from protein binding. EPR was used to monitor the mobility of the spin labels that were strategically placed in position 6, 9, or 11 with respect to the dyad axis of the 26mer. These positions are located within the flanking region on either side of the EcoRI hexamer binding site. This allows the monitoring of potential distal structural changes in the DNA helix caused by protein binding. The spectral line shapes indicate that the spin label closest to the EcoRI endonuclease binding site, i.e., in position 6, is most influenced by the binding event. The EPR data are analyzed according to a model that distinguishes between spectral effects due to a change in the hydrodynamic shape of the complex and those resulting from local variations in the spin-label mobility as characterized by a local order parameter S. S reflecting the motional restriction of the spin-labeled base is 0.20 ± 0.01 for all three oligomers as well as for the two complexes with the label in position 9 or 11, while the position 6 labeled complex yields S=0.25. To further evaluate the origin of the slightly larger EPR effect observed with position 6 labeled material, molecular dynamics (MD) simulations were used to explore the space accessible to the probes in positions 6, 9, and 11. MD results gave similar nitroxide trajectories for all three labeled 26mers in the absence or presence of EcoRI. Thus, the small position 6 effect is attributed to a structural distortion in the major groove of the DNA at this location possibly corresponding to a bend induced by protein binding. The observation that the spectral changes are small indicates the absence of any significant structural disruption being propagated along the helix as a result of protein binding. Also, the fact that the line shape of the 26mers did not change as expected from hydrodynamic theory in view of the significant increase in molecular volume upon protein binding suggests that there are additional relaxation processes involving the protein and nucleic acid.     

Assignment of proton resonances, identification of secondary structural elements, and analysis of backbone chemical shifts for the C102T variant of yeast iso-1-cytochrome c and horse cytochrome c

Resonance assignments for the main-chain, side-chain, exchangeable side chain, and heme protons of the C102T variant of Saccharomyces cerevisiae iso-1-cytochrome c in both oxidation states (with the exception of Gly-83) are reported. (We have also independently assigned horse cytochrome c.) Some additional assignments for the horse protein extend those of Wand and co-workers [Wand, A. J., Di Stefano, D. L., Feng, Y., Roder, H., & Englander, S. W. (1989) Biochemistry 28, 186-194; Feng, Y., Roder, H., Englander, S. W., Wand, A. J., & Di Stefano, D. L. (1989) Biochemistry 28, 195-203]. Qualitative interpretation of nuclear Overhauser enhancement data allows the secondary structure of these two proteins to be described relative to crystal structures. Comparison of the chemical shift of the backbone protons of the C102T variant and horse protein reveals significant differences resulting from amino acid substitution at positions 56 and 57 and further substitutions between residue 60 and residue 69. Although the overall folding of yeast iso-1-cytochrome c and horse cytochrome c is very similar, there can be large differences in chemical shift for structurally equivalent residues. Chemical shift differences of amide protons (and to a lesser extent alpha protons) represent minute changes in hydrogen bonding. Therefore, great care must be taken in the use of differences in chemical shift as evidence for structural changes even between highly homologous proteins.     

Hydrogen exchange behavior of [U-15N]-labeled oxidized and reduced iso-1-cytochrome c

Heteronuclear NMR spectroscopy was used to measure the hydrogen-deuterium exchange rates of backbone amide hydrogens in both oxidized and reduced [U-15N]iso-1-cytochrome c from the yeast Saccharomyces cerevisiae. The exchange data confirm previously reported data [Marmorino et al. (1993) Protein Sci. 2, 1966-1974], resolve several inconsistencies, and provide more thorough coverage of exchange rates throughout the cytochrome c protein in both oxidation states. Combining the data previously collected on unlabeled C102T with the current data collected on [U-15N]C102T, exchange rates for 53 protons in the oxidized state and 52 protons in the reduced state can now be reported. Most significantly, hydrogen exchange measurements on [U-15N]iso-1-cytochrome c allowed the observation of exchange behavior of the secondary structures, such as large loops, that are not extensively hydrogen-bonded. For the helices, the most slowly exchanging protons are found in the middle of the helix, with more rapidly exchanging protons at the helix ends. The observation for the Omega-loops in cytochrome c is just the opposite. In the loops, the ends contain the most slowly exchanging protons and the loop middles allow more rapid exchange. This is found to be true in cytochrome c loops, even though the loop ends are not attached to any regular secondary structures. Some of the exchange data are strikingly inconsistent with data collected on the C102S variant at a different pH, which suggests pH-dependent dynamic differences in the protein structure. This new hydrogen exchange data for loop residues could have implications for the substructure model of eukaryotic cytochrome c folding. Isotopic labeling of variant forms of cytochrome c can now be used to answer many questions about the structure and folding of this model protein.