Search

The association of complement component C3d with B or T cell complement receptor 2 (CR2 or CD21) is a link between innate and adaptive immunity. It has been recognized in experimental studies that the C3d-CR2 association is pH- and ionic strength-dependent. This led us to perform electrostatic calculations to obtain a theoretical understanding of the mechanism of C3d-CR2 association. We used the crystallographic structures of human free C3d, free CR2 (short consensus repeat (SCR)1-2), and the C3d-CR2(SCR1-2) complex, and continuum solvent representation, to obtain a detailed atomic-level picture of the components of the two molecules that contribute to association. Based on the calculation of electrostatic potentials for the free and bound species and apparent pK(a) values for each ionizable residue, we show that C3d-CR2(SCR1-2) recognition is electrostatic in nature and involves not only the association interface, but also the whole molecules. Our results are in qualitative agreement with experimental data that measured the ionic strength and pH dependence of C3d-CR2 association. Also, our results for the native molecules and a number of theoretical mutants of C3d explain experimental mutagenesis studies of amino acid replacements away from the association interface that modulate binding of iC3b with full-length CR2. Finally, we discuss the packing of the two SCR domains. Overall, our data provide global and site-specific explanations of the physical causes that underlie the ionic strength dependence of C3d-CR2 association in a unified model that accounts for all experimental data, some of which were previously thought to be contradictory.

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR)., (Copyright 2003 Wiley-Liss, Inc.)

We have predicted the free energy of unfolding for the pH-dependent helix-coil transition of the activation helix of GART using continuum electrostatic calculations and structural modeling. We have assigned the contributions of each element of secondary structure and of each ionizable residue, within and in the vicinity of the activation helix, to the stability of several fragments of GART that participate in the formation of the catalytic site. We demonstrate that the interaction of His121-His132 contributes 2.2 kcal/mol to the ionization free energy between pH 0 and approximately 6. We also show that the ionization state of a network of five histidines, His108, His119, His121, His132 and His137, and two aspartic acids Asp141 and Asp144, contributes approximately 12 kcal/mol to the stability of the catalytic site of GART, out of a total stability of 16 kcal/mol of the whole enzyme. These interactions are important for the formation of the catalytic site of GART.

Compstatin, a 13-mer cyclic peptide, is a novel and promising inhibitor of the activation of the complement system. In our search for a more active analog and better understanding of structure-functions relations, we designed a phage-displayed random peptide library based on previous knowledge of structure activity relations, in which seven amino acids deemed necessary for structure and activity were kept fixed while the remaining six were optimized. Screening of this library against C3 identified four binding clones. Synthetic peptides corresponding to these clones revealed one analog, called acetylated Ile(1)Leu/His(9)Trp/Thr(13)Gly triple replacement analog of compstatin corresponding to clone 640 (Ac-I1L/H9W/T13G), which was more active than compstatin. This newly identified peptide had 4-fold higher activity when compared with the originally isolated form of compstatin and 1.6-fold higher activity when compared with acetylated compstatin (Ac-compstatin). The structures of Ac-I1L/H9W/T13G and Ac-compstatin were studied by nuclear magnetic resonance, compared with the structure of compstatin, and found to be very similar. The binding of Ac-I1L/H9W/T13G and the equally active acetylated analog with His(9)Ala replacement (Ac-H9A) to C3 was evaluated by surface plasmon resonance, which suggested similarity in their binding mechanism but difference when compared with Ac-compstatin. Compensatory effects of flexibility outside the beta-turn and tryptophan ring stacking may be responsible for the measured activity increase in Ac-I1L/H9W/T13G and acetylated analog with His(9)Ala replacement and the variability in binding mechanism compared with Ac-compstatin. These data demonstrate that tryptophan is a key amino acid for activity. Finally, the significance of the N-terminal acetylation was examined and it was found that the hydrophobic cluster at the linked termini of compstatin is essential for binding to C3 and for activity.

