Photooxidation and carbethoxylation of a minor ribonuclease from Aspergillus saitoi.

In order to investigate the nature of amino acid residues involved in the active in the active site of a ribonuclease from Aspergillus saitoi, the pH dependence of the rates of inactivation of RNase Ms by photooxidation and modification with diethylpyrocarbonate were studied. (1) RNase Ms was inactivated by illumination in the presence of methylene blue at various pH's. The pH dependence of the rate of photooxidative inactivation of RNase Ms indicated that at least one functional group having pKa 7.2 was involved in the active site. (2) Amino acid analyses of photooxidized RNase Ms at various stages of photooxidative inactivation at pH's 4.0 and 6.0 indicated that one histidine residue was related to the activity of RNase Ms, but that no tryptophan residue was involved in the active site. (3) 2',(3')-AMP prevented the photooxidative inactivation of RNase Ms. The results also indicated the presence of a histidine residue in the active site. (4) Modification of RNase Ms with diethylpyrocarbonate was studied at various pH's. The results indicated that a functional group having pKa 7.1 was involved in the active site of RNase Ms.     

Ionization properties of titratable groups in ribonuclease T1. I. pKa values in the native state determined by two-dimensional heteronuclear NMR spectroscopy

pKa values of amino acid side chains of ribonuclease T1 have been determined from the pH dependence of 13C and 15N resonances. It was possible to derive pKa values of single protonation or deprotonation sites of carboxylate and imidazole groups. Deviations from pKa values of free amino acids could be interpreted with electrostatic interactions of corresponding side chains with the protein environment. In particular, the interaction between H27 and E82 led to an increase of the H27 pKa and a decrease of the E82 pKa. The pKa of E28 at the C-terminal end of the alpha-helix was increased because of the dipolar character of the alpha-helix. D76 did not titrate in the investigated pH range of about 2-9. From the chemical shift value this buried side chain seems to be protonated. The pKa values of side chains in the active site deviate from a normal behaviour. The lower pKa value of E58 may be interpreted with the close proximity of this side chain with positively charged H40 and R77. A novel two-dimensional 1H(13Cdelta)13Cgamma correlation experiment was developed to observe the pH dependence of the chemical shifts of the Cgamma resonances of histidine residues. From the inspection of the Cgamma chemical shift-pH profiles it was possible to determine the predominant tautomeric form for the histidine residues at higher pH values.     

Structural perturbation of human hemoglobin on glutathionylation probed by hydrogen-deuterium exchange and MALDI mass spectrometry

Glutathionyl hemoglobin, an example of post-translationally modified hemoglobin, has been studied as a marker of oxidative stress in various diseased conditions. Compared to normal hemoglobin, glutathionyl hemoglobin has been found to have increased oxygen affinity and reduced cooperativity. However, detailed information concerning the structural perturbation of hemoglobin associated with glutathionylation is lacking. In the present study, we report structural changes associated with glutathionylation of deoxyhemoglobin by hydrogen/deuterium (H/D) exchange coupled to matrix assisted laser desorption ionization (MALDI) mass spectrometry. We analyzed isotope exchange kinetics of backbone amide hydrogen of eleven peptic peptides in the deoxy state of both hemoglobin and glutathionyl hemoglobin molecules. Analysis of the deuterium incorporation kinetics for both molecules showed structural changes associated with the following peptides: α34-46, α1-29, β32-41, β86-102, β115-129, and β130-146. H/D exchange experiments suggest that glutathionylation of hemoglobin results in a change in conformation located at the above-mentioned regions of the hemoglobin molecule. MALDI mass spectrometry based H/D exchange experiment might be a simple way of monitoring structural changes associated with post-translational modification of protein.     

Assignment of the histidine proton magnetic resonance peaks of soybean trypsin inhibitor (Kunitz) by a differertial deuterium exchange technique.

