Role of Heme Pocket Water in Allosteric Regulation of Ligand Reactivity in Human Hemoglobin.

Water molecules can enter the heme pockets of unliganded myoglobins and hemoglobins, hydrogen bond with the distal histidine, and introduce steric barriers to ligand binding. The spectrokinetics of photodissociated CO complexes of human hemoglobin and its isolated α and ß chains were analyzed for the effect of heme hydration on ligand rebinding. A strong coupling was observed between heme hydration and quaternary state. This coupling may contribute significantly to the 20-60-fold difference between the R- and T-state bimolecular CO binding rate constants and thus to the modulation of ligand reactivity that is the hallmark of hemoglobin allostery. Heme hydration proceeded over the course of several kinetic phases in the tetramer, including the R to T quaternary transition. An initial 150 ns hydration phase increased the R-state distal pocket water occupancy, nw(R), to a level similar to that of the isolated α (∼60%) and ß (∼10%) chains, resulting in a modest barrier to ligand binding. A subsequent phase, concurrent with the first step of the R → T transition, further increased the level of heme hydration, increasing the barrier. The final phase, concurrent with the final step of the allosteric transition, brought the water occupancy of the T-state tetramer, nw(T), even higher and close to full occupancy in both the α and ß subunits (∼90%). This hydration level could present an even larger barrier to ligand binding and contribute significantly to the lower iron reactivity of the T state toward CO.

Probing kinetic mechanisms of protein function and folding with time-resolved natural and magnetic chiroptical spectroscopies.

Recent and ongoing developments in time-resolved spectroscopy have made it possible to monitor circular dichroism, magnetic circular dichroism, optical rotatory dispersion, and magnetic optical rotatory dispersion with nanosecond time resolution. These techniques have been applied to determine structural changes associated with the function of several proteins as well as to determine the nature of early events in protein folding. These studies have required new approaches in triggering protein reactions as well as the development of time-resolved techniques for polarization spectroscopies with sufficient time resolution and sensitivity to probe protein structural changes.

Kinetic spectroscopy of heme hydration and ligand binding in myoglobin and isolated hemoglobin chains: an optical window into heme pocket water dynamics.

The entry of a water molecule into the distal heme pocket of pentacoordinate heme proteins such as myoglobin and the alpha,beta chains of hemoglobin can be detected by time-resolved spectroscopy in the heme visible bands after photolysis of the CO complex. Reviewing the evidence from spectrokinetic studies of Mb variants, we find that this optical method measures the occupancy of non(heme)coordinated water in the distal pocket, n(w), with high fidelity. This evidence further suggests that perturbation of the kinetic barrier presented by distal pocket water is often the dominant mechanism by which active site mutations affect the bimolecular rate constant for CO binding. Water entry into the heme pockets of isolated hemoglobin subunits was detected by optical methods. Internal hydration is higher in the native alpha chains than in the beta chains, in agreement with previous crystallographic results for the subunits within Hb tetramers. The kinetic parameters obtained from modeling of the water entry and ligand rebinding in Mb mutants and native Hb chains are consistent with an inverse dependence of the bimolecular association rate constant on the water occupancy factor. This correlation suggests that water and ligand mutually exclude one another from the distal pockets of both types of hemoglobin chains and myoglobin.

Nanosecond time-resolved polarization spectroscopies: tools for probing protein reaction mechanisms.

