Homodimeric Hemoglobin Research Articles

Cooperative Protein Dynamics of Heterotetrameric Hemoglobin from Scapharca inaequivalvis.

Hemoglobins achieve cooperative oxygen binding by diverse strategies based on different assemblies of globin subunits. Heterotetrameric hemoglobin from Scapharca inaequivalvis (HbII) consists of two AB-dimers, whose structure closely resembles that of homodimeric hemoglobin from the same organism (HbI). Herein, we investigated the structural dynamics of HbII following carbon monoxide (CO) dissociation using time-resolved resonance Raman (RR) spectroscopy. The observed spectra showed that the heme structure of the transient dissociated form of HbII was similar to that of HbI; however, the transition from the transient dissociated form to the equilibrium unligated form was faster for HbII than for HbI. Furthermore, the dependence of the time-resolved spectra on the yield of CO dissociation revealed that the transition became faster as the number of dissociated ligands increased from one to four. The positive correlation between the rate constants and number of dissociated ligands indicates that the structural transition of HbII following CO dissociation is cooperative.

Ultrafast coherent motion and helix rearrangement of homodimeric hemoglobin visualized with femtosecond X-ray solution scattering

Ultrafast motion of molecules, particularly the coherent motion, has been intensively investigated as a key factor guiding the reaction pathways. Recently, X-ray free-electron lasers (XFELs) have been utilized to elucidate the ultrafast motion of molecules. However, the studies on proteins using XFELs have been typically limited to the crystalline phase, and proteins in solution have rarely been investigated. Here we applied femtosecond time-resolved X-ray solution scattering (fs-TRXSS) and a structure refinement method to visualize the ultrafast motion of a protein. We succeeded in revealing detailed ultrafast structural changes of homodimeric hemoglobin involving the coherent motion. In addition to the motion of the protein itself, the time-dependent change of electron density of the hydration shell was tracked. Besides, the analysis on the fs-TRXSS data of myoglobin allows for observing the effect of the oligomeric state on the ultrafast coherent motion.

Change in vibrational entropy with change in protein volume estimated with mode Grüneisen parameters.

For a small adjustment in average volume, due to a change in state of a protein or other macromolecule at constant temperature, the change in vibrational entropy is related to the mode Grüneisen parameters, which relate shifts in frequency to a small volume change. We report here values of mode Grüneisen parameters computed for two hydrated proteins, cytochrome c and myoglobin, which exhibit trends with mode frequency resembling those of glassy systems. We use the mode Grüneisen parameters to relate volumetric thermal expansion to previously computed values of the isothermal compressibility for several proteins. We also estimate changes in vibrational entropy resulting from the change in volume upon ligand bonding of myoglobin and the homodimeric hemoglobin from Scapharca inaequivalvis (HbI). We compare estimates of the change in entropy upon ligation obtained in terms of mode Grüneisen parameters with the results of normal mode analysis for myoglobin and earlier molecular dynamics simulations of HbI. The results illustrate how small changes in average volume can yield changes in entropy that contribute to ligand binding and allostery.

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.

Locating and Navigating Energy Transport Networks in Proteins.

We review computational methods to locate energy transport networks in proteins that are based on the calculation of local energy diffusion in nanoscale systems. As an illustrative example, we discuss energy transport networks computed for the homodimeric hemoglobin from Scapharca inaequivalvis, where channels for facile energy transport, which include the cluster of water molecules at the interface of the globules, have been found to lie along pathways that experiments reveal are important in allosteric processes. We also review recent work on master equation simulations to model energy transport dynamics, including efforts to relate rate constants in the master equation to protein structural dynamics. Results for apomyoglobin involving relations between fluctuations in the length of hydrogen bonds and the energy flux between them are presented.

Water-mediated biomolecular dynamics and allostery.

