Ligand binding to heme proteins: III. FTIR studies of His-E7 and Val-E11 mutants of carbonmonoxymyoglobin

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

Infrared spectroscopy from electrostatic embedding QM/MM: local normal mode analysis of infrared spectra of arabidopsis thaliana plant cryptochrome

Infrared (IR) spectroscopy is an undoubtedly valuable tool for analyzing vibrations, conformational changes, and chemical reactions of biological macromolecules. Currently, there is a lack of theoretical methods to create a model successfully and efficiently simulate and interpret the origin of the spectral signatures, which are often complex to analyze. Here, we develop a new method for IR vibrational spectroscopy based on analytic second derivatives of electrostatic embedding QM/MM energy, the computation of electric dipole moments with respect to nuclear perturbations and the localization of normal modes. In addition to the IR spectrum, the method can provide the origin of each peak from clearly identified molecular motions of constituent fragments. As a proof of concept, we analyze the IR spectra of flavin adenine dinucleotides in water and in Arabidopsis thaliana cryptochrome proteins for four redox forms, in addition to the difference IR spectra before and after illumination with blue light. We show that the main peaks in the difference spectrum are due to N-H hydrogen out-of-plane motions and hydrogen bendings.

Different Mechanisms of Catalytic Complex Formation in Two L-Tryptophan Processing Dioxygenases

The human heme enzymes tryptophan 2,3-dioxygenase (hTDO) and indoleamine 2,3 dioxygenase (hIDO) catalyze the initial step in L-tryptophan (L-Trp) catabolism, the insertion of dioxygen into L-Trp. Overexpression of these enzymes causes depletion of L-Trp and accumulation of metabolic products, and thereby contributes to tumor immune tolerance and immune dysregulation in a variety of disease pathologies. Understanding the assembly of the catalytically active, ternary enzyme-substrate-ligand complexes is not yet fully resolved, but an essential prerequisite for designing efficient and selective de novo inhibitors. Evidence is mounting that the ternary complex forms by sequential binding of ligand and substrate in a specific order. In hTDO, the apolar L-Trp binds first, decreasing active-site solvation and, as a result, reducing non-productive oxidation of the heme iron by the dioxygen ligand, which may leave the substrate bound to a ferric heme iron. In hIDO, by contrast, dioxygen must first coordinate to the heme iron because a bound substrate would occlude ligand access to the heme iron, so the ternary complex can no longer form. Consequently, faster association of L-Trp at high concentrations results in substrate inhibition. Here, we summarize our present knowledge of ternary complex formation in hTDO and hIDO and relate these findings to structural peculiarities of their active sites.

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.

Rebinding dynamics of NO to microperoxidase-8 probed by time-resolved vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding kinetics of NO to microperoxidase-8 (Mp), an ideal model system for the active site of ligand-binding heme proteins, including myoglobin and hemoglobin, after the photodeligation of MpNO in glycerol/water (G/W) solutions at 294 K. The geminate rebinding (GR) of NO to Mp in viscous solutions was highly efficient and ultrafast and negligibly dependent on the solution viscosity, which was adjusted by changing the glycerol content from 65% to 90% by volume in G/W mixtures. The kinetics of the GR of NO to Mp in viscous solutions was well represented by an exponential function with a time constant of ca. 11 ps. Although the kinetic traces of the GR of NO to Mp in solutions with three different viscosities (18, 81, and 252 cP) almost overlap, they show a slight difference early in the decay process. The kinetic traces were also described by the diffusion-controlled reaction theory with a Coulomb potential. Since the ligand is deligated in a neutral form, an ionic pair of NO(-) and Mp(+) may be produced before forming the Mp-NO bond by an electron transfer from Mp to NO as the deligated NO is sufficiently near to the Fe atom of Mp. The strong reactivity between NO and ferrous heme may arise from the Coulomb interaction between the reacting pair, which is consistent with the harpooning mechanism for NO binding to heme.

An engineered heme-copper center in myoglobin: CO migration and binding