The subcellular localization of cAMP-dependent protein kinase (PKA) occurs through interaction with A-Kinase Anchoring Proteins (AKAPs). AKAPs bind to the PKA regulatory subunit dimer of both type Ialpha and type IIalpha (RIalpha and RIIalpha). RIalpha and RIIalpha display characteristic localization within different cell types, which is maintained by interaction of AKAPs with the N-terminal dimerization and docking domain (D/D) of the respective regulatory subunit. Previously, we reported the solution structure of RIIa D/D module, both free and bound to AKAPs. We have now solved the solution structure of the dimerization and docking domain of the type Ialpha regulatory dimer subunit (RIalpha D/D). RIalpha D/D is a compact docking module, with unusual interchain disulfide bonds that help maintain the AKAP interaction surface. In contrast to the shallow hydrophobic groove for AKAP binding across the surface of the RIIalpha D/D dimeric interface, the RIalpha D/D module presents a deep cleft for proposed AKAP binding. RIalpha and RIIalpha D/D interaction modules present drastically differing dimeric topographies, despite a conserved X-type four-helix bundle structure.

A novel structure-activity-based combinatorial computational optimization methodology for the design of peptides that are candidates to become therapeutics is presented. This methodology has been successfully applied in the design of a 7-fold more active analogue, among other active analogues, in the case of the complement inhibitor compstatin. The main steps of the approach involve the availability of NMR-derived structural templates, combinatorial selection of sequences based on optimization of parametrized pairwise residue interaction potentials, prediction of fold stabilities using deterministic global optimization, and experimental validation with immunological activity measurements. This work is direct evidence that an integrated experimental and theoretical approach can make the engineering of compounds with enhanced immunological properties possible.

Although activation of complement protein C3 is essential for the generation of normal inflammatory responses against pathogens, its unregulated activation during various pathological conditions leads to host cell damage. Previously we have identified a 13-residue cyclic peptide, Compstatin, that inhibits C3 activation. In this study, we have examined the species-specificity of Compstatin. Bimolecular interaction analysis using a real-time surface plasmon resonance-based assay showed that Compstatin exhibits exclusive specificity for primate C3s and does not bind either to C3s from lower mammalian species or to two structural homologs of C3, human C4 and C5. Furthermore, it showed that although the kinetics of binding of Compstatin to non-human primate C3s were distinctly different from those to human C3, like human C3 its mechanism of binding to non-human primate C3 was biphasic and did not follow a simple 1:1 interaction, suggesting that this binding mechanism could be important for its inhibitory activity. Analysis of Ala substitution analogs of Compstatin for their inhibitory activities against mouse and rat complement suggested that the lack of binding of Compstatin to mouse and rat C3s was not a result of sterically hindered access to the binding pocket due to individual bulky side chains or the presence of charge on the Compstatin molecule. These results suggest that Compstatin's exclusive specificity for primate C3s could be exploited for the development of species-specific complement inhibitors.

The complement (C') system has long been recognized as an important mediator of innate immune defense and inflammation. In recent years there is increasing evidence suggesting that complement components may also participate in non-inflammatory and developmental processes. Here we review our current work on the structural-functional aspects of C3-ligand interactions and the rational design of small-sized complement inhibitors. We present a novel, proteomics-based, approach to studying protein-protein interactions within the C' system and discuss our progress in the study of viral immune evasion strategies. Furthermore we discuss the involvement of complement proteins in organ regeneration and hematopoietic development.

We have previously identified compstatin, a 13-residue cyclic peptide, that inhibits complement activation by binding to C3 and preventing C3 cleavage to C3a and C3b. The structure of compstatin consists of a disulfide bridge and a type I beta-turn located at opposite sides to each other. The disulfide bridge is part of a hydrophobic cluster, and the beta-turn is part of a polar surface. We present the design of compstatin analogs in which we have introduced a series of perturbations in key structural elements of their parent peptide, compstatin. We have examined the consistency of the structures of the designed analogs compared with compstatin using NMR, and we have used the resulting structural information to make structure-complement inhibitory activity correlations. We propose the following. 1) Even in the absence of the disulfide bridge, a linear analog has a propensity for structure formation consistent with a turn of a 3(10)-helix or a beta-turn. 2) The type I beta-turn is a necessary but not a sufficient condition for activity. 3) Our substitutions outside the type I beta-turn of compstatin have altered the turn population but not the turn structure. 4) Flexibility of the beta-turn is essential for activity. 5) The type I beta-turn introduces reversibility and sufficiently separates the two sides of the peptide, whereas the disulfide bridge prevents the termini from drifting apart, thus aiding in the formation of the hydrophobic cluster. 6) The hydrophobic cluster at the linked termini is involved in binding to C3 and activity but alone is not sufficient for activity. 7) beta-Turn residues Gln(5) (Asn(5))-Asp(6)-Trp(7)(Phe(7))-Gly(8) are specific for the turn formation, but only Gln(5)(Asn(5))-Asp(6)-Trp(7)-Gly(8) residues are specific for activity. 8) Trp(7) is likely to be involved in direct interaction with C3, possibly through the formation of a hydrogen bond. Finally we propose a binding model for the C3-compstatin complex.

