Dynamics of Myoglobin−CO with the Proximal Histidine Removed: Vibrational Echo Experiments

Computational infrared and two-dimensional infrared photon echo spectroscopy of both wild-type and double mutant myoglobin-CO proteins

The CO stretching mode of both wild-type and double mutant ( T67R / S92D ) MbCO (carbonmonoxymyoglobin) proteins is an ideal infrared (IR) probe for studying the local electrostatic environment inside the myoglobin heme pocket. Recently, to elucidate the conformational switching dynamics between two distinguishable states, extensive IR absorption, IR pump-probe, and two-dimensional (2D) IR spectroscopic studies for various mutant MbCO's have been performed by the Fayer group. They showed that the 2D IR spectroscopy of the double mutant, which has a peroxidase enzyme activity, reveals a rapid chemical exchange between two distinct states, whereas that of the wild-type does not. Despite the fact that a few simulation studies on these systems were already performed and reported, such complicated experimental results have not been fully reproduced nor described in terms of conformational state-to-state transition processes. Here, we first develop a distributed vibrational solvatochromic charge model for describing the CO stretch frequency shift reflecting local electric potential changes. Then, by carrying out molecular dynamic simulations of the two MbCO's and examining their CO frequency trajectories, it becomes possible to identify a proper reaction coordinate consisting of His64 imidazole ring rotation and its distance to the CO ligand. From the 2D surfaces of the resulting potential of mean forces, the spectroscopically distinguished A1 and A3 states of the wild-type as well as two more substates of the double mutant are identified and their vibrational frequencies and distributions are separately examined. Our simulated IR absorption and 2D IR spectra of the two MbCO's are directly compared with the previous experimental results reported by the Fayer group. The chemical exchange rate constants extracted from the two-state kinetic analyses of the simulated 2D IR spectra are in excellent agreement with the experimental values. On the basis of the quantitative agreement between the simulated spectra and experimental ones, we further examine the conformational differences in the heme pockets of the two proteins and show that the double mutation, T67R / S92D , suppresses the A1 population, restricts the imidazole ring rotation, and increases hydrogen-bond strength between the imidazole Nε-H and the oxygen atom of the CO ligand. It is believed that such delicate change of distal His64 imidazole ring dynamics induced by the double mutation may be responsible for its enhanced peroxidase catalytic activity as compared to the wild-type myoglobin.

Cytochrome c552 Mutants: Structure and Dynamics at the Active Site Probed by Multidimensional NMR and Vibration Echo Spectroscopy

Spectrally resolved infrared stimulated vibrational echo experiments are used to measure the vibrational dephasing of a CO ligand bound to the heme cofactor in two mutated forms of the cytochrome c552 from Hydrogenobacter thermophilus. The first mutant (Ht-M61A) is characterized by a single mutation of Met61 to an Ala (Ht-M61A), while the second variant is doubly modified to have Gln64 replaced by an Asn in addition to the M61A mutation (Ht-M61A/Q64N). Multidimensional NMR experiments determined that the geometry of residue 64 in the two mutants is consistent with a non-hydrogen-bonding and hydrogen-bonding interaction with the CO ligand for Ht-M61A and Ht-M61A/Q64N, respectively. The vibrational echo experiments reveal that the shortest time scale vibrational dephasing of the CO is faster in the Ht-M61A/Q64N mutant than that in Ht-M61A. Longer time scale dynamics, measured as spectral diffusion, are unchanged by the Q64N modification. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to confirm that the dynamical difference induced by the Q64N mutation is primarily an increase in the fast (hundreds of femtoseconds) frequency fluctuations, while the slower (tens of picoseconds) dynamics are nearly unaffected. We conclude that the faster dynamics in Ht-M61A/Q64N are due to the location of Asn64, which is a hydrogen bond donor, above the heme-bound CO. A similar difference in CO ligand dynamics has been observed in the comparison of the CO derivative of myoglobin (MbCO) and its H64V variant, which is caused by the difference in axial residue interactions with the CO ligand. The results suggest a general trend for rapid ligand vibrational dynamics in the presence of a hydrogen bond donor.

Dynamics of proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.

Redox-Coupled Dynamics and Folding in Cytochrome c

