PubMed Health⌕ Search

Biomedical subjects

E R Henry

Publications and source records attributed to E R Henry.

At least 19 recordsLinked to original sources

Oxygen binding by alpha(Fe2+)2beta(Ni2+)2 hemoglobin crystals.

Oxygen binding by hemoglobin fixed in the T state either by crystallization or by encapsulation in silica gels is apparently noncooperative. However, cooperativity might be masked by different oxygen affinities of alpha and beta subunits. Metal hybrid hemoglobins, where the noniron metal does not bind oxygen, provide the opportunity to determine the oxygen affinities of alpha and beta hemes separately. Previous studies have characterized the oxygen binding by alpha(Ni2+)2beta(Fe2+)2 crystals. Here, we have determined the three-dimensional (3D) structure and oxygen binding of alpha(Fe2+)2beta(Ni2+)2 crystals grown from polyethylene glycol solutions. Polarized absorption spectra were recorded at different oxygen pressures with light polarized parallel either to the b or c crystal axis by single crystal microspectrophotometry. The oxygen pressures at 50% saturation (p50s) are 95 +/- 3 and 87 +/- 4 Torr along the b and c crystal axes, respectively, and the corresponding Hill coefficients are 0.96 +/- 0.06 and 0.90 +/- 0.03. Analysis of the binding curves, taking into account the different projections of the alpha hemes along the optical directions, indicates that the oxygen affinity of alpha1 hemes is 1.3-fold lower than alpha2 hemes. Inspection of the 3D structure suggests that this inequivalence may arise from packing interactions of the Hb tetramer within the monoclinic crystal lattice. A similar inequivalence was found for the beta subunits of alpha(Ni2+)2beta(Fe2+)2 crystals. The average oxygen affinity of the alpha subunits (p50 = 91 Torr) is about 1.2-fold higher than the beta subunits (p50 = 110 Torr). In the absence of cooperativity, this heterogeneity yields an oxygen binding curve of Hb A with a Hill coefficient of 0.999. Since the binding curves of Hb A crystals exhibit a Hill coefficient very close to unity, these findings indicate that oxygen binding by T-state hemoglobin is noncooperative, in keeping with the Monod, Wyman, and Changeux model.

Crystallization↗

Fast kinetics and mechanisms in protein folding.

This review describes how kinetic experiments using techniques with dramatically improved time resolution have contributed to understanding mechanisms in protein folding. Optical triggering with nanosecond laser pulses has made it possible to study the fastest-folding proteins as well as fundamental processes in folding for the first time. These include formation of alpha-helices, beta-sheets, and contacts between residues distant in sequence, as well as overall collapse of the polypeptide chain. Improvements in the time resolution of mixing experiments and the use of dynamic nuclear magnetic resonance methods have also allowed kinetic studies of proteins that fold too fast (greater than approximately 10(3) s-1) to be observed by conventional methods. Simple statistical mechanical models have been extremely useful in interpreting the experimental results. One of the surprises is that models originally developed for explaining the fast kinetics of secondary structure formation in isolated peptides are also successful in calculating folding rates of single domain proteins from their native three-dimensional structure.

Animals↗

Is cooperative oxygen binding by hemoglobin really understood?

The enormous success of structural biology challenges the physical scientist. Can biophysical studies provide a truly deeper understanding of how a protein works than can be obtained from static structures and qualitative analysis of biochemical data? We address this question in a case study by presenting the key concepts and experimental results that have led to our current understanding of cooperative oxygen binding by hemoglobin, the paradigm of structure function relations in multisubunit proteins. We conclude that the underlying simplicity of the two-state allosteric mechanism could not have been demonstrated without novel physical experiments and a rigorous quantitative analysis.

Chemistry, Physical↗

A statistical mechanical model for beta-hairpin kinetics.

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a beta-hairpin, the minimal structural element of the antiparallel beta-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, beta-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 2(15) possible conformations and 15 x 2(15) possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the beta-hairpin peptide.

Animals↗

Can a two-state MWC allosteric model explain hemoglobin kinetics?

