Dynamics of dioxygen and carbon monoxide binding to soybean leghemoglobin

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

Distal heme pocket regulation of ligand binding and stability in soybean leghemoglobin

Leghemoglobins facilitate diffusion of oxygen through root tissue to a bacterial terminal oxidase in much the same way that myoglobin transports oxygen from blood to muscle cell mitochondria. Leghemoglobin serves an additional role as an oxygen scavenger to prevent inhibition of nitrogen fixation. For this purpose, the oxygen affinity of soybean leghemoglobin is 20-fold greater than myoglobin, resulting from an 8-fold faster association rate constant combined with a 3-fold slower dissociation rate constant. Although the biochemical mechanism used by myoglobin to bind oxygen has been described in elegant detail, an explanation for the difference in affinity between these two structurally similar proteins is not obvious. The present work demonstrates that, despite their similar structures, leghemoglobin uses methods different from myoglobin to regulate ligand affinity. Oxygen and carbon monoxide binding to a comprehensive set of leghemoglobin distal heme pocket mutant proteins in comparison to their myoglobin counterparts has revealed some of these mechanisms. The "distal histidine" provides a crucial hydrogen bond to stabilize oxygen in myoglobin but has little effect on bound oxygen in leghemoglobin and is retained mainly for reasons of protein stability and prevention of heme loss. Furthermore, soybean leghemoglobin uses an unusual combination of HisE7 and TyrB10 to sustain a weak stabilizing interaction with bound oxygen. Thus, the leghemoglobin distal heme pocket provides a much lower barrier to oxygen association than occurs in myoglobin and oxygen dissociation is regulated from the proximal heme pocket.

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.

Cytochrome P-450-CO and substrates: lessons from ligand binding under high pressure

An overview of the application of high-pressure studies on the carbon monoxide complex of cytochrome P-450 is given. Different approaches to characterize ligand binding steps, the conformational states and substates and the compressibility of the ligand-bound complex are reviewed. A particular focus is the effect of substrates on these properties. It is shown that substrate mobility, compressibility and water accessibility are interrelated and may have functional meaning.

Insight into protein structure and protein-ligand recognition by Fourier transform infrared spectroscopy

An overview of the application of Fourier transform infrared spectroscopy for the analysis of the structure of proteins and protein-ligand recognition is given. The principle of the technique and of the spectra analysis is demonstrated. Spectral signal assignments to vibrational modes of the peptide chromophore, amino acid side chains, cofactors and metal ligands are summarized. Several examples for protein-ligand recognition are discussed. A particular focus is heme proteins and, as an example, studies of cytochrome P450 are reviewed. Fourier transform infrared spectroscopy in combination with the various techniques such as time-resolved and low-temperature methods, site-directed mutagenesis and isotope labeling is a helpful approach to studying protein-ligand recognition.

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.

Characterization of recombinant soybean leghemoglobin a and apolar distal histidine mutants

The cDNA for soybean leghemoglobin a (Lba) was cloned from a root nodule cDNA library and expressed in Escherichia coli. The crystal structure of the ferric acetate complex of recombinant wild-type Lba was determined at a resolution of 2.2 A. Rate constants for O2, CO and NO binding to recombinant Lba are identical with those of native soybean Lba. Rate constants for hemin dissociation and auto-oxidation of wild-type Lba were compared with those of sperm whale myoglobin. At 37 degrees C and pH 7, soybean Lba is much less stable than sperm whale myoglobin due both to a fourfold higher rate of auto-oxidation and to a approximately 600-fold lower affinity for hemin. The role of His61(E7) in regulating oxygen binding was examined by site-directed mutagenesis. Replacement of His(E7) with Ala, Val or Leu causes little change in the equilibrium constant for O2 binding to soybean Lba, whereas the same mutations in sperm whale myoglobin cause 50 to 100-fold decreases in K(O2). These results show that, at neutral pH, hydrogen bonding with His(E7) is much less important in regulating O2 binding to the soybean protein. The His(E7) to Phe mutation does cause a significant decrease in K(O2) for Lba, apparently due to steric hindrance of the bound ligand. The rate constants for O2 dissociation from wild-type and native Lba decrease significantly with decreasing pH. In contrast, the O2 dissociation rate constants for mutants with apolar E7 residues are independent of pH, suggesting that hydrogen bonding to the distal histidine residue in the native protein is enhanced under acid conditions. All of these results support the hypothesis that the high affinity of Lba for oxygen and other ligands is determined primarily by enhanced accessibility and reactivity of the heme group.

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.

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.

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.

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.

Computational studies of ligand diffusion in globins: I. Leghemoglobin