Polarization methods, introduced in the 1800s, offered one of the earliest ways to examine protein structure. Since then, many other structure-sensitive probes have been developed, but circular dichroism (CD) remains a powerful technique because of its versatility and the specificity of protein structural information that can be explored. With improvements in time resolution, from millisecond to picosecond CD measurements, it has proven to be an important tool for studying the mechanism of folding and function in many biomolecules. For example, nanosecond time-resolved CD (TRCD) studies of the sub-microsecond events of reduced cytochrome c folding have provided direct experimental evidence of kinetic heterogeneity, which is an inherent property of the diffusional nature of early folding dynamics on the energy landscape. In addition, TRCD has been applied to the study of many biochemical processes, such as ligand rebinding in hemoglobin and myoglobin and signaling state formation in photoactive yellow protein and prototropin 1 LOV2. The basic approach to TRCD has also been extended to include a repertoire of nanosecond polarization spectroscopies: optical rotatory dispersion (ORD), magnetic CD and ORD, and linear dichroism. This article will discuss the details of the polarization methods used in this laboratory, as well as the coupling of time-resolved ORD with the temperature-jump trigger so that protein folding can be studied in a larger number of proteins.

Optical detection of disordered water within a protein cavity.

Internal water molecules are important to protein structure and function, but positional disorder and low occupancies can obscure their detection by X-ray crystallography. Here, we show that water can be detected within the distal cavities of myoglobin mutants by subtle changes in the absorbance spectrum of pentacoordinate heme, even when the presence of solvent is not readily observed in the corresponding crystal structures. A well-defined, noncoordinated water molecule hydrogen bonded to the distal histidine (His64) is seen within the distal heme pocket in the crystal structure of wild type (wt) deoxymyoglobin. Displacement of this water decreases the rate of ligand entry into wt Mb, and we have shown previously that the entry of this water is readily detected optically after laser photolysis of MbCO complexes. However, for L29F and V68L Mb no discrete positions for solvent molecules are seen in the electron density maps of the crystal structures even though His64 is still present and slow rates of ligand binding indicative of internal water are observed. In contrast, time-resolved perturbations of the visible absorption bands of L29F and V68L deoxyMb generated after laser photolysis detect the entry and significant occupancy of water within the distal pockets of these variants. Thus, the spectral perturbation of pentacoordinate heme offers a potentially robust system for measuring nonspecific hydration of the active sites of heme proteins.

Probing early events in ferrous cytochrome c folding with time-resolved natural and magnetic circular dichroism spectroscopies.

In a 1998 collaboration with Tony Fink, we coupled nanosecond circular dichroism methods (TRCD) with a CO-photolysis system for quickly triggering folding in cytochrome c (cyt c) in order to make the first time-resolved far-UV CD measurement of early secondary structure formation in a protein. The small signal observed in that initial study, approximately 10% of native helicity, became the seed for increasingly robust results from subsequent studies bringing additional natural and magnetic circular polarization dichroism and optical rotatory dispersion detection methods (e.g., TRORD, TRMCD, and TRMORD), coupled to fast photolysis and photoreduction triggers, to the study of early folding events. Nanosecond polarization methods are reviewed here in the context of the range of initiation methods and structure-sensitive probes currently available for fast folding studies. We also review the impact of experimental results from fast polarization studies on questions in folding dynamics such as the possibility of multiple folding pathways implied by energy landscape models, the sequence dependence of ultrafast helix formation, and the simultaneity of chain collapse and secondary structure formation implicit in molten globule models for kinetic folding intermediates.

Early events, kinetic intermediates and the mechanism of protein folding in cytochrome C.

Kinetic studies of the early events in cytochrome c folding are reviewed with a focus on the evidence for folding intermediates on the submillisecond timescale. Evidence from time-resolved absorption, circular dichroism, magnetic circular dichroism, fluorescence energy and electron transfer, small-angle X-ray scattering and amide hydrogen exchange studies on the t < or = 1 ms timescale reveals a picture of cytochrome c folding that starts with the approximately 1-micros conformational diffusion dynamics of the unfolded chains. A fractional population of the unfolded chains collapses on the 1 - 100 micros timescale to a compact intermediate I(C) containing some native-like secondary structure. Although the existence and nature of I(C) as a discrete folding intermediate remains controversial, there is extensive high time-resolution kinetic evidence for the rapid formation of I(C) as a true intermediate, i.e., a metastable state separated from the unfolded state by a discrete free energy barrier. Final folding to the native state takes place on millisecond and longer timescales, depending on the presence of kinetic traps such as heme misligation and proline mis-isomerization. The high folding rates observed in equilibrium molten globule models suggest that I(C) may be a productive folding intermediate. Whether it is an obligatory step on the pathway to the high free energy barrier associated with millisecond timescale folding to the native state, however, remains to be determined.

