New compstatin variants through two de novo protein design frameworks.

Two de novo protein design frameworks are applied to the discovery of new compstatin variants. One is based on sequence selection and fold specificity, whereas the other approach is based on sequence selection and approximate binding affinity calculations. The proposed frameworks were applied to a complex of C3c with compstatin variant E1 and new variants with improved binding affinities are predicted and experimentally validated. The computational studies elucidated key positions in the sequence of compstatin that greatly affect the binding affinity. Positions 4 and 13 were found to favor Trp, whereas positions 1, 9, and 10 are dominated by Asn, and position 11 consists mainly of Gln. A structural analysis of the C3c-bound peptide analogs is presented.

From atoms to systems: a cross-disciplinary approach to complement-mediated functions.

With an ever-increasing wealth of information made available to researchers from expanding genomic sequence and protein structure databases, traditional experimentation and research are being drastically revisited. The unidirectional study of single molecules and pathways is being replaced by a combinatorial and cross-disciplinary platform that investigates interactive biological systems and dynamic networks. The complement system constitutes an ideal paradigm of how this concept is being applied in the field of contemporary immunology. Our laboratory has adopted such a cross-disciplinary approach in elucidating key aspects of complement functions and determining the role of several complement proteins in both inflammatory and developmental processes. Here we discuss recent findings pertaining to the rational development of complement inhibitors, our studies on protein-protein interactions and our progress in the study of viral immune evasion and complement evolution. Furthermore, we present recent studies implicating complement components in complex developmental processes, such as organ regeneration, hematopoietic development, and stem cell engraftment.

Improvement of the anti-C3 activity of compstatin using rational and combinatorial approaches.

Compstatin is a 13-residue cyclic peptide that has the ability to inhibit the cleavage of C3 to C3a and C3b. The effects of targeting C3 cleavage are threefold, and result in hindrance of: (i) the generation of the pro-inflammatory peptide C3a, (ii) the generation of opsonin C3b (or its fragment C3d), and (iii) further complement activation of the common pathway (beyond C3) with the end result of the generation of the membrane attack complex. We will report on our progress on: (i) rational design of more active compstatin analogues based on the three-dimensional structure of compstatin, (ii) experimental combinatorial design based on the generation of a phage-displayed peptide library partially randomized with the implementation of structure-induced restraints, and (iii) theoretical combinatorial design based on a novel computational optimization method, structure-induced restraints and flexible structural templates. All three approaches have resulted in analogues with improved activities. Currently, the lead analogue has the sequence acetyl-I[CVYQDWGAHRC]T-NH(2) (where the brackets denote cyclization), and is 16-fold more active than the parent peptide. We will also report on our progress towards understanding the dynamic character of compstatin using molecular dynamics simulations. The identification of an ensemble of interconverting conformers of compstatin with variable populations is a first step towards the incorporation of dynamic elements in the design of new analogues using dynamics-activity relationships in addition to structure-activity relationships.

Structural aspects and design of low-molecular-mass complement inhibitors.

We present a mini-review on the structure-based design of three promising complement inhibitors. Firstly, we review compstatin, a 13-residue cyclic peptide that binds to C3 and inhibits the cleavage of C3 to C3a and C3b. Secondly, we review a six-residue cyclic peptide that binds to C5aR and antagonizes the binding of C5a to its receptor C5aR. Finally, we review three small molecules that bind to Factor D and inhibit the enzymic action of Factor D, during which Factor D proteolytically cleaves Factor B in complex with C3 or C3b.

Native-state conformational dynamics of GART: a regulatory pH-dependent coil-helix transition examined by electrostatic calculations.

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.

Proton transfer dynamics of GART: the pH-dependent catalytic mechanism examined by electrostatic calculations.

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.

A novel mechanism of PKA anchoring revealed by solution structures of anchoring complexes.

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.

Binding kinetics, structure-activity relationship, and biotransformation of the complement inhibitor compstatin.

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.

The molecular basis for protein kinase A anchoring revealed by solution NMR.

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.

Solution structure of Compstatin, a potent complement inhibitor.

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 exchange in the carbon monoxide complex of soybean leghemoglobin.

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.

Hydration of the partially folded peptide RN-24 studied by multidimensional NMR.

