Picosecond Dynamics of Ligand Interconversion in the Primary Docking Site of Heme Proteins

Enhancing Entangled Two-Photon Absorption for Picosecond Quantum Spectroscopy

Entangled two-photon absorption (ETPA) is known to create photoinduced transitions with extremely low light intensity, reducing the risk of phototoxicity compared to classical two-photon absorption. Previous works have predicted the ETPA cross-section, σ_e, to vary inversely with the product of entanglement time (T_e) and entanglement area (A_e), i.e., σ_e ∼ 1/A_eT_e. The decreasing σ_e with increasing T_e has limited ETPA to fs-scale T_e, while ETPA applications for ps-scale spectroscopy have been unexplored. However, we show that spectral-spatial coupling, which reduces A_e as the SPDC bandwidth (σ_f) decreases, plays a significant role in determining σ_e when T_e > ∼100 fs. We experimentally measured σ_e for zinc tetraphenylporphyrin at several σ_f values. For type-I ETPA, σ_e increases as σ_f decreases down to 0.1 ps^-1. For type-II SPDC, σ_e is constant for a wide range of σ_f. With a theoretical analysis of the data, the maximum type-I σ_e would occur at σ_f = 0.1 ps^-1 (T_e = 10 ps). At this maximum, σ_e is 1 order of magnitude larger than fs-scale σ_e and 3 orders of magnitude larger than previous predictions of ps-scale σ_e. By utilizing this spectral-spatial coupling, narrowband type-I ETPA provides a new opportunity to increase the efficiency of measuring nonlinear optical signals and to control photochemical reactions requiring ps temporal precision.

The time scale of the quaternary structural changes in hemoglobin revealed using the transient grating technique

The quaternary structural transition between the R and T states of human hemoglobin was investigated using the transient grating technique. The results presented herein reveal that the quaternary structural change accompanied by the R-T transition occurs within a few microseconds.

Dynamics of geminate rebinding of CO to cytochrome c in guanidine HCl probed by femtosecond vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding dynamics of CO to cytochrome c (Cytc) in 1.8 and 7 M guanidine HCl (GdnHCl) after photodeligation of the corresponding CO-bound protein in D2O buffer (pD = 7.4) at 283 K. Geminate rebinding (GR) dynamics of CO to the folded Cytc in 1.8 M GdnHCl (nCytc) is similar to that to chemically modified cytochrome c (cCytc), suggesting that the overall conformations of nCytcCO and cCytcCO are similar. About 86% of the dissociated CO molecules were geminately rebound to nCytc nonexponentially within 1 ns. The efficient GR of CO to the folded Cytc can be attributed to the organized protein matrix near the active site of nCytc that provides an efficient trap for the diffusing CO ligand after photodissociation. Although the concentration of nCytc did not affect its GR yield of CO, GR yield of CO to the unfolded Cytc in 7 M GdnHCl (uCytc) increased from 5 to 30% as the protein concentration increased from 0.3 to 9 mM. Time-resolved spectra of the (13)CO dissociated from both 9 mM nCytc(13)CO and 9 mM uCytc(13)CO showed a growing band with a peak at 2090 cm(-1) on the picosecond time scale, which was assigned to (13)CO in D2O solvent. At 1 ns, the fraction of the CO band in the solvent was about 10% of the nascent photodeligated protein in nCytc and more than 50% in the concentrated uCytc. Whereas a small opening in the active site of nCytc is responsible for the ultrafast escape of CO to solution in the folded protein, a large fraction of the CO escape to the solvent in uCytc results from the denatured structure of the active site in the unfolded protein. The spectrum of the CO dissociated from the concentrated uCytcCO contained a band that decayed as efficiently as that for the folded protein, suggesting that some fraction of uCytcCO might form aggregates even in 7 M denaturant, such that the aggregate acts as an efficient trap for the diffusing CO after deligation. No hint of precipitate in the concentrated uCytcCO and protein refolding upon dilution of the GdnHCl indicate that the aggregate does not grow continuously but remains as a soluble oligomer. The delayed appearance of the solvated CO and the inefficient GR of CO in uCytcCO suggest that the monomeric unfolded CytcCO so loosely arranged that the protein matrix cannot trap CO efficiently but the bound CO is still buried within hydrophobic residues even under the harsh denaturation condition.

Geminate rebinding dynamics of nitric oxide to ferric hemoglobin in D2O solution

Femtosecond mid-infrared (mid-IR) spectroscopy was used to probe geminate rebinding (GR) dynamics of photo-released nitric oxide (NO) to ferric hemoglobin (Hb(III)) in D2O solution at room temperature. Time-resolved vibrational spectra exhibit two overlapping NO bands for NO-bound Hb(III) (Hb(III)NO), a major band at 1925 cm(-1) (89%) and a minor one at 1905 cm(-1) (11%), suggesting that Hb(III)NO has at least two conformational substates. Both bands decay nonexponentially, each with a different time scale, and the decays are described by a stretched exponential function; the major band's decay is described by 0.96 exp(-t/40 ps)(0.86) + 0.04 and the minor band's decay is described by exp(-t/85 ps)(0.75). These decays arise mainly from the GR of the photo-released NO to Hb(III), indicating that the bound state's conformer influences the NO binding. In particular, the His64 residue, known to have inward conformation in the major band and outward conformation in the minor band, plays a significant role in controlling the binding of NO to Hb(III). The GR of NO to ferric Hb is slower than that to ferrous Hb, which shows fast and efficient GR due to the high reactivity of NO to the heme Fe(ii). The slower GR of NO to Hb(III) may be caused by the lower reactivity of NO to the heme Fe(iii).