The pH dependence of heme pocket hydration and ligand rebinding kinetics in photodissociated carbonmonoxymyoglobin.

We monitored the occupancy of a functionally important non-coordinated water molecule in the distal heme pocket of sperm whale myoglobin over the pH range 4.3-9.4. Water occupancy was assessed by using time-resolved spectroscopy to detect the perturbation of the heme visible band absorption spectrum caused by water entry after CO photodissociation ( Goldbeck, R. A., Bhaskaran, S., Ortega, C., Mendoza, J. L., Olson, J. S., Soman, J., Kliger, D. S., and Esquerra, R. M. (2006) Proc. Natl. Acad. Sci. U. S. A. 103, 1254-1259 ). We found that the water occupancy observed during the time interval between ligand photolysis and diffusive recombination decreased by nearly 20% as the pH was lowered below 6. This decrease accounted for most of the concomitant increase in the observed CO bimolecular recombination rate constant, as the lower water occupancy presented a smaller kinetic barrier to CO entry into the pocket at lower pH. These results were consistent with a model in which the distal histidine, which stabilizes the water molecule within the distal pocket by accepting a hydrogen bond, tends to swing out of the pocket upon protonation and destabilize the water occupancy at low pH. Extrapolation of this model to lower pH suggests that the additional increase in ligand association rate constant observed previously in stopped-flow studies at pH 3 may also be due in part to reduced distal water occupancy concomitant with further His64 protonation and coupled protein conformational change.

Running, swimming and diving modifies neuroprotecting globins in the mammalian brain.

The vulnerability of the human brain to injury following just a few minutes of oxygen deprivation with submergence contrasts markedly with diving mammals, such as Weddell seals (Leptonychotes weddellii), which can remain underwater for more than 90 min while exhibiting no neurological or behavioural impairment. This response occurs despite exposure to blood oxygen levels concomitant with human unconsciousness. To determine whether such aquatic lifestyles result in unique adaptations for avoiding ischaemic-hypoxic neural damage, we measured the presence of circulating (haemoglobin) and resident (neuroglobin and cytoglobin) oxygen-carrying globins in the cerebral cortex of 16 mammalian species considered terrestrial, swimming or diving specialists. Here we report a striking difference in globin levels depending on activity lifestyle. A nearly 9.5-fold range in haemoglobin concentration (0.17-1.62 g Hb 100 g brain wet wt(-1)) occurred between terrestrial and deep-diving mammals; a threefold range in resident globins was evident between terrestrial and swimming specialists. Together, these two globin groups provide complementary mechanisms for facilitating oxygen transfer into neural tissues and the potential for protection against reactive oxygen and nitrogen groups. This enables marine mammals to maintain sensory and locomotor neural functions during prolonged submergence, and suggests new avenues for averting oxygen-mediated neural injury in the mammalian brain.

Measurement of the heme affinity for yeast dap1p, and its importance in cellular function.

