Ultrafast dynamics of myoglobin probed by time-resolved resonance Raman spectroscopy

Influence of pump laser fluence on ultrafast myoglobin structural dynamics

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump-probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer1,2. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore3-17. As multiphoton absorption may force the protein response into non-physiological pathways, it is of great concern18,19 whether this experimental approach20 allows valid conclusions to be drawn vis-à-vis biologically relevant single-photon-induced reactions18,19. Here we describe ultrafast pump-probe SFX experiments on the photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe-CO bond distance (predicted by recent quantum wavepacket dynamics21) are seen to depend strongly on pump laser energy, in line with quantum chemical analysis. Our results confirm both the feasibility and necessity of performing ultrafast TR-SFX pump-probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing both the design and the interpretation of ultrafast TR-SFX pump-probe experiments20 such that mechanistically relevant insight emerges.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Dependence of Vibrational Energy Transfer on Distance in a Four-Helix Bundle Protein: Equidistant Increments with the Periodicity of α Helices

Vibrational energy exchanges between various degrees of freedom are critical to barrier-crossing processes in proteins. Heme proteins are highly suitable for studies of the vibrational energy exchanges in proteins. The migration of excess energy released by heme in a protein moiety can be observed using time-resolved anti-Stokes ultraviolet resonance Raman spectroscopy. The anti-Stokes resonance Raman intensity of a tryptophan residue is an excellent probe for the excess energy and the spatial resolution of a single amino acid residue can be achieved. Here, we studied dependence of vibrational energy transfer on the distance in cytochrome b₅₆₂, which is a heme-containing, four-helix bundle protein. The vibrational energy transfer from the heme group to a single tryptophan residue introduced by site-directed mutagenesis was examined for different heme-tryptophan distances by a quasi-constant length with the periodicity of α helices. Taken together with structural data obtained by molecular dynamics simulations, the energy transfer could be well described by the model of classical thermal diffusion, which suggests that continuum media provide a good approximation of the protein interior, of which the atomic packing density is very high.

A computational methodology to diagnose sequence-variant dynamic perturbations by comparing atomic protein structures

MOTIVATION

The objective is to diagnose dynamics perturbations caused by amino-acid mutations as prerequisite to assess protein functional health or drug failure, simply using network models of protein X-ray structures.

RESULTS

We find that the differences in the allocation of the atomic interactions of each amino acid to 1D, 2D, 3D, 4D structural levels between variants structurally robust, recover experimental dynamic perturbations. The allocation measure validated on two B-pentamers variants of AB5 toxins having 17 mutations, also distinguishes dynamic perturbations of pathogenic and non-pathogenic Transthyretin single-mutants. Finally, the main proteases of the coronaviruses SARS-CoV and SARS-CoV-2 exhibit changes in the allocation measure, raising the possibility of drug failure despite the main proteases structural similarity.

AVAILABILITY AND IMPLEMENTATION

The Python code used for the production of the results is available at github.com/lorpac/protein_partitioning_atomic_contacts. The authors will run the analysis on any PDB structures of protein variants upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Primary Electronic and Vibrational Dynamics of Cytochrome c Observed by Sub-10 fs NUV Laser Pulses

The primary reaction mechanism of cytochrome c (Cyt c) was elucidated for two redox forms of ferric (oxidized) and ferrous (reduced) Cyt c by measuring their transient absorption (TA) spectra using a homemade sub-10 fs broadband NUV laser pulses system. The TA traces measured in the broad probe wavelength region were analyzed by the global analysis method to study the electronic dynamics. The difference of relaxation dynamics dependent on the excitation bandwidth enabled us to elucidate that the 2.5 ps component in ferrous Cyt c can be assigned to intramolecular vibration energy redistribution and not to vibrational cooling, which was not clear until this work. The temporal resolution of 10 fs observes TA signal modulation caused by the molecular vibration in the time domain, which can be used to calculate the instantaneous frequency of the molecular vibration mode. The observed vibrational dynamics has visualized that the heme structure changes in 0.8 ps for ferric Cyt c and in >1.0 ps for ferrous Cyt c. These estimated lifetimes of vibrational dynamics reflect vibrational relaxation in the ground state of ferric Cyt c and electronic transition from the S2 state to the S1 state in ferrous Cyt c, respectively.

