Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

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.

Neuroglobin, clues to function and mechanism

Neuroglobin is expressed in vertebrate brain and belongs to a branch of the globin family that diverged early in evolution. Sequence conservation and presence in nervous cells of several taxa suggests a relevant role in the nervous system, with tight structural restraints. Twenty years after its discovery, a rich scientific literature provides convincing evidence of the involvement of neuroglobin in sustaining neuron viability in physiological and pathological conditions however, a full and conclusive picture of its specific function, or set of functions is still lacking. The difficulty of unambiguously assigning a precise mechanism and biochemical role to neuroglobin might arise from the participation to one or more cell mechanism that redundantly guarantee the functioning of the highly specialized and metabolically demanding central nervous system of vertebrates. Here we collect findings and hypotheses arising from recent biochemical, biophysical, structural, in cell and in vivo experimental work on neuroglobin, aiming at providing an overview of the most recent literature. Proteins are said to have jobs and hobbies, it is possible that, in the case of neuroglobin, evolution has selected for it more than one job, and support to cover for its occasional failings. Disentangling the mechanisms and roles of neuroglobin is thus a challenging task that might be achieved by considering data from different disciplines and experimental approaches.

Structural and (Pseudo-)Enzymatic Properties of Neuroglobin: Its Possible Role in Neuroprotection

Neuroglobin (Ngb), the third member of the globin family, was discovered in human and murine brains in 2000. This monomeric globin is structurally similar to myoglobin (Mb) and hemoglobin (Hb) α and β subunits, but it hosts a bis-histidyl six-coordinated heme-Fe atom. Therefore, the heme-based reactivity of Ngb is modulated by the dissociation of the distal HisE7-heme-Fe bond, which reflects in turn the redox state of the cell. The high Ngb levels (~100–200 μM) present in the retinal ganglion cell layer and in the optic nerve facilitate the O₂ buffer and delivery. In contrast, the very low levels of Ngb (~1 μM) in most tissues and organs support (pseudo-)enzymatic properties including NO/O₂ metabolism, peroxynitrite and free radical scavenging, nitrite, hydroxylamine, hydrogen sulfide reduction, and the nitration of aromatic compounds. Here, structural and (pseudo-)enzymatic properties of Ngb, which are at the root of tissue and organ protection, are reviewed, envisaging a possible role in the protection from neuronal degeneration of the retina and the optic nerve.

Ligand pathways in neuroglobin revealed by low-temperature photodissociation and docking experiments

A combined biophysical approach was applied to map gas-docking sites within murine neuroglobin (Ngb), revealing snapshots of events that might govern activity and dynamics in this unique hexacoordinate globin, which is most likely to be involved in gas-sensing in the central nervous system and for which a precise mechanism of action remains to be elucidated. The application of UV-visible microspectroscopy in crystallo, solution X-ray absorption near-edge spectroscopy and X-ray diffraction experiments at 15-40 K provided the structural characterization of an Ngb photolytic intermediate by cryo-trapping and allowed direct observation of the relocation of carbon monoxide within the distal heme pocket after photodissociation. Moreover, X-ray diffraction at 100 K under a high pressure of dioxygen, a physiological ligand of Ngb, unravelled the existence of a storage site for O2 in Ngb which coincides with Xe-III, a previously described docking site for xenon or krypton. Notably, no other secondary sites were observed under our experimental conditions.

Proximal and distal control for ligand binding in neuroglobin: role of the CD loop and evidence for His64 gating

Neuroglobin (Ngb) is predominantly expressed in neurons of the central and peripheral nervous systems and it clearly seems to be involved in neuroprotection. Engineering Ngb to observe structural and dynamic alterations associated with perturbation in ligand binding might reveal important structural determinants, and could shed light on key features related to its mechanism of action. Our results highlight the relevance of the CD loop and of Phe106 as distal and proximal controls involved in ligand binding in murine neuroglobin. We observed the effects of individual and combined mutations of the CD loop and Phe106 that conferred to Ngb higher CO binding velocities, which we correlate with the following structural observations: the mutant F106A shows, upon CO binding, a reduced heme sliding hindrance, with the heme present in a peculiar double conformation, whereas in the CD loop mutant “Gly-loop”, the original network of interactions between the loop and the heme was abolished, enhancing binding via facilitated gating out of the distal His64. Finally, the double mutant, combining both mutations, showed a synergistic effect on CO binding rates. Resonance Raman spectroscopy and MD simulations support our findings on structural dynamics and heme interactions in wild type and mutated Ngbs.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Engineering the internal cavity of neuroglobin demonstrates the role of the haem-sliding mechanism