Deuterium exchange at the C(2)-H position of the two histidine residues of native soybean trypsin inhibitor (Kunitz) in 2-H2O was followed by 1-H nuclear magnetic resonance (NMR) spectroscopy. The two histidine residues of soybean trypsin inhibitor exchange at significantly different rates at pH* 5.00, 40 degrees. Half-times observed were: peak H1, t1/2=61 plus or minus 2 days; peak H2, T1/2=24 plus or minus 2 days. Differentially deuterated soybean trypsin inhibitor was cleaved by cyanogen bromide into two fragments each containing one histidine residue. The deuterium content of the histidine residue of each separated fragment was analyzed by 1H NMR spectroscopy. Hisidine-71 in fragment 1-114 showed approximately twice the deuterium content of His-157 in fragment 115-181. These results lead to the assignment of 1H NMR peak H1 to His-157 and peak H2 to His-71. These assignments were extended to the histidine peaks of trypsin-modified soybean trypsin inhibitor by converting the differentially deuterated virgin soybean trypsin inhibitor to the modified form. The correlation of histidine peaks in virgin amd modified soybean trypsin inhibitors was the same as proposed earlier on the basis of pK arguments. The results demonstrate that His-71 is the residue whose pK value is raised from 5.27 to 5.91 on trypsin modification of soybean trypsin inhibitor [Markley, J. L., (1973), Biochemistry 12, 2245].     

Copper and the oxidation of hemoglobin: a comparison of horse and human hemoglobins.

Oxidation studies of hemoglobin by Cu(II) indicate that for horse hemoglobin, up to a Cu(II)/heme molar ratio of 0.5, all of the Cu(II) added is used to rapidly oxidize the heme. On the other hand, most of the Cu(II) added to human hemoglobin at low Cu(II)/heme molar ratios is unable to oxidize the heme. Only at Cu(II)/heme molar ratios greater than 0.5 does the amount of oxidation per added Cu(II) approach that of horse hemoglobin. At the same time, binding studies indicate that human hemoglobin has an additional binding site involving one copper for every two hemes, which has a higher copper affinity than the single horse hemoglobin binding site. The Cu(II) oxidation of human hemoglobin is explained utilizing this additional binding site by a mechanism where a transfer of electrons cannot occur between the heme and the Cu(II) bound to the high affinity human binding site. The electron transfer must involve the Cu(II) bound to the lower affinity human hemoglobin binding site, which is similar to the only horse hemoglobin site. The involvement of beta-2 histidine in the binding of this additional copper is indicated by a comparison of the amino acid sequences of various hemoglobins which possess the additional site, with the amino acid sequences of hemoglobins which do not possess the additional site. Zn(II), Hg(II), and N-ethylmaleimide (NEM) are found to decrease the Cu(II) oxidation of hemoglobin. The sulfhydryl reagents, Hg(II) and NEM, produce a very dramatic decrease in the rate of oxidation, which can only be explained by an effect on the rate for the actual transfer of electrons between the Cu(II) and the Fe(II). The effect of Zn(II) is much smaller and can, for the most part, be explained by the increased oxygen affinity, which affects the ligand dissociation process that must precede the electron transfer process.     

Ionizable groups linked to the reaction of 2,2'-dithiobispyridine with hemoglobin.

The pH-dependence of the second-order rate-constant for the reaction of 2,2'-dithiobispyridine with the CysF9(93) beta sulphydryl group of hemoglobin in the R quaternary structure is analyzed in terms of a tentative model based on the observation that this sulphydryl exists as a mixture of two tertiary conformations in dynamic equilibrium. For the four aquomethemoglobins studied (human A and S, dog and rabbit), the equation derived from this model gives a better fit than a simpler equation based on the assumption of only one tertiary conformation. For the corresponding carbonmonoxyhemoglobins the simpler equation gives a better fit. The dog and rabbit oxy and azidomet data are better fitted by the model equation, whereas the data for the corresponding human A and S derivatives are better fitted by the simpler equation. From the analysis pKa values of 6.1 and 8.7 are obtained for the ionization of groups coupled to the presumed conformational transition. The pKa of 6.1 is assigned to HisHC3(146) beta; the pKa of 8.7 is assigned to the CysF9(93) beta sulphydryl group in its external conformation. It is estimated that the pKa of this sulphydryl may be as high as 12.9 in its internal conformation.     

