The contribution of heme propionate groups to the conformational dynamics associated with CO photodissociation from horse heart myoglobin

Photoacoustic calorimetry and transient absorption spectroscopy were used to study conformational dynamics associated with CO photodissociation from horse heart myoglobin (Mb) reconstituted with either Fe protoporphyrin IX dimethylester (FePPDME), Fe octaethylporphyrin (FeOEP), or with native Fe protoporphyrin IX (FePPIX). The volume and enthalpy changes associated with the Fe-CO bond dissociation and formation of a transient deoxyMb intermediate for the reconstituted Mbs were found to be similar to those determined for native Mb (DeltaV1 = -2.5+/-0.6 ml mol(-1) and DeltaH1 = 8.1+/-3.0 kcal mol(-1)). The replacement of FePPIX by FeOEP significantly alters the conformational dynamics associated with CO release from protein. Ligand escape from FeOEP reconstituted Mb was determined to be roughly a factor of two faster (tau=330 ns) relative to native protein (tau=700 ns) and accompanying reaction volume and enthalpy changes were also found to be smaller (DeltaV2 = 5.4+/-2.5 ml mol(-1) and DeltaH2 = 0.7+/-2.2 kcal mol(-1)) than those for native Mb (DeltaV2 = 14.3+/-0.8 ml mol(-1) and DeltaH2 = 7.8+/-3.5 kcal mol(-1)). On the other hand, volume and enthalpy changes for CO release from FePPIX or FePPDME reconstituted Mb were nearly identical to those of the native protein. These results suggest that the hydrogen bonding network between heme propionate groups and nearby amino acid residues likely play an important role in regulating ligand diffusion through protein matrix. Disruption of this network leads to a partially open conformation of protein with less restricted ligand access to the heme binding pocket.     

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3-7 cm(-1) (compared to the experimental value of 10 cm(-1)) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe...CO and Fe...OC) is 0.3 kcal mol(-1) and agrees well with the experimentally estimated value of 0.29 kcal mol(-1). A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.     

The effect of ligand dynamics on heme electronic transition band III in myoglobin.

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.     

Evidence for fast conformational change upon ligand dissociation in the HemAT class of bacterial oxygen sensors

Here we report the results of transient absorption and photoacoustic calorimetry studies of CO photodissociation from the heme domain of the bacterial oxygen sensor HemAT-Bs. The results indicate that CO photolysis is accompanied by an overall DeltaH of -19 kcal mol(-1) and DeltaV of +4 ml mol(-1) as well as a red-shifted kinetic difference spectrum all occurring in <50 ns. Analysis of the DeltaH/DeltaV reveals that a conformational change takes place with a DeltaH(conf) of -40 kcal mol(-1) and DeltaV(conf) of -22 ml mol(-1). These thermodynamic changes are consistent with an increase in the solvent accessible surface area of the protein upon ligand dissociation, as observed in the X-ray structure of the ferric CN-bound and CN free forms of HemAT-Bs.     

Ligand selectivity of soluble guanylyl cyclase: effect of the hydrogen-bonding tyrosine in the distal heme pocket on binding of oxygen, nitric oxide, and carbon monoxide

Although soluble guanylyl cyclase (sGC) functions in an environment in which O(2), NO, and CO are potential ligands for its heme moiety, the enzyme displays a high affinity for only its physiological ligand, NO, but has a limited affinity for CO and no affinity for O(2). Recent studies of a truncated version of the sGC beta(1)-subunit containing the heme-binding domain (Boon, E. M., Huang, S H., and Marletta, M. A. (2005) Nat. Chem. Biol., 1, 53-59) showed that introduction of the hydrogen-bonding tyrosine into the distal heme pocket changes the ligand specificity of the heme moiety and results in an oxygen-binding sGC. The hypothesis that the absence of hydrogen-bonding residues in the distal heme pocket is sufficient to provide oxygen discrimination by sGC was put forward. We tested this hypothesis in a context of a complete sGC heterodimer containing both the intact alpha(1)- and beta(1)-subunits. We found that the I145Y substitution in the full-length beta-subunit of the sGC heterodimer did not produce an oxygen-binding enzyme. However, this substitution impeded the association of NO and destabilized the NO.heme complex. The tyrosine in the distal heme pocket also impeded both the binding and dissociation of the CO ligand. We propose that the mechanism of oxygen exclusion by sGC not only involves the lack of hydrogen bonding in the distal heme pocket, but also depends on structural elements from other domains of sGC.     