Current studies on the Saccharomyces cerevisiae protein Dap1p have demonstrated a heme-related function within the ergosterol biosynthetic pathway. Here we present data to further the understanding of the role of heme in the proper biological functioning of Dap1p in cellular processes. First, we examined the role of Dap1p in stabilizing the P450 enzyme, Erg11p, a key regulatory protein in ergosterol biosynthesis. Our data indicate that the absence of Dap1p does not affect Erg11p mRNA, protein expression levels, or the protein degradation rates in S. Cerevisaie. Second, in order to probe the role of heme in the biological functioning of Dap1p, we measured ferric and ferrous heme binding affinities for Dap1p and the mutant Dap1pY138F, as well as equilibrium midpoint reduction potentials of the Fe(III)/Fe(II) couples. Our results show that both wild-type and mutant proteins bind heme in a 1:1 fashion, possessing tight ferric heme affinities, KD values of 400 pM and 200 nM, respectively, but exhibiting weak ferrous affinities, 2 and 10 microM, respectively. Additionally, the measured reduction potential of Dap1p, which was found to be -307 mV, is similar to that of other monotyrosinate hemoproteins. Although previous reports show the weaker affinity of Dap1pY138F for ferric heme lowers the production of ergosterol with respect to wild-type Dap1p in S. pombe, we find that Dap1pY138F expression is still sufficient to rescue the growth sensitivity of dap1Delta to fluconazole and methyl methanesulfonate in S. cerevisiae. Various interpretations of these results are discussed with respect to Dap1p function in the cell.

Non-native heme-histidine ligation promotes microsecond time scale secondary structure formation in reduced horse heart cytochrome c.

Previous far-UV time-resolved optical rotatory dispersion (TRORD) studies of the sub-millisecond (burst) phase of secondary structure formation in horse and tuna cytochromes c after photoreduction in denaturant suggested that the non-native His18-Fe-His33 heme ligation dominant in the unfolded horse protein facilitated this fast folding better than did the His18-Fe-His26 coordination dominant in tuna [Chen, E., Goldbeck, R.A., and Kliger, D.S. (2003) J. Phys. Chem. A 107, 8149-8155; Chen, E., Goldbeck, R.A., and Kliger, D.S. (2004) J. Am. Chem. Soc. 126, 11175-11181]. Whether His18-Fe-His33 coordination actually facilitates fast secondary structure formation or just slows folding less than His18-Fe-His26 coordination is probed by examining the double histidine mutant H26QH33N of horse heart cytochrome c. The fast folding phase is absent in H26QH33N, indicating that His18-Fe-His33 misligation does promote fast secondary structure formation, as does His18-Fe-His26 to a lesser extent. His33 may be better able to facilitate folding because it is not as constrained by hydrogen bonding interactions in the denatured state as is His26. A greater flexibility, not only because of weakened or disrupted Van der Waals interactions in the presence of guanidine hydrochloride (GuHCl) but also because of its position relative to His18, may allow His33 to ligate to the heme group more easily than His26. These results are discussed along with the results of far-UV CD and Soret and visible region MCD measurements, which were performed to probe heme ligation in H26QH33N and to understand how GuHCl affects its folding stability and cooperativity.

Conformational equilibration time of unfolded protein chains and the folding speed limit.

The speed with which the conformers of unfolded protein chains interconvert is a fundamental question in the study of protein folding. Kinetic evidence is presented here for the time constant for interconversion of disparate unfolded chain conformations of a small globular protein, cytochrome c, in the presence of guanidine hydrochloride denaturant. The axial binding reactions of histidine and methionine residues with the Fe(II) heme cofactor were monitored with time-resolved magnetic circular dichroism spectroscopy after photodissociation of the CO complexes of unfolded protein obtained from horse and tuna and from several histidine mutants of the horse protein. A kinetic model fitting both the reaction rate constants and spectra of the intermediates was used to obtain a quantitative estimate of the conformational diffusion time. The latter parameter was approximated as a first-order time constant for exchange between conformational subensembles presenting either a methionine or a histidine residue to the heme iron for facile binding. The mean diffusional time constant of the wild type and variants was 3 +/- 2 mus, close to the folding "speed limit". The implications of the relatively rapid conformational equilibration time observed are discussed in terms of the energy landscape and classical pathway time regimes of folding, for which the conformational diffusion time can be considered a pivot point.

Water and ligand entry in myoglobin: assessing the speed and extent of heme pocket hydration after CO photodissociation.

