Binding of carbon monoxide to isolated hemoglobin chains

Evolution of the concept of conformational dynamics of enzyme functions over half of a century: A personal view

To most physicists, it was always evident that conformational fluctuation is an inherent property of all molecules. Its existence in proteins was mentioned first by Linderström-Lang and Schellman in 1959 based on their hydrogen-deuterium exchange experiments. The "induced fit" mechanism to explain ligand-induced conformational changes was suggested by Koshland in 1958. Straub combined these two concepts in his "fluctuation fit" theory in 1964. The era of protein X-ray crystallography imposed a static view of protein structures. With proteins becoming accessible to NMR analysis, conformational dynamics could be mapped, and a new wave of dynamic interpretation of enzymatic catalysis and molecular recognition appeared. Energy landscapes, energy funnels, conformational selection, conformational distribution shifts are now frequent terms in interpreting biomolecular recognition and enzymatic catalysis. All these interpretations are based on the concept that evolution uses the conformational fluctuations of enzymes to develop efficient and dynamic catalytic machines. In a resurrection of the original "fluctuation fit" concept, it is generally recognized now that spatial and temporal events of catalysis are equally important to describe its mechanism. This special issue, dedicated to the memory of Henryk Eisenberg, prompted us to look back at the last 50 years of development of a concept that-like other important concepts-appeared, evolved and became accepted during the period covered by the scientific lifespan of Henryk.

Proteins and Glasses

The structure, the energy landscape, and the dynamics of proteins and glasses are similar. Both types of systems display characteristic nonexponential time dependencies of relaxation phenomena. Experiments suggest that both, proteins and glasses, are heterogeneous and that this fact causes the observed time dependence. This result is discussed in terms of the rough energy landscape characteristic of complex systems.

Room-temperature magnetic properties of oxy- and carbonmonoxyhemoglobin

The magnetic susceptibility and the density of human oxy-(HbO(2)) and carbonmonoxyhemoglobin (HbCO) solutions of various concentrations have been measured at room temperature, with pure water used as a calibrant. Solutions of unstripped and stripped HbO(2) at pH 7.2 in unbuffered water solvent were always found to be less diamagnetic than pure water, whereas solutions of HbCO in identical conditions were always found to be more diamagnetic than pure water. After correcting for concentration-dependent density changes and assuming the HbCO samples to be fully diamagnetic, the paramagnetic reduction of the diamagnetic susceptibility of HbO(2) corresponds to a molar susceptibility per heme (chi(M) (heme)) of 2460 +/- 600 x 10(-6) cgs/mol.

Probing static disorder in Arrhenius kinetics by single-molecule force spectroscopy

The widely used Arrhenius equation describes the kinetics of simple two-state reactions, with the implicit assumption of a single transition state with a well-defined activation energy barrier DeltaE, as the rate-limiting step. However, it has become increasingly clear that the saddle point of the free-energy surface in most reactions is populated by ensembles of conformations, leading to nonexponential kinetics. Here we present a theory that generalizes the Arrhenius equation to include static disorder of conformational degrees of freedom as a function of an external perturbation to fully account for a diverse set of transition states. The effect of a perturbation on static disorder is best examined at the single-molecule level. Here we use force-clamp spectroscopy to study the nonexponential kinetics of single ubiquitin proteins unfolding under force. We find that the measured variance in DeltaE shows both force-dependent and independent components, where the force-dependent component scales with F(2), in excellent agreement with our theory. Our study illustrates a novel adaptation of the classical Arrhenius equation that accounts for the microscopic origins of nonexponential kinetics, which are essential in understanding the rapidly growing body of single-molecule data.

Exogenous carbon monoxide protects the bystander Chinese hamster ovary cells in mixed coculture system after alpha-particle irradiation