pH and kinetic isotope effects on the reductive half-reaction of D-amino acid oxidase.

Primary deuterium kinetic isotope and pH effects on the reduction of D-amino acid oxidase by amino acid substrates were determined using steady-state and rapid reaction methods. With D-serine as substrate, reduction of the enzyme-bound FAD requires that a group with a pKa value of 8.7 be unprotonated and that a group with a pKa value of 10.7 be protonated. The DV/Kser value of 4.5 is pH-independent, establishing that these pKa values are intrinsic. The limiting rate of reduction of the enzyme shows a kinetic isotope effect of 4.75, consistent with this as the intrinsic value. At high enzyme concentration (approximately 15 microM) at pH 9,D-serine is slightly sticky (k3/k2 = 0.8), consistent with a decrease in the rate of substrate dissociation. With D-alanine as substrate, the pKa values are perturbed to 8.1 and 11.5. The DV/Kala value increases from 1.3 at pH 9.5 to 5.1 at pH 4, establishing that D-alanine is sticky with a forward commitment of approximately 10. The effect of pH on the DV/Kala value is consistent with a model in which exchange with solvent of the proton from the group with pKa 8.7 is hindered and is catalyzed by H2O and OH- above pH 7 and by H3O+ and H2O below pH 7. With glycine, the pH optimum is shifted to a more basic value, 10.3. The DV/Kgly value increases from 1.26 at pH 6.5 to 3.1 at pH 10.7, consistent with fully reversible CH bond cleavage followed by a pH-dependent step. At pH 10.5, the kinetic isotope effect on the limiting rate of reduction is 3.4.     

1H electron-nuclear double resonance of cobalt hemoglobin

Proton 1H electron-nuclear double resonance spectra were obtained from cobalt-substituted hemoglobin A. For oxy cobalt-substituted hemoglobin (CoHb), two major couplings were found. One is exchangeable by deuterium and has a coupling of approximately 5.6 MHz which is assigned to the N epsilon 2-proton of the distal histidine (E7), interacting with the ligand. The other splitting (approximately 2.4-3.2 MHz) is attributed to a methyl-proton on Val (E11). No coupling of protons on the proximal side is resolved. Deoxy CoHb exhibits only one prominent interaction (approximately 1.3 MHz) which is assigned to one of the protons on the proximal histidine (F8). The influence of the distal amino acids His (E7) and Val (E11) on the ligand affinity in oxy CoHb is discussed.     

Crystallographic analysis of the interaction of nitric oxide with quaternary-T human hemoglobin