A previously undescribed spectrokinetic assay for the entry of water into the distal heme pocket of wild-type and mutant myoglobins is presented. Nanosecond photolysis difference spectra were measured in the visible bands of sperm whale myoglobin as a function of distal pocket mutation and temperature. A small blue shift in the 560-nm deoxy absorption peak marked water entry several hundred nanoseconds after CO photodissociation. The observed rate suggests that water entry is rate-limited by the escape of internal dissociated CO. The heme pocket hydration and geminate recombination yields were found to be the primary factors controlling the overall bimolecular association rate constants for CO binding to the mutants studied. The kinetic analysis provides estimates of 84%, 60%, 40%, 0%, and 99% for the steady-state hydrations of wild-type, H64Q, H64A, H64L, and V68F deoxymyoglobin, respectively. The second-order rate constants for CO and H(2)O entry into the empty distal pocket of myoglobin are markedly different, 8 x 10(7) and 2 x 10(5) M(-1).s(-1), respectively, suggesting that hydrophobic partitioning of the apolar gas from the aqueous phase into the relatively apolar protein interior lowers the free energy barrier for CO entry.

Spectroscopic and biochemical characterization of heme binding to yeast Dap1p and mouse PGRMC1p.

Yeast damage-associated response protein (Dap1p) and mouse progesterone receptor membrane component-1 protein (mPGRMC1p) belong to a highly conserved class of putative membrane-associated progesterone binding proteins (MAPR), with Dap1p and inner zone antigen (IZA), the rat homologue of mPGRMC1p, recently being reported to bind heme. While primary structure analysis reveals similarities to the cytochrome b(5) motif, neither of the two axial histidines responsible for ligation to the heme is present in any of the MAPR proteins. In this paper, EPR, MCD, CD, UV-vis, and general biochemical methods have been used to characterize the nature of heme binding in both Dap1p and a His-tagged, membrane anchor-truncated mPGRMC1p. As isolated, Dap1p is a tetramer which can be converted to a dimer upon addition of 150 mM salt. The heme is noncovalently attached, with a maximal, in vitro, heme loading of approximately 30%, for both proteins. CD and fluorescence spectroscopies indicate a well-ordered structure, suggesting the low level of heme loading is probably not due to improperly folded protein. EPR confirmed a five-coordinate, high-spin, ferric resting state for both proteins, indicating one axial amino acid ligand, in contrast to the six-coordinate, low-spin, ferric state of cytochrome b(5). The MCD spectrum confirmed this conclusion for Dap1p and indicated the axial ligand is most likely a tyrosine and not a histidine, or a cysteine; however, an aspartic acid residue could not be conclusively ruled out. Potential axial ligands, which are conserved in all MAPRs, were mutated (Y78F, D118A, and Y138F) and purified to homogeneity. The Y78F and D118A mutants were found to bind heme; however, Y138F did not. This result is consistent with the MCD data and indicates that Tyr138 is most likely the axial ligand to the heme in Dap1p.

The molecular code for hemoglobin allostery revealed by linking the thermodynamics and kinetics of quaternary structural change. 1. Microstate linear free energy relations.