Neuroglobin is a member of the globin family involved in neuroprotection; it is primarily expressed in the brain and retina of vertebrates. Neuroglobin belongs to the heterogeneous group of hexacoordinate globins that have evolved in animals, plants and bacteria, endowed with the capability of reversible intramolecular coordination, allowing the binding of small gaseous ligands (O2, NO and CO). In a unique fashion among haemoproteins, ligand-binding events in neuroglobin are dependent on the sliding of the haem itself within a preformed internal cavity, as revealed by the crystal structure of its CO-bound derivative. Point mutants of the neuroglobin internal cavity have been engineered and their functional and structural characterization shows that hindering the haem displacement leads to a decrease in CO affinity, whereas reducing the cavity volume without interfering with haem sliding has negligible functional effects.

The crystal structure of wild-type human brain neuroglobin reveals flexibility of the disulfide bond that regulates oxygen affinity

Neuroglobin plays an important function in the supply of oxygen in nervous tissues. In human neuroglobin, a cysteine at position 46 in the loop connecting the C and D helices of the globin fold is presumed to form an intramolecular disulfide bond with Cys55. Rupture of this disulfide bridge stabilizes bi-histidyl haem hexacoordination, causing an overall decrease in the affinity for oxygen. Here, the first X-ray structure of wild-type human neuroglobin is reported at 1.74 Å resolution. This structure provides a direct observation of two distinct conformations of the CD region containing the intramolecular disulfide link and highlights internal cavities that could be involved in ligand migration and/or are necessary to enable the conformational transition between the low and high oxygen-affinity states following S-S bond formation.

Crucial roles of Glu60 in human neuroglobin as a guanine nucleotide dissociation inhibitor and neuroprotective agent

Mammalian neuroglobin (Ngb) protects neuronal cells under conditions of oxidative stress. We previously showed that human Ngb acts as a guanine nucleotide dissociation inhibitor (GDI) for the α-subunits of heterotrimeric Gi/o proteins and inhibits reductions in cAMP concentration, leading to protection against cell death. In the present study, we created human E60Q Ngb mutant and clarified that Glu60 of human Ngb is a crucial residue for its GDI and neuroprotective activities. Moreover, we investigated structural and functional properties of several human Ngb mutants and demonstrated that the neuroprotective effect of human Ngb is due to its GDI activity and not due to its scavenging activity against reactive oxygen species.

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

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

Human neuroglobin functions as an oxidative stress-responsive sensor for neuroprotection

Mammalian neuroglobin (Ngb) protects neuronal cells under conditions of oxidative stress. The mechanism underlying this function is only partly understood. Here, we report that human Ngb exists in lipid rafts only during oxidative stress and that lipid rafts are crucial for neuroprotection by Ngb. The ferrous oxygen-bound form of Ngb, which exists under normoxia, is converted to the ferric bis-His conformation during oxidative stress, inducing large tertiary structural changes. We clarified that ferric bis-His Ngb, but not ferrous ligand-bound Ngb, specifically binds to flotillin-1, a lipid raft microdomain-associated protein, as well as to α-subunits of heterotrimeric G proteins (Gα(i/o)). Moreover, we found that human ferric bis-His Ngb acts as a guanine nucleotide dissociation inhibitor for Gα(i/o) that has been modified by oxidative stress. In addition, our data shows that Ngb inhibits the decrease in cAMP concentration that occurs under oxidative stress, leading to protection against cell death. Furthermore, by using a mutated Ngb protein that cannot form the bis-His conformation, we demonstrate that the oxidative stress-induced structural changes of human Ngb are essential for its neuroprotective activity.

Carbon Dioxide Migration Pathways in Proteins

Some of the most important biological processes, such as carbon fixation, are dependent on protein-gas interactions. The motion of CO2 through the enzyme phosphoenolpyruvate carboxykinase was investigated using extensive all-atom molecular dynamics simulations. Three discrete migration pathways were located, suggesting the protein directs the movement of CO2. The chemical nature of these pathways is discussed, as are their biotechnological ramifications.

Structure and dynamics of Antarctic fish neuroglobin assessed by computer simulations