Protein Collective Motions Coupled to Ligand Migration in Myoglobin

Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.     

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.     

On CO Egress from,and Re‐Uptake by,the Enzyme MauG,as a Mimic of the Acquisition of Oxidizing Agents by the pre‐MADH_MauG System. A Molecular Mechanics Approach

In this work, by applying a non‐deterministic, randomly‐oriented minimal force to the dissociated CO ligand of the MauG‐CO system, the molecular‐dynamics (MD) behavior of this system could be quickly unraveled. It turned out that CO has no marked directional egress from the high‐spin c‐heme iron distal pocket. Rather, CO is able to exploit all interstices created during the protein fluctuations. Nonetheless, no steady route toward the surrounding solvent was ever observed: CO jumped first into other binding pockets before being able to escape the protein. In a few cases, on hitting the surrounding H₂O molecules, CO was observed to reverse direction, re‐entering the protein. A contention that conformational inversion of the P107 ring provides a gate to the iron ion is not supported by the present simulations.     

Molecular dynamics simulation of sperm whale myoglobin: effects of mutations and trapped CO on the structure and dynamics of cavities

The results of extended (80-ns) molecular dynamics simulations of wild-type and YQR triple mutant of sperm whale deoxy myoglobin in water are reported and compared with the results of the simulation of the intermediate(s) obtained by photodissociation of CO in the wild-type protein. The opening/closure of pathways between preexistent cavities is different in the three systems. For the photodissociated state, we previously reported a clear-cut correlation between the opening probability and the presence of the photolyzed CO in the proximity of the passage; here we show that in wild-type deoxy myoglobin, opening is almost random. In wild-type deoxy myoglobin, the passage between the distal pocket and the solvent is strictly correlated to the presence/absence of a water molecule that simultaneously interacts with the distal histidine side chain and the heme iron; conversely, in the photodissociated myoglobin, the connection with the bulk solvent is always open when CO is in the vicinity of the A pyrrole ring. In YQR deoxy myoglobin, the mutated Gln(E7)64 is stably H-bonded with the mutated Tyr(B10)29. The essential dynamics analysis unveils a different behavior for the three systems. The motion amplitude is progressively restricted in going from wild-type to YQR deoxy myoglobin and to wild-type myoglobin photoproduct. In all cases, the principal motions involve mainly the same regions, but their directions are different. Analysis of the dynamics of the preexisting cavities indicates large fluctuations and frequent connections with the solvent, in agreement with the earlier hypothesis that some of the ligand may escape from the protein through these pathways.     

Proton-coupled structural changes upon binding of carbon monoxide to cytochrome cd1: a combined flash photolysis and X-ray crystallography study

We have investigated dynamic events after flash photolysis of CO from reduced cytochrome cd(1) nitrite reductase (NiR) from Paracoccus pantotrophus (formerly Thiosphaera pantotropha). Upon pulsed illumination of the cytochrome cd(1)-CO complex, at 460 nm, a rapid (<50 ns) absorbance change, attributed to dissociation of CO, was observed. This was followed by a biphasic rearrangement with rate constants of 1.7 x 10(4) and 2.5 x 10(3) s(-1) at pH 8.0. Both parts of the biphasic rearrangement phases displayed the same kinetic difference spectrum in the region of 400-660 nm. The slower of the two processes was accompanied by proton uptake from solution (0.5 proton per active site at pH 7.5-8.5). After photodissociation, the CO ligand recombined at a rate of 12 s(-1) (at 1 mM CO and pH 8.0), accompanied by proton release. The crystal structure of reduced cytochrome cd(1) in complex with CO was determined to a resolution of 1.57 A. The structure shows that CO binds to the iron of the d(1) heme in the active site. The ligation of the c heme is unchanged in the complex. A comparison of the structures of the reduced, unligated NiR and the NiR-CO complex indicates changes in the puckering of the d(1) heme as well as rearrangements in the hydrogen-bonding network and solvent organization in the substrate binding pocket at the d(1) heme. Since the CO ligand binds to heme d(1) and there are structural changes in the d(1) pocket upon CO binding, it is likely that the proton uptake or release observed after flash-induced CO dissociation is due to changes of the protonation state of groups in the active site. Such proton-coupled structural changes associated with ligand binding are likely to affect the redox potential of heme d(1) and may regulate the internal electron transfer from heme c to heme d(1).     