Cytochrome c functions as an electron carrier in the mitochondrial electron-transport chain using the Fe(II)-Fe(III) redox couple of a covalently attached heme prosthetic group, and it has served as a paradigm for both biological redox activity and protein folding. On the basis of a wide variety of biophysical techniques, it has been suggested that the protein is more flexible in the oxidized state than in the reduced state, which has led to speculation that it is the dynamics of the protein that has been evolved to control the cofactor's redox properties. To test this hypothesis, we incorporated carbon-deuterium bonds throughout cytochrome c and characterized their absorption frequencies and line widths using IR spectroscopy. The absorption frequencies of several residues on the proximal side of the heme show redox-dependent changes, but none show changes in line width, implying that the flexibility of the oxidized and reduced proteins is not different. However, the spectra demonstrate that folded protein is in equilibrium with a surprisingly large amount of locally unfolded protein, which increases with oxidation for residues localized to the proximal side of the heme. The data suggest that while the oxidized protein is not more flexible than the reduced protein, it is more locally unfolded. Local unfolding of cytochrome c might be one mechanism whereby the protein evolved to control electron transfer.

The Influence of Aqueous versus Glassy Solvents on Protein Dynamics: Vibrational Echo Experiments and Molecular Dynamics Simulations

Spectrally resolved infrared stimulated vibrational echo measurements are used to measure the vibrational dephasing of the CO stretching mode of carbonmonoxy-hemoglobin (HbCO), a myoglobin mutant (H64V), and a bacterial cytochrome c(552) mutant (Ht-M61A) in aqueous solution and trehalose glasses. The vibrational dephasing of the heme-bound CO is significantly slower for all three proteins embedded in trehalose glasses compared to that of aqueous protein solutions. All three proteins exhibit persistent but notably slower spectral diffusion when the protein surface is fixed by the glassy solvent. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to reveal that the structural dynamics, as sensed by the CO, of the three proteins in trehalose and aqueous solution are dominated by fast (tens of femtoseconds), motionally narrowed fluctuations. MD simulations of H64V in dynamic and "static" water are presented as models of the aqueous and glassy environments. FFCFs are calculated from the H64V simulations and qualitatively reproduce the important features of the experimentally extracted FFCFs. The suppression of long time scale (picoseconds to tens of picoseconds) frequency fluctuations (spectral diffusion) in the glassy solvent is the result of a damping of atomic displacements throughout the protein structure and is not limited to structural dynamics that occur only at the protein surface. The analysis provides evidence that some dynamics are coupled to the hydration shell of water, supporting the idea that the bioprotection offered by trehalose is due to its ability to immobilize the protein surface through a thin, constrained layer of water.

Adsorption of Atomic Hydrogen at a Nanostructured Electrode of Polyacrylate-Capped Pt Nanoparticles in Polyelectrolyte

Atomic hydrogen electrosorption is reported at crystallite sites of polyacrylate-capped Pt nanoparticles (d = 2.5 +/- 0.6 nm), by assembling nanostructured electrodes of polyacrylate-Pt nanocrystallites layer-by-layer in a cationic polyelectrolyte, poly(diallyldimethylammonium chloride). Cyclic voltammetry in 1 M H2SO4 revealed a strongly adsorbed hydrogen state and a weakly adsorbed hydrogen state assigned to adsorption at (100) and (110) sites of the modified nanocrystallites, respectively. Resolving hydrogen adsorption states signifies that surface capping by the carboxylate groups is not irreversibly blocking hydrogen adsorption sites at the modified Pt nanoparticle surface. Adsorption peak currents increased with increasing the number of layers up to 16 bilayers, indicating the feasibility of nanoparticle charging via interparticle charge hopping and the accessibility of adsorption states within the thickness of the nanoparticle/polyelectrolyte multilayers. Despite similarity in hydrogen adsorption in the cyclic voltammorgrams in 1 M H2SO4, negative shifts in adsorption potentials were measured at the nanocrystallite Pt-polyelectrolyte multilayers relative to a polycrystalline bulk Pt surface. This potential shift is attributed to a kinetic limitation in the reductive hydrogen adsorption as a result of the Pt nanoparticle surface modification and the polyelectrolyte environment.

Human myoglobin recognition of oxygen: dynamics of the energy landscape

Femtosecond to nanosecond dynamics of O(2) rebinding to human WT myoglobin and its mutants, V68F and I107F, have been studied by using transient absorption. The results are compared with NO rebinding. Even though the immediate environment around the heme binding site is changed by the mutations, the picosecond geminate rebinding of oxygen is at most minimally affected. On the other hand, the V68F (E11) mutation causes drastic differences in rebinding on the nanosecond time scale, whereas the effect of the I107F (G8) mutation remains relatively small within our 10-ns time window. Unlike traditional homogeneous kinetics and molecular dynamics collisional simulations, we propose a "bifurcation model" for populations of directed and undirected dynamics on the ultrafast time scale, reflecting the distribution of initial protein conformations. The major mutation effect occurs on the time scale on which global protein conformational change is possible, consistent with transitions between the conformations of directed and undirected population playing a role in the O(2) binding. We discuss the relevance of these findings to the bimolecular function of the protein.

