Structural factors controlling ligand binding to myoglobin: a kinetic hole-burning study

Synthesis of carbon nanotubes-chitosan nanocomposite and immunosensor fabrication for myoglobin detection: An acute myocardial infarction biomarker

The use of single-walled carbon nanotubes (SWCNTs) in biomedical applications is limited due to their inability to disperse in aqueous solutions. In this study, dispersed -COOH functionalized CNTs with N-succinylated chitosan (CS), greatly increasing the water solubility of CNTs and forming a uniformly dispersed nanocomposite solution of CNTs@CS. Coupling reagent EDC/NHS was used as a linker with the -COOH groups present on the N-succinylated chitosan which significantly improved the affinity of the CNTs for biomolecules. Myoglobin (Mb) is a promising biomarker for the precise assessment of cardiovascular risk, type 2 diabetes, metabolic syndrome, hypertension and several types of cancer. A high level of Mb can be used to diagnose the mentioned pathogenic diseases. The CNTs@CS-FET demonstrates superior sensing performance for Mb antigen fortified in buffer, with a wide linear range of 1 to 4000 ng/mL. The detection limit of the developed Mb immunosensor was estimated to be 4.2 ng/mL. The novel CNTs@CS-FET immunosensor demonstrates remarkable capability in detecting Mb without being affected by interferences from nonspecific antigens. Mb spiked serum showed a recovery rate of 100.262 to 118.55 % indicating great promise for Mb detection in clinical samples. The experimental results confirmed that the CNTs@CS-FET immunosensor had excellent selectivity, reproducibility and storage stability.

Substrate Binding Primes Human Tryptophan 2,3-Dioxygenase for Ligand Binding

The human heme enzyme tryptophan 2,3-dioxygenase (hTDO) catalyzes the insertion of dioxygen into its cognate substrate, l-tryptophan (l-Trp). Its active site structure is highly dynamic, and the mechanism of enzyme-substrate-ligand complex formation and the ensuing enzymatic reaction is not yet understood. Here we have studied complex formation in hTDO by using time-resolved optical and infrared spectroscopy with carbon monoxide (CO) as a ligand. We have observed that both substrate-free and substrate-bound hTDO coexist in two discrete conformations with greatly different ligand binding rates. In the fast rebinding hTDO conformation, there is facile ligand access to the heme iron, but it is greatly hindered in the slowly rebinding conformation. Spectroscopic evidence implicates active site solvation as playing a crucial role for the observed kinetic differences. Substrate binding shifts the conformational equilibrium markedly toward the fast species and thus primes the active site for subsequent ligand binding, ensuring that formation of the ternary complex occurs predominantly by first binding l-Trp and then the ligand. Consequently, the efficiency of catalysis is enhanced because O₂ binding prior to substrate binding, resulting in nonproductive oxidation of the heme iron, is greatly suppressed.

Hydrogen bonding plays a significant role in the binding of coomassie brilliant blue-R to hemoglobin: FT-IR, fluorescence and molecular dynamics studies

Coomassie brilliant blue-R (CBB-R) specifically binds to bovine hemoglobin with a stoichiometric ratio of 1 : 1.

A spectroscopic study of structural heterogeneity and carbon monoxide binding in neuroglobin

Neuroglobin (Ngb) is a small globular protein that binds diatomic ligands like oxygen, carbon monoxide (CO) and nitric oxide at a heme prosthetic group. We have performed FTIR spectroscopy in the infrared stretching bands of CO and flash photolysis with monitoring in the electronic heme absorption bands to investigate structural heterogeneity at the active site of Ngb and its effects on CO binding and migration at cryogenic temperatures. Four CO stretching bands were identified; they correspond to discrete conformations that differ in structural details and CO binding properties. Based on a comparison of bound-state and photoproduct IR spectra of the wild-type protein, Ngb distal pocket mutants and myoglobin, we have provided structural interpretations of the conformations associated with the different CO bands. We have also studied ligand migration to the primary docking site, B. Rebinding from this site is governed by very low enthalpy barriers (∼1 kJ/mol), indicating an extremely reactive heme iron. Moreover, we have observed ligand migration to a secondary docking site, C, from which CO rebinding involves higher enthalpy barriers.

Hole-burning spectroscopy as a probe of nano-environments and processes in biomolecules: a review