The structure of the N-terminal docking and dimerization domain of the type IIalpha regulatory subunit (RIIalpha D/D) of protein kinase A (PKA) forms a noncovalent stand-alone X-type four-helix bundle structural motif, consisting of two helix-loop-helix monomers. RIIalpha D/D possesses a strong hydrophobic core and two distinct, exposed faces. A hydrophobic face with a groove is the site of protein-protein interactions necessary for subcellular localization. A highly charged face, opposite to the former, may be involved in regulation of protein-protein interactions as a result of changes in phosphorylation state of the regulatory subunit. Although recent studies have addressed the hydrophobic character of packing of RIIalpha D/D and revealed the function of the hydrophobic face as the binding site to A-kinase anchoring proteins (AKAPs), little attention has been paid to the charges involved in structure and function. To examine the electrostatic character of the structure of RIIalpha D/D we have predicted mean apparent pKa values, based on Poisson-Boltzmann electrostatic calculations, using an ensemble of calculated dimer structures. We propose that the helix promoting sequence Glu34-X-X-X-Arg38 stabilizes the second helix of each monomer, through the formation of a (i, i +4) side chain salt bridge. We show that a weak inter-helical hydrogen bond between Tyr35-Glu19 of each monomer contributes to tertiary packing and may be responsible for discriminating from alternative quaternary packing of the two monomers. We also show that an inter-monomer hydrogen bond between Asp30-Arg40 contributes to quaternary packing. We propose that the charged face comprising of Asp27-Asp30-Glu34-Arg38-Arg40-Glu41-Arg43-Arg44 may be necessary to provide flexibility or stability in the region between the C-terminus and the interdomain/autoinhibitory sequence of RIIalpha, depending on the activation state of PKA. We also discuss the structural requirements necessary for the formation of a stacked (rather than intertwined) dimer, which has consequences for the orientation of the functionally important and distinct faces.