Higher-order optical correlation spectroscopy in liquids

Linear optical spectroscopies have long been used to study the behavior of liquids. Laser technology has progressed to the point that it has become possible to perform nonlinear optical experiments that probe higher-order correlation functions in liquids, opening a new window into our understanding of the microscopic details of solution-phase processes. Here we review advances that have been made in recent years in employing higher-order electronic and vibrational spectroscopies to study liquid-state dynamics and structure.

Myoglobin-CO conformational substate dynamics: 2D vibrational echoes and MD simulations

Two-dimensional (2D) infrared vibrational echoes were performed on horse heart carbonmonoxymyoglobin (MbCO) in water over a range of temperatures. The A(1) and A(3) conformational substates of MbCO are found to have different dephasing rates with different temperature dependences. A frequency-frequency correlation function derived from molecular dynamics simulations on MbCO at 298 K is used to calculate the vibrational echo decay. The calculated decay shows substantial agreement with the experimentally measured decays. The 2D vibrational echo probes protein dynamics and provides an observable that can be used to test structural assignments for the MbCO conformational substates.

Fast protein dynamics probed with infrared vibrational echo experiments

IR vibrational echo experiments are used to study dynamics in myoglobin (Mb) by investigating the dephasing of the CO-stretching mode of CO bound at the active site of the protein (Mb-CO). The temperature dependence and the viscosity dependence of Mb-CO pure dephasing have been measured in several solvents. In low-temperature, glassy solvents, the pure dephasing has a power law temperature dependence, T(1.3), that reflects glasslike protein dynamics. In liquids, the temperature dependence is much steeper and arises from a combination of pure temperature dependence and the influence of decreasing solvent viscosity with increasing temperature. As the solvent viscosity decreases, the ability of the protein's surface to undergo topological fluctuations increases, which in turn increases the internal protein-structural fluctuations. The protein-structural motions are coupled to the CO bound at the active site by electric field fluctuations that accompany movements of polar residues. The dynamic electric field-coupling mechanism is tested by observing differences in the temperature dependence of the pure dephasing of Mb-CO mutations.

Replacement of the axial histidine ligand with imidazole in cytochrome c peroxidase. 2. Effects on heme coordination and function

The inability of imidazole to complement function in the axial histidine deletion mutant, H175G, of yeast cytochrome c peroxidase has been an intriguing but unresolved issue that impacts our understanding of the role of axial ligands in heme catalysis. Here we report the functional and spectroscopic properties of H175G and of its complexes with imidazole. Combined with the crystal structures for these complexes, the data provide a detailed and consistent account of the modes of Im binding in the H175G cavity and their dependence on buffer and pH. UV--vis, EPR, and resonance Raman spectra reveal multiple coordination states for H175G/Im which can be correlated with the crystal structures to assign the following heme environments: H175G/H(2)O/H(2)O, H175G/Im(d)/phosphate(c), H175G/Im(d)/H(2)O(c), H175G/Im(c)/H(2)O(d), and H175G/Im(c)/OH(-)(c), where H175G/X/Y defines the proximal species as X and the distal species as Y and c and d subscripts refer, where known, to the coordinated and dissociated states, respectively. Resonance Raman data for reduced H175G/Im show two substates for heme-coordinated Im differing in the strength of their hydrogen bond to Asp-235, in a fashion similar to WT CCP. NO binding to ferrous H175G/Im results in dissociation of Im from the heme but not from the cavity, while no dissociation is observed for WT CCP, indicating that steric tethering may, in part, control NO-induced dissociation of trans ligands. H175G/Im forms an oxidized compound I state with two distinct radical species, each with a dramatically different anisotropy and spin relaxation from that of the Trp-191 radical of WT CCP. It is suggested that these signals arise from alternate conformations of Trp191 having different degrees of exchange coupling to the ferryl heme, possibly mediated by the conformational heterogeneity of Im within the H175G cavity. The kinetics of the reaction of H175G/Im with H(2)O(2) are multiphasic, also reflecting the multiple coordination states of Im. The rate of the fastest phase is essentially identical to that of WT CCP, indicating that the H175G/Im(c)/H(2)O(d) state is fully reactive with peroxide. However, the overall rate of enzyme turnover using cytochrome c as a substrate is <5% of WT and is unaffected by Im coordination. In summary, Im coordination to H175G results in a number of conformers, one of which is structurally and spectroscopically very similar to WT CCP. However, while this form is fully reactive with peroxide, the reaction with cytochrome c remains inefficient, perhaps implicating the altered Trp-191 radical species.

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the approximately 10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.