Temperature dependence of the backbone dynamics of ribonuclease A in the ground state and bound to the inhibitor 5'-phosphothymidine (3'-5')pyrophosphate adenosine 3'-phosphate

The interaction of the dinucleotide inhibitor 5'-phosphothymidine(3',5')pyrophosphate adenosine 3'-phosphate (pTppAp) with bovine pancreatic ribonuclease A (RNase A) was characterized by calorimetry and solution NMR spectroscopy. Calorimetric data show that binding of pTppAp to RNase A is exothermic (DeltaH = -60.1 +/- 4.1 kJ/mol) with a dissociation constant of 16 nM at 298 K. At this temperature, the binding results in an entropy loss (TDeltaS = -16.8 +/- 7.3 kJ/mol) that is more favorable than that with the product analogue, 2'-CMP (TDeltaS = -31.3 +/- 0.9 kJ/mol). Temperature-dependent calorimetric experiments give a DeltaC(p) for ligand binding of -230 +/- 100 J/mol K. Binding of pTppAp results in noticeable effects on the backbone amide chemical shifts and dynamics. Amide backbone (15)N NMR spin-relaxation studies were performed on both apo RNase A and RNase A/pTppAp as a function of temperature. At each temperature, the model-free-determined order parameters, S(2), were significantly higher for RNase A/pTppAp than for the apo enzyme indicating a decrease in the conformational entropy of the protein upon ligand binding. Furthermore, the magnitude of this difference varies along the amino acid sequence specifically locating the entropic changes. The temperature dependence of S(2) at each residue enabled assessment of the local heat capacity changes (DeltaC(p)) from ligand binding. In an overall, average sense, DeltaC(p) for the protein backbone, determined from the NMR dynamics measurements, did not differ between apo RNase A and RNase A/pTppAp indicating that backbone dynamics contribute little to DeltaC(p) for protein-ligand interactions in this system. However, residue-by-residue comparison of the temperature-dependent change in entropy (DeltaS(B)) between free and bound forms reveals nonzero contributions to DeltaC(p) at individual sites. The balance of positive and negative changes reveals a redistribution of energetics upon binding. Furthermore, experiment and semiempirical estimates suggest that a large negative DeltaC(p) should accompany binding of pTppAp, and we conclude that this contribution must arise from factors other than amide backbone dynamics.     

Picosecond phase grating spectroscopy of hemoglobin and myoglobin: energetics and dynamics of global protein motion.

Phase grating spectroscopy has been used to follow the optically triggered tertiary structural changes of carboxymyoglobin (MbCO) and carboxyhemoglobin (HbCO). Probe wavelength and temperature dependencies have shown that the grating signal arises from nonthermal density changes induced by the protein structural changes. The material displaced through the protein structural changes leads to the excitation of coherent acoustic modes of the surrounding water. The coupling of the structural changes to the fluid hydrodynamics demonstrates that a global change in the protein structure is occurring in less than 30 ps. The global relaxation is on the same time scale as the local changes in structure in the vicinity of the heme pocket. The observed dynamics for global relaxation and correspondence between the local and global structural changes provides evidence for the involvement of collective modes in the propagation of the initial tertiary conformational changes. The energetics can also be derived from the acoustic signal. For MbCO, the photodissociation process is endothermic by 21 +/- 2 kcal/mol, which corresponds closely to the expected Fe-CO bond enthalpy. In contrast, HbCO dissipates approximately 10 kcal/mol more energy relative to myoglobin during its initial tertiary structural relaxation. The difference in energetics indicates that significantly more energy is stored in the hemoglobin structure and is believed to be related to the quaternary structure of hemoglobin not present in the monomeric form of myoglobin. These findings provide new insight into the biomechanics of conformational changes in proteins and lend support to theoretical models invoking stored strain energy as the driving force for large amplitude correlated motions.     