Hole-burning spectroscopy, a high-resolution spectroscopic technique, allows details of heterogeneous nano-environments in biological systems to be obtained from broad absorption bands. Recently, this technique has been applied to proteins, nucleic acids, cells, and substructures of water to probe the electrostatic conditions created by macromolecules and the surrounding solvent. Starting with the factors that obscure the homogeneous linewidth of a chromophore within an inhomogeneously broadened absorption or emission band, we describe properties and processes in biological systems that are reflected in the measured hole spectra. The technique also lends itself to the resolution of perturbation experiments, such as temperature cycling to elucidate energy landscape barriers, applied external electric fields (Stark effect) to measure net internal electric fields, and applied hydrostatic pressure to find the volume compressibility of proteins.

Ligand migration and binding in nonsymbiotic hemoglobins of Arabidopsis thaliana

We have studied carbon monoxide (CO) migration and binding in the nonsymbiotic hemoglobins AHb1 and AHb2 of Arabidopsis thaliana using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures. Both proteins have similar amino acid sequences but display pronounced differences in ligand binding properties, at both physiological and cryogenic temperatures. Near neutral pH, the distal HisE7 side chain is close to the heme-bound ligand in the majority of AHb1-CO molecules, as indicated by a low CO stretching frequency at 1921 cm(-1). In this fraction, two CO docking sites can be populated, the primary site B and the secondary site C. When the pH is lowered, a high-frequency stretching band at approximately 1964 cm(-1) grows at the expense of the low-frequency band, indicating that HisE7 protonates and, concomitantly, moves away from the bound ligand. Geminate rebinding barriers are markedly different for the two conformations, and docking site C is not accessible in the low-pH conformation. Rebinding of NO ligands was observed only from site B of AHb1, regardless of conformation. In AHb2, the HisE7 side chain is removed from the bound ligand; rebinding barriers are low, and CO molecules can populate only primary docking site B. These results are interpreted in terms of differences in the active site structures and physiological functions.

Ligand dynamics in heme proteins observed by Fourier transform infrared-temperature derivative spectroscopy

Fourier transform infrared (FTIR) spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. Nitric oxide and carbon monoxide are attractive physiologically relevant ligands because their bond stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands display changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300K) using specific temperature protocols for sample photodissociation can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy (TDS) has proven to be a particularly powerful technique to study protein-ligand interactions. The FTIR-TDS technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. Here we describe infrared cryo-spectroscopy and present a variety of applications to the study of protein-ligand interactions in heme proteins. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

Influence of distal residue B10 on CO dynamics in myoglobin and neuroglobin

For many years, myoglobin has served as a paradigm for structure-function studies in proteins. Ligand binding and migration within myoglobin has been studied in great detail by crystallography and spectroscopy, showing that gaseous ligands such as O(2), CO, and NO not only bind to the heme iron but may also reside transiently in three internal ligand docking sites, the primary docking site B and secondary sites C and D. These sites affect ligand association and dissociation in specific ways. Neuroglobin is another vertebrate heme protein that also binds small ligands. Ligand migration pathways in neuroglobin have not yet been elucidated. Here, we have used Fourier transform infrared temperature derivative spectroscopy at cryogenic temperatures to compare the influence of the side chain volume of amino acid residue B10 on ligand migration to and rebinding from docking sites in myoglobin and neuroglobin.

Ligand dynamics in heme proteins observed by Fourier transform infrared spectroscopy at cryogenic temperatures

Fourier transform infrared spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. From the variety of ligands that bind to the heme iron, nitric oxide and carbon monoxide are particularly attractive, as their bond-stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, the ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand-stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands reveal changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites, and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300 K) using specific temperature protocols for sample photodissociation thus can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy has proven to be a particularly powerful technique to study protein-ligand interactions. This technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. This chapter describes infrared cryospectroscopy techniques and presents examples that demonstrate their applicability to nitric oxide binding to heme proteins.

Asymmetric band profile of the Soret band of deoxymyoglobin is caused by electronic and vibronic perturbations of the heme group rather than by a doming deformation