Analysis of Nanoconfined Protein Dielectric Signals Using Charged Amino Acid Network Models

Protein slow motions involving collective molecular fluctuations on the timescale of microseconds to seconds are difficult to measure and not well understood despite being essential to sustain protein folding and protein function. Broadband dielectric spectroscopy (BDS) is one of the most powerful experimental techniques to monitor, over a broad frequency and temperature range, the molecular dynamics of soft matter through the orientational polarisation of permanent dipole moments that are generated by the chemical structure and morphological organisation of matter. Its typical frequency range goes from 107 Hz down to 10−3 Hz, being thus suitable for investigations on slow motions in proteins. Moreover, BDS has the advantage of providing direct experimental access to molecular fluctuations taking place on different length-scales, from local to cooperative dipolar motions. The unfolding of the cholera toxin B pentamer (CtxB5) after thermal treatment for 3h at 80°C is investigated by BDS under nanoconfined and dehydrated conditions. From the X-ray structure of the toxin pentamer, network-based models are used to infer the toxin dipoles present in the native state and to compute their stability and dielectric properties. Network analyses highlight three domains with distinct dielectric and stability properties that support a model where the toxin unfolds into three conformations after the treatment at 80°C. This novel integrative approach offers some perspective into the investigation of the relation between local perturbations (e.g. mutation, thermal treatment) and larger scale protein conformational changes. It might help ranking protein sequence variants according to their respective scale of dynamics perturbations.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Experimental Protein Molecular Dynamics: Broadband Dielectric Spectroscopy coupled with nanoconfinement

Protein dynamics covers multiple spatiotemporal scale processes, among which slow motions, not much understood even though they are underlying protein folding and protein functions. Protein slow motions are associated with structural heterogeneity, short-lived and poorly populated conformations, hard to detect individually. In addition, they involve collective motions of many atoms, not easily tracked by simulation and experimental devices. Here we propose a biophysical approach, coupling geometrical nanoconfinement and broadband dielectric spectroscopy (BDS), which distinguishes protein conformations by their respective molecular dynamics. In particular, protein-unfolding intermediates, usually poorly populated in macroscopic solutions are detected. The protein dynamics is observed under unusual conditions (sample nanoconfinement and dehydration) highlighting the robustness of protein structure and protein dynamics to a variety of conditions consistent with protein sustainability. The protein dielectric signals evolve with the temperature of thermal treatments indicating sensitivity to atomic and molecular interaction changes triggered by the protein thermal unfolding. As dipole fluctuations depend on both collective large-scale motions and local motions, the approach offers a prospect to track in-depth unfolding events.

Vibrational relaxation and energy migration of N-methylacetamide in water: the role of nonbonded interactions

Nonequilibrium molecular dynamics (MD) simulations together with physics-based force fields are used to follow energy flow between vibrationally excited N-methylacetamide (NMA) and water. The simulations are carried out with a previously validated force field for NMA, based on a multipolar representation of the electrostatics, and with a new fluctuating point charge model. For the water solvent, a flexible and a rigid model was employed to distinguish between the role of inter- and intramolecular degrees of freedom. On a 10 ps time scale about 90% of the available energy goes into the solvent. The remaining energy resides within internal NMA-degrees of freedom from where energy flow takes place on longer time scales. The total amount of energy transferred to the solvent on the 10 ps time scale does not depend on whether the water molecules are rigid or flexible during the simulations. Vibrational energy relaxation time scales include two regimes: one on the several 100 fs time scale and a longer one, ranging from 6 to 10 ps. This longer time scale agrees with previous simulations but overestimates the experimentally determined relaxation time by a factor of 2, which can be explained by the classical treatment of the vibrations. Including a previously determined quantum correction factor brings the long time scale into quite favorable agreement with experiment. Coupling to the bending vibration of the water molecules in H-bonding contact with the excited C═O chromophore is substantial. The equilibrium and nonequilibrium distribution of the bending angles of the water molecules in contact with the local oscillator are non-Gaussian, and one approaches the other on the subpicosecond time scale. Analysis of the water velocity distribution suggests that the C═O vibrational energy relaxes into the solvent water shells in an impulsive fashion on a picosecond time scale.