We have investigated CO migration and binding in CuBMb, a copper-binding myoglobin double mutant (L29H-F43H), by using Fourier transform infrared spectroscopy and flash photolysis over a wide temperature range. This mutant was originally engineered with the aim to mimic the catalytic site of heme-copper oxidases. Comparison of the wild-type protein Mb and CuBMb shows that the copper ion in the distal pocket gives rise to significant effects on ligand binding to the heme iron. In Mb and copper-free CuBMb, primary and secondary ligand docking sites are accessible upon photodissociation. In copper-bound CuBMb, ligands do not migrate to secondary docking sites but rather coordinate to the copper ion. Ligands entering the heme pocket from the outside normally would not be captured efficiently by the tight distal pocket housing the two additional large imidazole rings. Binding at the Cu ion, however, ensures efficient trapping in CuBMb. The Cu ion also restricts the motions of the His64 side chain, which is the entry/exit door for ligand movement into the active site, and this restriction results in enhanced geminate and slow bimolecular CO rebinding. These results support current mechanistic views of ligand binding in hemoglobins and the role of the CuB in the active of heme-copper oxidases. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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.

Ligand binding to heme proteins: a comparison of cytochrome c variants with globins

We have studied the binding of carbon monoxide (CO) in mutants of Cyt c having its methionine at position 80 replaced by alanine, aspartate, and arginine, so that the sixth coordination is available for ligand binding. We have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy to examine interactions of the heme-bound and photolyzed CO (and also nitric oxide, NO) in the small heme pocket created by the mutations. By using FTIR temperature derivative spectroscopy (TDS) and nanosecond flash photolysis, the enthalpy barrier distributions for CO rebinding were determined. In flash photolysis experiments, the majority of ligands rebind to the heme iron on picosecond time scales so that only the high-barrier tail of the distributions is visible on the nanosecond scale. By continuous wave excitation prior to TDS characterization of the barriers, however, each Cyt c molecule is photoexcited multiple times and complete photodissociation can be achieved, which likely arises from a rotation of the CO within the heme pocket so that the oxygen faces the heme iron. Apparently, reorientation prior to rebinding constitutes an additional and significant contribution to the rebinding barrier. Our experiments reveal that the compact, rigid structure of Cyt c offers no alternative binding sites for photodissociated ligands in the protein matrix. A comparison of ligand binding in these Cyt c mutants and hemoglobins underscores the importance of internal ligand docking sites and ligand migration routes for conveying a ligand binding function to heme proteins.

NMR, IR, Mössbauer and quantum chemical investigations of metalloporphyrins and metalloproteins

We review contributions made towards the elucidation of CO and O2binding geometries in respiratory proteins. Nuclear magnetic resonance, infrared spectroscopy, Mössbauer spectroscopy, X-ray crystallography and quantum chemistry have all been used to investigate the Fe –ligand interactions. Early experimental results showed linear correlations between17O chemical shifts and the infrared stretching frequency (νCO) of the CO ligand in carbonmonoxyheme proteins and between the17O chemical shift and the13CO shift. These correlations led to early theoretical investigations of the vibrational frequency of carbon monoxide and of the13C and17O NMR chemical shifts in the presence of uniform and non-uniform electric fields. Early success in modeling these spectroscopic observables then led to the use of computational methods, in conjunction with experiment, to evaluate ligand-binding geometries in heme proteins. Density functional theory results are described which predict57Fe chemical shifts and Mössbauer electric field gradient tensors,17O NMR isotropic chemical shifts, chemical shift tensors and nuclear quadrupole coupling constants (e2qQ/h) as well as13C isotropic chemical shifts and chemical shift tensors in organometallic clusters, heme model metalloporphyrins and in metalloproteins. A principal result is that CO in most heme proteins has an essentially linear and untilted geometry (τ = 4 °, β = 7 °) which is in extremely good agreement with a recently published X-ray synchrotron structure. CO / O2discrimination is thus attributable to polar interactions with the distal histidine residue, rather than major Fe–C–O geometric distortions.