Neuroglobin (Ngb) is a heme protein, highly conserved along evolution, predominantly found in the nervous system. It is upregulated by hypoxia and ischemia and may have a neuroprotective role under hypoxic stress. Although many other roles have been proposed, the physiological function is still unclear. Antarctic icefishes lack hemoglobin and some species also lack myoglobin, but all have Ngb and thus may help the elucidation of Ngb function. We present the first theoretically derived structure of fish Ngb and describe its behavior using molecular dynamics simulations. Specifically, we sequenced and analyzed Ngbs from a colorless-blooded Antarctic icefish species Chaenocephalus aceratus and a related red-blooded species (Dissostichus mawsoni). Both fish Ngbs are 6-coordinated but have some peculiarities that differentiate them from mammalian counterparts: they have extensions in the N and C termini that can interact with the EF loop, and a gap in the alignment that changes the CD-region structure/dynamics that has been found to play a key role in human neuroglobin. Our results suggest that a single mutation between both fish Ngbs is responsible for significant difference in the behavior of the proteins. The functional role of these characteristics is discussed.

Kinetics of carbon monoxide migration and binding in solvated neuroglobin as revealed by molecular dynamics simulations and quantum mechanical calculations

Neuroglobin (Ngb) is a globular protein that reversibly binds small ligands at the six coordination position of the heme. With respect to other globins similar to myoglobin, Ngb displays some peculiarities as the topological reorganization of the internal cavities coupled to the sliding of the heme, or the binding of the endogenous distal histidine to the heme in the absence of an exogenous ligand. In this Article, by using multiple (independent) molecular dynamics trajectories (about 500 ns in total), the migration pathways of photolized carbon monoxide (CO) within solvated Ngb were analyzed, and a quantitative description of CO migration and corresponding kinetics was obtained. MD results, combined with quantum mechanical calculations on the CO-heme binding-unbinding reaction step in Ngb, allowed construction of a quantitative model representing the relevant steps of CO migration and rebinding.

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

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.

Probing the role of the internal disulfide bond in regulating conformational dynamics in neuroglobin

In this report, we demonstrate that the internal disulfide bridge in human neuroglobin modulates structural changes associated with ligand photo-dissociation from the heme active site. This is evident from time-resolved photothermal studies of CO photo-dissociation, which reveal a 13.4+/-0.9 mL mol(-1) volume expansion upon ligand photo-release from human neuroglobin, whereas the CO dissociation from rat neuroglobin leads to a significantly smaller volume change (DeltaV=4.6+/-0.3 mL mol(-1)). Reduction of the internal disulfide bond in human neuroglobin leads to conformational changes (reflected by DeltaV) nearly identical to those observed for rat Ngb. Our data favor the hypothesis that the disulfide bond between Cys46 and Cys55 modulates the functioning of human neuroglobin.

Ligand migration through the internal hydrophobic cavities in human neuroglobin

Neuroglobin (Ngb), a member of the globin superfamily, was found in the brain of vertebrates and is suggested to play a neuroprotective function under hypoxic conditions by scavenging nitrogen monoxide (NO) through a dioxygenase activity. In order for such a reaction to efficiently take place and to minimize the release of reactive intermediates in the cytosol, the cosubstrates O(2) and NO and other unstable reaction intermediates should bind sequentially to docking sites in the protein matrix. We have characterized the accessibility of these sites by analyzing the geminate CO rebinding kinetics to the heme moiety observed upon nanosecond flash photolysis of the Ngb-CO complex encapsulated in silica gels. The geminate rebinding phase showed a remarkable complexity, revealing the presence of a system of secondary docking sites where ligands are stored for hundreds of microseconds. Most kinetics steps display little temperature dependence, demonstrating that ligands can easily migrate through the cavities, except for the slowest reaction intermediate, possibly reflecting a structural conformational change reshaping the system of cavities. This conformational change is unrelated with distal His E7 binding to the heme, as it persists for the HE7L mutant. Overall, data are consistent with the presence of a discrete system of docking sites, possibly acting as reservoirs for the putative cosubstrates and for other reactive species involved in the physiologically relevant reaction.

The structure of neuroglobin at high Xe and Kr pressure reveals partial conservation of globin internal cavities

Neuroglobin (Ngb) is a hexacoordinate globin expressed in the brain of vertebrates. Ferrous Ngb binds dioxygen with high affinity and the O(2) adduct is able to scavenge NO. Convincing in vitro and in vivo data indicate that Ngb is involved in neuroprotection during hypoxia and ischemia. The 3D structure of Ngb reveals the presence of a wide internal cavity connecting its heme active site with the bulk. To explore the role of this "tunnel" in the control of ligand binding, we determined the structure of metNgb and NgbCO equilibrated with Xe or Kr. We show four docking sites for Xe (only two for Kr); two of the four Xe sites are within the large cavity. They are only partially conserved in globins, since the two proximal Xe sites identified in myoglobin (Xe1 and Xe2) are absent in Ngb, as well as in cytoglobin. The Xe docking sites in Ngb map a pathway within the protein matrix, leading to the heme, which becomes more accessible in the ligand-bound species. This may be of significance in connection with the redox chemistry that may be the primary function of this hexacoordinate globin.