Observing Vibrational Energy Flow in a Protein with the Spatial Resolution of a Single Amino Acid Residue

One of the challenges in physical chemistry has been understanding how energy flows in a condensed phase from the microscopic viewpoint. To address this, space-resolved information at the molecular scale is required but has been lacking due to experimental difficulties. We succeeded in the real-time mapping of the vibrational energy flow in a protein with the spatial resolution of a single amino acid residue by combining time-resolved resonance Raman spectroscopy and site-directed single-Trp mutagenesis. Anti-Stokes Raman intensities of the Trp residues at different sites exhibited different temporal evolutions, reflecting propagation of the energy released by the heme group. A classical heat transport model was not able to reproduce the entire experimental data set, showing that we need a molecular-level description to explain the energy flow in a protein. The systematic application of our general methodology to proteins with different structural motifs may provide a greatly increased understanding of the energy flow in proteins.

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.

Resonance Raman spectroscopy

The main topics in resonance Raman spectroscopy presented at ICPP-2 in Kyoto are briefly discussed. These include: (i) coherent spectroscopy and low frequency vibrations of ligand-photodissociated heme proteins, (ii) vibrational relaxation revealed by time-resolved anti-Stokes Raman spectroscopy, (iii) electron transfer in porphyrin arrays, (iv) vibrational assignments of tetraazaporphyrins and (v) resonance Raman spectra of an NO storing protein, nitrophorin.

On the role of nonbonded interactions in vibrational energy relaxation of cyanide in water

The vibrationally excited cyanide ion (CN(-)) in H2O or D2O relaxes back to the ground state within several tens of picoseconds. Pump-probe infrared spectroscopy has determined relaxation times of T1 = 28 ± 7 and 71 ± 3 ps in H2O and D2O, respectively. Atomistic simulations of this process using nonequilibrium molecular dynamics simulations allow determination of whether it is possible at all to describe such a process, what level of accuracy in the force fields is required, and whether the information can be used to understand the molecular mechanisms underlying vibrational relaxation. It is found that, by using the best electrostatic models investigated, absolute relaxation times can be described rather more qualitatively (T1(H2O) = 19 ps and T1(D2O) = 34 ps) whereas the relative change in going from water to deuterated water is more quantitatively captured (factor of 2 vs 2.5 from experiment). However, moderate adjustment of the van der Waals ranges by less than 20% (for NVT) and 7.5% (for NVE), respectively, leads to almost quantitative agreement with experiment. Analysis of the energy redistribution establishes that the major pathway for CN(-) relaxation in H2O or D2O proceeds through coupling to the water-bending plus libration mode.

Direct evidence for mode-specific vibrational energy relaxation from quantum time-dependent perturbation theory. III. The nu(4) and nu(7) modes of nonplanar nickel porphyrin models

The time scales and pathways of vibrational energy relaxation (VER) of the nu(4) and nu(7) modes of three nickel porphyrin models, nickel porphine (NiP), nickel protoporphyrin IX (Ni-heme), and nickel octaethylporphyrin (NiOEP), were studied using a non-Markovian time-dependent perturbation theory at the B3LYP/6-31G(d) level. When NiP is calculated with D(4h) symmetry, it has the planar structure and the same VER properties as ferrous iron porphine (FeP). The porphine cores of both Ni-heme and NiOEP were distorted from a planar geometry, assuming a nonplanar structure, similar to that of the heme structure in cytochrome c. The VER time scales of Ni-heme are found to be similar to those predicted for a planar iron heme, but the derived pathways have distinctly different features. In particular, the strong coupling between the nu(7) mode and the overtone of the approximately 350 cm(-1) gamma(7) mode, observed for planar porphyrins, is absent in both nonplanar nickel porphyrins. Direct energy exchange between the nu(4) and nu(7) modes is not observed in NiOEP, but is found to play an essential role in the VER of the nu(4) mode in Ni-heme. The Ni-heme isopropionate groups are involved in the dominant VER pathways of both the nu(4) and nu(7) modes of Ni-heme. However, in contrast with VER pathways derived in planar iron heme, the isopropionate groups are not observed to play an essential role relative to other side chains in spatially directing the vibrational energy flow.