On CO Egress from, and Re-Uptake by, the Enzyme MauG, as a Mimic of the Acquisition of Oxidizing Agents by the pre-MADH_MauG System. A Molecular Mechanics Approach

In this work, by applying a non-deterministic, randomly-oriented minimal force to the dissociated CO ligand of the MauG-CO system, the molecular-dynamics (MD) behavior of this system could be quickly unraveled. It turned out that CO has no marked directional egress from the high-spin c-heme iron distal pocket. Rather, CO is able to exploit all interstices created during the protein fluctuations. Nonetheless, no steady route toward the surrounding solvent was ever observed: CO jumped first into other binding pockets before being able to escape the protein. In a few cases, on hitting the surrounding H(2) O molecules, CO was observed to reverse direction, re-entering the protein. A contention that conformational inversion of the P107 ring provides a gate to the iron ion is not supported by the present simulations.     

A comparative resonance Raman analysis of heme-binding PAS domains: heme iron coordination structures of the BjFixL, AxPDEA1, EcDos, and MtDos proteins

The heme-PAS is a specialized domain with which a broad class of signal-transducing heme proteins detect physiological heme ligands. Such domains exhibit a wide range of ligand binding parameters, yet they are all expected to feature an alpha-beta heme binding fold and a predominantly hydrophobic heme distal pocket without a distal histidine. We have compared, for the first time, the resonance Raman spectra of several heme-PASs: the heme-binding domains of Bradyrhizobium japonicum FixL, Escherichia coli Dos, Acetobacter xylinum PDEA1, and Methanobacterium thermoautotrophicum Dos. In all cases, the nu(Fe)-(CO) and nu(C-O) values of the carbonmonoxy forms were consistent with coordination of the heme iron to histidine on the proximal side and binding of the CO without electrostatic interaction with the heme distal pocket. EcDos was unusual in having predominantly hexacoordinate heme iron in the deoxy and met forms. Despite an evident lack of CO interaction with the EcDos heme pocket, relatively low Fe-O(2) (562 cm(-1)) and N-O (1576 cm(-1)) stretching frequencies indicated that strong polar interactions with that heme distal pocket are possible for highly bent ligands such as O(2) or NO. None of the newly studied NO adducts exhibited evidence of the Fe-His rupture and pentacoordination previously noted for Sinorhizobium meliloti FixL. A low Fe-His stretching frequency, formerly interpreted as a strained Fe-His bond, and the slow association of O(2) with S. meliloti FixL failed to correlate with the newly studied proteins having low association rate or low equilibrium association constants for binding of O(2). We conclude that although heme-PASs share some features, they represent distinct signal transduction mechanisms.     

The position 68(E11) side chain in myoglobin regulates ligand capture, bond formation with heme iron, and internal movement into the xenon cavities