A novel model linking the thermodynamics and kinetics of hemoglobin's allosteric (R --> T) and ligand binding reactions is applied to photolysis data for human HbCO. To describe hemoglobin's kinetics at the microscopic level of structural transitions and ligand-binding events for individual [ij]-ligation microstates ((ij)R --> (ij)T, (ij)R + CO --> ((i)(+1))(k)R, and (ij)T + CO --> ((i)(+1))(k)T), the model calculates activation energies, (ij)DeltaG(++), from previously measured cooperative free energies of the equilibrium microstates (Huang, Y., and Ackers, G. K. (1996) Biochemistry 35, 704-718) by using linear free energy relations ((ij)DeltaG(++) - (01)DeltaG(++) = alpha[(ij)DeltaG - (01)DeltaG], where the parameter alpha, describing the variation of activation energy with reaction energy perturbation, can depend on the natures of both the reaction and the perturbation). The alpha value measured here for the allosteric dynamics, 0.21 +/- 0.03, corresponds closely to values observed previously, strongly suggesting that the thermodynamic microstate energies directly underlie the allosteric kinetics (as opposed to the alpha((ij)DeltaG(RT)) serving merely as arbitrary fitting parameters). Besides systematizing the study of hemoglobin kinetics, the utility of the microstate linear free energy model lies in the ability to test microscopic aspects of allosteric dynamics such as the "symmetry rule" for quaternary change deduced previously from thermodynamic evidence (Ackers, G. K., et al. (1992) Science 255, 54-63). Reflecting a remarkably detailed correspondence between thermodynamics and kinetics, we find that a kinetic model that includes the large free energy splitting between doubly ligated T microstates implied by the symmetry rule fits the data significantly better than one that does not.

The molecular code for hemoglobin allostery revealed by linking the thermodynamics and kinetics of quaternary structural change. 2. Cooperative free energies of (alphaFeCObetaFe)2 and (alphaFebetaFeCO)2 T-state tetramers.

Ligand photodissociation experiments are used to measure the prephotolysis equilibria between doubly liganded R and T quaternary conformers of the symmetric Fe-Co HbCO hybrids, (alpha(FeCO)beta(Co))(2) and (alpha(Co)beta(FeCO))(2). The free energies obtained from these data are used to calculate the cooperative free energies of the (alpha(FeCO)beta(Fe))(2) and (alpha(Fe)beta(FeCO))(2) intermediate CO-ligation states of normal hemoglobin in the T conformation, quantities important to the evaluation of current models of cooperativity. The symmetry rule model, incorporating sequential cooperativity of T-state ligand binding within an alphabeta dimer in addition to the traditional two-state cooperativity of the tetramer, predicts a larger free energy penalty for disturbing both dimers in a doubly liganded T tetramer than would be expected in the two-state model as currently formulated. (Cooperative energy penalties are simply proportional to the number of tetramer-bound ligands in the traditional two-state model.) The value found here for the energies of doubly liganded T microstates in which both dimers are perturbed, 7.9 +/- 0.3 kcal/mol, is consistent with the symmetry rule model but significantly higher than that expected (5-6 kcal/mol) in the two-state model of cooperativity.

The earliest events in protein folding: a structural requirement for ultrafast folding in cytochrome C.

The folding dynamics of reduced cytochrome c (redcyt c) obtained from tuna heart, which contains a tryptophan residue at the site occupied by His33 in horse heart cytochrome c, was studied using nanosecond time-resolved optical rotatory dispersion spectroscopy. As observed previously for horse heart redcyt c, two time regimes were observed for secondary structure formation in tuna redcyt c: a fast (microseconds) and a slow (milliseconds) phase. However, the fast phase of tuna redcyt c folding was much slower and smaller in amplitude than the same phase in horse. The differences in the fast folding phases suggest that for horse heart redcyt c, the conformers that undergo the fastest observed folding have the His18-Fe-His33 heme configuration, which appears to be necessary, but not sufficient, to poise an unfolded chain conformation for fastest folding in redcyt c.

Hydrogen bonding to Trp beta37 is the first step in a compound pathway for hemoglobin allostery.

Human hemoglobin is widely thought to change from the R to the T quaternary structure in a single rate process requiring tens of microseconds. Here we present kinetic evidence that the R --> T allosteric pathway in hemoglobin requires more than one step. We use magnetic circular dichroism (MCD) spectroscopy of the aromatic amino acid bands to show that formation of a tryptophan-aspartate hydrogen bond in the hinge region of the dimer-dimer interface is part of an obligatory R --> T step proceeding more than a factor of 10 faster than the kinetic step previously identified in heme-band absorption studies.