Diversity of solvent dependent energy transfer pathways in heme proteins

The time scales and pathways of heme cooling in both reduced cytochrome c and oxidized cytochrome c following heme photoexcitation were studied using molecular dynamics simulation. Five different solvent models, including normal water, heavy water, normal glycerol, deuterated glycerol, and a nonpolar solvent, were used in the simulation. Single exponential decay of the excess kinetic energy of the heme following photoexcitation was observed in all systems studied. The simulated time scale for heme cooling in normal water agrees with recent experimental results. In contrast to heme cooling in myoglobin, no solvent dependence was observed for the time scale for heme cooling in cytochrome c. The diversity of solvent dependence results from the different local heme environments in the two proteins. In myoglobin, it has been established that the dominant mechanism for heme cooling is direct energy transfer from the heme to the solvent. In cytochrome c, direct interaction between heme and protein residues forms the dominant energy transfer pathway. This distinction is dictated by protein topology and linked to function.

The impact of altered protein-heme interactions on the resonance Raman spectra of heme proteins. Studies of heme rotational disorder

Heme proteins that have been reconstituted with certain hemins may contain substantial fractions of a minor component in which the orientation of the heme in the folded pocket differs from the major ("native") conformation by a 180 degrees rotation about the alpha-gamma meso axis. In fact, this minor component has also been shown to exist in some native proteins, including several mammalian globins. While resonance Raman spectroscopy has emerged as a powerful probe of active site structure of heme proteins, no systematic study has yet been undertaken to elucidate the specific spectral changes associated with this disorder. In the present work, combined analyses of the temporal behavior of both NMR and RR data sets have been completed to permit the extraction of a unique RR spectrum for the disoriented form, documenting rather dramatic changes associated with this rotational disorder. In addition, the use of protohemes bearing selectively deuterated peripheral methyl groups has permitted the association of the observed modes with specific fragments of the heme residing in the reversed orientation. The studies conducted here clearly illustrate the exquisite sensitivity of low frequency heme deformation modes to altered protein-heme interactions.

Investigations of vibrational coherence in the low-frequency region of ferric heme proteins

Femtosecond coherence spectroscopy is applied to a series of ferric heme protein samples. The low-frequency vibrational spectra that are revealed show dominant oscillations near 40 cm(-1). MbCN is taken as a typical example of a histidine-ligated, six-coordinate, ferric heme and a comprehensive spectroscopic analysis is carried out. The results of this analysis reveal a new heme photoproduct species, absorbing near 418 nm, which is consistent with the photolysis of the His(93) axial ligand. The photoproduct undergoes subsequent rebinding/recovery with a time constant of approximately 4 ps. The photoproduct lineshapes are consistent with a photolysis quantum yield of 75-100%, although the observation of a relatively strong six-coordinate heme coherence near 252 cm(-1) (assigned to nu(9) in the MbCN Raman spectrum) suggests that the 75% lower limit is much more likely. The phase and amplitude excitation profiles of the low-frequency mode at 40 cm(-1) suggest that this mode is strongly coupled to the MbCN photoproduct species and it is assigned to the doming mode of the transient penta-coordinated material. The absolute phase of the 40 cm(-1) mode is found to be pi/2 on the red side of 418 nm and it jumps to 3pi/2 as excitation is tuned to the blue side of 418 nm. The absolute phase of the 40 cm(-1) signal is not explained by the standard theory for resonant impulsive stimulated Raman scattering. New mechanisms that give a dominant momentum impulse to the resonant wavepacket, rather than a coordinate displacement, are discussed. The possibilities of heme iron atom recoil after photolysis, as well as ultrafast nonradiative decay, are explored as potential ways to generate the strong momentum impulse needed to understand the phase properties of the 40 cm(-1) mode.

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.

Low frequency spectral density of ferrous heme: perturbations induced by axial ligation and protein insertion