In addition to interacting with hemoglobin as a heme ligand to form nitrosylhemoglobin, NO can react with cysteine sulfhydryl groups to form S-nitrosocysteine or cysteine oxides such as cysteinesulfenic acid. Both modes of interaction are very sensitive to the quaternary structure of hemoglobin. To directly view the interaction of NO with quaternary-T deoxyhemoglobin, crystallographic studies were carried out on crystals of deoxyhemoglobin that were exposed to gaseous NO under a variety of conditions. Consistent with previous spectroscopic studies in solution, these crystallographic studies show that the binding of NO to the heme groups of crystalline wild-type deoxyhemoglobin ruptures the Fe-proximal histidine bonds of the alpha-subunits but not the beta-subunits. This finding supports Perutz's theory that ligand binding induces tension in the alpha Fe-proximal histidine bond. To test Perutz's theory, deoxy crystals of the mutant hemoglobin betaW37E were exposed to NO. This experiment was carried out because previous studies have shown that this mutation greatly reduces the quaternary constraints that oppose the ligand-induced movement of the alpha-heme Fe atom into the plane of the porphyrin ring. As hypothesized, the Fe-proximal histidine bonds in both the beta- and the alpha-subunits remain intact in crystalline betaW37E after exposure to NO. With regard to S-nitrosocysteine or cysteine oxide formation, no evidence for the reaction of NO with any cysteine residues was detected under anaerobic conditions. However, when deoxyhemoglobin crystals are first exposed to air and then to NO, the appearance of additional electron density indicates that Cys93(F9)beta has been modified, most likely to cysteinesulfenic acid. This modification of Cys93(F9)beta disrupts the intrasubunit salt bridge between His146(HC3)beta and Asp94(FG1)beta, a key feature of the quaternary-T hemoglobin structure. Also presented is a reanalysis of our previous crystallographic studies [Chan, N.-L., et al. (1998) Biochemistry 37, 16459-16464] of the interaction of NO with liganded hemoglobin in the quaternary-R2 structure. These studies showed additional electron density at Cys93(F9)beta that was consistent with an NO adduct. However, for reasons discussed in this paper, we now believe that this adduct may be the Hb-S-N.-O-H radical intermediate and not Hb-S-N=O as previously suggested.     

Biochemical and genetic studies on rabbit hemoglobin. I. Electrophoretic polymorphism of theβ chain

A genetic polymorphism of beta-chain rabbit (Oryctolagus cuniculus) hemoglobin is demonstrated by means of acid starch gel electrophoresis. The biochemical evidence presented suggests that a previously reported substitution of a neutral amino acid for a histidine is responsible for the detected genetic variation. Segregation analysis was performed in a sample of 15 matings with 49 offspring and confirmed the genetic hypothesis: two common alleles at an autosomal locus. The calculated gene frequencies in a random sample of 125 individuals are HBB*1 = 0.48 and HBB*2 = 0.52.     

Hydrogen-tritium exchange titration of the histidine residues in bovine heart cytochrome c and analysis of their microenvironment.

Microenvironments of the three histidine residues located at the positions 18, 26, and 33 from the amino terminus in bovine heart cytochrome c were analysed in solution by the hydrogen-tritium exchange titration method, which has been developed in this laboratory. Histidine-18, which is liganded to the heme iron, and histidine-26 did not incorporate tritium in native state, indicating that the two are located in solvent inaccessible hydrophobic regions. Histidine-33 was labeled with tritium to an appreciable extent and seemed to be partially buried in the molecule. The pKa value estimated for histidine-33 was 6.1 at 37 degrees by the tritium exchange titration, suggesting that the residue interacts very weakly with a neighboring cationic group. These results seem to be compatible with the tertiary structure of the protein deduced from the X-ray crystallographic analysis.     

Heme-linked ionization of horseradish peroxidase compound II monitored by the resonance Raman Fe(IV)=O stretching vibration

Fe(IV)=O resonance Raman stretching vibrations were recently identified by this laboratory for horseradish peroxidase compound II and ferryl myoglobin. In the present report it is shown that Fe(IV)=O stretching frequency for horseradish peroxidase compound II will switch between two values depending on pH, with pKa values corresponding to the previously reported compound II heme-linked ionizations of pKa = 6.9 for isoenzyme A-2 and pKa = 8.5 for isoenzyme C. Similar pH-dependent shifts of the Fe(IV)=O frequency of ferryl myoglobin were not detected above pH 6. The Fe(IV)=O stretching frequencies of compound II of the horseradish peroxidase isoenzymes at pH values above the transition points were at a high value approaching the Fe(IV)=O stretching frequency of ferryl myoglobin. Below the transition points the horseradish peroxidase frequencies were found to be 10 cm-1 lower. Frequencies of the Fe(IV)=O stretching vibrations of horseradish peroxidase compound II for one set of isoenzymes were found to be sensitive to deuterium exchange below the transition point but not above. These results were interpreted to be indicative of an alkaline deprotonation of a distal amino acid group, probably histidine, which is hydrogen bonded to the oxyferryl group below the transition point. Deprotonation of this group at pH values above the pKa disrupts hydrogen bonding, raising the Fe(IV)=O stretching frequency, and is proposed to account for the lowering of compound II reactivity at alkaline pH. The high value of the Fe(IV)=O vibration of compound II above the transition point appears to be identical in frequency to what is believed to be the Fe(IV)=O vibration of compound X.     