In the present work, the inhibitory effect of carbon monoxide (CO), generated by tricarbonyldichlororuthenium (II) dimer [CO-releasing molecule (CORM-2)], on the toxicity of radiation-induced bystander effect (RIBE) after alpha-particle irradiation was studied in a mixed coculture system. CO (CORM-2) treatment showed a significant inhibitory effect to the formation of p53 binding protein 1 (BP1) and micronuclei (MN) induced by RIBE in a concentration-dependent manner, but in the directly irradiated cell population no distinct decreases of BP1 and MN formation were observed. In this mixed coculture system, nitric oxide (NO) or superoxide anion (O2(*-)) was also proved to mediate the transduction of RIBE by using a NO synthase inhibitor or NADPH-oxidase-specific inhibitor treatment. The elevated O2(*-) was attenuated by CO (CORM-2) treatment in the bystander cells as measured by hydroethidine staining and fluorescence assessment. The exogenous NO (sper) or O2(*-) (H2O2) was used to mimic NO/O(2)-mediated RIBE, and CO (CORM-2) treatment also showed a protective effect to cells against the toxicity of these exogenous factors. Considering the inhibitory effect of CO on RIBE and the wide use of CO in therapy of diseases, it is hoped that a low concentration of CO can protect normal tissues against RIBE during radiotherapy.

Slaving: solvent fluctuations dominate protein dynamics and functions

Protein motions are essential for function. Comparing protein processes with the dielectric fluctuations of the surrounding solvent shows that they fall into two classes: nonslaved and slaved. Nonslaved processes are independent of the solvent motions; their rates are determined by the protein conformation and vibrational dynamics. Slaved processes are tightly coupled to the solvent; their rates have approximately the same temperature dependence as the rate of the solvent fluctuations, but they are smaller. Because the temperature dependence is determined by the activation enthalpy, we propose that the solvent is responsible for the activation enthalpy, whereas the protein and the hydration shell control the activation entropy through the energy landscape. Bond formation is the prototype of nonslaved processes; opening and closing of channels are quintessential slaved motions. The prevalence of slaved motions highlights the importance of the environment in cells and membranes for the function of proteins.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Multiple geminate ligand recombinations in human hemoglobin

The geminate ligand recombination reactions of photolyzed carbonmonoxyhemoglobin were studied in a nanosecond double-excitation-pulse time-resolved absorption experiment. The second laser pulse, delayed by intervals as long as 400 ns after the first, provided a measure of the geminate kinetics by rephotolyzing ligands that have recombined during the delay time. The peak-to-trough magnitude of the Soret band photolysis difference spectrum measured as a function of the delay between excitation pulses showed that the room temperature kinetics of geminate recombination in adult human hemoglobin are best described by two exponential processes, with lifetimes of 36 and 162 ns. The relative amounts of bimolecular recombination to T- and R-state hemoglobins and the temperature dependence of the submicrosecond kinetics between 283 and 323 K are also consistent with biexponential kinetics for geminate recombination. These results are discussed in terms of two models: geminate recombination kinetics modulated by concurrent protein relaxation and heterogeneous kinetics arising from alpha and beta chain differences.

Dynamic properties of monomeric insect erythrocruorin III from Chironomus thummi-thummi: relationships between structural flexibility and functional complexity

We have investigated the kinetics of geminate carbon monoxide binding to the monomeric component III of Chironomus thummi-thummi erythrocruorin, a protein that undergoes pH-induced conformational changes linked to a pronounced Bohr effect. Measurements were performed from cryogenic temperatures to room temperature in 75% glycerol and either 0.1 M potassium phosphate (pH 7) or 0.1 potassium borate (pH 9) after nanosecond laser photolysis. The distributions of the low temperature activation enthalpy g(H) for geminate ligand binding derived from the kinetic traces are quite narrow and are influenced by temperature both below and above approximately 170 K, the glass transition temperature. The thermal evolution of the CO binding kinetics between approximately 50 K and approximately 170 K indicates the presence of some degree of structural relaxation, even in this temperature range. Above approximately 220 K the width of the g(H) progressively decreases, and at 280 K geminate CO binding becomes exponential in time. Based on a comparison with analogous investigations of the homodimeric hemoglobin from Scapharca inaequivalvis, we propose a link between dynamic properties and functional complexity.

Effects of the intramolecular disulfide bond on ligand binding dynamics in myoglobin

