Direct observation of photolysis-induced tertiary structural changes in hemoglobin

Role of hemoglobin structural-functional relationships in oxygen transport

The molecular mechanism of O₂ binding to hemoglobin (Hb) has been critically reviewed on the basis of the information built up in the last decades. It allows to describe in detail from the kinetic and thermodynamic viewpoint the process of O₂ uptake in the lungs and release to the tissues, casting some light on the physiological and pathological aspects of this process. The relevance of structural-functional relationships for O₂ binding is particularly outlined in the case of poorly vascularized tissues, such as retina, briefly discussing of strategies employed for optimization of oxygen supply to this type of tissues.

Effect of the abolition of intersubunit salt bridges on allosteric protein structural dynamics

A salt bridge, one of the representative structural factors established by non-covalent interactions, plays a crucial role in stabilizing the structure and regulating the protein function, but its role in dynamic processes has been elusive. Here, to scrutinize the structural and functional roles of the salt bridge in the process of performing the protein function, we investigated the effects of salt bridges on the allosteric structural transition of homodimeric hemoglobin (HbI) by applying time-resolved X-ray solution scattering (TRXSS) to the K30D mutant, in which the interfacial salt bridges of the wild type (WT) are abolished. The TRXSS data of K30D are consistent with the kinetic model that requires one monomer intermediate in addition to three structurally distinct dimer intermediates (I1, I2, and I3) observed in WT and other mutants. The kinetic and structural analyses show that K30D has an accelerated biphasic transition from I2 to I3 by more than nine times compared to WT and lacks significant structural changes in the transition from R-like I2 to T-like I3 observed in WT, unveiling that the loss of the salt bridges interrupts the R-T allosteric transition of HbI. Besides, the correlation between the bimolecular CO recombination rates in K30D, WT, and other mutants reveals that the bimolecular CO recombination is abnormally decelerated in K30D, indicating that the salt bridges also affect the cooperative ligand binding in HbI. These comparisons of the structural dynamics and kinetics of K30D and WT show that the interfacial salt bridges not only assist the physical connection of two subunits but also play a critical role in the global structural signal transduction of one subunit to the other subunit via a series of well-organized structural transitions.

Haemoglobin(βK120C)-albumin trimer as an artificial O2 carrier with sufficient haemoglobin allostery

The allosteric O₂ release of haemoglobin (Hb) allows for efficient O₂ delivery from the lungs to the tissues. However, allostery is weakened in Hb-based O₂ carriers because the chemical modifications of the Lys- and Cys-β93 residues prevent the quaternary transition of Hb. In this paper, we describe the synthesis and O₂ binding properties of a recombinant Hb [rHb(βK120C)]–albumin heterotrimer that maintains sufficient Hb allostery. The rHb(βK120C) core, with two additional cysteine residues at the symmetrical positions on its protein surface, was expressed using yeast cells. The mutations did not influence either the O₂ binding characteristics or the quaternary transition of Hb. Maleimide-activated human serum albumins (HSAs) were coupled with rHb(βK120C) at the two Cys-β120 positions, yielding the rHb(βK120C)–HSA₂ trimer, in which the Cys-β93 residues were unreacted. Molecular dynamics simulation demonstrated that the HSA moiety does not interact with the amino acid residues around the haem pockets and the α₁β₂ surfaces of the rHb(βK120C) core, the alteration of which retards Hb allostery. Circular dichroism spectroscopy demonstrated that the quaternary transition between the relaxed (R) state and the tense (T) state of the Hb core occurred upon both the association and dissociation of O₂. In phosphate-buffered saline solution (pH 7.4) at 37 °C, the rHb(βK120C)–HSA₂ trimer exhibited a sigmoidal O₂ equilibrium curve with the O₂ affinity and cooperativity identical to those of native Hb (p₅₀ = 12 Torr, n = 2.4). Moreover, we observed an equal Bohr effect and 2,3-diphosphoglycerate response in the rHb(βK120C)–HSA₂ trimer compared with naked Hb.

Recombinant haemoglobin [rHb(βK120C)] was coupled with two human serum albumins (HSAs), yielding a rHb(βK120C)–HSA₂ heterotrimer, which shows a sigmoidal O₂ equilibrium curve and sufficient Hb allostery identical to those of native Hb.

Effect of Occluded Ligand Migration on the Kinetics and Structural Dynamics of Homodimeric Hemoglobin

Small molecules such as molecular oxygen, nitric oxide, and carbon monoxide play important roles in life, and many proteins require the transport of small molecules to and from the bulk solvent for their function. Ligand migration within a protein molecule is expected to be closely related to the overall structural changes of the protein, but the detailed and quantitative connection remains elusive. For example, despite numerous studies, how occluded ligand migration affects the kinetics and structural dynamics of the R-T transition remains unclear. To shed light on this issue, we chose homodimeric hemoglobin (HbI) with the I114F mutation (I114F), which is known to interfere with ligand migration between the primary and secondary docking sites, and studied its kinetics and structural dynamics using time-resolved X-ray solution scattering. The kinetic analysis shows that I114F has three structurally distinct intermediates (I₁, I₂, and I₃) as in the wild type (WT), but its geminate CO recombination occurs directly from I₁ without the path via I₂ observed in WT. Moreover, the structural transitions, which involve ligand migration (the transitions from I₁ to I₂ and from I₃ to the initial state), are decelerated compared to WT. The structural analysis revealed that I114F involves generally smaller structural changes in all three intermediates compared to WT.