The hemoglobin of the sea lamprey,Petromyzon marinus

The blood hemoglobin of the sea lamprey presents a curious mixture of primitive and highly specialized properties. Like muscle hemoglobin, it has a molecular weight of about 17,000, and apparently contains a single heme. Its isoelectric point is like that of a typical invertebrate hemoglobin. Its amino acid composition is partly characteristic of invertebrate) partly of vertebrate hemoglobins (Pedersen; Roche and Fontaine). In the present experiments, the oxygen equilibrium curve of this pigment was measured at several pH's. As expected, it is a rectangular hyperbola, the first such function to be observed in a vertebrate blood hemoglobin. Other hemoglobins known to possess this type of oxygen dissociation curve—those of vertebrate muscle, the worm Nippostrongylus, and the bot-fly larva—appear to serve primarily the function of oxygen storage rather than transport. Lamprey hemoglobin on the contrary is an efficient oxygen-transporting agent. It achieves this status by having, unlike muscle hemoglobin, a relatively low oxygen affinity, and a very large Bohr effect. In these properties it rivals the most effective vertebrate blood hemoglobins.     

1H-n.m.r. studies of the isolated activation segment from pig procarboxypeptidase A.

The isolated activation segment (asA) from pig pancreatic procarboxypeptidase A was studied by 1H-n.m.r. spectroscopy over a wide range of solution conditions. Isolated asA shows many characteristics of compactly folded globular proteins, such as the observation of perturbed positions for resonances from methyl groups, alpha-carbon atoms, histidine residues and the tyrosine residue. The single tyrosine residue (Tyr-70) exhibits a very high pKa, and both histidine and tyrosine residues show slow chemical modification (deuteration and iodination). In contrast, asA shows rapid NH exchange. Analysis of the spectra by pH titration and nuclear Overhauser effects revealed several residue interactions. Quantitative analysis of deuterium and tritium exchange allowed the assignment of the histidine C-2-H resonances to their respective residues in the sequence. His-66, the closest to the sites of proteolytic attack in the proenzyme, is shown to be the most accessible to solvent in procarboxypeptidase A. It was also shown that asA is thermally very stable ['melting' temperature (Tm) 88 degrees C] and requires a high urea concentration for denaturation (6.25 M, at pH 7.5). Evidence is presented for some degree of conformational flexibility in the premelting range, a feature that could be ascribed to the preponderance of helical secondary structure and to the lack of disulphide bridges. The free solution structure of asA is probably unchanged when it binds to carboxypeptidase A.     

A competitive labeling method for the determination of the chemical properties of solitary functional groups in proteins.