In order to investigate the effects of an intramolecular disulfide bond on protein structure and ligand binding dynamics in myoglobin, we prepared a mutant myoglobin having a disulfide bond at the EF corner by introducing two cysteine at the position of Ile 75 and Glu 85. On the basis of the spectral features of the mutant, the formation of the disulfide bond only affected minor structural deviations of the heme environmental structure in the carbonmonoxy form, whereas more substantial structural alterations were induced in the deoxygenated form. Laser photolysis experiments for carbon monoxide rebinding clearly showed that the artificial S-S bond accelerates the bimolecular rebinding rate from 1.0 to 1.8 microM-1 s-1 and increases the geminate yield from 0.072 to 0.092. The ligand migration rate from the solvent to the heme pocket and the bond formation rate from the heme pocket to the heme iron also increased. The free energy diagram for the mutants indicates that the energy barrier for the bond formation was raised as well as that for the ligand migration by introduction of the disulfide bond. However, the effects of the disulfide linkage at the EF corner on the kinetic parameter is much smaller than those of the amino acid substitutions located in the heme cavity. We can conclude that the perturbation of the protein fluctuations by formation of the disulfide bond would be localized at the mutation site or the contributions from other regions and motions might be more important for the ligand binding dynamics.

Picosecond geminate recombination of CO to the complexes calmodulin*heme-CO and calmodulin*heme-CO*melittin

Picosecond CO recombination kinetics have been measured after photodissociation of the artificial complexes calmodulin*heme-CO and calmodulin*heme-CO*melittin. These systems show an enhancement of the geminate fraction of kinetics relative to unbound heme-CO, due in part to fast geminate kinetics (tau=50ps for the initial phase), as well as a decrease in the rate of migration of CO away from the binding site. This indicates that calmodulin provides a complete pocket around the heme group. Rather than competing with the hemes for binding to calmodulin, the melittin seems to act as a cap to further enclose the hemes; melittin increases the affinity of calmodulin for heme-CO, but only weakly affects the CO recombination kinetics.

Two-dimensional distributions of activation enthalpy and entropy from kinetics by the maximum entropy method

The maximum entropy method (MEM) is used to numerically invert the kinetics of ligand rebinding at low temperatures to obtain the underlying two-dimensional distribution of activation enthalpies and entropies, g(H,S). A global analysis of the rebinding of carbon monoxide (CO) to myoglobin (Mb), monitored in the Soret band at temperatures from 60 to 150 K, is performed using a Newton-Raphson optimization algorithm. The MEM approach describes the data much better than traditional least-squares analyses, reducing chi 2 by an order of magnitude. The MEM resolves two barrier distributions suggestive of rebinding to different bound conformations of MbCO, the so-called A1 and A3 substates, whose activation barriers have been independently estimated from kinetics monitored in the infrared. The distribution corresponding to A3 possesses higher activation entropies, also consistent with infrared measurements. Within an A substate, correlations of S and H are recovered qualitatively from simulated data but can be difficult to obtain from experimental data. When the rebinding measured at 60 K is excluded from the inversion, two peaks are no longer clearly resolved. Thus, data of very high quality are required to unambiguously determine the kinetic resolvability of subpopulations and the shape of the barrier distribution for a single A substate.

Conformational substates in azurin

Azurin is a small blue copper protein in the electron transfer chain of denitrifying bacteria. It forms a photolabile complex with nitric oxide (NO) at low temperatures. We studied the temperature dependence of the ligand binding equilibrium and the kinetics of the association reaction after photodissociation over a wide range of temperature (80-280 K) and time (10(-6)-10(2) s). The nonexponential rebinding below 200 K is independent of the NO concentration and is interpreted as internal recombination. The rebinding can be modeled with the Arrhenius law by using a single preexponential factor of 6.3 x 10(8) s-1 and a Gaussian distribution of enthalpy barriers centered at 23 kJ/mol with a width of 11 kJ/mol. Above 200 K, a slower, exponential rebinding process appears. The dependence of the kinetics on the NO concentration characterizes this reaction as bimolecular rebinding. The binding kinetics of NO to azurin show impressive analogies to the binding of carbon monoxide to myoglobin. We conclude that conformational substates occur not only in heme proteins but also in proteins with different active sites and secondary structures.