After photodissociation, ligand rebinding to myoglobin exhibits complex kinetic patterns associated with multiple first-order geminate recombination processes occurring within the protein and a simpler bimolecular phase representing second-order ligand rebinding from the solvent. A smooth transition from cryogenic-like to solution phase properties can be obtained by using a combination of sol-gel encapsulation, addition of glycerol as a bathing medium, and temperature tuning (-15 --> 65 degrees C). This approach was applied to a series of double mutants, myoglobin CO (H64L/V68X, where X = Ala, Val, Leu, Asn, and Phe), which were designed to examine the contributions of the position 68(E11) side chain to the appearance and disappearance of internal rebinding phases in the absence of steric and polar interactions with the distal histidine. Based on the effects of viscosity, temperature, and the stereochemistry of the E11 side chain, the three major phases, B --> A, C --> A, and D --> A, can be assigned, respectively, to ligand rebinding from the following: (i) the distal heme pocket, (ii) the xenon cavities prior to large amplitude side chain conformational relaxation, and (iii) the xenon cavities after significant conformational relaxation of the position 68(E11) side chain. The relative amplitudes of the B --> A and C --> A phases depend markedly on the size and shape of the E11 side chain, which regulates sterically both ligand return to the heme iron atom and ligand migration to the xenon cavities. The internal xenon cavities provide a transient docking site that allows side chain relaxations and the entry of water into the vacated distal pocket, which in turn slows ligand recombination markedly.     

Structural dynamics of distal histidine replaced mutants of myoglobin accompanied with the photodissociation reaction of the ligand

Protein dynamics observed by the transient grating (TG) method are studied for some site-directed mutants at the distal histidine of myoglobin (H64L, H64Q, H64V). The time profiles of the TG signals are very sensitive to the amino acid residue of the 64 position. It was found that the sensitivity is mostly caused by the different rates of the ligand escape from the protein to solvent and the magnitude of the molecular volume change. Several molecular origins of the volume difference between MbCO and Mb, such as the electrostatic interaction in the distal pocket, movement of helices, and distal water, are proposed. Interestingly, the volume difference between the CO-trapped Mb inside the protein interior and Mb is similar to that of the partial molar volume of CO in organic solvent. The effect of mutation on the nature of the CO trapped site is discussed.     

Structural dynamics of myoglobin: spectroscopic and structural characterization of ligand docking sites in myoglobin mutant L29W

We have studied CO binding to the heme and CO migration among protein internal cavities after photodissociation in sperm whale carbonmonoxy myoglobin (MbCO) mutant L29W using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) and kinetic experiments at cryogenic temperatures. Photoproduct intermediates, characterized by CO at particular locations in the protein, were selectively enhanced by applying special laser illumination protocols. These studies were performed on the L29W mutant protein and a series of double mutants constructed so that bulky amino acid side chains block passageways between cavities or fill these sites. Binding of xenon was also employed as an alternative means of occluding cavities. All mutants exhibit two conformations, A(I) and A(II), with distinctly different photoproduct states and ligand binding properties. These differences arise mainly from different positions of the W29 and H64 side chains in the distal heme pocket [Ostermann, A., et al. (2000) Nature 404, 205-208]. The detailed knowledge of the interplay between protein structure, protein dynamics, and ligand migration at cryogenic temperatures allowed us to develop a dynamic model that explains the slow CO and O(2) bimolecular association observed after flash photolysis at ambient temperature.     

Volume and enthalpy profiles of CO rebinding to horse heart myoglobin

Carbon monoxide binding to myoglobin was characterized using the photothermal beam deflection method. The volume and enthalpy changes coupled to CO dissociation were found to be 9.3+/-0.8 mL x mol(-1) and 7.4+/-2.8 kcal x mol(-1), respectively. The corresponding values observed for CO rebinding have the same magnitude but opposite sign: Delta V=-8.6+/-0.9 mL x mol(-1) and Delta H=-5.8+/-2.9 kcal x mol(-1). Ligand rebinding occurs as a single conformational step with a rate constant of 5 x 10(5) M(-1) s(-1) and with activation enthalpy of 7.1+/-0.8 kcal x mol(-1) and activation entropy of -22.4+/-2.8 cal x mol(-1) K(-1). Activation parameters for the ligand binding correspond to the activation parameters previously obtained using the transient absorption methods. Hence, at room temperature the CO binding to Mb can be described as a two-state model and the observed volume contraction occurs during CO-Fe bond formation. Comparing these results with CO dissociation reactions, for which two discrete intermediates were characterized, indicates differences in mechanism by which the protein modulates ligand association and dissociation.