Femtosecond coherence spectroscopy is used to probe low frequency (20-400 cm(-1)) modes of the ferrous heme group in solution, with and without 2-methyl imidazole (2MeIm) as an axial ligand. The results are compared to heme proteins (CPO, P450(cam), HRP, Mb) where insertion of the heme into the protein results in redistribution of the low frequency spectral density and in (approximately 60%) longer damping times for the coherent signals. The major effect of imidazole ligation to the ferrous heme is the "softening" of the low frequency force constants by a factor of approximately 0.6 +/- 0.1. The functional consequences of imidazole ligation are assessed and it is found that the enthalpic CO rebinding barrier is increased significantly when imidazole is bound. The force constant softening analysis, combined with the kinetics results, indicates that the iron is displaced by only approximately 0.2 A from the heme plane in the absence of the imidazole ligand, whereas it is displaced by approximately 0.4 A when imidazole (histidine) is present. This suggests that binding of imidazole (histidine) as an axial ligand, and the concomitant softening of the force constants, leads to an anharmonic distortion of the heme group that has significant functional consequences.

Crystallographic and Spectroscopic Studies of Peroxide-derived Myoglobin Compound II and Occurrence of Protonated FeIV–O

High resolution crystal structures of myoglobin in the pH range 5.2-8.7 have been used as models for the peroxide-derived compound II intermediates in heme peroxidases and oxygenases. The observed Fe-O bond length (1.86-1.90 A) is consistent with that of a single bond. The compound II state of myoglobin in crystals was controlled by single-crystal microspectrophotometry before and after synchrotron data collection. We observe some radiation-induced changes in both compound II (resulting in intermediate H) and in the resting ferric state of myoglobin. These radiation-induced states are quite unstable, and compound II and ferric myoglobin are immediately regenerated through a short heating above the glass transition temperature (<1 s) of the crystals. It is unclear how this influences our compound II structures compared with the unaffected compound II, but some crystallographic data suggest that the influence on the Fe-O bond distance is minimal. Based on our crystallographic and spectroscopic data we suggest that for myoglobin the compound II intermediate consists of an Fe(IV)-O species with a single bond. The presence of Fe(IV) is indicated by a small isomer shift of delta = 0.07 mm/s from Mössbauer spectroscopy. Earlier quantum refinements (crystallographic refinement where the molecular-mechanics potential is replaced by a quantum chemical calculation) and density functional theory calculations suggest that this intermediate H species is protonated.

Ultrafast heme dynamics in ferrous versus ferric cytochrome c studied by time-resolved resonance Raman and transient absorption spectroscopy

Cytochrome c (Cyt c) is a heme protein involved in electron transfer and also in apoptosis. Its heme iron is bisaxially ligated to histidine and methionine side chains and both ferric and ferrous redox states are physiologically relevant, as well as a ligand exchange between internal residue and external diatomic molecule. The photodissociation of internal axial ligand was observed for several ferrous heme proteins including Cyt c, but no time-resolved studies have been reported on ferric Cyt c. To investigate how the oxidation state of the heme influences the primary photoprocesses, we performed a comprehensive comparative study on horse heart Cyt c by subpicosecond time-resolved resonance Raman and femtosecond transient absorption spectroscopy. We found that in ferric Cyt c, in contrast to ferrous Cyt c, the photodissociation of an internal ligand does not take place, and relaxation dynamics is dominated by vibrational cooling in the ground electronic state of the heme. The intermolecular vibrational energy transfer was found to proceed in a single phase with a temperature decay of approximately 7 ps in both ferric and ferrous Cyt c. For ferrous Cyt c, the instantaneous photodissociation of the methionine side chain from the heme iron is the dominant event, and its rebinding proceeds in two phases, with time constants of approximately 5 and approximately 16 ps. A mechanism of this process is discussed, and the difference in photoinduced coordination behavior between ferric and ferrous Cyt c is explained by an involvement of the excited electronic state coupled with conformational relaxation of the heme.

Primary protein response after ligand photodissociation in carbonmonoxy myoglobin