Direct observation of ligand migration within human hemoglobin at work

Human hemoglobin is the textbook example of the stereochemistry of an allosteric protein and of the exquisite control that a protein can exert over ligand binding. However, the fundamental basis by which the protein facilitates the ligand movement remains unknown. In this study, we used cryogenic X-ray crystallography and a high-repetition pulsed laser irradiation technique to elucidate the atomic details of ligand migration processes in hemoglobin after photolysis of the bound CO. Our data clarify the distinct CO migration pathways in the individual subunits of hemoglobin and unravel the functional roles of the internal cavities and neighboring amino acid residues in ligand exit and entry. Our results also demonstrate the high gas permeability and porosity of hemoglobin, facilitating O₂ delivery.

Hemoglobin is one of the best-characterized proteins with respect to structure and function, but the internal ligand diffusion pathways remain obscure and controversial. Here we captured the CO migration processes in the tense (T), relaxed (R), and second relaxed (R2) quaternary structures of human hemoglobin by crystallography using a high-repetition pulsed laser technique at cryogenic temperatures. We found that in each quaternary structure, the photodissociated CO molecules migrate along distinct pathways in the α and β subunits by hopping between the internal cavities with correlated side chain motions of large nonpolar residues, such as α14Trp(A12), α105Leu(G12), β15Trp(A12), and β71Phe(E15). We also observe electron density evidence for the distal histidine [α58/β63His(E7)] swing-out motion regardless of the quaternary structure, although less evident in α subunits than in β subunits, suggesting that some CO molecules have escaped directly through the E7 gate. Remarkably, in T-state Fe(II)-Ni(II) hybrid hemoglobins in which either the α or β subunits contain Ni(II) heme that cannot bind CO, the photodissociated CO molecules not only dock at the cavities in the original Fe(II) subunit, but also escape from the protein matrix and enter the cavities in the adjacent Ni(II) subunit even at 95 K, demonstrating the high gas permeability and porosity of the hemoglobin molecule. Our results provide a comprehensive picture of ligand movements in hemoglobin and highlight the relevance of cavities, nonpolar residues, and distal histidines in facilitating the ligand migration.

Ligand pathways in neuroglobin revealed by low-temperature photodissociation and docking experiments

A combined biophysical approach was applied to map gas-docking sites within murine neuroglobin (Ngb), revealing snapshots of events that might govern activity and dynamics in this unique hexacoordinate globin, which is most likely to be involved in gas-sensing in the central nervous system and for which a precise mechanism of action remains to be elucidated. The application of UV-visible microspectroscopy in crystallo, solution X-ray absorption near-edge spectroscopy and X-ray diffraction experiments at 15-40 K provided the structural characterization of an Ngb photolytic intermediate by cryo-trapping and allowed direct observation of the relocation of carbon monoxide within the distal heme pocket after photodissociation. Moreover, X-ray diffraction at 100 K under a high pressure of dioxygen, a physiological ligand of Ngb, unravelled the existence of a storage site for O2 in Ngb which coincides with Xe-III, a previously described docking site for xenon or krypton. Notably, no other secondary sites were observed under our experimental conditions.

Protein Structural Dynamics of Wild-Type and Mutant Homodimeric Hemoglobin Studied by Time-Resolved X-Ray Solution Scattering

The quaternary transition between the relaxed (R) and tense (T) states of heme-binding proteins is a textbook example for the allosteric structural transition. Homodimeric hemoglobin (HbI) from Scapharca inaequivalvis is a useful model system for investigating the allosteric behavior because of the relatively simple quaternary structure. To understand the cooperative transition of HbI, wild-type and mutants of HbI have been studied by using time-resolved X-ray solution scattering (TRXSS), which is sensitive to the conformational changes. Herein, we review the structural dynamics of HbI investigated by TRXSS and compare the results of TRXSS with those of other techniques.

Pressure Effect on Zero-Field Splitting Parameter of Hemin: Model Case of Hemoproteins under Pressure

We experimentally studied the pressure dependence of the zero-field splitting (ZFS) parameter of hemin (iron(III) protoporphyrin IX chloride), which is a model complex of hemoproteins, via high-frequency and high-field electron paramagnetic resonance (HFEPR) under pressure. Owing to the large ZFS, the pressure effect on the electronic structure of iron-porphyrin complexes has not yet been explored using EPR. Therefore, we systematically studied this effect using our newly developed sub-terahertz EPR spectroscopy system in the frequency range of 80-515 GHz, under magnetic fields up to 10 T and pressure up to 2 GPa. We observed a systematic shift of the resonance fields of hemin upon pressure application, from which the axial component of the ZFS parameter was found to increase from D = 6.9 to 7.9 cm^-1 at 2 GPa. In contrast to the previous methods used to study proteins under pressure, which mainly focused on conformational changes, our HFEPR technique can obtain more microscopic insights into the electronic structures of metal ions under pressure. In this sense, our technique provides novel opportunities to study the pressure effects on biofunctional active centers of versatile metalloproteins.

Size and Shape Controlled Crystallization of Hemoglobin for Advanced Crystallography