Dynamic coupling with water contributes to regulating the functional dynamics of a biomolecule. We discuss protein-water dynamics, with emphasis on water that is partially confined, and the role of protein-confined water dynamics in allosteric regulation. These properties are illustrated with two systems, a homodimeric hemoglobin from Scapharca inaequivalvis (HbI) and an A2A adenosine receptor (A2AAR). For HbI, water-protein interactions, long known to contribute to the thermodynamics of cooperativity, are seen to influence the dynamics of the protein not only around the protein-water interface but also into the core of each globule, where dynamic and entropic changes upon ligand binding are coupled to protein-water contact dynamics. Similarly, hydration waters trapped deep inside the core region of A2AAR enable the formation of an allosteric network made of water-mediated inter-residue contacts. Extending from the ligand binding pocket to the G-protein binding site, this allosteric network plays key roles in regulating the activity of the receptor.

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 (I1, I2, and I3) as in the wild type (WT), but its geminate CO recombination occurs directly from I1 without the path via I2 observed in WT. Moreover, the structural transitions, which involve ligand migration (the transitions from I1 to I2 and from I3 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.

Variation of Energy Transfer Rates across Protein-Water Contacts with Equilibrium Structural Fluctuations of a Homodimeric Hemoglobin.

Molecular dynamics simulations of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) have been carried out to examine relations between rates of vibrational energy transfer across nonbonded contacts and equilibrium structural fluctuations, with emphasis on protein-water contacts. The scaling of rates of energy transfer with equilibrium fluctuations of the contact length is found to hold up well for contacts between residues and hemes at the interface and the cluster of 17 interface water molecules in the unliganded state of HbI, as well as for the liganded state, for which the cluster contains on average 11 water molecules. In both states, the rate of energy transfer is also found to satisfy a diffusion relation. Within each globule, the scaling for polar contacts is similar to that found in an earlier analysis of myoglobin. Entropy associated with dynamics of polar contacts within each globule and with contacts between the hemes and water cluster is found to increase upon ligation. Energy exchange networks (EENs) for liganded and unliganded states obtained from the simulations are also presented and discussed. Energy transport networks through which nonbonded contacts transport energy in HbI, referred to as nonbonded networks (NBNs), are determined from the EENs and compared for the two states.

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.

Water-Mediated Energy Dynamics in a Homodimeric Hemoglobin.

We examine energy dynamics in the unliganded and liganded states of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI), which exhibits cooperativity mediated by the cluster of water molecules at the interface upon ligand binding and dissociation. We construct and analyze a dynamic network in which nodes representing the residues, hemes, and water cluster are connected by edges that represent energy transport times, as well as a nonbonded network (NBN) indicating regions that respond rapidly to local strain within the protein via nonbonded interactions. One of the two largest NBNs includes the Lys30-Asp89 salt bridge critical for stabilizing the dimer. The other includes the hemes and surrounding residues, as well as, in the unliganded state, the cluster of water molecules between the globules. Energy transport in the protein appears to be controlled by the Lys30-Asp89 salt bridge critical for stabilizing the dimer, as well as the interface water cluster in the unliganded state. Possible connections between energy transport dynamics in response to local strain identified here and allosteric transitions in HbI are discussed.

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.

Ligand and interfacial dynamics in a homodimeric hemoglobin.

The structural dynamics of dimeric hemoglobin (HbI) from Scapharca inaequivalvis in different ligand-binding states is studied from atomistic simulations on the μs time scale. The intermediates are between the fully ligand-bound (R) and ligand-free (T) states. Tertiary structural changes, such as rotation of the side chain of Phe97, breaking of the Lys96–heme salt bridge, and the Fe–Fe separation, are characterized and the water dynamics along the R-T transition is analyzed. All these properties for the intermediates are bracketed by those determined experimentally for the fully ligand-bound and ligand-free proteins, respectively. The dynamics of the two monomers is asymmetric on the 100 ns timescale. Several spontaneous rotations of the Phe97 side chain are observed which suggest a typical time scale of 50–100 ns for this process. Ligand migration pathways include regions between the B/G and C/G helices and, if observed, take place in the 100 ns time scale.