The properties of the functional groups in a protein can be used as built-in-probes of the structure of the protein. We have developed a general procedure whereby the ionization constant and chemical reactivity of solitary functional groups in proteins may be determined. The method may be applied to the side chain of histidine, tyrosine, lysine, and cysteine, as well as to the amino terminus of the protein. The method, which is an extension of the competitive labeling technique using [3H]- and [14C]1-fluoro-2,4-dinitrobenzene (N2ph-F) in a double-labeling procedure, is rapid and sensitive. Advantage is taken of the fact that after acid hydrolysis of a dinitrophenylated protein, a derivative is obtained which must be derived from a unique position in the protein. The method has been applied to the solitary histidine residue of lysozyme, alpha-lytic protease, and Streptomyces griseus (S.G.) trypsin, as well as to the amino terminus of the latter protein. The following parameters were obtained for reaction with N2ph-F at 20 degrees C in 0.1 N KCl: the histidine of hen egg-white lysozyme, pKa of 6.4 and second-order velocity constant of 0.188 M-1 min-1; the histidine of alpha-lytic protease, pKa of 6.5 and second-order velocity constant of 0.0235 M-1 min-1; the histidine of S.G. trypsin, pKa of 6.5 and second-order velocity constant of 0.0328 M-1 min-1; the valyl amino terminus of S.G. trypsin, pKa of 8.1 and second-order velocity constant of 0.403 M-1 min-1. In addition, the results obtained provide clues as to the microenvironments of these functional groups, and indicate that the proteins studied undergo pH-dependent conformational changes which affect the microenvironment, and hence the chemical reactivity of these groups.     

Nuclear magnetic resonance titration curves of histidine ring protons. A direct assignment of the resonances of the active site histidine residues of ribonuclease.

One of the four titrating histidine ring C-2 proton resonances of bovine pancreatic ribonuclease has been assigned to histidine residue 12. This was accomplished by a direct comparison of the rate of tritium incorporation into position C-2 of histidine 12 of S-peptide (residues 1 to 20) derived from ribonuclease S, with the rates of deuterium exchange of the four histidine C-2 proton resonances of ribonuclease S under the same experimental conditions. The same assignment was obtained by a comparison of the NMR titration curves of ribonuclease S, the noncovalent complex of S-peptide and S-protein (residues 21 to 124) with the results for the recombined complex in which position C-2 of histidine 12 was fully deuterated. The second active site histidine resonance was assigned to histidine residue 119 by consideration of the NMR titration results fro carboxymethylated histidines and 1-carboxymethylhistidine 119 ribonuclease. This assignment is a reversal of that originally reported, and has important implications for the interpretation of NMR titration data of ribonuclease.     

The amino acid sequence and oxygen-binding properties of the single hemoglobin of the cold-adapted Antarctic teleost Gymnodraco acuticeps.

The complete amino acid sequence of the single hemoglobin of the Antarctic teleost Gymnodraco acuticeps has been determined. The alpha chain contains 142 amino acid residues; an acetylated seryl residue is at the amino terminal. The beta chain contains 146 residues. A very high degree of sequence identity has been found with hemoglobins of other Antarctic fishes. Oxygen binding is not modulated by pH and allosteric effectors. The Bohr and Root effects are absent, although specific amino acid residues, considered responsible of most of these functions, are conserved in the sequence, thus posing new questions about the molecular basis of these mechanisms. The low heat of oxygenation may be interpreted as one of the mechanisms involved in the process of cold adaptation.     

Properties of carboxymethylated cross-linked hemoglobin A

The selective carboxymethylation of the N-terminal amino groups of hemoglobin A with glyoxylic acid and sodium cyanoborohydride has been studied as a function of the state of ligation of hemoglobin. The N-terminal residues have been established as the primary sites of reaction by peptide mapping of the tryptic digest of each chain and subsequent amino acid analysis of the modified peptides. With oxyhemoglobin, the desired derivatives with a carboxymethyl group at the N-terminal of either or both chains amounted to 55% [Di Donato, A., Fantl, W. J., Acharya, A. S., & Manning, J. M. (1983) J. Biol. Chem. 258, 11890-11895]. In the present study it is shown that with deoxyhemoglobin the amount of the desired derivative is increased to 75%. The oxygen equilibrium curve of hemoglobin A carboxymethylated on its four N-terminal residues [0.5 mM as tetramer in 50 mM [bis(2-hydroxyethyl)amino]tris(hydroxymethyl)methane (Bis-Tris), pH 7.5, 37 degrees C] had a P50 value of 30 mmHg (Hill coefficient n = 2.8, alkaline Bohr value = 0.4) compared to a P50 of 9 mmHg for unmodified hemoglobin under the same conditions (n = 2.5, alkaline Bohr value = 0.5). In carboxymethylated oxyhemoglobin A, cross-linked with the mild agent glycolaldehyde for 3.5 h, there was 85% of Mr 64,000 species and 15% of Mr 128,000 or higher species. For the former, the extent of cross-linking between two subunits was 19%. For the latter, there was 29% of two cross-linked subunits and 13% of three cross-linked subunits. Termination of cross-linking, which may be desirable in some circumstances, can be successfully achieved with isonicotinic acid hydrazide.(ABSTRACT TRUNCATED AT 250 WORDS)     