While high-throughput screening for protein crystallization conditions have rapidly evolved in the last few decades, it is also becoming increasingly necessary for the control of crystal size and shape as increasing diversity of protein crystallographic experiments. For example, X-ray crystallography (XRC) combined with photoexcitation and/or spectrophotometry requires optically thin but well diffracting crystals. By contrast, large-volume crystals are needed for weak signal experiments, such as neutron crystallography (NC) or recently developed X-ray fluorescent holography (XFH). In this article, we present, using hemoglobin as an example protein, some techniques for obtaining the crystals of controlled size, shape, and adequate quality. Furthermore, we describe a few case studies of applications of the optimized hemoglobin crystals for implementing the above mentioned crystallographic experiments, providing some hints and tips for the further progress of advanced protein crystallography.

Alteration of the α1β2/α2β1 subunit interface contributes to the increased hemoglobin-oxygen affinity of high-altitude deer mice

BACKGROUND

Deer mice (Peromyscus maniculatus) that are native to high altitudes in the Rocky Mountains have evolved hemoglobins with an increased oxygen-binding affinity relative to those of lowland conspecifics. To elucidate the molecular mechanisms responsible for the evolved increase in hemoglobin-oxygen affinity, the crystal structure of the highland hemoglobin variant was solved and compared with the previously reported structure for the lowland variant.

RESULTS

Highland hemoglobin yielded at least two crystal types, in which the longest axes were 507 and 230 Å. Using the smaller unit cell crystal, the structure was solved at 2.2 Å resolution. The asymmetric unit contained two tetrameric hemoglobin molecules.

CONCLUSIONS

The analyses revealed that αPro50 in the highland hemoglobin variant promoted a stable interaction between αHis45 and heme that was not seen in the αHis50 lowland variant. The αPro50 mutation also altered the nature of atomic contacts at the α1β2/α2β1 intersubunit interfaces. These results demonstrate how affinity-altering changes in intersubunit interactions can be produced by mutations at structurally remote sites.

Cooperative protein structural dynamics of homodimeric hemoglobin linked to water cluster at subunit interface revealed by time-resolved X-ray solution scattering

Homodimeric hemoglobin (HbI) consisting of two subunits is a good model system for investigating the allosteric structural transition as it exhibits cooperativity in ligand binding. In this work, as an effort to extend our previous study on wild-type and F97Y mutant HbI, we investigate structural dynamics of a mutant HbI in solution to examine the role of well-organized interfacial water cluster, which has been known to mediate intersubunit communication in HbI. In the T72V mutant of HbI, the interfacial water cluster in the T state is perturbed due to the lack of Thr72, resulting in two less interfacial water molecules than in wild-type HbI. By performing picosecond time-resolved X-ray solution scattering experiment and kinetic analysis on the T72V mutant, we identify three structurally distinct intermediates (I1, I2, and I3) and show that the kinetics of the T72V mutant are well described by the same kinetic model used for wild-type and F97Y HbI, which involves biphasic kinetics, geminate recombination, and bimolecular CO recombination. The optimized kinetic model shows that the R-T transition and bimolecular CO recombination are faster in the T72V mutant than in the wild type. From structural analysis using species-associated difference scattering curves for the intermediates, we find that the T-like deoxy I3 intermediate in solution has a different structure from deoxy HbI in crystal. In addition, we extract detailed structural parameters of the intermediates such as E-F distance, intersubunit rotation angle, and heme-heme distance. By comparing the structures of protein intermediates in wild-type HbI and the T72V mutant, we reveal how the perturbation in the interfacial water cluster affects the kinetics and structures of reaction intermediates of HbI.

O2 and Water Migration Pathways between the Solvent and Heme Pockets of Hemoglobin with Open and Closed Conformations of the Distal HisE7

Hemoglobin transports O2 by binding the gas at its four hemes. Hydrogen bonding between the distal histidine (HisE7) and heme-bound O2 significantly increases the affinity of human hemoglobin (HbA) for this ligand. HisE7 is also proposed to regulate the release of O2 to the solvent via a transient E7 channel. To reveal the O2 escape routes controlled by HisE7 and to evaluate its role in gating heme access, we compare simulations of O2 diffusion from the distal heme pockets of the T and R states of HbA performed with HisE7 in its open (protonated) and closed (neutral) conformations. Irrespective of HisE7's conformation, we observe the same four or five escape routes leading directly from the α- or β-distal heme pockets to the solvent. Only 21-53% of O2 escapes occur via these routes, with the remainder escaping through routes that encompass multiple internal cavities in HbA. The conformation of the distal HisE7 controls the escape of O2 from the heme by altering the distal pocket architecture in a pH-dependent manner, not by gating the E7 channel. Removal of the HisE7 side chain in the GlyE7 variant exposes the distal pockets to the solvent, and the percentage of O2 escapes to the solvent directly from the α- or β-distal pockets of the mutant increases to 70-88%. In contrast to O2, the dominant water route from the bulk solvent is gated by HisE7 because protonation and opening of this residue dramatically increase the rate of influx of water into the empty distal heme pockets. The occupancy of the distal heme site by a water molecule, which functions as an additional nonprotein barrier to binding of the ligand to the heme, is also controlled by HisE7. Overall, analysis of gas and water diffusion routes in the subunits of HbA and its GlyE7 variant sheds light on the contribution of distal HisE7 in controlling polar and nonpolar ligand movement between the solvent and the hemes.

Reverse engineering the cooperative machinery of human hemoglobin