We measured the Soret band of deoxymyoglobin (deoxyMb), myoglobin cyanide (MbCN), and aquo-metmyoglobin (all from horse heart) with absorption and circular dichroism (CD) spectroscopies. A clear non-coincidence was observed between the absorption and CD profiles of deoxyMb and MbCN, with the CD profiles red- and blueshifted with respect to the absorption band position, respectively. On the contrary, the CD and absorption profiles of aquametMb were nearly identical. The observed noncoincidence indicates a splitting of the excited B state due to heme-protein interactions. CD and absorption profiles of deoxyMb and MbCN were self-consistently analyzed by employing a perturbation approach for weak vibronic coupling as well as the relative intensities and depolarization ratios of seven bands in the respective resonance Raman spectra measured with B-band excitation. The respective B(y) component was found to dominate the observed Cotton effect of both myoglobin derivatives. The different signs of the noncoincidences between CD and absorption bands observed for deoxyMb and MbCN are due to different signs of the respective matrix elements of A(1g) electronic interstate coupling, which reflects an imbalance of Gouterman's 50:50 states. The splitting of the B band reflects contributions from electronic and vibronic perturbations of B(1g) symmetry. The results of our analysis suggest that the broad and asymmetric absorption band of deoxyMb results from this band splitting rather than from its dependence on heme doming. Thus, we are able to explain recent findings that the temperature dependences of CO rebinding to myoglobin and the Soret band profile are uncorrelated[Ormos et al., Proc. Natl. Acad. Sci U.S.A. 95, 6762 (1998)].

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.

Spectroscopic study of substrate binding to the carbonmonoxy form of dehaloperoxidase from Amphitrite ornata

Dehaloperoxidase (DHP) is a globular heme enzyme found in the marine worm Amphitrite ornata that can catalyze the dehalogenation of halophenols to the corresponding quinones by using hydrogen peroxide as a cosubstrate. Its three-dimensional fold is surprisingly similar to that of the oxygen storage protein myoglobin (Mb). A key structural feature common to both DHP and Mb is the existence of multiple conformations of the distal histidine. In DHP, the conformational flexibility may be involved in promotion of substrate and cosubstrate entry and exit. Here we have explored the dynamics of substrate binding in DHP using Fourier transform infrared spectroscopy and flash photolysis. A number of discrete conformations at the active site were identified from the appearance of multiple CO absorbance bands in the infrared region of the spectrum. Upon photolysis at cryogenic temperatures, the CO molecules are trapped at docking sites within the protein matrix, as inferred from the appearance of several photoproduct bands characteristic of each site. Substrate binding stabilizes the protein by approximately 20 kJ/mol. The low yield of substrate-bound DHP at ambient temperature points toward a steric inhibition of substrate binding by carbon monoxide.

Spectroscopic markers of the T<-->R quaternary transition in human hemoglobin

In this work, we use a sol-gel protocol to trap and compare the R and T quaternary states of both the deoxygenated (deoxyHb) and carbonmonoxide (HbCO) derivatives of human hemoglobin. The near infrared optical absorption band III and the infrared CO stretching band are used to detect the effect of quaternary structure on the spectral properties of deoxyHb and HbCO; comparison with myoglobin allows for an assessment of tertiary and quaternary contributions to the measured band shifts. The R<-->T transition is shown to cause a blue shift of the band III by approximately 35 cm(-1) for deoxyHb and a red shift of the CO stretching band by only approximately 0.3 cm(-1) for HbCO. This clearly shows that quaternary structure changes are transmitted to the heme pocket and that effects on deoxyHb are much larger than on HbCO, at least as far as the band energies are concerned. Experiments performed in the ample temperature interval of 300-10K show that the above quaternary structure effects are "static" and do not influence the dynamic properties of the heme pocket, at least as probed by the temperature dependence of band III and of the CO stretching band. The availability of quaternary structure sensitive spectroscopic markers and the quantitative measurement of the quaternary structure contribution to band shifts will be of considerable help in the analysis of flash-photolysis experiments on hemoglobin. Moreover, it will enable one to characterize the dynamic properties of functionally relevant hemoglobin intermediates and to study the kinetics of both the T-->R and R-->T quaternary transitions through time-resolved spectroscopy.

Probing electric fields in protein cavities by using the vibrational stark effect of carbon monoxide

To determine the magnitude and direction of the internal electric field in the Xe4 cavity of myoglobin mutant L29W-S108L, we have studied the vibrational Stark effect of carbon monoxide (CO) using infrared spectroscopy at cryogenic temperatures. CO was photodissociated from the heme iron and deposited selectively in Xe4. Its infrared spectrum exhibits Stark splitting into two bands associated with CO in opposite orientations. Two different photoproduct states can be distinguished, C' and C'', with markedly different properties. For C', characteristic temperature-dependent changes of the area, shift, and width were analyzed, based on a dynamic model in which the CO performs fast librations within a double-well model potential. For the barrier between the wells, a height of approximately 1.8 kJ/mol was obtained, in which the CO performs oscillations at an angular frequency of approximately 25 cm(-1). The magnitude of the electric field in the C' conformation was determined as 11.1 MV/cm; it is tilted by an angle of 29 degrees to the symmetry axis of the potential. Above 140 K, a protein relaxation leads to a significantly altered photoproduct, C'', with a smaller Stark splitting and a more confining potential (barrier >4 kJ/mol) governing the CO librations.