Local structure and dynamics of hemeproteins by X-ray absorption near edge structure spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a synchrotron radiation technique sensitive to the local structure and dynamics around the metal site of a heme containing protein. Advances in detection techniques and theoretical/computational platforms in the last 15 years allowed the use of XANES as a quantitative probe of the key structural determinants driving functional changes, both in a concerted way with protein crystallography and EXAFS (extended X-ray absorption fine structure), or as a stand-alone method to apply in the crystal state as well as in solution. Moreover, the local dynamics of the heme site has been deeply investigated, on one hand, coupling XANES to classical photolysis experiments at cryogenic temperatures; on the other hand, the intrinsic property of the synchrotron radiation to induce radiolysis events, has been exploited to investigate specific cryotrapped intermediates, using X-rays both as a pump and a probe. Insights on the XANES method and some specific examples are presented to illustrate these topics.

Ligand migration in human indoleamine-2,3 dioxygenase

Human indoleamine 2,3-dioxygenase (hIDO), a monomeric heme enzyme, catalyzes the oxidative degradation of L-tryptophan (L-Trp) and other indoleamine derivatives. Its activity follows typical Michaelis-Menten behavior only for L-Trp concentrations up to 50 μM; a further increase in the concentration of L-Trp causes a decrease in the activity. This substrate inhibition of hIDO is a result of the binding of a second L-Trp molecule in an inhibitory substrate binding site of the enzyme. The molecular details of the reaction and the inhibition are not yet known. In the following, we summarize the present knowledge about this heme enzyme.

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.

CO binding and ligand discrimination in human myeloperoxidase

Despite the fact that ferrous myeloperoxidase (MPO) can bind both O(2) and NO, its ability to bind CO has been questioned. UV/visible spectroscopy was used to confirm that CO induces small spectral shifts in ferrous MPO, and Fourier transform infrared difference spectroscopy showed definitively that these arose from formation of a heme ferrous-CO compound. Recombination rates after CO photolysis were monitored at 618 and 645 nm as a function of CO concentration and pH. At pH 6.3, k(on) and k(off) were 0.14 mM(-1) x s(-1) and 0.23 s(-1), respectively, yielding an unusually high K(D) of 1.6 mM. This affinity of MPO for CO is 10 times weaker than its affinity for O(2). The observed rate constant for CO binding increased with increasing pH and was governed by a single protonatable group with a pK(a) of 7.8. Fourier transform infrared spectroscopy revealed two different conformations of bound CO with frequencies at 1927 and 1942 cm(-1). Their recombination rate constants were identical, indicative of two forms of bound CO that are in rapid thermal equilibrium rather than two distinct protein populations with different binding sites. The ratio of bound states was pH-dependent (pK(a) approximately 7.4) with the 1927 cm(-1) form favored at high pH. Structural factors that account for the ligand-binding properties of MPO are identified by comparisons with published data on a range of other ligand-binding heme proteins, and support is given to the recent suggestion that the proximal His336 in MPO is in a true imidazolate state.

Ligand and substrate migration in human indoleamine 2,3-dioxygenase

Human indoleamine 2,3-dioxygenase (hIDO), a monomeric heme enzyme, catalyzes the oxidative degradation of L-Trp and other indoleamine derivatives. Using Fourier transform infrared and optical absorption spectroscopy, we have investigated the interplay between ferrous hIDO, the ligand analog CO, and the physiological substrate L-Trp. These data provide the long sought evidence for two distinct L-Trp binding sites. Upon photodissociation from the heme iron at T > 200 K, CO escapes into the solvent. Concomitantly, L-Trp exits the active site and, depending on the l-Trp concentration, migrates to a secondary binding site or into the solvent. Although L-Trp is spectroscopically silent at this site, it is still noticeable due to its pronounced effect on the CO association kinetics, which are significantly slower than those of L-Trp-free hIDO. L-Trp returns to its initial site only after CO has rebound to the heme iron.

Structural assignment of spectra by characterization of conformational substates in bound MbCO