Hemoglobin transports molecular oxygen from the lungs to all human tissues for cellular respiration. Its α₂β₂ tetrameric assembly undergoes cooperative binding and releasing of oxygen for superior efficiency and responsiveness. Over past decades, hundreds of hemoglobin structures were determined under a wide range of conditions for investigation of molecular mechanism of cooperativity. Based on a joint analysis of hemoglobin structures in the Protein Data Bank (Ren, companion article), here I present a reverse engineering approach to elucidate how two subunits within each dimer reciprocate identical motions that achieves intradimer cooperativity, how ligand-induced structural signals from two subunits are integrated to drive quaternary rotation, and how the structural environment at the oxygen binding sites alter their binding affinity. This mechanical model reveals the intricate design that achieves the cooperative mechanism and has previously been masked by inconsistent structural fluctuations. A number of competing theories on hemoglobin cooperativity and broader protein allostery are reconciled and unified.

Reaction trajectory revealed by a joint analysis of protein data bank

Structural motions along a reaction pathway hold the secret about how a biological macromolecule functions. If each static structure were considered as a snapshot of the protein molecule in action, a large collection of structures would constitute a multidimensional conformational space of an enormous size. Here I present a joint analysis of hundreds of known structures of human hemoglobin in the Protein Data Bank. By applying singular value decomposition to distance matrices of these structures, I demonstrate that this large collection of structural snapshots, derived under a wide range of experimental conditions, arrange orderly along a reaction pathway. The structural motions along this extensive trajectory, including several helical transformations, arrive at a reverse engineered mechanism of the cooperative machinery (Ren, companion article), and shed light on pathological properties of the abnormal homotetrameric hemoglobins from α-thalassemia. This method of meta-analysis provides a general approach to structural dynamics based on static protein structures in this post genomics era.

Real-time tracking of CO migration and binding in the α and β subunits of human hemoglobin via 150-ps time-resolved Laue crystallography

We have developed the method of picosecond Laue crystallography and used this capability to probe ligand dynamics in tetrameric R-state hemoglobin (Hb). Time-resolved, 2 Å-resolution electron density maps of photolyzed HbCO reveal the time-dependent population of CO in the binding (A) and primary docking (B) sites of both α and β subunits from 100 ps to 10 μs. The proximity of the B site in the β subunit is about 0.25 Å closer to its A binding site, and its k_BA rebinding rate (~300 μs^-1) is six times faster, suggesting distal control of the rebinding dynamics. Geminate rebinding in the β subunit exhibits both prompt and delayed geminate phases. We developed a microscopic model to quantitatively explain the observed kinetics, with three states for the α subunit and four states for the β subunit. This model provides a consistent framework for interpreting rebinding kinetics reported in prior studies of both HbCO and HbO₂.

Deer mouse hemoglobin exhibits a lowered oxygen affinity owing to mobility of the E helix

The deer mouse, Peromyscus maniculatus, exhibits altitude-associated variation in hemoglobin oxygen affinity. To examine the structural basis of this functional variation, the structure of the hemoglobin was solved. Recombinant hemoglobin was expressed in Escherichia coli and was purified by ion-exchange chromatography. Recombinant hemoglobin was crystallized by the hanging-drop vapor-diffusion method using polyethylene glycol as a precipitant. The obtained orthorhombic crystal contained two subunits in the asymmetric unit. The refined structure was interpreted as the aquo-met form. Structural comparisons were performed among hemoglobins from deer mouse, house mouse and human. In contrast to human hemoglobin, deer mouse hemoglobin lacks the hydrogen bond between α1Trp14 in the A helix and α1Thr67 in the E helix owing to the Thr67Ala substitution. In addition, deer mouse hemoglobin has a unique hydrogen bond at the α1β1 interface between residues α1Cys34 and β1Ser128.

Coarse-grained and all-atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen-binding protein composed of four subunits (α1, α2, β1, and β2), is a well-known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse-grained (CG) modeling, all-atom simulation, and structural data from X-ray crystallography and wide-angle X-ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot-spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all-atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all-atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling.

Discrimination between CO and O(2) in heme oxygenase: comparison of static structures and dynamic conformation changes following CO photolysis

Heme oxygenase (HO) catalyzes heme degradation, one of its products being carbon monoxide (CO). It is well known that CO has a higher affinity for heme iron than does molecular oxygen (O(2)); therefore, CO is potentially toxic. Because O(2) is required for the HO reaction, HO must discriminate effectively between CO and O(2) and thus escape product inhibition. Previously, we demonstrated large conformational changes in the heme-HO-1 complex upon CO binding that arise from steric hindrance between CO bound to the heme iron and Gly-139. However, we have not yet identified those changes that are specific to CO binding and do not occur upon O(2) binding. Here we determine the crystal structure of the O(2)-bound form at 1.8 Å resolution and reveal the structural changes that are specific to CO binding. Moreover, difference Fourier maps comparing the structures before and after CO photolysis at <160 K clearly show structural changes such as movement of the distal F-helix upon CO photolysis. No such changes are observed upon O(2) photolysis, consistent with the structures of the ligand-free, O(2)-bound, and CO-bound forms. Protein motions even at cryogenic temperatures imply that the CO-bound heme-HO-1 complex is severely constrained (as in ligand binding to the T-state of hemoglobin), indicating that CO binding to the heme-HO-1 complex is specifically inhibited by steric hindrance. The difference Fourier maps also suggest new routes for CO migration.

Direct observation of cooperative protein structural dynamics of homodimeric hemoglobin from 100 ps to 10 ms with pump-probe X-ray solution scattering