The enzyme glycinamide ribonucleotide transformylase (GART) catalyzes the transfer of a formyl group from formyl tetrahydrofolate (fTHF) to glycinamide ribonucleotide (GAR), a process that is pH-dependent with pK(a) of approximately 8. Experimental studies of pH-rate profiles of wild-type and site-directed mutants of GART have led to the proposal that His108, Asp144, and GAR are involved in catalysis, with His108 being an acid catalyst, while forming a salt bridge with Asp144, and GAR being a nucleophile to attack the formyl group of fTHF. This model implied a protonated histidine with pK(a) of 9.7 and a neutral GAR with pK(a) of 6.8. These proposed unusual pK(a)s have led us to investigate the electrostatic environment of the active site of GART. We have used Poisson-Boltzmann-based electrostatic methods to calculate the pK(a)s of all ionizable groups, using the crystallographic structure of a ternary complex of GART involving the pseudosubstrate 5-deaza-5,6,7,8-THF (5dTHF) and substrate GAR. Theoretical mutation and deletion analogs have been constructed to elucidate pairwise electrostatic interactions between key ionizable sites within the catalytic site. Also, a construct of a more realistic catalytic site including a reconstructed pseudocofactor with an attached formyl group, in an environment with optimal local van der Waals interactions (locally minimized) that imitates closely the catalytic reactants, has been used for pK(a) calculations. Strong electrostatic coupling among catalytic residues His108, Asp144, and substrate GAR was observed, which is extremely sensitive to the initial protonation and imidazole ring flip state of His108 and small structural changes. We show that a proton can be exchanged between GAR and His108, depending on their relative geometry and their distance to Asp144, and when the proton is attached on His108, catalysis could be possible. Using the formylated locally minimized construct of GART, a high pK(a) for His108 was calculated, indicating a protonated histidine, and a low pK(a) for GAR(NH(2)) was calculated, indicating that GAR is in neutral form. Our results are in qualitative agreement with the current mechanistic picture of the catalytic process of GART deduced from the experimental data, but they do not reproduce the absolute magnitude of the pK(a)s extracted from fits of k(cat)-pH profiles, possibly because the static time-averaged crystallographic structure does not describe adequately the dynamic nature of the catalytic site during binding and catalysis. In addition, a strong effect on the pK(a) of GAR(NH(2)) is produced by the theoretical mutations of His108Ala and Asp144Ala, which is not in agreement with the observed insensitivity of the pK(a) of GAR(NH(2)) modeled from the experimental data using similar mutations. Finally, we show that important three-way electrostatic interactions between highly conserved His137, with His108 and Asp144, are responsible for stabilizing the electrostatic microenvironment of the catalytic site. In conclusion, our data suggest that further detailed computational and experimental work is necessary.

Glycinamide ribonucleotide transformylase (GART) undergoes a pH-dependent coil-helix transition with pK(a) approximately 7. An alpha-helix is formed at high pH spanning 8 residues of a 21-residue-long loop, comprising the segment Thr120-His121-Arg122-Gln123-Ala124-Leu125-Glu126-Asn127. To understand the electrostatic nature of this loop-helix, called the activation loop-helix, which leads to the formation and stability of the alpha-helix, pK(a) values of all ionizable residues of GART have been calculated, using Poisson-Boltzmann electrostatic calculations and crystallographic data. Crystallographic structures of high and low pH E70A GART have been used in our analysis. Low pK(a) values of 5.3, 5.3, 3.9, 1.7, and 4.7 have been calculated for five functionally important histidines, His108, His119, His121, His132, and His137, respectively, using the high pH E70A GART structure. Ten theoretical single and double mutants of the high pH E70A structure have been constructed to identify pairwise interactions of ionizable residues, which have aided in elucidating the multiplicity of electrostatic interactions of the activation loop-helix, and the impact of the activation helix on the catalytic site. Based on our pK(a) calculations and structural data, we propose that: (1) His121 forms a molecular switch for the coil-helix transition of the activation helix, depending on its protonation state; (2) a strong electrostatic interaction between His132 and His121 is observed, which can be of stabilizing or destabilizing nature for the activation helix, depending on the relative orientation and protonation states of the rings of His121 and His132; (3) electrostatic interactions involving His119 and Arg122 play a role in the stability of the activation helix; and (4) the activation helix contains the helix-promoting sequence Arg122-Gln123-Ala124-Leu125-Glu126, but its alignment relative to the N and C termini of the helix is not optimal, and is possibly of a destabilizing nature. Finally, we provide electrostatic evidence that the formation and closure of the activation helix create a hydrophobic environment for catalytic-site residue His108, to facilitate catalysis.

The specificity of intracellular signaling events is controlled, in part, by compartmentalization of protein kinases and phosphatases. The subcellular localization of these enzymes is often maintained by protein- protein interactions. A prototypic example is the compartmentalization of the cAMP-dependent protein kinase (PKA) through its association with A-kinase anchoring proteins (AKAPs). A docking and dimerization domain (D/D) located within the first 45 residues of each regulatory (R) subunit protomer forms a high affinity binding site for its anchoring partner. We now report the structures of two D/D-AKAP peptide complexes obtained by solution NMR methods, one with Ht31(493-515) and the other with AKAP79(392-413). We present the first direct structural data demonstrating the helical nature of the peptides. The structures reveal conserved hydrophobic interaction surfaces on the helical AKAP peptides and the PKA R subunit, which are responsible for mediating the high affinity association in the complexes. In a departure from the dimer-dimer interactions seen in other X-type four-helix bundle dimeric proteins, our structures reveal a novel hydrophobic groove that accommodates one AKAP per RIIalpha D/D.