Protein dynamics. An overview on flash-photolysis over broad temperature ranges

Ligand binding kinetics to heme-proteins between 40 and 300 K point to a regulatory role of protein dynamics. A protein-specific susceptibility of the heme-iron reactivity to dynamic fluctuations emerges from the distribution of reaction enthalpies derived from flash-photolysis measurements below ca. 180 K; we quantify it in terms of 'intramolecular viscosity', postulating that narrow low-temperature enthalpy distributions correspond to low internal viscosity and vice versa. The thermal evolution of ligand binding kinetics suggests, with other results, an interplay between high-frequency transitions of the amino acid side chains and low-frequency collective motions as a possible regulatory mechanism of protein dynamics.

Inhomogeneous broadening and kinetic carbon monoxide isotope effects in low-temperature carbon monoxide recombination with myoglobin and hemoglobin

We have analyzed the non-exponential kinetics, the temperature variation, and the CO isotope effects of the CO recombination reactions with myoglobin and single-chain hemoglobin. The analysis rests on multiphonon quantum-mechanical chemical-rate theory combined with static inhomogeneous broadening of either the reorganization free energy or the reaction Gibbs free energy. The simplest specific model which can account for all the data contains an inhomogeneous distribution function of width 0.2-0.3 eV, independent of temperature down to the tunnel transition at about 20 K, two discrete nuclear coordinates of low vibrational frequency (60-150 cm-1) representing iron-heme and CO bending motion, the CO stretching motion of frequency about 2000 cm-1, and additional inhomogeneous broadening of the protein and CO bending configuration below the tunnel transition temperature. The model appears somewhat involved but in return provides corresponding insight in the dynamics of this important class of processes.

Development of an extended simulated annealing method: application to the modeling of complementary determining regions of immunoglobulins

An extended simulated annealing process (ESAP) has been developed in order to obtain an ensemble of conformations of a peptide segment from a protein fluctuating at a given temperature. The annealing process was performed with a fast Monte Carlo method using the scaled collective variables developed by Noguti and Go. The system was divided into two parts: one consists of one or more peptide segments and is flexible around the main-chain and side-chain torsional angles; the other represents the rest of the molecule and was maintained fixed at the atomic positions determined by x-ray experiments. The target function included the nonbonding atomic interactions and a distance function to anchor the N and C terminal ends of each segment to the fixed part. Three systems of complementary determining regions (CDR) of antibodies were tested and compared to x-ray data: L2 loop (7 residues) of the light chain of lambda-type Bence-Jones protein, H1 and the H2 loops (14 residues) of McPC603, and H1 and H2 loops (12 residues) of HyHEL-5. Each state of CDR conformations was characterized at room temperature by the average of their coordinates (average conformation) and the internal energy. With a limited number of annealing processes (10), starting from the extended conformation, we have obtained states with conformations close to the observed x-ray structures, from 1.1 to 1.7 A root mean square deviation (rmsd) of main-chain atoms depending on the system. These states were identical or within 0.25 A rmsd of those of lowest internal energy. For unknown CDR structures the criteria of lowest internal energies from ESAP can be used to predict hypervariable loop structures in antibodies with an accuracy comparable to other methods.

Biophysical aspects of neutron scattering from vibrational modes of proteins

This review describes a major portion of the published work on neutron scattering experiments aimed at measuring large scale motions in proteins. The importance of these motions for enzyme function and oxygen transport is indicated. The theory applicable to each type of neuron scattering measurement is given and results are discussed with a view to biological relevance. New experiments are suggested and a comparison of neutron scattering data is made with results from other techniques such as raman scattering, infrared absorption, photolysis and molecular dynamics simulations.

Biological transport processes and space dimension