Proteins serve as molecular machines in performing their biological functions, but the detailed structural transitions are difficult to observe in their native aqueous environments in real time. For example, despite extensive studies, the solution-phase structures of the intermediates along the allosteric pathways for the transitions between the relaxed (R) and tense (T) forms have been elusive. In this work, we employed picosecond X-ray solution scattering and novel structural analysis to track the details of the structural dynamics of wild-type homodimeric hemoglobin (HbI) from the clam Scapharca inaequivalvis and its F97Y mutant over a wide time range from 100 ps to 56.2 ms. From kinetic analysis of the measured time-resolved X-ray solution scattering data, we identified three structurally distinct intermediates (I(1), I(2), and I(3)) and their kinetic pathways common for both the wild type and the mutant. The data revealed that the singly liganded and unliganded forms of each intermediate share the same structure, providing direct evidence that the ligand photolysis of only a single subunit induces the same structural change as the complete photolysis of both subunits does. In addition, by applying novel structural analysis to the scattering data, we elucidated the detailed structural changes in the protein, including changes in the heme-heme distance, the quaternary rotation angle of subunits, and interfacial water gain/loss. The earliest, R-like I(1) intermediate is generated within 100 ps and transforms to the R-like I(2) intermediate with a time constant of 3.2 ± 0.2 ns. Subsequently, the late, T-like I(3) intermediate is formed via subunit rotation, a decrease in the heme-heme distance, and substantial gain of interfacial water and exhibits ligation-dependent formation kinetics with time constants of 730 ± 120 ns for the fully photolyzed form and 5.6 ± 0.8 μs for the partially photolyzed form. For the mutant, the overall kinetics are accelerated, and the formation of the T-like I(3) intermediate involves interfacial water loss (instead of water entry) and lacks the contraction of the heme-heme distance, thus underscoring the dramatic effect of the F97Y mutation. The ability to keep track of the detailed movements of the protein in aqueous solution in real time provides new insights into the protein structural dynamics.

Heme reactivity is uncoupled from quaternary structure in gel-encapsulated hemoglobin: a resonance Raman spectroscopic study

Encapsulation of hemoglobin (Hb) in silica gel preserves structure and function but greatly slows protein motion, thereby providing access to intermediates along the allosteric pathway that are inaccessible in solution. Resonance Raman (RR) spectroscopy with visible and ultraviolet laser excitation provides probes of heme reactivity and of key tertiary and quaternary contacts. These probes were monitored in gels after deoxygenation of oxyHb and after CO binding to deoxyHb, which initiate conformational change in the R-T and T-R directions, respectively. The spectra establish that quaternary structure change in the gel takes a week or more but that the evolution of heme reactivity, as monitored by the Fe-histidine stretching vibration, ν(FeHis), is completed within two days, and is therefore uncoupled from the quaternary structure. Within each quaternary structure, the evolving ν(FeHis) frequencies span the full range of values between those previously associated with the high- and low-affinity end states, R and T. This result supports the tertiary two-state (TTS) model, in which the Hb subunits can adopt high- and low-affinity tertiary structures, r and t, within each quaternary state. The spectra also reveal different tertiary pathways, involving the breaking and reformation of E and F interhelical contacts in the R-T direction but not the T-R direction. In the latter, tertiary motions are restricted by the T quaternary contacts.

An atomistic view on human hemoglobin carbon monoxide migration processes

A significant amount of work has been devoted to obtaining a detailed atomistic knowledge of the human hemoglobin mechanism. Despite this impressive research, to date, the ligand diffusion processes remain unclear and controversial. Using recently developed computational techniques, PELE, we are capable of addressing the ligand migration processes. First, the methodology was tested on myoglobin's CO migration, and the results were compared with the wealth of theoretical and experimental studies. Then, we explored both hemoglobin tense and relaxed states and identified the differences between the α-and β-subunits. Our results indicate that the proximal site, equivalent to the Xe1 cavity in myoglobin, is never visited. Furthermore, strategically positioned residues alter the diffusion processes within hemoglobin's subunits and suggest that multiple pathways exist, especially diversified in the α-globins. A significant dependency of the ligand dynamics on the tertiary structure is also observed.

Modulating distal cavities in the α and β subunits of human HbA reveals the primary ligand migration pathway

The free volume in the active site of human HbA plays a crucial role in governing the bimolecular rates of O(2), CO, and NO binding, the fraction of geminate ligand recombination, and the rate of NO dioxygenation by the oxygenated complex. We have decreased the size of the distal pocket by mutating Leu(B10), Val(E11), and Leu(G8) to Phe and Trp and that of other more internal cavities by filling them with Xe at high gas pressures. Increasing the size of the B10 side chain reduces bimolecular rates of ligand binding nearly 5000-fold and inhibits CO geminate recombination due to both reduction of the capture volume in the distal pocket and direct steric hindrance of Fe-ligand bond formation. Phe and Trp(E11) mutations also cause a decrease in distal pocket volume but, at the same time, increase access to the Fe atom because of the loss of the γ2 CH(3) group of the native Val(E11) side chain. The net result of these E11 substitutions is a dramatic increase in the rate of geminate recombination because dissociated CO is sequestered close to the Fe atom and can rapidly rebind without steric resistance. However, the bimolecular rate constants for binding of ligand to the Phe and Trp(E11) mutants are decreased 5-30-fold, because of a smaller capture volume. Geminate and bimolecular kinetic parameters for Phe and Trp(G8) mutants are similar to those for the native HbA subunits because the aromatic rings at this position cause little change in distal pocket volume and because ligands do not move past this position into the globin interior of wild-type HbA subunits. The latter conclusion is verified by the observation that Xe binding to the α and β Hb subunits has little effect on either geminate or bimolecular ligand rebinding. All of these experimental results argue strongly against alternative ligand migration pathways that involve movements through the protein interior in HbA. Instead, ligands appear to enter through the His(E7) gate and are captured directly in the distal cavity.