Molecular dynamics simulation of the neuroglobin crystal: comparison with the simulation in solution

Neuroglobin (Ngb) is a monomeric protein that, despite the small sequence similarity with other globins, displays the typical globin fold. In the absence of exogenous ligands, the ferric and the ferrous forms of Ngb are both hexacoordinated to the distal and proximal histidines. In the ferrous form, oxygen, nitric oxide or carbon monoxide can displace the distal histidine, yielding a reversible adduct. Crystallographic data show that the binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities. Molecular dynamics (MD) simulations in solution show that the heme oscillates between two positions, much as the ones observed in the crystal structure, although the occupancy is different. The simulations also suggest that ligand binding in solution can affect the flexibility and conformation of residues connecting the C and D helices, referred to as the CD corner, which is coupled to the configuration adopted by the distal histidine. In this study, we report the results of 30 ns MD simulations of CO-bound Ngb in the crystal. Our goal was to compare the protein dynamical behavior in the crystal with the results supplied by the previous MD simulation of CO-bound Ngb in solution and the x-ray experimental data. The results show that the different environments (crystal or solution) affect the dynamics of the heme group and of the CD corner.

Exploring the molecular basis of heme coordination in human neuroglobin

Neuroglobin (Ngb), a recently discovered ancient heme protein, presents the typical globin fold and is around 20% identical to myoglobin (Mb). In contrast with Mb, however, its heme is hexacoordinated (6c). It is expressed in the nervous system and has been the subject of numerous investigations in the last years, but its function is still unclear. The proposed roles include oxygen transport, reactive oxygen species (ROS) detoxification, hypoxia protection, and redox state sensing. All proposed functions require distal histidine dissociation from the heme to yield a reactive iron. With the aim of understanding the 6c to 5c transition, we have performed molecular dynamics simulations for ferrous Ngb in the 6c, 5c, and oxy states. We also computed free energy profiles associated with the transition employing an advanced sampling technique. Finally, we studied the effect of the redox state of CysCD7 and CysD5, which are known to form a disulfide bridge. Our results show that protein oxidation promotes a stabilization of the pentacoordinated species, thus favoring the protein to adopt the more reactive state and supporting the existence of a molecular mechanism whereby O2 would be released under hypoxic conditions, thereby suggesting an O(2) storage function for Ngb. Taken together, our results provide structural information not available experimentally which may shed light on the protein proposed functions, particularly as a redox sensor.

Disulfide bond influence on protein structural dynamics probed with 2D-IR vibrational echo spectroscopy

Intramolecular disulfide bonds are understood to play a role in regulating protein stability and activity. Because disulfide bonds covalently link different components of a protein, they influence protein structure. However, the effects of disulfide bonds on fast (subpicosecond to approximately 100 ps) protein equilibrium structural fluctuations have not been characterized experimentally. Here, ultrafast 2D-IR vibrational echo spectroscopy is used to examine the constraints an intramolecular disulfide bond places on the structural fluctuations of the protein neuroglobin (Ngb). Ngb is a globin family protein found in vertebrate brains that binds oxygen reversibly. Like myoglobin (Mb), Ngb has the classical globin fold and key residues around the heme are conserved. Furthermore, the heme-ligated CO vibrational spectra of Mb (Mb-CO) and Ngb (Ngb-CO) are virtually identical. However, in contrast to Mb, human Ngb has an intramolecular disulfide bond that affects its oxygen affinity and protein stability. By using 2D-IR vibrational echo spectroscopy, we investigated the equilibrium protein dynamics of Ngb-CO by observing the CO spectral diffusion (time dependence of the 2D-IR line shapes) with and without the disulfide bond. Despite the similarity of the linear FTIR spectra of Ngb-CO with and without the disulfide bond, 2D-IR measurements reveal that the equilibrium sampling of different protein configurations is accelerated by disruption of the disulfide bond. The observations indicate that the intramolecular disulfide bond in Ngb acts as an inhibitor of fast protein dynamics even though eliminating it does not produce significant conformational change in the protein's structure.