Protein structural dynamics revealed by time-resolved X-ray solution scattering.

One of the most important questions in biological science is how a protein functions. When a protein performs its function, it undergoes regulated structural transitions. In this regard, to better understand the underlying principle of a protein function, it is desirable to monitor the dynamic evolution of the protein structure in real time. To probe fast and subtle motions of a protein in physiological conditions demands an experimental tool that is not only equipped with superb spatiotemporal resolution but also applicable to samples in solution phase. Time-resolved X-ray solution scattering (TRXSS), discussed in this Account, fits all of those requirements needed for probing the movements of proteins in aqueous solution. The technique utilizes a pump-probe scheme employing an optical pump pulse to initiate photoreactions of proteins and an X-ray probe pulse to monitor ensuing structural changes. The technical advances in ultrafast lasers and X-ray sources allow us to achieve superb temporal resolution down to femtoseconds. Because X-rays scatter off all atomic pairs in a protein, an X-ray scattering pattern provides information on the global structure of the protein with subangstrom spatial resolution. Importantly, TRXSS is readily applicable to aqueous solution samples of proteins with the aid of theoretical models and therefore is well suited for investigating structural dynamics of protein transitions in physiological conditions. In this Account, we demonstrate that TRXSS can be used to probe real-time structural dynamics of proteins in solution ranging from subtle helix movement to global conformational change. Specifically, we discuss the photoreactions of photoactive yellow protein (PYP) and homodimeric hemoglobin (HbI). For PYP, we revealed the kinetics of structural transitions among four transient intermediates comprising a photocycle and, by applying structural analysis based on ab initio shape reconstruction, showed that the signaling of PYP involves the protrusion of the N-terminus with significant increase of the overall protein size. For HbI, we elucidated the dynamics of complex allosteric transitions among transient intermediates. In particular, by applying structural refinement analysis based on rigid-body modeling, we found that the allosteric transition of HbI accompanies the rotation of quaternary structure and the contraction between two heme domains. By making use of the experimental and analysis methods presented in this Account, we envision that the TRXSS can be used to probe the structural dynamics of various proteins, allowing us to decipher the working mechanisms of their functions. Furthermore, when combined with femtosecond X-ray pulses generated from X-ray free electron lasers, TRXSS will gain access to ultrafast protein dynamics on sub-picosecond time scales.

Communication Maps: Exploring Energy Transport through Proteins and Water

AbstractFrequency‐resolved communication maps provide a coarse‐grained, global mapping of energy transport channels in a protein as a function of frequency of modes that carry energy. We illustrate the approach with a study of the homodimeric hemoglobin of Scapharca inaequivalvis, which exhibits cooperativity during ligand binding. We compare energy transport between the two hemes of the unliganded and oxygenated protein, which is mediated by water as well as residues forming a hydrogen‐bonding network at the interface between the globules, and lies along the pathway for allosteric transitions observed in time‐resolved X‐ray studies. Non‐equilibrium molecular simulations on energy transport from the heme corroborate the energy transport pathways identified by the communication maps.

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 (I1, I2, and I3) 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 I1 intermediate is generated within 100 ps and transforms to the R-like I2 intermediate with a time constant of 3.2 ± 0.2 ns. Subsequently, the late, T-like I3 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 I3 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.

Cooperative macromolecular device revealed by meta-analysis of static and time-resolved structures

Here we present a meta-analysis of a large collection of static structures of a protein in the Protein Data Bank in order to extract the progression of structural events during protein function. We apply this strategy to the homodimeric hemoglobin HbI from Scapharca inaequivalvis. We derive a simple dynamic model describing how binding of the first ligand in one of the two chemically identical subunits facilitates a second binding event in the other partner subunit. The results of our ultrafast time-resolved crystallographic studies support this model. We demonstrate that HbI functions like a homodimeric mechanical device, such as pliers or scissors. Ligand-induced motion originating in one subunit is transmitted to the other via conserved pivot points, where the E and F' helices from two partner subunits are "bolted" together to form a stable dimer interface permitting slight relative rotation but preventing sliding.