Time-resolved UV resonance Raman (UVRR) spectroscopic studies of WT and mutant myoglobin were performed to reveal the dynamics of protein motion after ligand dissociation. After dissociation of carbon monoxide (CO) from the heme, UVRR bands of Tyr showed a decrease in intensity with a time constant of 2 ps. The intensity decrease was followed by intensity recovery with a time constant of 8 ps. On the other hand, UVRR bands of Trp residues located in the A helix showed an intensity decrease that was completed within the instrument response time. The intensity decrease was followed by an intensity recovery with a time constant of approximately 50 ps and lasted up to 1 ns. The time-resolved UVRR study of the myoglobin mutants demonstrated that the hydrophobicity of environments around Trp-14 decreased, whereas that around Trp-7 barely changed in the primary protein response. The present data indicate that displacement of the E helix toward the heme occurs within the instrument response time and that movement of the FG corner takes place with a time constant of 2 ps. The finding that the instantaneous motion of the E helix strongly suggests a mechanism in which protein structural changes are propagated from the heme to the A helix through the E helix motion.

The formation of hydrogen bond in the proximal heme pocket of HemAT-Bs upon ligand binding

HemAT-Bs is the heme-based O(2) sensor responsible for aerotaxis control in Bacillus subtilis. In this study, we measured the time-resolved resonance Raman spectra of full-length HemAT-Bs wild-type (WT) and Y133F in the deoxy form and the photoproduct after photolysis of CO-bound form. In WT, the nu(Fe-His) band for the 10 ps photoproduct was observed at higher frequency by about 2 cm(-1) compared with that of the deoxy form. This frequency difference is relaxed in hundreds of picoseconds. This time-dependent frequency shift would reflect the conformational change of the protein matrix. On the other hand, Y133F mutant did not show such a substantial nu(Fe-His) frequency shift after photolysis. Since a hydrogen bond to the proximal His induces an up-shift of the nu(Fe-His) frequency, these results indicate that Tyr133 forms a hydrogen bond to the proximal His residue upon the ligand binding. We discuss a functional role of this hydrogen bond formation for the signal transduction in HemAT-Bs.

Molecular dynamics study on the solvent dependent heme cooling following ligand photolysis in carbonmonoxy myoglobin

The time scale and mechanism of vibrational energy relaxation of the heme moiety in myoglobin was studied using molecular dynamics simulation. Five different solvent models, including normal water, heavy water, normal glycerol, deuterated glycerol and a nonpolar solvent, and two forms of the heme, one native and one lacking acidic side chains, were studied. Structural alteration of the protein was observed in native myoglobin glycerol solution and native myoglobin water solution. The single-exponential decay of the excess kinetic energy of the heme following ligand photolysis was observed in all systems studied. The relaxation rate depends on the solvent used. However, this dependence cannot be explained using bulk transport properties of the solvent including macroscopic thermal diffusion. The rate and mechanism of heme cooling depends upon the detailed microscopic interaction between the heme and solvent. Three intermolecular energy transfer mechanisms were considered: (i) energy transfer mediated by hydrogen bonds, (ii) direct vibration-vibration energy transfer via resonant interaction, and (iii) energy transfer via vibration-translation or vibration-rotation interaction, or in other words, thermal collision. The hydrogen bond interaction and vibration-vibration interaction between the heme and solvent molecules dominates the energy transfer in native myoglobin aqueous solution and native myoglobin glycerol solutions. For modified myoglobin, the vibration-vibration interaction is also effective in glycerol solution, different from aqueous solution. Thermal collisions form the dominant energy transfer pathway for modified myoglobin in water solution, and for both native myoglobin and modified myoglobin in a nonpolar environment. For native myoglobin in a nonpolar solvent solution, hydrogen bonds between heme isopropionate side chains and nearby protein residues, absent in the modified myoglobin nonpolar solvent solution, are key interactions influencing the relaxation pathways.

State Preparation and Excited Electronic and Vibrational Behavior in Hemes