We discuss the generic time behavior of reaction-diffusion processes capable of modeling various types of biological transport processes, such as ligand migration in proteins and gating fluctuations in ion channel proteins. The main observable in these two cases, the fraction of unbound ligands and the probability of finding the channel in the closed state, respectively, exhibits an algebraic t-1/2 decay at intermediate times, followed by an exponential cutoff. We provide a simple framework for understanding these observations and explain their ubiquity by showing that these qualitative results are independent of space dimension. We also derive an experimental criterion to distinguish between a one-dimensional process and one whose effective dimension is higher.

Protein dynamics. Comparative investigation on heme-proteins with different physiological roles

We report the low temperature carbon monoxide recombination kinetics after photolysis and the temperature dependence of the visible absorption spectra of the isolated alpha SH-CO and beta SH-CO subunits from human hemoglobin A in ethylene glycol/water and in glycerol/water mixtures. Kinetic measurements on sperm whale (Physeter catodon) myoglobin and previously published optical spectroscopy data on the latter protein and on human hemoglobin A, in both solvents, (Cordone, L., A. Cupane, M. Leone, E. Vitrano, and D. Bulone. 1988. J. Mol. Biol. 199:312-218) are taken as reference. Low temperature flash photolysis data are analyzed within the multiple substates model proposed by Frauenfelder and co-workers (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355-5373). Within this model a distribution of activation enthalpies for ligand binding accounts for the structural heterogeneity of the protein, while the preexponential factor, containing also the entropic contribution to the free energy of the process, is considered to be constant for all conformational substates. Optical spectra are deconvoluted in gaussian components and the temperature dependence of the moments of the resulting bands is analyzed, within the harmonic Frank-Condon approximation, to obtain information on the stereodynamic properties of the heme pocket. The kinetic and spectral parameters thus obtained are found to be protein dependent also with respect to their sensitivity to changes in the composition of the external medium. A close correlation between the kinetic and spectral features is observed for the proteins examined under all experimental conditions studied. The results reported are discussed in terms of differences in the heme pocket structure and in the conformational heterogeneity among the various proteins, as related to their different capability to accommodate constraints imposed by the external medium.

Quaternary structure and the geminate recombination of carp hemoglobin with methylisocyanide

The kinetics of geminate recombination were studied for the methylisocyanide derivative of carp hemoglobin. Carp hemoglobin is of interest because it has been established that the fully liganded form switches between a high affinity R state at pH 9 and a low affinity T state at pH 6 in the presence of IHP. Geminate recombination was observed on both the picosecond and the nanosecond time scales under all conditions; however, only a small variation is observed in the rates and the yields of geminate recombination as the protein switches from the R to the T state. Taken together with overall "on" and "off" rates, the data indicate that the change from the R to the T configuration affects bond breaking most, but also influences subsequent escape from the protein as well as both entry into the protein and bond formation. There is some reason to postulate tertiary conformational change in the T state on the microsecond time scale following ligand escape from the protein.

Infrared spectra of carbonyl hemoglobins: characterization of dynamic heme pocket conformers

The infrared spectra for carbon monoxide complexed to hemoglobins were examined in the C-O stretch region. Deconvolution of the spectra requires four bands and supports the presence of four distinct conformers at the ligand binding site. Most typical hemoglobins exhibit only one predominant conformer for each subunit represented by a band at 1951 cm-1 in contrast to myoglobins, which typically exist in two major conformations. Several hemoglobins with an enlarged heme pocket are shown to shift the C-O frequency into the higher frequency conformer regions. Many factors, including pH, temperature, solvents, and divalent metals, are also shown to be capable of expanding the heme pocket. Only very specific structural changes that can reduce the size of the heme pocket will result in the lower frequency conformers. The weighted averages of the multiple CO vibrational frequencies are linearly related to the single 13CO NMR chemical shift values and to the exponential of fast CO on-rates. Conformer interconversion occurs at a rate greater than 10(4) s-1. The infrared C-O stretch spectra provide qualitative and quantitative information on the structural dynamics, stability, and ligand binding properties of hemoglobins.

O2 and CO reactions with heme proteins: quantum yields and geminate recombination on picosecond time scales

