Distal Interactions in the Cyanide Complex of Ferric Chlamydomonas Hemoglobin

Vibronic Coupling in Vitamin B12: A Combined Spectroscopic and Computational Study

Understanding the diverse reactivities of vitamin B12 and its derivatives, collectively called cobalamins, requires detailed knowledge of their geometric and electronic structures. Electronic absorption (Abs) and resonance Raman (rR) spectroscopies have proven invaluable in this area, particularly when used in concert with computational techniques such as density functional theory (DFT). There remain, however, lingering uncertainties in the computational description of electronic excited states of cobalamins, particularly surrounding the vibronic coupling that impacts the Abs bandshapes and gives rise to rR enhancement of vibrational modes. Past computational analyses of the vibrational spectra of cobalamins have either neglected rR enhancement or calculated rR enhancement for only a small number of modes. In the present study, we used the recently developed ORCA_ASA computational tool in conjunction with the popular B3LYP and BP86 functionals to predict Abs bandshapes and rR spectra for vitamin B12. The ORCA_ASA/B3LYP-computed Abs envelope in the visible spectral region and rR spectra of vitamin B12 agree remarkably well with our experimental data, while BP86 fails to reproduce both. This finding represents a significant advance in our understanding of how these two commonly used density functionals differently model the electronic properties of cobalamins. Guided by the computed frequencies for the Co-C stretching and Co-C-N bending modes, we identified, for the first time, isotope-sensitive features in our rR spectra of 12CNCbl and 13CNCbl that can be assigned to these modes. A normal coordinate analysis of the experimentally determined Co-C stretching and Co-C-N bending frequencies indicates that the Co-C force constant for vitamin B12 is 2.67 mdyn/Å, considerably larger than the Co-C force constants reported for alkylcobalamins.

Distinguishing Active Site Characteristics of Chlorite Dismutases with Their Cyanide Complexes

O2-evolving chlorite dismutases (Clds) efficiently convert chlorite (ClO2-) to O2 and Cl-. Dechloromonas aromatica Cld ( DaCld) is a highly active chlorite-decomposing homopentameric enzyme, typical of Clds found in perchlorate- and chlorate-respiring bacteria. The Gram-negative, human pathogen Klebsiella pneumoniae contains a homodimeric Cld ( KpCld) that also decomposes ClO2-, albeit with an activity 10-fold lower and a turnover number lower than those of DaCld. The interactions between the distal pocket and heme ligand of the DaCld and KpCld active sites have been probed via kinetic, thermodynamic, and spectroscopic behaviors of their cyanide complexes for insight into active site characteristics that are deterministic for chlorite decomposition. At 4.7 × 10-9 M, the KD for the KpCld-CN- complex is 2 orders of magnitude smaller than that of DaCld-CN- and indicates an affinity for CN- that is greater than that of most heme proteins. The difference in CN- affinity between Kp- and DaClds is predominantly due to differences in koff. The kinetics of binding of cyanide to DaCld, DaCld(R183Q), and KpCld between pH 4 and 8.5 corroborate the importance of distal Arg183 and a p Ka of ∼7 in stabilizing complexes of anionic ligands, including the substrate. The Fe-C stretching and FeCN bending modes of the DaCld-CN- (νFe-C, 441 cm-1; δFeCN, 396 cm-1) and KpCld-CN- (νFe-C, 441 cm-1; δFeCN, 356 cm-1) complexes reveal differences in their FeCN angle, which suggest different distal pocket interactions with their bound cyanide. Conformational differences in their catalytic sites are also reported by the single ferrous KpCld carbonyl complex, which is in contrast to the two conformers observed for DaCld-CO.

Investigations of ferric heme cyanide photodissociation in myoglobin and horseradish peroxidase