The temporally overlapping, ultrafast electronic and vibrational dynamics of a model five-coordinate, high-spin heme in a nominally isotropic solvent environment has been studied for the first time with three complementary ultrafast techniques: transient absorption, time-resolved resonance Raman Stokes, and time-resolved resonance Raman anti-Stokes spectroscopies. Vibrational dynamics associated with an evolving ground-state species dominate the observations. Excitation into the blue side of the Soret band led to very rapid S2 --> S1 decay (sub-100 fs), followed by somewhat slower (800 fs) S1 --> S0 nonradiative decay. The initial vibrationally excited, non-Boltzmann S0 state was modeled as shifted to lower energy by 300 cm(-1) and broadened by 20%. On a approximately 10 ps time scale, the S0 state evolved into its room-temperature, thermal distribution S0 profile largely through VER. Anti-Stokes signals disappear very rapidly, indicating that the vibrational energy redistributes internally in about 1-3 ps from the initial accepting modes associated with S1 --> S0 internal conversion to the rest of the macrocycle. Comparisons of anti-Stokes mode intensities and lifetimes from TRARRS studies in which the initial excited state was prepared by ligand photolysis [Mizutani, T.; Kitagawa, T. Science 1997, 278, 443, and Chem. Rec. 2001, 1, 258] suggest that, while transient absorption studies appear to be relatively insensitive to initial preparation of the electronic excited state, the subsequent vibrational dynamics are not. Direct, time-resolved evaluation of vibrational lifetimes provides insight into fast internal conversion in hemes and the pathways of subsequent vibrational energy flow in the ground state. The overall similarity of the model heme electronic dynamics to those of biological systems may be a sign that the protein's influence upon the dynamics of the heme active site is rather subtle.

Pathway of Information Transmission from Heme to Protein upon Ligand Binding/Dissociation in Myoglobin Revealed by UV Resonance Raman Spectroscopy

Gas sensory heme proteins respond to their environment by binding a specific gas molecule to heme and transmitting this primary binding signal to the protein. How the binding signal is transmitted from the heme to the protein remains to be clarified. Using UV resonance Raman (UVRR) spectroscopy, we investigated this pathway in sperm whale myoglobin as a model gas sensory heme protein. Based on the UVRR data and the effects of deleting one of three important pathways (His-93, 6-propionate, or 7-propionate), we determined the changes in the conformation of globin that occur upon binding of CO, nitric oxide (NO), or O(2) to heme and how they are transmitted from heme to globin. The UVRR results show that heme discriminates different ligands, resulting in different conformations in the globin protein. Specifically, NO induces changes in the spectrum of Trp residues in the A-helix that are significantly different from those induced by O(2) or CO binding. On the other hand, binding of O(2) to heme produces changes in the Tyr residues of the H-helix that are different from those induced by CO or NO binding. Furthermore, we found that cleavage of the Fe-His-93 covalent bond eliminates communication to the terminal region of the H-helix and that the 7-propionate hydrogen-bonding network is essential for transmitting the CO or NO binding signal to the N and C termini. Finally, the 6-propionate is important only for NO binding. Thus, the hydrogen-bonding network in the protein appears to be critical for intramolecular signal transduction in gas sensory heme proteins.

Evidence for displacements of the C-helix by CO ligation and DNA binding to CooA revealed by UV resonance Raman spectroscopy

The UV and visible resonance Raman spectra are reported for CooA from Rhodospirillum rubrum, which is a transcriptional regulator activated by growth in a CO atmosphere. CO binding to heme in its sensor domain causes rearrangement of its DNA-binding domain, allowing binding of DNA with a specific sequence. The sensor and DNA-binding domains are linked by a hinge region that follows a long C-helix. UV resonance Raman bands arising from Trp-110 in the C-helix revealed local movement around Trp-110 upon CO binding. The indole side chain of Trp-110, which is exposed to solvent in the CO-free ferrous state, becomes buried in the CO-bound state with a slight change in its orientation but maintains a hydrogen bond with a water molecule at the indole nitrogen. This is the first experimental data supporting a previously proposed model involving displacement of the C-helix and heme sliding. The UV resonance Raman spectra for the CooA-DNA complex indicated that binding of DNA to CooA induces a further displacement of the C-helix in the same direction during transition to the complete active conformation. The Fe-CO and C-O stretching bands showed frequency shifts upon DNA binding, but the Fe-His stretching band did not. Moreover, CO-geminate recombination was more efficient in the DNA-bound state. These results suggest that the C-helix displacement in the DNA-bound form causes the CO binding pocket to narrow and become more negative.