We have previously identified a 13-residue cyclic peptide, Compstatin, that binds to complement component C3 and inhibits complement activation. Herein, we describe the binding kinetics, structure-activity relationship, and biotransformation of Compstatin. Biomolecular interaction analysis using surface-plasmon resonance showed that Compstatin bound to native C3 and its fragments C3b and C3c, but not C3d. While binding of Compstatin to native C3 was biphasic, binding to C3b and C3c followed the 1:1 Langmuir binding model; the affinities of Compstatin for C3b and C3c were 22- and 74-fold lower, respectively, than that of native C3. Analysis of Compstatin analogs synthesized for structure-function studies indicated that 1) the 11-membered ring between disulfide-linked Cys2-Cys12 constitutes a minimal structure required for optimal activity; 2) retro-inverso isomerization results in loss of inhibitory activity; and 3) some residues of the type I beta-turn segment also interact with C3. In vitro studies of Compstatin in human blood indicated that a major pathway of biotransformation was the removal of Ile1, which could be blocked by N-acetylation of the peptide. These findings indicate that acetylated Compstatin is stable against enzymatic degradation and that the type I beta-turn segment is not only critical for preservation of the conformational stability, but also involved in intermolecular recognition.

Compartmentalization of signal transduction enzymes into signaling complexes is an important mechanism to ensure the specificity of intracellular events. Formation of these complexes is mediated by specialized protein motifs that participate in protein-protein interactions. The adenosine 3',5'-cyclic monophosphate (cAMP)-dependent protein kinase (PKA) is localized through interaction of the regulatory (R) subunit dimer with A-kinase-anchoring proteins (AKAPs). We now report the solution structure of the type II PKA R-subunit fragment RIIalpha(1-44), which encompasses both the AKAP-binding and dimerization interfaces. This structure incorporates an X-type four-helix bundle dimerization motif with an extended hydrophobic face that is necessary for high-affinity AKAP binding. NMR data on the complex between RIIalpha(1-44) and an AKAP fragment reveals extensive contacts between the two proteins. Interestingly, this same dimerization motif is present in other signaling molecules, the S100 family. Therefore, the X-type four-helix bundle may represent a conserved fold for protein-protein interactions in signal transduction.

The third component of complement, C3, plays a central role in activation of the classical, alternative, and lectin pathways of complement activation. Recently, we have identified a 13-residue cyclic peptide (named Compstatin) that specifically binds to C3 and inhibits complement activation. To investigate the topology and the contribution of each critical residue to the binding of Compstatin to C3, we have now determined the solution structure using 2D NMR techniques; we have also synthesized substitution analogues and used these to study the structure-function relationships involved. Finally, we have generated an ensemble of a family of solution structures of the peptide with a hybrid distance geometry-restrained simulated-annealing methodology, using distance, dihedral angle, and 3J(NH-Halpha)-coupling constant restraints. The Compstatin structure contained a type I beta-turn comprising the segment Gln5-Asp6-Trp7-Gly8. Preference for packing of the hydrophobic side chains of Val3, Val4, and Trp7 was observed. The generated structure was also analyzed for consistency using NMR parameters such as NOE connectivity patterns, 3J(NH-Halpha)-coupling constants, and chemical shifts. Analysis of Ala substitution analogues suggested that Val3, Gln5, Asp6, Trp7, and Gly8 contribute significantly to the inhibitory activity of the peptide. Substitution of Gly8 caused a 100-fold decrease in inhibitory potency. In contrast, substitution of Val4, His9, His10, and Arg11 resulted in minimal change in the activity. These findings indicate that specific side-chain interactions and the beta-turn are critical for preservation of the conformational stability of Compstatin and they might be significant for maintaining the functional activity of Compstatin.