We have analyzed the nanosecond-millisecond kinetics of ligand binding and conformational changes in hemoglobin. The kinetics were determined from measurements of precise time-resolved optical spectra following nanosecond photodissociation of the heme-carbon monoxide complex. To fit the data, it was necessary to extend the two-state allosteric model of Monod, Wyman, and Changeux (MWC) to include geminate ligand rebinding and nonexponential tertiary relaxation within the R quaternary structure. Considerable simplification of the model is obtained by using a linear free energy relation for the rates of quaternary transitions, and by incorporating concepts from recent studies on the physics of geminate rebinding and conformational changes in myoglobin. The model, described by 85 coupled differential equations, quantitatively explains a demanding set of complex kinetic data. Moreover, with the same set of kinetic parameters it simultaneously fits the equilibrium data on ligand binding and the distribution of ligation states. The present results, together with those from single-crystal oxygen binding studies, indicate that the two-state MWC allosteric model has survived its most critical tests.

Allosteric Site↗

Allosteric effectors do not alter the oxygen affinity of hemoglobin crystals.

In solution, the oxygen affinity of hemoglobin in the T quaternary structure is decreased in the presence of allosteric effectors such as protons and organic phosphates. To explain these effects, as well as the absence of the Bohr effect and the lower oxygen affinity of T-state hemoglobin in the crystal compared to solution, Rivetti C et al. (1993a, Biochemistry 32:2888-2906) suggested that there are high- and low-affinity subunit conformations of T, associated with broken and unbroken intersubunit salt bridges. In this model, the crystal of T-state hemoglobin has the lowest possible oxygen affinity because the salt bridges remain intact upon oxygenation. Binding of allosteric effectors in the crystal should therefore not influence the oxygen affinity. To test this hypothesis, we used polarized absorption spectroscopy to measure oxygen binding curves of single crystals of hemoglobin in the T quaternary structure in the presence of the "strong" allosteric effectors, inositol hexaphosphate and bezafibrate. In solution, these effectors reduce the oxygen affinity of the T state by 10-30-fold. We find no change in affinity (< 10%) of the crystal. The crystal binding curve, moreover, is noncooperative, which is consistent with the essential feature of the two-state allosteric model of Monod J, Wyman J, and Changeux JP (1965, J Mol Biol 12:88-118) that cooperative binding requires a change in quaternary structure. Noncooperative binding by the crystal is not caused by cooperative interactions being masked by fortuitous compensation from a difference in the affinity of the alpha and beta subunits. This was shown by calculating the separate alpha and beta subunit binding curves from the two sets of polarized optical spectra using geometric factors from the X-ray structures of deoxygenated and fully oxygenated T-state molecules determined by Paoli M et al. (1996, J Mol Biol 256:775-792).

Allosteric Regulation↗

The use of matrix methods in the modeling of spectroscopic data sets.

We describe a general approach to the model-based analysis of sets of spectroscopic data that is built upon the techniques of matrix analysis. A model hypothesis may often be expressed by writing a matrix of measured spectra as the product of a matrix of spectra of individual molecular species and a matrix of corresponding species populations as a function of experimental conditions. The modeling procedure then requires the simultaneous determination of a set of species spectra and a set of model parameters (from which the populations are derived), such that this product yields an optimal description of the measured spectra. This procedure may be implemented as an optimization problem in the space of the (possibly nonlinear) model parameters alone, coupled with the efficient solution of a corollary linear optimization problem using matrix decomposition methods to obtain a set of species spectra corresponding to any set of model parameters. Known species spectra, as well as other information and assumptions about spectral shapes, may be incorporated into this general framework, using parametrized analytical functional forms and basis-set techniques. The method by which assumed relationships between global features (e.g., peak positions) of different species spectra may be enforced in the modeling without otherwise specifying the shapes of the spectra will be shown. We also consider the effect of measurement errors on this approach and suggest extensions of the matrix-based least-squares procedures applicable to situations in which measurement errors may not be assumed to be normally distributed. A generalized analysis procedure is introduced for cases in which the species spectra vary with experimental conditions.

Ligands↗

Oxygen binding by single crystals of hemoglobin: the problem of cooperativity and inequivalence of alpha and beta subunits.