pH-dependent structural changes in haemoglobin component V from the midge larvaPropsilocerus akamusi(Orthocladiinae, Diptera)

Haemoglobin component V (Hb V) from the midge larvaPropsilocerus akamusiexhibits oxygen affinity despite the replacement of HisE7 and a pH-dependence of its functional properties. In order to understand the contribution of the distal residue to the ligand-binding properties and the pH-dependent structural changes in this insect Hb, the crystal structure of Hb V was determined under five different pH conditions. Structural comparisons of these Hb structures indicated that at neutral pH ArgE10 contributes to the stabilization of the haem-bound ligand molecule as a functional substitute for the nonpolar E7 residue. However, ArgE10 does not contribute to stabilization at acidic and alkaline pH because of the swinging movement of the Arg side chain under these conditions. This pH-dependent behaviour of Arg results in significant differences in the hydrogen-bond network on the distal side of the haem in the Hb V structures at different pH values. Furthermore, the change in pH results in a partial movement of the F helix, considering that coupled movements of ArgE10 and the F helix determine the haem location at each pH. These results suggested that Hb V retains its functional properties by adapting to the structural changes caused by amino-acid replacements.

Tracking ligand-migration pathways of carbonmonoxy myoglobin in crystals at cryogenic temperatures

In order to explore the ligand-migration dynamics in myoglobin induced by photodissociation, cryogenic X-ray crystallographic investigations of carbonmonoxy myoglobin crystals illuminated by continuous wave and pulsed lasers at 1–15 kHz repetition rate have been carried out. Here it is shown that this novel method, extended pulsed-laser pumping of carbonmonoxy myoglobin, promotes ligand migration in the protein matrix by crossing the glass transition temperature repeatedly, and enables the visualization of the migration pathway of the photodissociated ligands in native Mb at cryogenic temperatures. It has revealed that the migration of the CO molecule into each cavity induces structural changes of the amino-acid residues around the cavity which result in the expansion of the cavity. The sequential motion of the ligand and the cavity suggests a self-opening mechanism of the ligand-migration channel arising by induced fit.

Biological and biomedical applications of two-dimensional vibrational spectroscopy: proteomics, imaging, and structural analysis

In the last 10 years, several forms of two-dimensional infrared (2DIR) spectroscopy have been developed, such as IR pump-probe spectroscopy and photon-echo techniques. In this Account, we describe a doubly vibrationally enhanced four-wave mixing method, in which a third-order nonlinear signal is generated from the interaction of two independently tunable IR beams and an electron-polarizing visible beam at 790 nm. When the IR beams are independently in resonance with coupled vibrational transitions, the signal is enhanced and cross-peaks appear in the spectrum. This method is known as either DOVE (doubly vibrationally enhanced) four-wave mixing or EVV (electron-vibration-vibration) 2DIR spectroscopy. We begin by discussing the basis and properties of EVV 2DIR. We then discuss several biological and potential biomedical applications. These include protein identification and quantification, as well as the potential of this label-free spectroscopy for protein and peptide structural analysis. In proteomics, we also show how post-translational modifications in peptides (tyrosine phosphorylation) can be detected by EVV 2DIR spectroscopy. The feasibility of EVV 2DIR spectroscopy for tissue imaging is also evaluated. Preliminary results were obtained on a mouse kidney histological section that was stained with hematoxylin (a small organic molecule). We obtained images by setting the IR frequencies to a specific cross-peak (the strongest for hematoxylin was obtained from its analysis in isolation; a general CH(3) cross-peak for proteins was also used) and then spatially mapping as a function of the beam position relative to the sample. Protein and hematoxylin distribution in the tissue were measured and show differential contrast, which can be entirely explained by the different tissue structures and their functions. The possibility of triply resonant EVV 2DIR spectroscopy was investigated on the retinal chromophore at the centre of the photosynthetic protein bacteriorhodopsin (bR). By putting the visible third beam in resonance with an electronic transition, we were able to enhance the signal and increase the sensitivity of the method by several orders of magnitude. This increase in sensitivity is of great importance for biological applications, in which the number of proteins, metabolites, or drug molecules to be detected is low (typically pico- to femtomoles). Finally, we present theoretical investigations for using EVV 2DIR spectroscopy as a structural analysis tool for inter- and intramolecular interaction geometries.

Visualizing breathing motion of internal cavities in concert with ligand migration in myoglobin

Proteins harbor a number of cavities of relatively small volume. Although these packing defects are associated with the thermodynamic instability of the proteins, the cavities also play specific roles in controlling protein functions, e.g., ligand migration and binding. This issue has been extensively studied in a well-known protein, myoglobin (Mb). Mb reversibly binds gas ligands at the heme site buried in the protein matrix and possesses several internal cavities in which ligand molecules can reside. It is still an open question as to how a ligand finds its migration pathways between the internal cavities. Here, we report on the dynamic and sequential structural deformation of internal cavities during the ligand migration process in Mb. Our method, the continuous illumination of native carbonmonoxy Mb crystals with pulsed laser at cryogenic temperatures, has revealed that the migration of the CO molecule into each cavity induces structural changes of the amino acid residues around the cavity, which results in the expansion of the cavity with a breathing motion. The sequential motion of the ligand and the cavity suggests a self-opening mechanism of the ligand migration channel arising by induced fit, which is further supported by computational geometry analysis by the Delaunay tessellation method. This result suggests a crucial role of the breathing motion of internal cavities as a general mechanism of ligand migration in a protein matrix.