Hydrogen/deuterium exchange rates for individual amide protons have been measured for the carbon monoxide complex of soybean leghemoglobin. Fast two-dimensional NOESY experiments were performed, with 5.2-min data-collection time for each spectrum, which made possible the measurement of NOE cross-peaks of relatively rapidly exchanging amide protons at early time points. Exchange rates were measured for 61 backbone amides, the protection factors were calculated to provide information on the packing and local stability of the protein. The data are consistent with the presence of transient cooperative local unfolding of helical segments. The B-, E-, G- and H-helices have extensive regions of slow-, medium- and fast-exchanging amide protons. For each of these helices, there is a progressive decrease in protection on moving from the helix center to the termini. This is consistent with a stable helix center, with dynamic fraying at the ends. Amide exchange from the A-helix and C-helix is rapid except in small local regions. The F-helix, which is located on the proximal side of the heme pocket and is well formed in solution as demonstrated by characteristic medium range NOE connectivities [Morikis, D. Lepre, C.A. & Wright, P.E. (1994) Eur. J. Biochem. 219, 611-626], exhibits fast exchange for all amide protons. The implied flexibility and low stability of the F-helix may be functionally important in facilitating movement of the helix upon ligand binding. Fast exchange has also been observed for all amide protons in the CE-loop and in turns, as expected for flexible or solvent exposed regions. A strong tertiary contact has been established between the A-, G- and H-helices by the presence of a slowly exchanging indole N epsilon H of Trp129.

Homonuclear two-dimensional and three-dimensional 1H-NMR spectroscopy has been utilized to study the 15.9-kDa protein soybean leghemoglobin. NMR experiments were performed on the diamagnetic carbon monoxide complex at two temperatures and two pH values. Sequence-specific assignments have been made for 94% of the backbone and approximately 70% of the expected side-chain resonances. The secondary structure of leghemoglobin in solution has been determined on the basis of NOE connectivity patterns, hydrogen exchange and chemical-shift analyses. Leghemoglobin consists of seven helices and, unlike mammalian myoglobins, is missing the D helix. Instead an extended loop, the CE loop, is observed which might have importance for ligand entry into and exit from the protein interior. The hydrogen exchange behavior for the F helix and at the beginning of the A helix suggests different dynamic stability compared to other helical regions in leghemoglobin. Population of a second protein conformation, in which there is perturbation at the A-G-H helix interface, is observed at low pH.

A monomeric model peptide, acetyl-WEAQAREALAKEAAARA-amide, has been structurally characterized using the complementary techniques of 1H 2D NMR and circular dichroism. Temperature-dependent CD measurements are consistent with a two-state helix/coil transition model and indicate a 65% contribution of helical conformers at 5 degrees C. Homonuclear 2D NMR experiments allowed the assignment of all proton resonances. The analysis of NOE-type cross-relaxation data established a large number of specific short- and medium-range NOE connectivities throughout the peptide, confirming the highly helical character of the peptide. However, the observation of long-range NOEs between the methyl protons of leucine-9 and backbone and side-chain protons of amino acids located at the N-terminus, as well as other unusual NOEs, unambiguously reflects the existence of significantly populated nonhelical structured conformers, indicating a multiconformational equilibrium. Implications of these observations with regard to secondary structure quantitation and current method limitations are discussed.

We explore heme structure and ligation subsequent to a low-pH conformational transition in sperm whale myoglobin. Below pH 4.0, the iron-histidine bond breaks in metMb and deoxyMb. In MbCO, the majority of the iron-histidine bonds remain intact down to pH 2.6; however, the observation of a weak Fe-CO mode at 526 cm-1 indicates that a small fraction of the sample has the histidine replaced by a weak ligand, possibly water. The existence of a sterically hindered CO subpopulation in MbCO and the continued association of the four-coordinate heme with the protein in deoxyMb suggest that the heme pocket remains at least partially intact in the acid-induced conformation. The global pH-dependent conformational change described here is clearly distinguished from the local "closed" to "open" transition described previously in MbCO [Morikis et al. (1989) Biochemistry 28, 4791-4800]. Further observations of the four-coordinate heme state yield insights on the mechanism of heme photoreduction and the assignment of the 760-nm band in deoxyMb.