Residue motions of the distal heme pocket and bound CO ligand of carbonmonoxy Myoglobin are studied using a combination of molecular dynamics simulations and quantum chemical methods. Using mixed quantum mechanics/molecular mechanics calculations together with sampling from molecular dynamics simulations (QM/MM(MD)), the experimentally observed spectroscopic A(0) and A(1) substates of the bound CO ligand are assigned to the open and closed conformation of His(64) and the His(epsilon)(64) tautomer, respectively. Several previously proposed origins of the A(3) substate, including rotamers of the doubly protonated His(64)H(+) side chain, His(64)H(+) inside the distal pocket, and cooperative motions with Arg(45), are investigated with QM/MM(MD). However, the signatures of the calculated infrared spectra do not agree with the experimentally observed ones. For additional insight on this, extensive molecular dynamics simulations are used together with improved electrostatics for the bound ligand. A CO fluctuating charge model is developed to describe the ab initio dipole and quadrupole moments of the bound ligand. CO absorption spectra are then obtained directly from the dynamics simulations. Finally, the electrostatics of the heme pocket is examined in detail in an attempt to determine the structural origins of the observed spectroscopic A-states from MD simulations. However, contrary to related simulations for unbound CO in myoglobin, the shifts and splittings for carbonmonoxy Myoglobin are generally small and difficult to relate to structural change. This suggests that coupling of the CO motion to other degrees of freedom, such as the Fe-CO stretching and bending, is important to correctly describe the dynamics of bound CO in myoglobin.

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.

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.

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Molecular Dynamics Simulations of CN− Dynamics and Spectroscopy in Myoglobin

The vibrational dynamics of the cyanide anion and the heme group in MbCN (CN complexed to Myoglobin) are investigated using molecular dynamics simulations. A previously calculated quantum-chemical heme-ligand potential-energy surface together with a three-center charge model for the iron-ligand center that captures both polarization and ligand-to-metal charge transfer allows for a detailed description of the interactions around the active site. It is found that the CN binding orientation (Fe--CN or Fe--NC) to the heme affects the stretching frequency of the ligand, with a 25 cm-1 difference in the fundamental wavenumber between the two orientations as well as a change in bond length. The charge model also captures such crucial interactions as the possible hydrogen bond between the ligand and the His64 residue. This interaction is weakened when the ligand binds in the Fe--NC conformation but is also sensitive to the protonation state of His64. The structural changes around the active site, the observation of water penetration for the Fe--NC conformation, the computed IR spectrum, and the energetics suggest that the Fe--CN conformation with Hisepsilon64 is the most likely one. The water accessibility of the active site is also found to be related to the protonation state of His64. The presence of water in the active site could also affect the IR band of the C--N stretch mode. Thus, IR spectroscopy of the C--N stretch is a potentially useful reporter of ligand isomers and active-site structure.

Transient ligand docking sites in Cerebratulus lacteus mini-hemoglobin

The monomeric hemoglobin of the nemertean worm Cerebratulus lacteus functions as an oxygen storage protein to maintain neural activity under hypoxic conditions. It shares a large, apolar matrix tunnel with other small hemoglobins, which has been implicated as a potential ligand migration pathway. Here we explore ligand migration and binding within the distal heme pocket, to which the tunnel provides access to ligands from the outside. FTIR/TDS experiments performed at cryogenic temperatures reveal the presence of three transient ligand docking sites within the distal pocket, the primary docking site B on top of pyrrole C and secondary sites C and D. Site C is assigned to a cavity adjacent to the distal portion of the heme pocket, surrounded by the B and E helices. It has an opening to the apolar tunnel and is expected to be on the pathway for ligand entry and exit, whereas site D, circumscribed by TyrB10, GlnE7, and the CD corner, most likely is located on a side pathway of ligand migration. Flash photolysis experiments at ambient temperatures indicate that the rate-limiting step for ligand binding to CerHb is migration through the apolar channel to site C. Movement from C to B and iron-ligand bond formation involve low energy barriers and thus are very rapid processes in the wt protein.

Probing dynamics of complex molecular systems with ultrafast 2D IR vibrational echo spectroscopy