The photodissociation of cyanide from ferric myoglobin (MbCN) and horseradish peroxidase (HRPCN) has definitively been observed. This has implications for the interpretation of ultrafast IR (Helbing et al. Biophys. J. 2004, 87, 1881-1891) and optical (Gruia et al. Biophys. J. 2008, 94, 2252-2268) studies that had previously suggested the Fe-CN bond was photostable in MbCN. The photolysis of ferric MbCN takes place with a quantum yield of ~75%, and the resonance Raman spectrum of the photoproduct observed in steady-state experiments as a function of laser power and sample spinning rate is identical to that of ferric Mb (metMb). The data are quantitatively analyzed using a simple model where cyanide is photodissociated and, although geminate rebinding with a rate of kBA ≈ (3.6 ps)(-1) is the dominant process, some CN(-) exits from the distal heme pocket and is replaced by water. Using independently determined values for the CN(-) association rate, we find that the CN(-) escape rate from the ferric myoglobin pocket to the solution at 293 K is kout ≈ (1-2) × 10(7) s(-1). This value is very similar to, but slightly larger than, the histidine gated escape rate of CO from Mb (1.1 × 10(7) s(-1)) under the same conditions. The analysis leads to an escape probability kout/(kout + kBA) ~ 10(-4), which is unobservable in most time domain kinetic measurements. However, the photolysis is surprisingly easy to detect in Mb using cw resonance Raman measurements. This is due to the anomalously slow CN(-) bimolecular association rate (170 M(-1) s(-1)), which arises from the need for water to exchange at the ferric heme binding site of Mb. In contrast, ferric HRP does not have a heme bound water molecule and its CN(-) bimolecular association rate is larger by ~10(3), making the CN(-) photolysis more difficult to observe.

H-bonding networks of the distal residues and water molecules in the active site of Thermobifida fusca hemoglobin

The ferric form of truncated hemoglobin II from Thermobifida fusca (Tf-trHb) and its triple mutant WG8F-YB10F-YCD1F at neutral and alkaline pH, and in the presence of CN(-) have been characterized by resonance Raman spectroscopy, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations. Tf-trHb contains three polar residues in the distal site, namely TrpG8, TyrCD1 and TyrB10. Whereas TrpG8 can act as a potential hydrogen-bond donor, the tyrosines can act as donors or acceptors. Ligand binding in heme-containing proteins is determined by a number of factors, including the nature and conformation of the distal residues and their capability to stabilize the heme-bound ligand via hydrogen-bonding and electrostatic interactions. Since both the RR Fe-OH(-) and Fe-CN(-) frequencies are very sensitive to the distal environment, detailed information on structural variations has been obtained. The hydroxyl ligand binds only the WT protein giving rise to two different conformers. In form 1 the anion is stabilized by H-bonds with TrpG8, TyrCD1 and a water molecule, in turn H-bonded to TyrB10. In form 2, H-bonding with TyrCD1 is mediated by a water molecule. Unlike the OH(-) ligand, CN(-) binds both WT and the triple mutant giving rise to two forms with similar spectroscopic characteristics. The overall results clearly indicate that H-bonding interactions both with distal residues and water molecules are important structural determinants in the active site of Tf-trHb. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Structural and Functional Properties of a Single Domain Hemoglobin from the Food-borne Pathogen Campylobactor jejuni

Campylobacter jejuni contains two globins, a truncated hemoglobin, Ctb, and a single domain hemoglobin, Cgb. The physiological function of Ctb remains unclear, whereas Cgb has been linked to NO detoxification. With resonance Raman scattering, the iron-histidine stretching mode of Cgb was identified at 251 cm(-1). This frequency is unusually high, suggesting an imidazolate character of the proximal histidine as a result of the H-bonding network linking the catalytic triad involving the F8His, H23Glu, and G5Tyr residues. In the CO-complex, two conformers were identified with the nuC-O/nuFe-CO at 529/1914 cm(-1) and 492/1963 cm(-1). The former is assigned to a "closed" conformation, in which the heme-bound CO is stabilized by the H-bond(s) donated from the B10Tyr-E7Gln residues, whereas the latter is assigned to an "open" conformer, in which the H-bonding interaction is absent. The presence of the two alternative conformations demonstrates the plasticity of the protein matrix. In the O2-complex, the iron-O2 stretching frequency was identified at 554 cm(-1), which is unusually low, indicating that the heme-bound O2 is stabilized by strong H-bond(s) donated by the B10Tyr-E7Gln residues. This scenario is consistent with its low O2 off-rate (0.87 s(-1)). Taken together the data suggest that the NO-detoxifying activity of Cgb is facilitated by the imidazolate character of the proximal F8His and the distal positive polar environment provided by the B10Tyr-E7Gln. They may offer electronic "push" and "pull," respectively, for the O-O bond cleavage reaction required for the isomerization of the presumed peroxynitrite intermediate to the product, nitrate.