Peptide-water interactions of a ribonuclease C-peptide analogue, RN-24 (Suc-AETAAAKFLRAHANH2), which exhibits significant helicity, have been studied in solution using homonuclear 2D and 3D NMR cross-relaxation experiments. Dipolar peptide proton-water proton interactions are indicated by a large number of NOESY-type cross peaks at the H2O resonance frequency, most of them with opposite sign relative to the diagonal. Some cross peaks arise from intrapeptide cross relaxation to labile protons of histidine, threonine, lysine and arginine side chains. The observed peptide-water interactions are rather uniformly distributed, involving peptide backbone and side chains equally. The data are consistent with rapid fluctuations of the conformational ensemble and the absence of peptide regions that are highly shielded from bulk solvent, even in a peptide that exhibits high propensities for formation of helical secondary structure.

1H resonance assignments and secondary structure of the carbon monoxide complex of soybean leghemoglobin determined by homonuclear two-dimensional and three-dimensional NMR spectroscopy.

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.

NMR evidence for multiple conformations in a highly helical model peptide.

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.

Low pH myoglobin photoproducts.

Recently, there has been interest in determining the conditions under which the iron-histidine bond ruptures in myoglobin at low pH, so that the effect of proximal heme ligation can be studied. A 220-cm-1 Raman mode, assigned to iron-histidine stretching, is clearly visible after photolysis of aqueous MbCO samples below pH4 at room temperature (Sage et al. Biochemistry. 30:1237-1247). In contrast, Iben et al. (Biophys. J. 59:908-919) do not observe this mode upon photolysis of a pH3 MbCO sample in a glycerol/water glass at low temperature. In order to account for both the low temperature and the room temperature experiments, Iben et al. suggest a scheme involving an unusual protonation state of the proximal histidine. Here, we discuss some inconsistencies in their explanation of the room temperature results and offer instead a simple modification of an earlier model. In addition, circular dichroism data are presented that indicate partial unfolding of MbCO in aqueous solution below pH4, and raise questions about the claim of Iben et al. that MbCO remains folded in 75% glycerol at pH3.

Conformational interconversion in protein crystals.

We present evidence that the structure of carbonmonoxy myoglobin crystals can be altered by lowering the pH. This structural change is monitored by the characteristic Fe-CO Raman modes at 508 and 491 cm-1 and is thought to involve a localized distal pocket transition from a "closed" conformation at pH 7 to a more "open" conformation at pH 4. These changes take place in the crystal without loss of intensity of a conformationally sensitive Raman mode at 252 cm-1 that signals a partial unfolding of the globin structure in solution. Quantitative studies, which monitor the open and closed populations as a function of laser photolysis, demonstrate that the interconversion rates (k+/-) in solution at 298 K are fast compared to the photolysis and CO entry rates (i.e. k+/- much greater than 10(3) s-1), while in frozen samples the interconversion is much slower than the experimental time scale (minutes). Since the open conformation is a minority species at pH 7, rapid exchange in aqueous solution is a necessary condition for this species to play a functional role. In the crystal, the interconversion rates are slowed compared to solution and begin to approach the photolysis rate (i.e. k+/- approximately 10(3) to 10(4) s-1). This indicates that the barriers for conformational exchange are increased in the crystal environment, compared to the solution, apparently due to the packing forces of the surrounding molecules. X-ray and neutron diffraction studies of MbCO crystals at high and low pH are needed to characterize the details of the structural changes and to test the hypothesis that closed and open distal pocket structures are associated with the 508 and 491 cm-1 Fe-CO modes.

Spectroscopic studies of myoglobin at low pH: heme structure and ligation.

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.

Alteration of sperm whale myoglobin heme axial ligation by site-directed mutagenesis.

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)

Resonance Raman studies of iron spin and axial coordination in distal pocket mutants of ferric myoglobin.

We have used resonance Raman spectroscopy to study 11 distal pocket mutants and the "wild type" and native ferric sperm whale myoglobin. The characteristic Raman core-size markers v4, v3, v2, and v10 are utilized to assign the spin and coordination state of each sample. It is demonstrated that replacements of the distal and proximal histidines can discriminate against H2O as a sixth ligand and favor a pentacoordinate Fe3+ atom. Soret absorption band blueshifts are correlated with the pentacoordinate heme environment. One E7 replacement (Arg) leads to an iron spin state change and produces a low spin species. The Glu and Ala mutations at position E11 leave the protein's spin and coordination unaltered. A laser-induced photoreduction effect is observed in all pentacoordinate mutants and seems to be correlated with the loss of the heme-bound water molecule.

Resonance raman investigations of site-directed mutants of myoglobin: effects of distal histidine replacement.

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.