Ultrafast 2D IR vibrational echo spectroscopy is described and a number of experimental examples are given. Details of the experimental method including the pulse sequence, heterodyne detection, and determination of the absorptive component of the 2D spectrum are outlined. As an initial example, the 2D spectrum of the stretching mode of CO bound to the protein myoglobin (MbCO) is presented. The time dependence of the 2D spectrum of MbCO, which is caused by protein structural evolution, is presented and its relationship to the frequency-frequency correlation function is described and used to make protein structural assignments based on comparisons to molecular dynamics simulations. The 2D vibrational echo experiments on the protein horseradish peroxidase are presented. The time dependence of the 2D spectra of the enzyme in the free form and with a substrate bound at the active site are compared and used to examine the influence of substrate binding on the protein's structural dynamics. The application of 2D vibrational echo spectroscopy to the study of chemical exchange under thermal equilibrium conditions is described. 2D vibrational echo chemical exchange spectroscopy is applied to the study of formation and dissociation of organic solute-solvent complexes and to the isomerization around a carbon-carbon single bond of an ethane derivative.

Energetics and Dynamics in MbCN: CN--Vibrational Relaxation from Molecular Dynamics Simulations

The dynamics of the cyanide anion bound to sperm-whale myoglobin is investigated using atomistic simulations. With density-functional theory, a 2D potential energy surface for the cyanide-heme complex is calculated. Two deep minima with a stabilization energy of approximately 50 kcal/mol corresponding to two different binding orientations (Fe-CN and Fe-NC) of the ligand are found. The Fe-CN conformation is favored over Fe-NC by several kcal/mol. Mixed quantum mechanics/molecular mechanics calculations show that the binding orientation affects the bond strength of the ligand, with a significantly different bond length and a 25 cm-1 shift in the fundamental CN-frequency. For the molecular dynamics (MD) simulations, a 3-center fluctuating charge model for the Fe-CN unit is developed that captures polarization and ligand-metal charge transfer. Stability arguments based on the energetics around the active site and the CN- frequency shifts suggest that the Fe-CN conformation with epsilon-protonation of His epsilon 64 are most likely, which is in agreement with experiment. Both equilibrium and nonequilibrium MD simulations are carried out to investigate the relaxation time scale and possible relaxation pathways in bound MbCN. The nonequilibrium MD simulations with a vibrationally excited ligand reveal that vibrational relaxation takes place on a time scale of hundreds of picoseconds within the active site. This finding supports the hypothesis that the experimentally observed relaxation rate (3.6 ps) reflects the repopulation of the electronic ground state.

The structure of carbonmonoxy neuroglobin reveals a heme-sliding mechanism for control of ligand affinity

Neuroglobin (Ngb), a globular heme protein expressed in the brain of vertebrates, binds oxygen reversibly, with an affinity comparable to myoglobin (Mb). Despite low sequence identity, the overall 3D fold of Ngb and Mb is very similar. Unlike in Mb, in Ngb the sixth coordination position of the heme iron is occupied by the distal histidine, in the absence of an exogenous ligand. Endogenous ligation has been proposed as a unique mechanism for affinity regulation and ligand discrimination in heme proteins. This peculiarity might be related to the still-unknown physiological function of Ngb. Here, we present the x-ray structure of CO-bound ferrous murine Ngb at 1.7 A and a comparison with the 1.5-A structure of ferric bis-histidine Ngb. We have also used Fourier transform IR spectroscopy of WT and mutant CO-ligated Ngb to examine structural heterogeneity in the active site. Upon CO binding, the distal histidine retains (by and large) its position, whereas the heme group slides deeper into a preformed crevice, thereby reshaping the large cavity ( approximately 290 A(3)) connecting the distal and proximal heme sides with the bulk. The heme relocation is accompanied by a significant decrease of structural disorder, especially of the EF loop, which may be the signal whereby Ngb communicates hypoxic conditions. This unexpected structural change unveils a heme-sliding mechanism of affinity control that may be of significance to understanding Ngb's role in the pathophysiology of the brain.

Structural Dynamics Controls Nitric Oxide Affinity in Nitrophorin 4