Carbon monoxide recombination dynamics in truncated hemoglobins studied with visible-pump midIR-probe spectroscopy

Carbon monoxide recombination dynamics upon photodissociation with visible light has been characterized by means of ultrafast visible-pump/MidIR probe spectroscopy for the truncated hemoglobins from Thermobifida fusca and Bacillus subtilis. Photodissociation has been induced by exciting the sample at two different wavelengths: 400 nm, corresponding to the heme absorption in the B-band, and 550 nm, in the Q-bands. The bleached iron-CO coordination band located at 1850-1950 cm(-1) and the free CO absorption band in the region 2050-2200 cm(-1) have been observed by probe pulses tuned in the appropriate infrared region. The kinetic traces measured at 1850-1950 cm(-1) reveal multiexponential subnanosecond dynamics that have been interpreted as arising from fast geminate recombination of the photolyzed CO. A compared analysis of the crystal structure of the two proteins reveals a similar structure of their distal heme pocket, which contains conserved polar and aromatic amino acid residues closely interacting with the iron ligand. Although fast geminate recombination is observed in both proteins, several kinetic differences can be evidenced, which can be interpreted in terms of a different structural flexibility of the corresponding heme distal pockets. The analysis of the free CO band-shape and of its dynamic evolution brings out novel features about the nature of the docking site inside the protein cavity.

From structure to cellular mechanism with infrared microspectroscopy

Current efforts in structural biology aim to integrate structural information within the context of cellular organization and function. X-rays and infrared radiation stand at opposite ends of the electromagnetic spectrum and act as complementary probes for achieving this goal. Intense and bright beams are produced by synchrotron radiation, and are efficiently used in the wavelength domain extending from hard X-rays to the far-infrared (or THz) regime. While X-ray crystallography provides exquisite details on atomic structure, Fourier transform infrared microspectroscopy (FTIRM) is emerging as a spectroscopic probe and imaging tool for correlating molecular structure to biochemical dynamics and function. In this manuscript, the role of synchrotron FTIRM in bridging the gap towards 'functional biology' is discussed based upon recent achievements, with a critical assessment of the contributions to biological and biomedical research.

The role of higher CO-multipole moments in understanding the dynamics of photodissociated carbonmonoxide in myoglobin

The influence of electrostatic multipole moments up to hexadecapole on the dynamics of photodissociated carbon monoxide (CO) in myoglobin is investigated. The CO electrostatic potential is expressed as an expansion into atomic multipole moments of increasing order up to octopole which are obtained from a distributed multipole analysis. Three models with increasingly accurate molecular multipoles (accurate quadrupole, octopole, and hexadecapole moments, respectively) are developed and used in molecular dynamics simulations. All models with a fluctuating quadrupole moment correctly describe the location of the B-state whereas the sign of the octopole moment differentiates between the Fe...CO and Fe...OC orientation. For the infrared spectrum of photodissociated CO, considerable differences between the three electrostatic models are found. The most detailed electrostatic model correctly reproduces the splitting, shift, and width of the CO spectrum in the B-state. From an analysis of the trajectories, the spectroscopic B(1) and B(2) states are assigned to the Fe...CO and Fe...OC substates, respectively.

Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.

Force-Field Development and Molecular Dynamics Simulations of Ferrocene–Peptide Conjugates as a Scaffold for Hydrogenase Mimics

The increasing importance of hydrogenase enzymes in the new energy research field has led us to examine the structure and dynamics of potential hydrogenase mimics, based on a ferrocene-peptide scaffold, using molecular dynamics (MD) simulations. To enable this MD study, a molecular mechanics force field for ferrocene-bearing peptides was developed and implemented in the CHARMM simulation package, thus extending the usefulness of the package into peptide-bioorganometallic chemistry. Using the automated frequency-matching method (AFMM), optimized intramolecular force-field parameters were generated through quantum chemical reference normal modes. The partial charges for ferrocene were derived by fitting point charges to quantum-chemically computed electrostatic potentials. The force field was tested against experimental X-ray crystal structures of dipeptide derivatives of ferrocene-1,1'-dicarboxylic acid. The calculations reproduce accurately the molecular geometries, including the characteristic C2-symmetrical intramolecular hydrogen-bonding pattern, that were stable over 0.1 micros MD simulations. The crystal packing properties of ferrocene-1-(D)alanine-(D)proline-1'-(D)alanine-(D)proline were also accurately reproduced. The lattice parameters of this crystal were conserved during a 0.1 micros MD simulation and match the experimental values almost exactly. Simulations of the peptides in dichloromethane are also in good agreement with experimental NMR and circular dichroism (CD) data in solution. The developed force field was used to perform MD simulations on novel, as yet unsynthesized peptide fragments that surround the active site of [Ni-Fe] hydrogenase. The results of this simulation lead us to propose an improved design for synthetic peptide-based hydrogenase models. The presented MD simulation results of metallocenes thereby provide a convincing validation of our proposal to use ferrocene-peptides as minimal enzyme mimics.