The thermally assisted diffusion of a small ligand (carbon monoxide) through a protein matrix (lupine leghemoglobin) is investigated computationally. The diffusion paths are calculated by a variant of the time-dependent Hartree approximation which we call LES (locally enhanced sampling). The variant which was recently introduced by Elber and Karplus is based on the classical TDSCF approximation of Gerber et al. The simulation enables more significant search for diffusion pathways than was possible before. This is done by increasing the number of ligand trajectories using a single trajectory for the protein. We compare qualitatively diffusion rates in leghemoglobin and in myoglobin. The calculation shows that the diffusion in leghemoglobin is much faster than the diffusion in myoglobin, in agreement with experiment. The gate in leghemoglobin is opened by fluctuations at a close contact between the B/C and the G helices. The most relevant fluctuation is the rigid shift of the C helix with respect to the G helix. This path is not observed in a comparable calculation for myoglobin. This finding is rationalized by the lack of the D helix in leghemoglobin and a significantly more flexible CE loop. Supporting experimental evidence for the importance of the CE loop in leghemoglobin can be found in the kinetics studies of Gibson et al.

Contribution of protein conformation to heme stereochemistry and reactivity. Low-temperature magnetic circular dichroism data

Visible and near infrared magnetic circular dichroism (MCD) spectra of heme proteins and enzymes as well as those of a protein-free heme bound to 2-methylimidazole were recorded and compared at 4.2 K in unrelaxed metastable and relaxed equilibrium heme stereochemistry. The relaxed and unrelaxed stereochemistries of a 5-coordinate ferrous heme were generated by chemical reduction of iron at room temperature before freezing the sample and by photolysis of CO or O2 complexes at 4.2 K, respectively. The results are discussed in terms of a protein contribution into energies of the Fe-N epsilon(His) and Fe-N(pyrrols) bonds and their change on a ligand binding. We observed and analyzed cases of weak (myoglobin, hemoglobin) and strong (leghemoglobin, peroxidases) constraints imposed by the protein conformation on the proximal heme stereochemistry by comparing the bond energies in proteins with those in the protoheme-(2-methylimidazole) model compound. The role of a protein moiety in modulating the ligand binding properties of leghemoglobin and the heme reactivity of horseradish peroxidase is discussed.

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.

Recombination of carbon monoxide to ferrous horseradish peroxidase types A and C

The recombination of carbon monoxide to isoenzymes A2 and C of horseradish peroxidase (HRP) was studied as a function of temperature (2 to 320 K) and pH (5 to 8.3) with flash photolysis and infrared difference absorption. At low temperatures three geminate recombination processes are observed. One of these internal processes, denoted by I*, is exponential in time with a rate coefficient that deviates strongly from an Arrhenius behavior below 100 K, implying phonon-assisted tunneling. The two other processes, denoted by I, are non-exponential in time and related to different carbonyl isomers, as shown by the infrared difference spectra. The existence of three internal processes indicates that HRP differs considerably from myoglobin where only one internal process, I, is seen. Moreover, the internal processes in HRP are faster than process I in myoglobin. At 300 K, only one recombination process from the solvent is observed and it is very slow (lambda s approximately 1 s-1 at 1 atm CO (1 atm = 101,325 Pa)), much slower than the corresponding association process in myoglobin. Since process I is fast, but binding from the solvent is slow, the barrier at the heme cannot be responsible for the small association rate. The infrared absorption difference spectra of the amide I/II bands indicate that photolysis and recombination trigger a two-step structural change. The slow recombination rate at 300 K can thus be explained by the large Gibbs energy of the conformational transition that is necessary to let CO move into the heme pocket. The partition coefficient for the CO in the heme pocket and the solvent is extremely small, while bond formation with the heme iron occurs in less than 100 nanoseconds.

Protoheme-carbon monoxide geminate kinetics

Recombination kinetics of CO to protoheme after laser photolysis have been measured vs. temperature and viscosity. A 25-ns laser pulse was focused on the sample to produce an excitation rate of 10(9)/s per heme. This temporarily populates the heme-CO state of dissociated pairs which either separate or recombine on a picosecond time scale in viscous glycerol-water solutions. From the equilibrium amplitude of the fraction dissociated during the laser pulse, the geminate recombination rate constant is calculated to be 3 X 10(9)/s. This rate coefficient is only weakly dependent on temperature or viscosity. As previously observed the fraction that escapes depends on the solvent viscosity [Marden, M. C. (1983) Ph.D. Thesis, University of Illinois-Urbana]. A model consisting of a single barrier plus diffusive escape is used to simulate the kinetics during and just after the flash.