Oxygen binding by the human hemoglobin tetramer in the T quaternary structure is apparently noncooperative in the crystalline state (Hill n = 1.0), as predicted by the two-state allosteric model of Monod, Wyman, and Changeux (MWC) (Mozzarelli et al., Nature 351:416-419, 1991; Rivetti et al., Biochemistry 32:2888-2906, 1993). However, cooperativity within the tetramer can be masked by a difference in affinity between the alpha and beta subunits. Indeed, analysis of the binding curves derived from absorption of light polarized along two different crystal directions, for which the projections of the alpha and beta hemes are slightly different, revealed an inequivalence in the intrinsic oxygen affinity of the alpha and beta subunits (p50(alpha) approximately 80 torr, p50(beta) approximately 370 torr at 15 degrees C) that compensates a small amount of cooperativity (Rivetti et al., Biochemistry 32:2888-2906, 1993). To further investigate this problem, we have measured oxygen binding curves of single crystals of hemoglobin (in a different lattice) in which the iron in the alpha subunits has been replaced by the non-oxygen-binding nickel(II). The Hill n is 0.90 +/- 0.06, and the p50 is slightly different for light polarized parallel to different crystal directions, indicating a very small difference in affinity between the two crystallographically inequivalent beta subunits. The average crystal p50 is 110 +/- 20 torr at 15 degrees C, close to the p50 of 80 torr observed in solution, but about threefold less than the p50 calculated by Rivetti et al. (Biochemistry 32:2888-2906, 1993) for the beta subunits of the unsubstituted tetramer. These results suggest that Rivetti et al., if anything, overestimated the alpha/beta inequivalence. They therefore did not underestimate the cooperativity within the T quaternary structure, when they concluded that it represents a small deviation from the perfectly noncooperative binding of an MWC allosteric model. Our conclusion of nearly perfect MWC behavior for binding to the T state of unmodified hemoglobin raises the question of the relevance of the large T-state cooperativity inferred for cyanide binding to partially oxidized hemoglobin (Ackers et al., Science 255:54-63, 1992).

Allosteric Site↗

Conformational relaxation and ligand binding in myoglobin.

Absorption spectroscopy with nanosecond time resolution shows that myoglobin undergoes conformational relaxation on the same time scale as geminate rebinding of carbon monoxide. Ligand rebinding following photodissociation of the heme-CO complex was measured from the amplitude of the average difference spectrum, while conformational changes were measured from changes in the detailed shape of the Soret spectra of the deoxyhemes. Experiments in which the solvent viscosity was varied between 1 and 300 cP and the temperature between 268 and 308 K were analyzed by fitting the multiwavelength kinetic data with both empirical and molecular models. Novel numerical techniques were employed in fitting the data, including the use of singular value decomposition to remove the effects of temperature and solvent on the spectra and of a Monte Carlo method to overcome the multiple minimum problem in searching parameter space. The molecular model is the minimal model that incorporates all of the major features of myoglobin kinetics at ambient temperatures, including a fast and slow rebinding conformation and two geminate states for each conformation. The results of fitting the kinetic data with this model indicate that the geminate-rebinding rates for the two conformations differ by at least a factor of 100. The differences between the spectra of the two conformations generated from the fits are similar to the differences between those of the R and T conformations of hemoglobin. In modeling the data, the dependence of the rates on temperature and viscosity was parametrized using a modification of Kramers theory which includes the contributions of both protein and solvent to the friction. The rate of the transition from the fast to the slow rebinding conformation is found to be inversely proportional to the viscosity when the viscosity exceeds about 30 cP and nearly viscosity independent at low viscosity. The viscosity dependence at high viscosities suggests that the two conformations differ by the global displacement of protein atoms on the proximal side of the heme observed by X-ray crystallography. We suggest that the conformational change observed in our experiments corresponds to the final portion of the nonexponential conformational relaxation recently observed by Anfinrud and co-workers, which begins on a picosecond time scale. Furthermore, extrapolation of our data to temperatures near that of the solvent glass transition suggests that this conformational relaxation may very well be the one postulated by Frauenfelder and co-workers to explain the decrease in the rate of geminate rebinding with increasing temperature above 180 K.

Animals↗

Fast events in protein folding initiated by nanosecond laser photolysis.