Quantitative vibrational dynamics of iron in carbonyl porphyrins

We use nuclear resonance vibrational spectroscopy and computational predictions based on density functional theory (DFT) to explore the vibrational dynamics of (57)Fe in porphyrins that mimic the active sites of histidine-ligated heme proteins complexed with carbon monoxide. Nuclear resonance vibrational spectroscopy yields the complete vibrational spectrum of a Mössbauer isotope, and provides a valuable probe that is not only selective for protein active sites but quantifies the mean-squared amplitude and direction of the motion of the probe nucleus, in addition to vibrational frequencies. Quantitative comparison of the experimental results with DFT calculations provides a detailed, rigorous test of the vibrational predictions, which in turn provide a reliable description of the observed vibrational features. In addition to the well-studied stretching vibration of the Fe-CO bond, vibrations involving the Fe-imidazole bond, and the Fe-N(pyr) bonds to the pyrrole nitrogens of the porphyrin contribute prominently to the observed experimental signal. All of these frequencies show structural sensitivity to the corresponding bond lengths, but previous studies have failed to identify the latter vibrations, presumably because the coupling to the electronic excitation is too small in resonance Raman measurements. We also observe the FeCO bending vibrations, which are not Raman active for these unhindered model compounds. The observed Fe amplitude is strongly inconsistent with three-body oscillator descriptions of the FeCO fragment, but agrees quantitatively with DFT predictions. Over the past decade, quantum chemical calculations have suggested revised estimates of the importance of steric distortion of the bound CO in preventing poisoning of heme proteins by carbon monoxide. Quantitative agreement with the predicted frequency, amplitude, and direction of Fe motion for the FeCO bending vibrations provides direct experimental support for the quantum chemical description of the energetics of the FeCO unit.

Time-resolved visible and infrared study of the cyano complexes of myoglobin and of hemoglobin I from Lucina pectinata

The dynamics of the ferric CN complexes of the heme proteins Myoglobin and Hemoglobin I from the clam Lucina pectinata upon Soret band excitation is monitored using infrared and broad band visible pump-probe spectroscopy. The transient response in the UV-vis spectral region does not depend on the heme pocket environment and is very similar to that known for ferrous proteins. The main feature is an instantaneous, broad, short-lived absorption signal that develops into a narrower red-shifted Soret band. Significant transient absorption is also observed in the 360-390 nm range. At all probe wavelengths the signal decays to zero with a longest time constant of 3.6 ps. The infrared data on MbCN reveal a bleaching of the C triple bond N stretch vibration of the heme-bound ligand, and the formation of a five-times weaker transient absorption band, 28 cm(-1) lower in energy, within the time resolution of the experiment. The MbC triple bond N stretch vibration provides a direct measure for the return of population to the ligated electronic (and vibrational) ground state with a 3-4 ps time constant. In addition, the CN-stretch frequency is sensitive to the excitation of low frequency heme modes, and yields independent information about vibrational cooling, which occurs on the same timescale.

Viscosity-dependent Relaxation Significantly Modulates the Kinetics of CO Recombination in the Truncated Hemoglobin TrHbN from Mycobacterium tuberculosis

Kinetic traces were generated for the nanosecond and slower rebinding of photodissociated CO to trHbN in solution and in porous sol-gel matrices as a function of viscosity, conformation, and mutation. TrHbN is one of the two truncated hemoglobins from Mycobacterium tuberculosis. The kinetic traces were analyzed in terms of three distinct phases. These three phases are ascribed to rebinding: (i) from the distal heme pocket, (ii) from the adjacent apolar tunnel prior to conformational relaxation, and (iii) from the apolar tunnel subsequent to conformational relaxation. The fractional content of each of these phases was shown to be a function of the viscosity and, in the case of the sol-gel-encapsulated samples, sample preparation history. The observed kinetic patterns support a model consisting of the following elements: (i) the viscosity and conformation-sensitive dynamics of the Tyr(B10) side chain facilitate diffusion of the dissociated ligand from the distal heme pocket into the adjacent tunnel; (ii) the distal heme pocket architecture determines ligand access from the tunnel back to the heme iron; (iii) the distal heme pocket architecture is governed by a ligand-dependent hydrogen bonding network that limits the range of accessible side chain positions; and (iv) the apolar tunnel linking the heme site to the solvent biases the competition between water and ligand for occupancy of the vacated polar distal heme pocket greatly toward the nonpolar ligand. Implications of these finding with respect to biological function are discussed.