Communication maps computed for homodimeric hemoglobin: Computational study of water-mediated energy transport in proteins

Frequency-resolved communication maps provide a coarse-grained picture of energy transport in nanoscale systems. We calculate communication maps for homodimeric hemoglobin from Scapharca inaequivalvis and sample them to elucidate energy transfer pathways between the binding sites and other parts of the protein with focus on the role of the cluster of water molecules at the interface between the globules. We complement analysis of communication maps with molecular simulations of energy flow. Both approaches reveal that excess energy in one heme flows mainly to regions of the interface where early hydrogen bond rearrangements occur in the allosteric transition. In particular, energy is carried disproportionately by the water molecules, consistent with the larger thermal conductivity of water compared to proteins.

Time Resolved Crystallographic Analysis of Cooperative Ligand Binding and Ligand Migration in a Dimeric Hemoglobin

Invertebrate hemoglobins, which range in size from dimers to assemblies of hundreds of subunits, offer excellent models for the investigation of allosteric protein function. A recurring theme among cooperative invertebrate hemoglobins is a subunit pairing involving the heme-embedding E and F helices, despite a lack of sequence conservation in these helices. This is quite different from the assembly of vertebrate hemoglobins and suggests that such a pairing is amenable for regulating oxygen delivery. We have carried out extensive time-resolved crystallographic analysis on the simplest of these hemoglobins, the cooperative homodimeric hemoglobin (HbI) from the blood clam Scapharca inaequivalvis. Limited ligand-linked subunit rotation permits the full allosteric transition to be followed in crystals of HbI. These studies have been combined with mutant and functional studies along with conventional crystallographic analysis to provide a comprehensive understanding of the communication between subunits and the movement of ligands through the protein. Investigation of ligand migration through this protein suggests, surprisingly, that ligands primarily escape through a distal histidine gate, despite its involvement in the subunit interface. As a result, ligand escape may require transient subunit movement, which is inhibited in the protein crystal lattice. These findings emphasize the central role of dynamics in protein function. Our studies are also directed towards understanding cooperative protein function in larger invertebrate hemoglobin assemblies. Preliminary time-resolved crystallogaphic analysis on the tetrameric Scapharca HbII suggests similar structural transitions as HbI and offers the possibility of investigating how changes in one subunit directly impact a neighboring subunit in the heterodimeric halves of HbII. In addition, we will discuss new insights that may provide unifying themes for the basis of cooperativity among diverse invertebrate hemoglobins that form EF dimeric pairing.

Dynamics of Water Clusters Confined in Proteins: A Molecular Dynamics Simulation Study of Interfacial Waters in a Dimeric Hemoglobin

Water confined in proteins exhibits dynamics distinct from the dynamics of water in the bulk or near the surface of a biomolecule. We examine the water dynamics at the interface of the two globules of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) by molecular dynamics (MD) simulations, with focus on water-protein hydrogen bond lifetimes and rotational anisotropy of the interfacial waters. We find that relaxation of the waters at the interface of both deoxy- and oxy-HbI, which contain a cluster of 17 and 11 interfacial waters, respectively, is well described by stretched exponentials with exponents from 0.1 to 0.6 and relaxation times of tens to thousands of picoseconds. The interfacial water molecules of oxy-HbI exhibit slower rotational relaxation and hydrogen bond rearrangement than those of deoxy-HbI, consistent with an allosteric transition from unliganded to liganded conformers involving the expulsion of several water molecules from the interface. Though the interfacial waters are translationally and rotationally static on the picosecond time scale, they contribute to fast communication between the globules via vibrations. We find that the interfacial waters enhance vibrational energy transport across the interface by ≈10%.

Protein Conformational Dynamics of Homodimeric Hemoglobin Revealed by Combined Time‐Resolved Spectroscopic Probes

Protein Conformational Dynamics of Homodimeric Hemoglobin Revealed by Combined Time‐Resolved Spectroscopic Probes