Nitrophorin 4 (NP4) is one of seven nitric oxide (NO) transporting proteins in the blood-sucking insect Rhodnius prolixus. In its physiological function, NO binds to a ferric iron centered in a highly ruffled heme plane. Carbon monoxide (CO) also binds after reduction of the heme iron. Here we have used Fourier transform infrared spectroscopy at cryogenic temperatures to study CO and NO binding and migration in NP4, complemented by x-ray cryo-crystallography on xenon-containing NP4 crystals to identify cavities that may serve as ligand docking sites. Multiple infrared stretching bands of the heme-bound ligands indicate different active site conformations with varying degrees of hydrophobicity. Narrow infrared stretching bands are observed for photodissociated CO and NO; temperature-derivative spectroscopy shows that these bands are associated with ligand docking sites close to the extremely reactive heme iron. No rebinding from distinct secondary sites was detected, although two xenon binding cavities were observed in the x-ray structure. Photolysis studies at approximately 200 K show efficient NO photoproduct formation in the more hydrophilic, open NP4 conformation. This result suggests that ligand escape is facilitated in this conformation, and blockage of the active site by water hinders immediate reassociation of NO to the ferric iron. In the closed, low-pH conformation, ligand escape from the active site of NP4 is prevented by an extremely reactive heme iron and the absence of secondary ligand docking sites.

FTIR study of ATP-induced changes in Na+/K+-ATPase from duck supraorbital glands

The Na+/K+-ATPase uses energy from the hydrolysis of ATP to pump Na+ ions out of and K+ ions into the cell. ATP-induced conformational changes in the protein have been examined in the Na+/K+-ATPase isolated from duck supraorbital salt glands using Fourier transform infrared spectroscopy. Both standard transmission and attenuated total internal reflection sample geometries have been employed. Under transmission conditions, enzyme at 75 mg/ml was incubated with dimethoxybenzoin-caged ATP. ATP was released by flashing with a UV laser pulse at 355 nm, which resulted in a large change in the amide I band. The absorbance at 1659 cm(-1) decreased with a concomitant increase in the absorbance at 1620 cm(-1). These changes are consistent with a partial conversion of protein secondary structure from alpha-helix to beta-sheet. The changes were approximately 8% of the total absorbance, much larger than those seen with other P-type ATPases. Using attenuated total internal reflection Fourier transform infrared spectroscopy, the decrease in absorbance at approximately 1650 cm(-1) was titrated with ATP, and the titration midpoint K0.5 was determined under different ionic conditions. In the presence of metal ions (Na+, Na+ and K+, or Mg2+), K0.5 was on the order of a few microM. In the absence of these ions, K0.5 was an order of magnitude lower (0.1 microM), indicating a higher apparent affinity. This effect suggests that the equilibrium for the ATP-induced conformational changes is dependent on the presence of metal ions.

Enhanced sampling by multiple molecular dynamics trajectories: carbonmonoxy myoglobin 10 μs A0 → A1–3 transition from ten 400 picosecond simulations

The utility of multiple trajectories to extend the time scale of molecular dynamics simulations is reported for the spectroscopic A-states of carbonmonoxy myoglobin (MbCO). Experimentally, the A0-->A(1-3) transition has been observed to be 10 micros at 300 K, which is beyond the time scale of standard molecular dynamics simulations. To simulate this transition, 10 short (400 ps) and two longer time (1.2 ns) molecular dynamics trajectories, starting from five different crystallographic and solution phase structures with random initial velocities centered in a 37 A radius sphere of water, have been used to sample the native-fold of MbCO. Analysis of the ensemble of structures gathered over the cumulative 5.6 ns reveals two biomolecular motions involving the side chains of His64 and Arg45 to explain the spectroscopic states of MbCO. The 10 micros A0-->A(1-3) transition involves the motion of His64, where distance between His64 and CO is found to vary up to 8.8 +/- 1.0 A during the transition of His64 from the ligand (A(1-3)) to bulk solvent (A0). The His64 motion occurs within a single trajectory only once, however the multiple trajectories populate the spectroscopic A-states fully. Consequently, multiple independent molecular dynamics simulations have been found to extend biomolecular motion from 5 ns of total simulation to experimental phenomena on the microsecond time scale.