Picosecond time-resolved absorption spectroscopy and low-temperature studies have been undertaken in order to understand the nature of the intrinsic quantum yields and geminate recombination of carbon monoxide and oxygen to hemoglobin and myoglobin. We find that the photoproduct yields at 40 ps and long times (minutes) after photolysis at 8 K are similar; however, the yield of oxygen photoproducts is 0.4 +/- 0.1 while the yield of carbon monoxide photoproducts is 1.0 +/- 0.1 for both myoglobin and hemoglobin. Measurements in the Soret, near-infrared, and far-IR are used to quantitate the photoproduct yields. These results call into question previous cryogenic kinetic studies of O2 recombination. Significant subnanosecond geminate recombination is observed in oxyhemoglobin down to 150 K, while below 100 K this geminate recombination disappears. The lower photoproduct yields for oxyheme protein complexes can be attributed to both subnanosecond and subpicosecond recombination events which are ligand and protein dynamics dependent.

E.p.r. studies of photolysis of nitrosyl haemoglobin at low temperatures

Photolysis of HbNO has been studied from 6.2 K to 15.5 K by electron spin resonance during and after continuous illumination. Non-exponential kinetics of both dissociation and reassociation of NO was observed. The prolonged illumination separates the fast and slow ligands. This picture is consistent with NO tunnelling from two sites at different distances from the bound position. This result is obtained using a model of a sum of two exponentials or of conformational substates.

Orientation of carbon monoxide and structure-function relationship in carbonmonoxymyoglobin

Fourier transform infrared spectroscopy of the CO stretch bands in carbonmonoxymyoglobin (MbCO) reveals three major bands implying that MbCO exists in three major substates, A0, A1, and A3. After photolysis at low temperatures the CO is in the heme pocket, and the resulting CO stretch bands represent the B substates. Photoselection experiments determine the orientation of CO in the A (bound) and B (photolyzed) substates: Small fractions of MbCO are photolyzed at 10 K with linearly polarized light at 540 nm. The resulting linear dichroism in the A and B IR bands yields the tilt angle between the heme normal and CO. The average angles are as follows: alpha (A0) = 15 degrees +/- 3 degrees; alpha (A1) = 28 degrees +/- 2 degrees, and alpha (A3) = 33 degrees +/- 4 degrees. The A bands are inhomogeneously broadened; the angle alpha shows a wavenumber dependence within the A bands. The wavenumber dependence is interpreted as a distribution of the tilt angle within the individually inhomogeneous A substates, thus providing a structural parameter to characterize the distribution of the conformational substates. The B substates exhibit no induced linear dichroism; in the photolyzed substates the ligand is randomly oriented with respect to the heme plane. The present results together with earlier data on static and kinetic properties of CO binding to Mb establish relations among spectroscopic, structural, energetic, and functional parameters.

Assignment of proton resonances in the NMR spectrum of carbonmonoxy hemoglobin beta subunit tetramers

Isolated beta chains from human adult hemoglobin at millimolar concentration are mainly associated to form beta 4 tetramers. We were able to obtain relevant two-dimensional proton nuclear magnetic resonance (NMR) spectra of such supermolecular complexes (Mr approximately 66,000) in the carboxylated state. Analysis of the spectra enabled us to assign the major part of the proton resonances corresponding to the heme substituents. We also report assignments of proton resonances originating from 12 amino acid side chains mainly situated in the heme pocket. These results provide a basis for a comparative analysis of the tertiary heme structure in isolated beta(CO) chains in solution and in beta(CO) subunits of hemoglobin crystals. The two structures are generally similar. A significantly different position, closer to the heme center, is predicted by the NMR for Leu-141 (H19) in isolated beta chains. Comparison of the assigned resonances of conserved amino acids in alpha chains, beta chains and sperm whale myoglobin indicates a close similarity of the tertiary heme pocket structure in the three homologous proteins. Significant differences were noted on the distal heme side, at the position of Val-E11, and on Leu-H19 and Phe-G5 position on the proximal side.

One- and two-dimensional NMR investigations of the heme pocket in free alpha(CO) chains from human hemoglobin