Direct observation of ligand transfer and bond formation in cytochrome c oxidase by using mid-infrared chirped-pulse upconversion

We have implemented the recently demonstrated technique of chirped-pulse upconversion of midinfrared femtosecond pulses into the visible in a visible pump-midinfrared probe experiment for high-resolution, high-sensitivity measurements over a broad spectral range. We have succeeded in time-resolving the CO ligand transfer process from the heme Fe to the neighboring Cu(B) atom in the bimetallic active site of mammalian cytochrome c oxidase, which was known to proceed in <1 ps, using the full CO vibrational signature of Fe-CO bond breaking and Cu(B)-CO bond formation. Our differential transmission results show a delayed onset of the appearance of the Cu(B)-bound species (200 fs), followed by a 450-fs exponential rise. Trajectories calculated by using molecular-dynamics simulations with a Morse potential for the Cu(B)-C interaction display a similar behavior. Both experimental and calculated data strongly suggest a ballistic contribution to the transfer process.

Angular momentum dependent friction slows down rotational relaxation under nonequilibrium conditions

It has recently been shown that relaxation of the rotational energy of hot nonequilibrium photofragments (i) slows down significantly with the increase of their initial rotational temperature and (ii) differs dramatically from the relaxation of the equilibrium rotational energy correlation function, manifesting thereby the breakdown of the linear response description [A. C. Moskun et al., Science 311, 1907 (2006)]. We demonstrate that this phenomenon may be caused by the angular momentum dependence of rotational friction. We have developed the generalized Fokker-Planck equation whose rotational friction depends upon angular momentum algebraically. The calculated rotational correlation functions correspond well to their counterparts obtained via molecular dynamics simulations in a broad range of initial nonequilibrium conditions. It is suggested that the angular momentum dependence of friction should be taken into account while describing rotational relaxation far from equilibrium.

Theoretical characterization of carbon monoxide vibrational spectrum in sperm whale myoglobin distal pocket

In this article we use the perturbed matrix method and an extended molecular dynamics sampling of the carbon monoxide (CO) in the myoglobin distal pocket to characterize the CO vibrational spectrum and hence to relate its spectroscopic features with the atomic-molecular behavior. Results show the accuracy of the method employed and confirm the assignment of the spectroscopic B1 and B2 states proposed by Lim et al.

Protein ligand migration mapped by nonequilibrium 2D-IR exchange spectroscopy

2D-IR exchange spectroscopy has been introduced recently to map chemical exchange networks in equilibrium with subpicosecond time resolution. Here, we demonstrate the generalization of 2D-IR exchange spectroscopy to nonequilibrium systems and its application to map light-triggered migration of ligands between different sites in a protein. Within picoseconds after a photodissociating laser pulse, carbon monoxide ligands relocate from their binding site A at the heme prosthetic group of myoglobin to a primary docking site B in the distal heme pocket. Multiple CO stretching bands are observed for the CO ligand in sites A and B, indicating that several distinct conformational substates of the myoglobin:ligand complex coexist in solution. Exchange cross-peaks between the bands associated with substates of heme-bound CO and photodissociated CO in the primary docking site reveal the substate connectivity at physiological temperature.

Ligand Dynamics in an Electron Transfer Protein

Substitution of the heme coordination residue Met-80 of the electron transport protein yeast iso-1-cytochrome c allows external ligands like CO to bind and thus increase the effective redox potential. This mutation, in principle, turns the protein into a quasi-native photoactivable electron donor. We have studied the kinetic and spectral characteristics of geminate recombination of heme and CO in a series of single M80X (X = Ala, Ser, Asp, Arg) mutants, using femtosecond transient absorption spectroscopy. In these proteins, all geminate recombination occurs on the picosecond and early nanosecond time scale, in a multiphasic manner, in which heme relaxation takes place on the same time scale. The extent of geminate recombination varies from >99% (Ala, Ser) to approximately 70% (Arg), the latter value being in principle low enough for electron injection studies. The rates and extent of the CO geminate recombination phases are much higher than in functional ligand-binding proteins like myoglobin, presumably reflecting the rigid and hydrophobic properties of the heme environment, which are optimized for electron transfer. Thus, the dynamics of CO recombination in cytochrome c are a tool for studying the heme pocket, in a similar way as NO in myoglobin. We discuss the differences in the CO kinetics between the mutants in terms of the properties of the heme environment and strategies to enhance the CO escape yield. Experiments on double mutants in which Phe-82 is replaced by Asp or Gly as well as the M80D substitution indicate that such steric changes substantially increase the motional freedom-dissociated CO.