Structural dynamics of myoglobin: effect of internal cavities on ligand migration and binding

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 sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation. Photoproduct intermediates, characterized by CO in different locations, were selectively enhanced by laser illumination at specific temperatures. Measurements were performed on the wild-type protein and a series of mutants (L104W, I107W, I28F, and I28W) in which bulky amino acid side chains were introduced to block passageways between cavities or to fill these sites. Binding of xenon was also employed as an alternative means of filling cavities. In all samples, photolyzed CO ligands were observed to initially bind at primary docking site B in the vicinity of the heme iron, from where they migrate to the secondary docking sites, the Xe4 and/or Xe1 cavities. To examine the relevance of these internal docking sites for physiological ligand binding, we have performed room-temperature flash photolysis on the entire set of proteins in the CO- and O(2)-bound form. Together with the cryospectroscopic results, these data provide a clear picture of the role of the internal sites for ligand escape from and binding to myoglobin.

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.

Conformational substates and dynamic properties of carbonmonoxy hemoglobin

Heme pocket dynamics of human carbonmonoxy hemoglobin (HbCO) is studied by Fourier transform infrared spectroscopy. The CO stretching band at various temperatures in the interval 300-10 K is analyzed in terms of three taxonomic A substates; however, in HbCO the band attributed to the A(1) taxonomic substate accounts for approximately 90% of the total intensity in the pH range 8.8-4.5. Two different regimes as a function of temperature are observed: below 160 K, the peak frequency and the bandwidth of the A(1) band have constant values whereas, above this temperature, a linear temperature dependence is observed, suggesting the occurrence of transitions between statistical substates within the A(1) taxonomic substate in this protein. The relationship between the heme pocket dynamics (as monitored by the thermal behavior of the CO stretching band), the overall dynamic properties of the protein matrix (as monitored by the thermal behavior of Amide II and Amide I' bands) and the glass transition of the solvent (as monitored by the thermal behavior of the bending band of water) is also investigated. From this analysis, we derive the picture of a very soft heme pocket of hemoglobin characterized by rather large anharmonic terms and strongly coupled to the dynamic properties of the solvent.

Structural dynamics of myoglobin: ligand migration among protein cavities studied by Fourier transform infrared/temperature derivative spectroscopy

Fourier transform infrared (FTIR) spectroscopy in the CO stretch bands combined with temperature derivative spectroscopy (TDS) was used to characterize intermediate states obtained by photolysis of two sperm whale mutant myoglobins, YQR (L29(B10)Y, H64(E7)Q, T67(E10)R) and YQRF (with an additional I107(G8)F replacement). Both mutants assume two different bound-state conformations, A(0) and A(3), which can be distinguished by their different CO bands near 1965 and 1933 cm(-1). They most likely originate from different conformations of the Gln-64 side chain. Within each A substate, a number of photoproduct states have been characterized on the basis of the temperature dependence of recombination in TDS experiments. Different locations and orientations of the ligand within the protein can be distinguished by the infrared spectra of the photolyzed CO. Recombination from the primary docking site, B, near the heme dominates below 50 K. Above 60 K, ligand rebinding occurs predominantly from a secondary docking site, C', in which the CO is trapped in the Xe4 cavity on the distal side, as shown by crystallography of photolyzed YQR and L29W myoglobin CO. Another kinetic state (C") has been identified from which rebinding occurs around 130 K. Moreover, a population appearing above the solvent glass transition at approximately 180 K (D state) is assigned to rebinding from the Xe1 cavity, as suggested by the photoproduct structure of the L29W sperm whale myoglobin mutant. For both the YQR and YQRF mutants, rebinding from the B sites near the heme differs for the two A substates, supporting the view that the return of the ligand from the C', C", and D states is not governed by the recombination barrier at the heme iron but rather by migration to the active site. Comparison of YQR and YQRF shows that access to the Xe4 site (C') is severely restricted by introduction of the bulky Phe side chain at position 107.

Cytochrome P-450-CO and substrates: lessons from ligand binding under high pressure

An overview of the application of high-pressure studies on the carbon monoxide complex of cytochrome P-450 is given. Different approaches to characterize ligand binding steps, the conformational states and substates and the compressibility of the ligand-bound complex are reviewed. A particular focus is the effect of substrates on these properties. It is shown that substrate mobility, compressibility and water accessibility are interrelated and may have functional meaning.

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.

FTIR and resonance Raman studies of nitric oxide binding to H93G cavity mutants of myoglobin

Nitric oxide (NO) binds to the myoglobin (Mb) cavity mutant, H93G, forming either a five- or six-coordinate Fe-NO complex. The H93G mutation eliminates the covalent attachment between the protein and the proximal ligand, allowing NO to bind H93G possibly from the proximal side of the heme rather than the typical diatomic binding pocket on the distal side. The question of whether NO binds on the distal or proximal side was addressed by FTIR spectroscopy of the N-O vibrational frequency nuN(-O) for a set of Mb mutants that perturb the electrostatic environment of the heme pocket. Vibrational spectra of five- and six-coordinate MbNO complexes indicate that nu(N-O) shifts (by as much as 26 cm(-1)) to higher energies for the distal mutants H64V and H64V/H93G relative to the energies of wild-type and H93G MbNO, while nu(N-O) is not affected by the proximal side mutation S92A/H93G. This result suggests that NO binds on the distal side of heme in the five- and six-coordinate MbNO complexes of H93G. Additionally, values of the Fe-NO vibrational frequency nu(Fe-NO) as measured by resonance Raman spectroscopy are reported for the distal and proximal double mutants of H93G. These results suggest that nu(Fe-NO) is not very sensitive to mutations that perturb the electrostatic environment of the heme pocket, leading to the observation that nu(N-O) and nu(Fe-NO) are not quantitatively correlated for the MbNO complexes presented here. Furthermore, nu(N-O) and nu(Fe-NO) do not correlate well with equilibrium constants for imidazole binding to the five-coordinate MbNO complexes of the H93G double mutants. The data presented here do not appear to support the presence of pi-back-bonding or an inverse trans effect of NO binding in Mb mutants that alter the electrostatic environment of the heme pocket.

The Fe(2+)-His(F8) Raman band shape of deoxymyoglobin reveals taxonomic conformational substates of the proximal linkage

The band shape of the Raman line attributed to the Fe(2+)-N(epsilon)(His(F8)) stretching mode in deoxymyoglobin contains significant information on the nature of the Fe-His proximal linkage. Raman lines appearing close to it, however, obscure the true line profile. To isolate this from its accompanying lines we use its isotopic shift of approximately 1 cm(-1) when (56)Fe in natural-abundance deoxymyoglobin is substituted by (54)Fe. This enables us to isolate the true line shape. We have measured this line shape in sperm whale myoglobin dissolved in a 66% vol/vol glycerol/water solution for nine temperatures from 10 K to 300 K. The nu(Fe-His) band shows a complex temperature-dependent profile, with a shoulder on its high-frequency wing, which becomes more prominent with increasing temperature. Detailed analysis reveals that the band is composed of five distinct lines attributable to taxonomic conformational substates of the nu(Fe-His) linkage. These are in thermodynamic equilibrium above the glass transition temperature T(f) but freeze in into the thermodynamic distribution at T(f) for lower temperatures. Alternative models that try to explain the nu(Fe-His) band shape by either an anharmonic coupling of the nu(Fe-His) to a low-frequency heme doming mode or by conformational substates related to a Gaussian distribution of iron out-of-plane displacements are at variance with the distinct features observed experimentally.

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.

Influence of the heme pocket conformation on the structure and vibrations of the Fe-CO bond in myoglobin: a QM/MM density functional study

The influence of the distal pocket conformation on the structure and vibrations of the heme-CO bond in carbonmonoxy myoglobin (MbCO) is investigated by means of hybrid QM/MM calculations based on density functional theory combined with a classical force field. It is shown that the heme-CO structure (QM treated) is quite rigid and not influenced by the distal pocket conformation (MM treated). This excludes any relation between FeCO distortions and the different CO absorptions observed in the infrared spectra of MbCO (A states). In contrast, both the CO stretch frequency and the strength of the CO...His64 interaction are very dependent on the orientation and tautomerization state of His64. Our calculations indicate that the CO...N(epsilon) type of approach does not contribute to the A states, whereas the CO...H-N(epsilon) interaction is the origin of the A(1) and A(3) states, the His64 residue being protonated at N(epsilon). The strength of the CO...His64 interaction is quantified, in comparison with the analogous O(2)...His64 interaction and with the observed changes in the CO stretch frequency. Additional aspects of the CO...His64 interaction and its biological implications are discussed.