Quaternary structures of intermediately ligated human hemoglobin a and influences from strong allosteric effectors: resonance Raman investigation

The Fe-histidine stretching (nu(Fe-His)) frequency was determined for deoxy subunits of intermediately ligated human hemoglobin A in equilibrium and CO-photodissociated picosecond transient species in the presence and absence of strong allosteric effectors like inositol(hexakis)phosphate, bezafibrate, and 2,3-bisphosphoglycerate. The nu(Fe-His) frequency of deoxyHb A was unaltered by the effectors. The T-to-R transition occurred around m = 2-3 in the absence of effectors but m > 3.5 in their presence, where m is the average number of ligands bound to Hb and was determined from the intensity of the nu(4) band measured in the same experiment. The alpha1-beta2 subunit contacts revealed by ultraviolet resonance Raman spectra, which were distinctly different between the T and R states, remained unchanged by the effectors. This observation would solve the recent discrepancy that the strong effectors remove the cooperativity of oxygen binding in the low-affinity limit, whereas the (1)H NMR spectrum of fully ligated form exhibits the pattern of the R state.

Interactions of soluble guanylate cyclase with diatomics as probed by resonance Raman spectroscopy

Soluble guanylate cyclase (sGC, EC 4.6.1.2) acts as a sensor for nitric oxide (NO), but is also activated by carbon monoxide in the presence of an allosteric modulator. Resonance Raman studies on the structure-function relations of sGC are reviewed with a focus on the CO-adduct in the presence and absence of allosteric modulator, YC-1, and substrate analogues. It is demonstrated that the sGC isolated from bovine lung contains one species with a five-coordinate (5c) ferrous high-spin heme with the Fe-His stretching mode at 204 cm(-1), but its CO adduct yields two species with different conformations about the heme pocket with the Fe-CO stretching (nuFe-CO) mode at 473 and 489 cm(-1), both of which are His- and CO-coordinated 6c ferrous adducts. Addition of YC-1 to it changes their population and further addition of GTP yields one kind of 6c (nuFe-CO=489 cm(-1)) in addition to 5c CO-adduct (nuFe-CO=521 cm(-1)). Under this condition the enzymatic activity becomes nearly the same level as that of NO adduct. Addition of gamma-S-GTP yields the same effect as GTP does but cGMP and GDP gives much less effects. Unexpectedly, ATP cancels the effects of GTP. The structural meaning of these spectroscopic observations is discussed in detail.

Kinetic modulation in carbonmonoxy derivatives of truncated hemoglobins: the role of distal heme pocket residues and extended apolar tunnel

Truncated hemoglobins (trHbs), are a distinct and newly characterized class of small myoglobin-like proteins that are widely distributed in bacteria, unicellular eukaryotes, and higher plants. Notable and distinctive features associated with trHbs include a hydrogen-bonding network within the distal heme pocket and a long apolar tunnel linking the external solvent to the distal heme pocket. The present work compares the geminate and solvent phase rebinding kinetics from two trHbs, one from the ciliated protozoan Paramecium caudatum (P-trHb) and the other from the green alga Chlamydomonas eugametos (C-trHb). Unusual kinetic patterns are observed including indications of ultrafast (picosecond) geminate rebinding of CO to C-trHb, very fast solvent phase rebinding of CO for both trHbs, time-dependent biphasic CO rebinding kinetics for P-trHb at low CO partial pressures, and for P-trHb, an increase in the geminate yield from a few percent to nearly 100% under high viscosity conditions. Species-specific differences in both the 8-ns photodissociation quantum yield and the rebinding kinetics, point to a pivotal functional role for the E11 residue. The response of the rebinding kinetics to temperature, ligand concentration, and viscosity (glycerol, trehalose) and the viscosity-dependent changes in the resonance Raman spectrum of the liganded photoproduct, together implicate both the apolar tunnel and the static and dynamic properties of the hydrogen-bonding network within the distal heme pocket in generating the unusual kinetic patterns observed for these trHbs.