Ab initio determination of dark structures in radiationless transitions for aromatic carbonyl compounds

Mechanistic photodissociation of a polyatomic molecule has long been regarded as an intellectually challenging area of chemical physics, the results of which are relevant to atmospheric chemistry, biological systems, and many application fields. Carbonyl compounds play a unique role in the development of our understanding of the spectroscopy, photochemistry, and photophysics of polyatomic molecules and their photodissociation has been the subject of numerous studies over many decades. Upon irradiation, a molecule can undergo internal conversion (IC) and intersystem crossing (ISC) processes, besides photochemical and other photophysical processes. Transient intermediates formed in the IC and ISC radiationless processes, which are termed "dark", are not amenable to detection by conventional light absorption or emission. However, these dark intermediates play critical roles in IC and ISC processes and thus are essential to understanding mechanistic photochemistry of a polyatomic molecule. We have applied the multiconfiguration complete active space self-consistent field (CASSCF) method to determine the dark transient structures involved in radiationless processes for acetophenone and the related aromatic carbonyl compounds. The electronic and geometric structures predicted for the dark states are in a good agreement with those determined by ultrafast electron diffraction experiments. Intersection structure of different electronic states provides a very efficient "funnel" for the IC or ISC process. However, experimental determination of the intersection structure involved in radiationless transitions of a polyatomic molecule is impossible at present. We have discovered a minimum energy crossing point among the three potential energy surfaces (S1, T1, and T2) that appears to be common to a wide variety of aromatic carbonyl compounds with a constant structure. This new type of crossing point holds the key to understanding much about radiationless processes after photoexcitation of aromatic carbonyl compounds. The importance of ab initio determination of transient structures in the photodissociation dynamics has been demonstrated for the case of the aromatic carbonyl compounds. In addition, the detailed knowledge of mechanistic photochemistry for aromatic carbonyl compounds forms the basis for further investigating photodissociation dynamics of a polyatomic molecule.

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

A quantum-chemical picture of hemoglobin affinity

Understanding the molecular mechanism of hemoglobin cooperativity remains an enduring challenge. Protein forces that control ligand affinity are not directly accessible by experiment. We demonstrate that computational quantum mechanics/molecular mechanics methods can provide reasonable values of ligand binding energies in Hb, and of their dependence on allostery. About 40% of the binding energy differences between the relaxed state and tense state quaternary structures result from strain induced in the heme and its ligands, especially in one of the pyrrole rings. The proximal histidine also contributes significantly, in particular, in the alpha-chains. The remaining energy difference resides in protein contacts, involving residues responsible for locking the quaternary changes. In the alpha-chains, the most important contacts involve the FG corner, at the "hinge" region of the alpha(1)beta(2) quaternary interface. The energy differences are spread more evenly among the beta-chain residues, suggesting greater flexibility for the beta- than for the alpha-chains along the quaternary transition. Despite this chain differentiation, the chains contribute equally to the relaxed substitute state energy difference. Thus, nature has evolved a symmetric response to the quaternary structure change, which is a requirement for maximum cooperativity, via different mechanisms for the two kinds of chains.

Kappa-alpha plot derived structural alphabet and BLOSUM-like substitution matrix for rapid search of protein structure database

3D BLAST, a novel protein structure database search tool, is a useful tool for analysing novel structures, capable of returning a list of aligned structures ordered according to E-values.

We present a novel protein structure database search tool, 3D-BLAST, that is useful for analyzing novel structures and can return a ranked list of alignments. This tool has the features of BLAST (for example, robust statistical basis, and effective and reliable search capabilities) and employs a kappa-alpha (κ, α) plot derived structural alphabet and a new substitution matrix. 3D-BLAST searches more than 12,000 protein structures in 1.2 s and yields good results in zones with low sequence similarity.

Allosteric action in real time: time-resolved crystallographic studies of a cooperative dimeric hemoglobin

Protein allostery provides mechanisms for regulation of biological function at the molecular level. We present here an investigation of global, ligand-induced allosteric transition in a protein by time-resolved x-ray diffraction. The study provides a view of structural changes in single crystals of Scapharca dimeric hemoglobin as they proceed in real time, from 5 ns to 80 micros after ligand photodissociation. A tertiary intermediate structure forms rapidly (<5 ns) as the protein responds to the presence of an unliganded heme within each R-state protein subunit, with key structural changes observed in the heme groups, neighboring residues, and interface water molecules. This intermediate lays a foundation for the concerted tertiary and quaternary structural changes that occur on a microsecond time scale and are associated with the transition to a low-affinity T-state structure. Reversal of these changes shows a considerable lag as a T-like structure persists well after ligand rebinding, suggesting a slow T-to-R transition.

Advances in kinetic protein crystallography