Low-lying electronic states of FeNC and FeCN: A theoretical journey into isomerization and quartet/sextet competition

With several levels of multireference and restricted open-shell single-reference electronic structure theory, optimum structures, relative energetics, and spectroscopic properties of the low-lying (6)Delta, (6)Pi, (4)Delta, (4)Pi, and (4)Sigma(-) states of linear FeNC and FeCN have been investigated using five contracted Gaussian basis sets ranging from Fe[10s8p3d], C/N[4s2p1d] to Fe[6s8p6d3f2g1h], C/N[6s5p4d3f2g]. Based on multireference configuration interaction (MRCISD+Q) results with a correlation-consistent polarized valence quadruple-zeta (cc-pVQZ) basis set, appended with core correlation and relativistic corrections, we propose the relative energies: T(e)(FeNC), (6)Delta(0)<(6)Pi (2300 cm(-1))<(4)Delta (2700 cm(-1))<(4)Pi (4200 cm(-1))<(4)Sigma(-); and T(e)(FeCN), (6)Delta(0)<(6)Pi (1800 cm(-1))<(4)Delta (2500 cm(-1))<(4)Pi (2900 cm(-1))<(4)Sigma(-). The (4)Delta and (4)Pi states have massive multireference character, arising mostly from 11sigma-->12sigma promotions, whereas the sextet states are dominated by single electronic configurations. The single-reference CCSDT-3 (coupled cluster singles and doubles with iterative partial triples) method appears to significantly overshoot the stabilization of the quartet states provided by both static and dynamical correlation. The (4,6)Delta and (4,6)Pi states of both isomers are rather ionic, and all have dipole moments near 5 D. On the ground (6)Delta surface, FeNC is predicted to lie 0.6 kcal mol(-1) below FeCN, and the classical barrier for isocyanide/cyanide isomerization is about 6.5 kcal mol(-1). Our data support the recent spectroscopic characterization by Lei and Dagdigian [J. Chem. Phys. 114, 2137 (2000)] of linear (6)Delta FeNC as the first experimentally observed transition-metal monoisocyanide. Their assignments for the ground term symbol, isotopomeric rotational constants, and the Fe-N omega(3) stretching frequency are confirmed; however, we find rather different structural parameters for (6)Delta FeNC:r(e)(Fe-N)=1.940 A and r(N-C)=1.182 A at the cc-pVQZ MRCISD+Q level. Our results also reveal that the observed band of FeNC originating at 27 236 cm(-1) should have an analog in FeCN near 23 800 cm(-1) of almost equal intensity. Therefore, both thermodynamic stability and absorption intensity factors favor the eventual observation of FeCN via a (6)Pi<--(6)Delta transition in the near-UV.

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.

Simultaneous observation of the O---O and Fe---O2 stretching modes in oxyhemoglobins

Understanding of the chemical nature of the dioxygen moiety of oxyhemoglobin is crucial for elucidation of its physiological function. In the present work, direct Raman spectroscopic observation of both the FeO(2) and OO stretching modes unambiguously establishes the vibrational characteristics of the oxygen-bound heme moiety in the hemoglobins of Chlamydomonas eugametos and Synechocystis PCC6803. In addition to providing the resonance Raman assignment of the OO stretching mode (1136 cm(-1) for Chlamydomonas, 1133 cm(-1) for Synechocystis) in an oxyhemoglobin with an iron-porphyrin, this study also reports unusually low frequencies for the FeO(2) stretching modes (554 cm(-1)). The effect of strong hydrogen bonding to the bound oxygen is confirmed by changes in the frequency of the FeO(2) stretching mode on mutation of distal residues. These findings suggest an enzymatic function rather than an oxygen transport role for these hemoglobins.