Two-dimensional nuclear magnetic resonance techniques were used to assign resonances corresponding to heme pocket residues of the isolated alpha(CO) subunits of the human adult hemoglobin (HbA). The assignment procedure was based on the partial identification of the amino acid spin system from the J-correlated (COSY) spectrum and on the nuclear Overhauser effect connectivities (from NOSEY spectra) with the heme substituents. We present here partial assignments corresponding to five amino acid residues: Leu86, Leu-91, Val-93, Leu-101 and Leu-136. Starting from the known crystallographic structure of the alpha subunit in the hemoglobin tetramer, we applied a dipolar model to compute the ring-current shift of the protons from fifteen amino acid residues in the heme pocket. Comparison of the predicted and observed chemical shifts suggests that there is a very close similarity between the heme pocket tertiary structure of the alpha(CO) subunits in crystals of HbA(CO) and of the free alpha(CO) chains. The one-dimensional NMR spectra were used to monitor the pH-induced structural changes, the effects of chemical modification and of ligand substitution. Upon increasing the pH from 5.6 to 9.0 the structure of the heme environment appears to be invariant with the exception of some residues in the CD corner. The structure is also largely conserved when p-chloromercuribenzoate is bound to Cys-104. In contrast, the substitution of CO by O2 as ligand induces many large changes in the heme cavity which can be partially characterized by NMR spectroscopy.

Rebinding and relaxation in the myoglobin pocket

The infrared stretching bands of carboxymyoglobin (MbCO) and the rebinding of CO to Mb after photodissociation have been studied in the temperature range 10-300 K in a variety of solvents. Four stretching bands imply that MbCO can exist in four substates, A0-A3. The temperature dependences of the intensities of the four bands yield the relative binding enthalpies and and entropies. The integrated absorbances and pH dependences of the bands permit identification of the substates with the conformations observed in the X-ray data (Kuriyan et al., J. Mol. Biol. 192 (1986) 133). At low pH, A0 is hydrogen-bonded to His E7. The substates A0-A3 interconvert above about 180 K in a 75% glycerol/water solvent and above 270 K in buffered water. No major interconversion is seen at any temperature if MbCO is embedded in a solid polyvinyl alcohol matrix. The dependence of the transition on solvent characteristics is explained as a slaved glass transition. After photodissociation at low temperature the CO is in the heme pocket B. The resulting CO stretching bands which are identified as B substates are blue-shifted from those of the A substates. At 40 K, rebinding after flash photolysis has been studied in the Soret, the near-infrared, and the integrated A and B substates. All data lie on the same rebinding curve and demonstrate that rebinding is nonexponential in time from at least 100 ns to 100 ks. No evidence for discrete exponentials is found. Flash photolysis with monitoring in the infrared region shows four different pathways within the pocket B to the bound substates Ai. Rebinding in each of the four pathways B----A is nonexponential in time to at least 10 ks and the four pathways have different kinetics below 180 K. From the time and temperature dependence of the rebinding, activation enthalpy distributions g(HBA) and preexponentials ABA are extracted. No pumping from one A substate to another, or one B substate to another, is observed below the transition temperature of about 180 K. If MbCO is exposed to intense white light for 10-10(3) s before being fully photolyzed by a laser flash, the amplitude of the long-lived states increases. The effect is explained in terms of a hierarchy of substates and substate symmetry breaking. The characteristics of the CO stretching bands and of the rebinding processes in the heme pocket depend strongly on the external parameters of solvent, pH and pressure. This sensitivity suggests possible control mechanisms for protein reactions.

Energy barriers in binding of carbon monoxide and oxygen to heme model compounds

Binding of carbon monoxide and oxygen to sterically protected heme model compounds (basket-handle porphyrins) was investigated in liquid toluene at temperatures from 180 to 300 K by laser flash photolysis. Only a single exponential rebinding process from the solvent could be seen in the time range of 20 nsec to milliseconds. The fraction of ligands that initially escaped into the solvent decreased when the temperature was lowered, and the Arrhenius plots for the rebinding rate coefficients were found to deviate significantly from linearity. These findings suggest that protected heme model compounds might react according to a double energy-barrier scheme. In contrast, the reaction of an unprotected porphyrin can be described by a single energy barrier.