Initiation of protein folding by light can dramatically improve the time resolution of kinetic studies. Here we present an example of an optically triggered folding reaction by using nanosecond photodissociation of the heme-carbon monoxide complex of reduced cytochrome c. The optical trigger is based on the observation that under destabilizing conditions cytochrome c can be unfolded by preferential binding of carbon monoxide to the covalently attached heme group in the unfolded state. Photodissociation of the carbon monoxide thus triggers the folding reaction. We used time-resolved absorption spectroscopy to monitor binding at the heme. Before folding begins we observe transient binding of both nonnative and native ligands from the unfolded polypeptide on a microsecond time scale. Kinetic modeling suggests that the intramolecular binding of methionine-65 and -80 is faster than that of histidine-26 and -33, even though the histidines are closer to the heme. This optical trigger should provide a powerful method for studying chain collapse and secondary structure formation in cytochrome c without any limitations in time resolution.

Amino Acid Sequence↗

Simulation of the kinetics of ligand binding to a protein by molecular dynamics: geminate rebinding of nitric oxide to myoglobin.

We have begun to use molecular dynamics to simulate the kinetics of nitric oxide rebinding to myoglobin after photodissociation. Rebinding was simulated using a potential function that switches smoothly between a nonbinding potential and a binding potential as a function of the position and orientation of the ligand, with no barrier arising from the crossing of potential surfaces of different electron spin. In 96 of 100 trajectories, the ligand rebound in < 15 ps. The kinetic progress curve was obtained by determining the time in each trajectory at which the ligand rebound and then calculating the fraction of unbound ligands as a function of time. The curve can be well reproduced by a simple model based on the dynamics of a Langevin particle moving on a one-dimensional potential of mean force calculated from nonreactive protein trajectories. The rate of escape from the energy well adjacent to the heme is in good agreement with the value calculated from experimental data, suggesting that a multiple-well model provides a plausible explanation for the nonexponential rebinding kinetics. A transition-state analysis suggests that protein conformational relaxation coupled to the displacement of the iron from the heme plane is an unlikely cause for the nonexponential rebinding of nitric oxide.

Computer Simulation↗

Oxygen binding by single crystals of hemoglobin.

Reversible oxygen binding curves for single crystals of hemoglobin in the T quaternary structure have been measured using microspectrophotometry. Saturations were determined from complete visible spectra measured with light linearly polarized parallel to the a and c crystal axes. Striking differences were observed between the binding properties of hemoglobin in the crystal and those of hemoglobin in solution. Oxygen binding to the crystal is effectively noncooperative, the Bohr effect is absent, and there is no effect of chloride ion. Also, the oxygen affinity is lower than that of the T quaternary structure in solution. The absence of the Bohr effect supports Perutz's hypothesis on the key role of the salt bridges, which are known from X-ray crystallography to remain intact upon oxygenation. The low affinity and absence of the Bohr effect can be explained by a generalization of the MWC-PSK model (Monod, Wyman, & Changeux, 1965; Perutz, 1970; Szabo & Karplus, 1972) in which both high- and low-affinity tertiary conformations, with broken and unbroken salt bridges, respectively, are populated in the T quaternary structure. Because the alpha and beta hemes make different projections onto the two crystal axes, separate binding curves for the alpha and beta subunits could be calculated from the two measured binding curves. The approximately 5-fold difference between the oxygen affinities of the alpha and beta subunits is much smaller than that predicted from the crystallographic study of Dodson, Liddington, and co-workers, which suggested that oxygen binds only to the alpha hemes.(ABSTRACT TRUNCATED AT 250 WORDS)

Crystallization↗

Photoselection in polarized photolysis experiments on heme proteins.

Polarized photolysis experiments have been performed on the carbon monoxide complex of myoglobin to assess the effects of photoselection on the kinetics of ligand rebinding and to investigate the reorientational dynamics of the heme plane. The results are analyzed in terms of the optical theory developed in the preceding paper by Ansari and Szabo. Changes in optical density arising from rotational diffusion of the photoselected population produce large deviations from the true geminate ligand rebinding curves if measurements are made with only a single polarization. The apparent ligand rebinding curves are significantly distorted even at photolysis levels greater than 90%. These deviations are eliminated by obtaining isotropically-averaged optical densities from measurements using both parallel and perpendicular polarizations of the probe pulse. These experiments also yield the optical anisotropy, which gives a novel method for accurately determining the degree of photolysis, as well as important information on the reorientational dynamics of the heme plane. The correlation time for the overall rotational diffusion of the molecule is obtained from the decay of the anisotropy. The anisotropy prior to rotational diffusion is lower than that predicted for a rigidly attached, perfectly circular absorber, corresponding to an apparent order parameter of S = 0.95 +/- 0.02. Polarized absorption data on single crystals suggest that the decreased anisotropy results more from internal motions of the heme plane which take place on time scales shorter than the duration of the laser pulse (10 ns) than from out-of-plane polarized transitions.