Different relaxations in myoglobin after photolysis

To clarify the interplay of kinetic hole-burning (KHB), structural relaxation, and ligand migration in myoglobin (Mb), we measured time-resolved absorption spectra in the Soret region after photolysis of carbon monoxide Mb (MbCO) in the temperature interval 120-260 K and in the time window 350 ns to 200 ms. The spectral contributions of both photolyzed (Mb*) and liganded Mb (MbCO) have been analyzed by taking into account homogeneous bandwidth, coupling to vibrational modes, and static conformational heterogeneity. We succeeded in separating the "time-dependent" spectral changes, and this work provides possibilities to identify the events in the process of ligand rebinding. KHB is dominant at T <190 K in both the Mb* and the MbCO components. For MbCO, conformational substates interconversion at higher temperatures tends to average out the KHB effect. At 230-260 K, whereas almost no shift is observed in the MbCO spectrum, a shift of the order of approximately 80 cm(-1) is observed in Mb*. We attribute this shift to protein relaxation coupled to ligand migration. The time dependence of the Mb* spectral shift is interpreted with a model that enables us to calculate the highly nonexponential relaxation kinetics. Fits of stretched exponentials to this kinetics yield Kohlrausch parameter values of 0.25, confirming the analogy between proteins and glasses.

The structure of murine neuroglobin: Novel pathways for ligand migration and binding

Neuroglobin, a recently discovered globin predominantly expressed in neuronal tissue of vertebrates, binds small, gaseous ligands at the sixth coordination position of the heme iron. In the absence of an exogenous ligand, the distal histidine (His64) binds to the heme iron in the ferrous and ferric states. The crystal structure of murine ferric (met) neuroglobin at 1.5 A reveals interesting features relevant to the ligand binding mechanism. Only weak selectivity is observed for the two possible heme orientations, the occupancy ratio being 70:30. Two small internal cavities are present on the heme distal side, which enable the His64(E7) side chain to move out of the way upon exogenous ligand binding. Moreover, a third, huge cavity (volume approximately 290 A3) connecting both sides of the heme, is open towards the exterior and provides a potential passageway for ligands. The CD and EF corners exhibit substantial flexibility, which may assist ligands in entering the protein and accessing the active site. Based on this high-resolution structure, further structure-function studies can be planned to elucidate the role of neuroglobin in physiological responses to hypoxia.

Ligand dynamics in a protein internal cavity

We have studied the temperature dependence of the IR stretch bands of carbon monoxide (CO) in the Xe 4 internal cavity of myoglobin mutant L29W-S108L at cryogenic temperatures. Pronounced changes of band areas and positions were analyzed quantitatively by using a simple dynamic model in which CO rotation in the cavity is constrained by a static potential. The librational dynamics of the CO causes a decrease of the total spectral area. A strong local electric field splits the CO stretch absorption into a doublet, indicating that CO can assume opposite orientations in the cavity. With increasing temperature, the two peaks approach each other, because the average angle of the CO with respect to the electric field increases. A combined classical and quantum-mechanical analysis precisely reproduces the observed temperature dependencies of both spectral area and peak shifts. It yields the height of the energy barrier between the two wells associated with opposite CO orientations, V0 approximately 2 kJ/mol, and the frequency of oscillation within a well, omega approximately 25 cm(-1). The electric field in the protein cavity was estimated as 10 MV/cm.

Heme charge-transfer band III is vibronically coupled to the Soret band