Three mutant proteins of sperm whale myoglobin (Mb) that exhibit altered axial ligations were constructed by site-directed mutagenesis of a synthetic gene for sperm whale myoglobin. Substitution of distal pocket residues, histidine E7 and valine E11, with tyrosine and glutamic acid generated His(E7)Tyr Mb and Val(E11)Glu Mb. The normal axial ligand residue, histidine F8, was also replaced with tyrosine, resulting in His(F8)Tyr Mb. These proteins are analogous in their substitutions to the naturally occurring hemoglobin M mutants (HbM). Tyrosine coordination to the ferric heme iron of His(E7)Tyr Mb and His(F8)Tyr Mb is suggested by optical absorption and EPR spectra and is verified by similarities to resonance Raman spectral bands assigned for iron-tyrosine proteins. His(E7)Tyr Mb is high-spin, six-coordinate with the ferric heme iron coordinated to the distal tyrosine and the proximal histidine, resembling Hb M Saskatoon [His(beta E7)Tyr], while the ferrous iron of this Mb mutant is high-spin, five-coordinate with ligation provided by the proximal histidine. His(F8)Tyr Mb is high-spin, five-coordinate in both the oxidized and reduced states, with the ferric heme iron liganded to the proximal tyrosine, resembling Hb M Iwate [His(alpha F8)Tyr] and Hb M Hyde Park [His(beta F8)Tyr]. Val(E11)Glu Mb is high-spin, six-coordinate with the ferric heme iron liganded to the F8 histidine. Glutamate coordination to the ferric iron of this mutant is strongly suggested by the optical and EPR spectral features, which are consistent with those observed for Hb M Milwaukee [Val(beta E11)Glu]. The ferrous iron of Val(E11)Glu Mb exhibits a five-coordinate structure with the F8 histidine-iron bond intact.(ABSTRACT TRUNCATED AT 250 WORDS)

The resonance Raman spectra of met-, deoxy-, and (carbonmonoxy)myoglobin (MbCO) are studied as a function of amino acid replacement at the distal histidine-E7 position. The synthetic wild type is found to be spectroscopically identical with the native material. The methionine and glycine replacements do not affect the met or deoxy spectra but do lead to distinct changes in the nu Fe-CO region of the MbCO spectrum. The native MbCO displays a pH-dependent population redistribution of the nu Fe-CO modes, while the analogous population in the mutant systems is found to be pH independent. This indicates that histidine-E7 is the titratable group in native MbCO. Moreover, the pH dependence of the population dynamics is found to be inconsistent with a simple two-state Henderson-Hasselbalch analysis. Instead, we suggest a four-state model involving the coupling of histidine protonation and conformational change. Within this model, the pK of the distal histidine is found to be 6.0 in the "open" configuration and 3.8 in the "closed" conformation. This corresponds to a 3 kcal/mol destabilization of the positively charged distal histidine within the hydrophobic pocket and suggests how protonation can lead to a larger population of the "open" conformation. At pH 7, the pocket is found to be "open" approximately 3% of the time. Further work, involving both IR and Raman measurements, allows the electron-nuclear coupling strengths of the various nu Fe-CO and nu C-O Raman modes to be determined. The slowly rebinding conformational state, corresponding to nu Fe-CO = 518 cm-1 (nu C-O = 1932 cm-1), displays unusually weak coupling of the Fe-CO mode to the Soret transition. Studies of the nu Fe-CO region as a function of temperature reveal that the equilibria between the conformational states are quenched in both the native and glycine mutant below the freezing point of the solvent. Unusual line narrowing of the nu Fe-CO modes at the phase transition is also observed in all samples studied. This line narrowing stands in marked contrast to the other heme Raman modes and suggests that Fe-CO librational motion and/or distal pocket vibrational (or conformational) excitations are involved in the line broadening at room temperature.