Many proteins function in the crystalline state, making crystallography a tool that can address mechanism, as well as structure. By initiating biological turnover in the crystal, transient structural species form, which may be filmed by Laue diffraction or captured by freeze-trapping methods. Laue diffraction has now reached an unprecedented level of sophistication and has found a 'niche of excellence' in the study of cyclic, ultra-fast, light-triggered reactions. Trapping methods, on the other hand, are more generally applicable, but require care to avoid artifacts. New strategies have been developed and difficulties such as radiation damage have received particular attention. Complementary methods--mainly UV/visible single-crystal spectroscopy--have proven essential to design, interpret and validate kinetic crystallography experiments.

Structure of an extracellular giant hemoglobin of the gutless beard worm Oligobrachia mashikoi

Mouthless and gutless marine animals, pogonophorans and vestimentiferans, obtain their nutrition solely from their symbiotic chemoautotrophic sulfur-oxidizing bacteria. These animals have sulfide-binding 400-kDa and/or 3,500-kDa Hb, which transports oxygen and sulfide simultaneously. The symbiotic bacteria are supplied with sulfide by Hb of the host animal and use it to provide carbon compounds. Here, we report the crystal structure of a 400-kDa Hb from pogonophoran Oligobrachia mashikoi at 2.85-A resolution. The structure is hollow-spherical, composed of a total of 24 globins as a dimer of dodecamer. This dodecameric assemblage would be a fundamental structural unit of both 400-kDa and 3,500-kDa Hbs. The structure of the mercury derivative used for phasing provides insights into the sulfide-binding mechanism. The mercury compounds bound to all free Cys residues that have been expected as sulfide-binding sites. Some of the free Cys residues are surrounded by Phe aromatic rings, and mercury atoms come into contact with these residues in the derivative structure. It is strongly suggested that sulfur atoms bound to these sites could be stabilized by aromatic-electrostatic interactions by the surrounding Phe residues.

CO-trapping site in heme oxygenase revealed by photolysis of its co-bound heme complex: mechanism of escaping from product inhibition

Heme oxygenase (HO) catalyzes physiological heme degradation using O(2) and reducing equivalents to produce biliverdin, iron, and CO. Notably, the HO reaction proceeds without product inhibition by CO, which is generated in the conversion reaction of alpha-hydroxyheme to verdoheme, although CO is known to be a potent inhibitor of HO and other heme proteins. In order to probe how endogenous CO is released from the reaction site, we collected X-ray diffraction data from a crystal of the CO-bound form of the ferrous heme-HO complex in the dark and under illumination by a red laser at approximately 35 K. The difference Fourier map indicates that the CO ligand is partially photodissociated from the heme and that the photolyzed CO is trapped in a hydrophobic cavity adjacent to the heme pocket. This hydrophobic cavity was occupied also by xenon, which is similar to CO in terms of size and properties. Taking account of the affinity of CO for the ferrous verdoheme-HO complex being much weaker than that for the ferrous heme complex, the CO derived from alpha-hydroxyheme would be trapped preferentially in the hydrophobic cavity but not coordinated to the iron of verdoheme. This structural device would ensure the smooth progression of the subsequent reaction, from verdoheme to biliverdin, which requires O(2) binding to verdoheme.

The structural dynamics of myoglobin

Conformational fluctuations in proteins were initially invoked to explain the observation that diffusion of small ligands through the matrix is a global phenomenon. Small globular proteins contain internal cavities that play a role not only in matrix dynamics but also in controlling function, tracing a pathway for the diffusion of the ligand to and from the active site. This is the main point addressed in this Review, which presents pertinent information obtained on myoglobin (Mb). Mb, a simple globular heme protein which binds reversibly oxygen and other ligands. The bond between the heme Fe(II) and gaseous ligands can be photodissociated by a laser pulse, generating a non-equilibrium population of protein structures that relaxes on a picosecond to millisecond time range. This process is associated with migration of the ligand to internal cavities of the protein, which are known to bind xenon. Some of the results obtained by laser photolysis, molecular dynamics simulations, and X-ray diffraction of intermediate states of wild-type and mutant myoglobins are summarized. The extended relaxation of the globin moiety directly observed by Laue crystallography reflects re-equilibration among conformational substates known to play an essential role in controlling protein function.

A High-Throughput Screen for Porphyrin Metal Chelatases: Application to the Directed Evolution of Ferrochelatases for Metalloporphyrin Biosynthesis

Porphyrins are of particular interest in a variety of applications ranging from biocatalysis and chemical synthesis to biosensor and electronic technologies as well as cancer treatment. Recently, we have developed a versatile system for the high-level production of porphyrins in engineered E. coli cells with the aim of diversifying substitution patterns and accessing porphyrin systems not readily available through chemical synthesis. However, this approach failed to produce significant amounts of the metalloporphyrin in vivo from overproduced protoporphyrin due to insufficient metal insertion. Therefore, we systematically assessed the activity of the B. subtilis ferrochelatase in vivo and in vitro. A true high-throughput-screening approach based on catalytic in vivo ferrochelatase activity was developed by using fluorescence-activated cell sorting (FACS). This assay was used to screen a library of 2.4 x 10(6) ferrochelatase mutants expressed in protoporphyrin-overproducing recombinant E. coli cells. Several selected protein variants were purified, and their improved catalytic activity was confirmed in vitro. In addition to ferrochelatase activity, metal transport into E. coli was identified as another limitation for in vivo heme overproduction. Overexpression of the metal transporter zupT as part of the assembled pathway increased the overall metalloporphyrin production twofold. This report represents the most exhaustive in vitro evolution study of a ferrochelatase and demonstrates the effectiveness of our novel high-throughput-screening system for directed evolution of ferrochelatases based on their catalytic activity.