Histidine-aromatic interactions in barnase. Elevation of histidine pKa and contribution to protein stability.

Aromatic side-chains are found in the vicinity of histidine residues in many proteins and protein complexes. We have studied the interaction between a histidine residue (His18) and aromatic residues at position 94 in barnase. Three different techniques have been applied to show that Trp94 interacts more strongly with the protonated form of His18. The aromatic-histidine interaction stabilizes the protonated form of histidine by 0.8 to 1 kcal mol-1 relative to the unprotonated and, thereby, increases its pKa value. This was shown indirectly from the pH dependence of the stability of the wild-type protein and the mutant Trp94----Leu; and directly from the difference in pKa of His18 between wild-type barnase and the same mutant protein, and from double-mutant cycles that measure the total interaction energy of Trp94 with His18 at both low and high pH. When Trp94 is replaced by other aromatic amino acids, the strength of the interaction decreases in the series His-Trp greater than His-Tyr greater than His-Phe. The interaction is not masked by high salt concentrations. The raising of the pKa value of His18 by interaction with Trp94 is shown to be consistent with solution studies with model compounds. The histidine-aromatic interaction could have implications in binding and catalysis for modulation of the histidine pKa value.     

Modification of human hemoglobin by glutathione. III. Perturbations of hemoglobin conformation analyzed by computer modeling

The perturbations of the conformation of human deoxyhemoglobin induced by the covalent attachment of glutathione at cysteine beta 93 have been investigated by computer simulation in conjunction with molecular graphics. In the first phase of the analysis, a systematic search was carried out of the conformational space of glutathione attached to deoxyhemoglobin. In this search, the conformation of the hemoglobin molecule was held constant, while the relative energies of a series of 186,624 glutathione conformations involving systematic variation of six dihedral angels were calculated. From this search, the most favorable conformation was selected as the starting conformation for energy minimization of the glutathionyl hemoglobin molecule as a function of all Cartesian coordinates. In order to provide a reference state, an independent minimization by the same procedures was carried out for deoxyhemoglobin in the absence of glutathione. Comparison of the minimized structures with and without glutathione attached revealed a number of significant differences. The most conspicuous difference in the protein moiety concerned the salt bridge between aspartate beta 94 and histidine beta 146 which is destabilized upon minimization of the glutathionyl-hemoglobin complex due to interactions of the aspartate residue with the glycyl NH group of glutathione. Other observed differences in the minimized structures are located at the alpha 1-beta 2 interface and include displacement of the carboxyl group of aspartate beta 99. In the minimized complex, the glutathione portion assumes a quasi-cyclic conformation stabilized through interactions between the free (gamma-glutamyl) amino and (glycyl) carboxyl ends of the tripeptide and between this carboxyl end and the epsilon amino group of lysine alpha 40. In a parallel conformational study of glutathione alone, a similar structure was found as the lowest energy form. These quasi-cyclic conformations contrast with the extended structures reported by Wright (Wright, W.B. (1955) Acta Crystallogr. 11, 632-642) for crystals of glutathione where interactions between molecules play a major role. The conclusions of our analysis are in agreement with the experimental investigations reported in the two preceding papers and permit, moreover, a coherent interpretation of the observed functional and structural changes in deoxyhemoglobin induced by glutathione.