Animals↗

Molecular dynamics simulations of heme reorientational motions in myoglobin.

Molecular dynamics simulations of 2-ns duration were performed on carbonmonoxymyoglobin and deoxymyoglobin in vacuo to study the reorientational dynamics of the heme group. The heme in both simulations undergoes reorientations of approximately 5 degrees amplitude on a subpicosecond time scale, which produce a rapid initial decay in the reorientational correlation function to about 0.99. The heme also experiences infrequent changes in average orientation of approximately 10 degrees amplitude, which lead to a larger slow decay of the reorientational correlation function over a period of hundreds of picoseconds. The simulations have not converged with respect to these infrequent transitions. However, an estimate of the order parameter for rapid internal motions of the heme from those orientations which are sampled by the simulations suggests that the subnanosecond orientational dynamics of the heme accounts for at least 30% of the unresolved initial anisotropy decay observed in the nanosecond time-resolved optical absorption experiments on myoglobin reported by Ansari et al. in a companion paper (Ansari, A., C.M. Jones, E.R. Henry, J. Hofrichter, and W.A. Eaton. 1992. Biophys. J. 64:852-868.). A more complete sampling of the accessible heme orientations would most likely increase this fraction further. The simulation of the liganded molecule also suggests that the conformational dynamics of the CO ligand may contribute significantly to discrepancies between the ligand conformation as probed by x-ray diffraction and by infrared-optical photoselection experiments. The protein back-bone explores multiple conformations during the simulations, with the largest structural changes appearing in the E and F helices, which are in contact with the heme. The variations in the heme orientation correlate with the conformational dynamics of the protein on a time scale of hundreds of picoseconds, suggesting that the heme orientation may provide a useful probe of dynamical processes in the protein.

Animals↗

Speed of intersubunit communication in proteins.

To determine the speed of communication between protein subunits, time-resolved absorption spectra were measured following partial photodissociation of the carbon monoxide complex of hemoglobin. The experiments were carried out using linearly polarized, 10-ns laser pulses, with the polarization of the excitation pulse both parallel and perpendicular to the polarization of the probe pulse. The substantial contribution to the observed spectra from photoselection effects was eliminated by isotropically averaging the polarized spectra, allowing a detailed comparison of the kinetics as a function of the degree of photolysis. These results show that prior to 1 microsecond both geminate ligand rebinding and conformational relaxation are independent of the number of ligands dissociated from the hemoglobin tetramer, as expected for a two-state allosteric model. After this time the kinetics depend on the ligation state of the tetramer. The conformational relaxation at 10 microseconds can be interpreted in terms of the two-state allosteric model as arising from the R to T quaternary conformational change of both unliganded and singly liganded molecules. These results suggest that communication between subunits requires about 1 microsecond and that the mechanism of the communication which occurs after this time is via the R to T conformational change. The optical anisotropy provides a novel means of accurately determining the extinction coefficients of the transient photoproduct. The decay in the optical anisotropy, moreover, provides an accurate determination of the rotational correlation time of 36 +/- 3 ns.

Carboxyhemoglobin↗

The role of solvent viscosity in the dynamics of protein conformational changes.

Nanosecond lasers were used to measure the rate of conformational changes in myoglobin after ligand dissociation at ambient temperatures. At low solvent viscosities the rate is independent of viscosity, but at high viscosities it depends on approximately the inverse first power of the viscosity. Kramers theory for unimolecular rate processes can be used to explain this result if the friction term is modified to include protein as well as solvent friction. The theory and experiment suggest that the dominant factor in markedly reducing the rate of conformational changes in myoglobin at low temperatures (less than 200 K) is the very high viscosity (greater than 10(7) centipoise) of the glycerol-water solvent. That is, at low temperatures conformational substates may not be "frozen" so much as "stuck."

Carbon Monoxide↗