A complete resonance Raman excitation profile of the heme charge-transfer band known as band III is presented. The data obtained throughout the near-infrared region show preresonance with the Q-band, but the data also clearly show the enhancement of a number of modes in the spectral region of band III. Only nontotally symmetric modes are observed to have resonance enhancement in the band III region. The observed resonance enhancements in modes of B(1g) symmetry are compared with the enhancements of those same modes in the excitation profiles of the Q-band of deoxy myoglobin, also presented here for this first time. The Q-band data agree well with the theory of vibronic coupling in metalloporphyrins (Shelnutt, J. A. J. Chem. Phys. 1981, 74, 6644-6657). The strong vibronic coupling of the Q-band of the deoxy form of hemes is discussed in terms of the enhancement of modes with both B(1g) and A(2g) symmetry. The comparison between the Q-band and band III reveals that, consistent with the theory, only modes of B(1g) symmetry are enhanced in the vicinity of band III. These results show that band III is vibronically coupled to the Soret band. The coupling of band III to modes with strong rhombic distortion of the heme macrocycle calls into question the hypothesis that the axial iron out-of-plane displacement is primarily responsible for the structure-dynamics correlations observed in myoglobin.

Ligand binding and protein dynamics in neuroglobin

Neuroglobin (Ngb) is a recently discovered protein in vertebrate brain tissue that belongs to the globin family of proteins. It has been implicated in the neuronal response to hypoxia or ischemia, although its physiological role has been hitherto unknown. Ngb is hexacoordinate in the ferrous deoxy form under physiological conditions. To bind exogenous ligands like O(2) and CO, the His E7 endogenous ligand is displaced from the sixth coordination. By using infrared spectroscopy and nanosecond time-resolved visible spectroscopy, we have investigated the ligand-binding reaction over a wide temperature range (3-353 K). Multiple, intrinsically heterogeneous distal heme pocket conformations exist in NgbCO. Photolysis at cryogenic temperatures creates a five-coordinate deoxy species with very low geminate-rebinding barriers. The photodissociated CO is observed to migrate within the distal heme pocket even at 20 K. Flash photolysis near physiological temperature (275-353 K) exhibits four sequential kinetic features: (i) geminate rebinding (t < 1 micros); (ii) extremely fast bimolecular exogenous ligand binding (10 micros < t < 1 ms) with a nontrivial temperature dependence; (iii) endogenous ligand binding (100 micros < t < 10 ms), which can be studied by using flash photolysis on deoxy Ngb; and (iv) displacement of the endogenous by the exogenous ligand (10 ms < t < 10 ks). All four processes are markedly nonexponential, suggesting that Ngb fluctuates among different conformations on surprisingly long time scales.

Protein dynamics in an intermediate state of myoglobin: optical absorption, resonance Raman spectroscopy, and x-ray structure analysis

A metastable state of myoglobin is produced by reduction of metmyoglobin at low temperatures. This is done either by irradiation with x-rays at 80 K or by electron transfer from photoexcited tris(2, 2'-bipyridine)-ruthenium(II) at 20 K. At temperatures above 150 K, the conformational transition toward the equilibrium deoxymyoglobin is observed. X-ray crystallography, Raman spectroscopy, and temperature-dependent optical absorption spectroscopy show that the metastable state has a six-ligated iron low-spin center. The x-ray structure at 115K proves the similarity of the metastable state with metmyoglobin. The Raman spectra yield the high-frequency vibronic modes and give additional information about the distortion of the heme. Analysis of the temperature dependence of the line shape of the Soret band reveals that a relaxation within the metastable state starts at approximately 120 K. Parameters representative of static properties of the intermediate state are close to those of CO-ligated myoglobin, while parameters representative of dynamics are close to deoxymyoglobin. Thus within the metastable state the relaxation to the equilibrium is initiated by changes in the dynamic properties of the active site.

Harmonic and anharmonic dynamics of Fe-CO and Fe-O(2) in heme models

We present density-functional molecular dynamics simulations of FeP(Im)(AB) heme models (AB = CO, O(2), Im = imidazole) as a way of sketching the dynamic motion of the axial ligands at room temperature. The FeP(Im)(CO) model is characterized by an essentially upright FeCO unit, undergoing small deviations with respect to its linear equilibrium structure (bending and tilting up to 10 degrees and 7 degrees, often occur). The motion of the carbon monoxide ligand is found to be quite complex and fast, its projection on the porphyrin plane sampling all the porphyrin quadrants in a short time ( approximately 0.5 ps). Simultaneously, the imidazole ligand rotates slowly around the Fe-N(epsilon) bond. In contrast to carbon monoxide, the oxygen ligand in FeP(Im)(O(2)) prefers a conformation where the projection of the O-O axis on the porphyrin plane bisects one of the porphyrin quadrants. A transition to other quadrants takes place through an O-O/Fe-N(p) overlapping conformation, within 4-6 ps. Further details of these mechanisms and their implications are discussed.