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.     

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₂, 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 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.     

Viscosity-dependent relaxation significantly modulates the kinetics of CO recombination in the truncated hemoglobin TrHbN from Mycobacterium tuberculosis

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

Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants

We have combined Fourier transform infrared/temperature derivative (FTIR-TDS) spectroscopy at cryogenic temperatures and flash photolysis at ambient temperature to examine the effects of polar and bulky amino acid replacements of the highly conserved distal valine 68 in sperm whale myoglobin. In FTIR-TDS experiments, the CO ligand can serve as an internal voltmeter that monitors the local electrostatic field not only at the active site but also at intermediate ligand docking sites. Mutations of residue 68 alter size, shape, and electric field of the distal pocket, especially in the vicinity of the primary docking site (state B). As a consequence, the infrared bands associated with the ligand at site B are shifted. The effect is most pronounced in mutants with large aromatic side chains. Polar side chains (threonine or serine) have only little effect on the peak frequencies. Ligands that migrate toward more remote sites C and D give rise to IR bands with altered frequencies. TDS experiments separate the photoproducts according to their recombination temperatures. The rates and extent of ligand migration among internal cavities at cryogenic temperatures can be used to interpret geminate and bimolecular O2 and CO recombination at room temperature. The kinetics of geminate recombination can be explained by steric arguments alone, whereas both the polarity and size of the position 68 side chain play major roles in regulating bimolecular ligand binding from the solvent.     

Structural heterogeneity and ligand gating in ferric Methanosarcina acetivorans protoglobin mutants

Protoglobin from Methanosarcina acetivorans C2A (MaPgb), a strictly anaerobic methanogenic Archaea, displays peculiar structural and functional properties within members of the hemoglobin superfamily. In fact, MaPgb-specific loops and a N-terminal extension (20 amino acid residues) completely bury the heme within the protein matrix. Therefore, the access of diatomic gaseous molecules to the heme is granted by two apolar tunnels reaching the heme distal site from locations at the B/G and B/E helix interfaces. The presence of two tunnels within the protein matrix could be partly responsible for the slightly biphasic ligand binding behavior. Unusually, MaPgb oxygenation is favored with respect to carbonylation. Here, the crucial role of Tyr(B10)61 and Ile(G11)149 residues, located in the heme distal site and lining the protein matrix tunnels 1 and 2, respectively, on ligand binding to the heme-Fe-atom and on distal site structural organization is reported. In particular, tunnel 1 accessibility is modulated by a complex reorganization of the Trp(B9)60 and Phe(E11)93 side-chains, triggered by mutations of the Tyr(B10)61 and Ile(G11)149 residues, and affected by the presence and type of the distal heme-bound ligand.     

Protein conformational relaxation and ligand migration in myoglobin: a nanosecond to millisecond molecular movie from time-resolved Laue X-ray diffraction.

A time-resolved Laue X-ray diffraction technique has been used to explore protein relaxation and ligand migration at room temperature following photolysis of a single crystal of carbon monoxymyoglobin. The CO ligand is photodissociated by a 7.5 ns laser pulse, and the subsequent structural changes are probed by 150 ps or 1 micros X-ray pulses at 14 laser/X-ray delay times, ranging from 1 ns to 1.9 ms. Very fast heme and protein relaxation involving the E and F helices is evident from the data at a 1 ns time delay. The photodissociated CO molecules are detected at two locations: at a distal pocket docking site and at the Xe 1 binding site in the proximal pocket. The population by CO of the primary, distal site peaks at a 1 ns time delay and decays to half the peak value in 70 ns. The secondary, proximal docking site reaches its highest occupancy of 20% at approximately 100 ns and has a half-life of approximately 10 micros. At approximately 100 ns, all CO molecules are accounted for within the protein: in one of these two docking sites or bound to the heme. Thereafter, the CO molecules migrate to the solvent from which they rebind to deoxymyoglobin in a bimolecular process with a second-order rate coefficient of 4.5 x 10(5) M(-1) s(-1). Our results also demonstrate that structural changes as small as 0.2 A and populations of CO docking sites of 10% can be detected by time-resolved X-ray diffraction.     

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

We have investigated CO migration and binding in Cu_BMb, 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 Cu_BMb 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 Cu_BMb, primary and secondary ligand docking sites are accessible upon photodissociation. In copper-bound Cu_BMb, 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 Cu_BMb. 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 Cu_B in the active of heme–copper oxidases. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non‐O2‐binding H‐NOX domain

The Nostoc sp (Ns) H‐NOX (heme‐nitric oxide or OXygen‐binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y‐shaped tunnel system hosted in the Ns H‐NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881–889] recently using x‐ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H‐NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O₂ migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H‐NOX was identified. On the contrary, CO and O₂ would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H‐NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group. Proteins 2013; 81:1363–1376. © 2013 Wiley Periodicals, Inc.     

Mapping protein matrix cavities in human cytoglobin through Xe atom binding

Cytoglobin is the fourth recognized globin type, almost ubiquitously distributed in human tissues; its function is still poorly understood. Cytoglobin displays a core region of about 150 residues, structurally related to hemoglobin and myoglobin, and two extra segments, about 20 residues each, at the N- and C-termini. The core region hosts a large apolar cavity, held to provide a ligand diffusion pathway to/from the heme, and/or ligand temporary docking sites. Here we report the crystal structure (2.4A resolution, R-factor 19.1%) of a human cytoglobin mutant bearing the CysB2(38) --> Ser and CysE9(83) --> Ser substitutions (CYGB*), treated under pressurized xenon. Three Xe atoms bind to the heme distal site region of CYGB* mapping the protein matrix apolar cavity. Despite the conserved globin fold, the cavity found in CYGB* is structured differently from those recognized to play a functional role in myoglobin, neuroglobin, truncated hemoglobins, and Cerebratulus lacteus mini-hemoglobin.     

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.     

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

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

Unbinding pathways of an agonist and an antagonist from the 5-HT3 receptor

The binding sites of 5-HT3 and other Cys-loop receptors have been extensively studied, but there are no data on the entry and exit routes of ligands for these sites. Here we have used molecular dynamics simulations to predict the pathway for agonists and antagonists exiting from the 5-HT3 receptor binding site. The data suggest that the unbinding pathway follows a tunnel at the interface of two subunits, which is approximately 8 A long and terminates approximately 20 A above the membrane. The exit routes for an agonist (5-HT) and an antagonist (granisetron) were similar, with trajectories toward the membrane and outward from the ligand binding site. 5-HT appears to form many hydrogen bonds with residues in the unbinding pathway, and experiments show that mutating these residues significantly affects function. The location of the pathway is also supported by docking studies of granisetron, which show a potential binding site for granisetron on the unbinding route. We propose that leaving the binding pocket along this tunnel places the ligands close to the membrane and prevents their immediate reentry into the binding pocket. We anticipate similar exit pathways for other members of the Cys-loop receptor family.     

Crystal structure of cytoglobin: the fourth globin type discovered in man displays heme hexa-coordination

Cytoglobin is a recently discovered hemeprotein belonging to the globin superfamily together with hemoglobin, myoglobin and neuroglobin. Although distributed in almost all human tissues, cytoglobin has not been ascribed a specific function. Human cytoglobin is composed of 190 amino acid residues. Sequence alignments show that a protein core region (about 150 residues) is structurally related to hemoglobin and myoglobin, being complemented by about 20 extra residues both on the N and C termini. In the absence of exogenous ligands (e.g. O2), the cytoglobin distal HisE7 residue is coordinated to the heme Fe atom, thus decreasing the ligand affinity. The crystal structure of human cytoglobin (2.1 A resolution, 21.3% R-factor) highlights a three-over-three alpha-helical globin fold, covering residues 18-171; the 1-17 N-terminal, and the 172-190 C-terminal residue segments are disordered in both molecules of the crystal asymmetric unit. Heme hexa-coordination is evident in one of the two cytoglobin chains, whereas alternate conformation for the heme distal region, achieving partial heme penta-coordination, is observed in the other. Human cytoglobin displays a large apolar protein matrix cavity, next to the heme, not related to the myoglobin cavities recognized as temporary ligand docking stations. The cavity, which may provide a heme ligand diffusion pathway, is connected to the external space through a narrow tunnel nestled between the globin G and H helices.     

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.     

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.     

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.     

Heme-ligand tunneling in group I truncated hemoglobins

Truncated hemoglobins (trHbs) are small hemoproteins forming a separate cluster within the hemoglobin superfamily; their functional roles in bacteria, plants, and unicellular eukaryotes are marginally understood. Crystallographic investigations have shown that the trHb fold (a two-on-two alpha-helical sandwich related to the globin fold) hosts a protein matrix tunnel system offering a potential path for ligand diffusion to the heme distal site. The tunnel topology is conserved in group I trHbs, although with modulation of its size/structure. Here, we present a crystallographic investigation on trHbs from Mycobacterium tuberculosis, Chlamydomonas eugametos, and Paramecium caudatum, showing that treatment of trHb crystals under xenon pressure leads to binding of xenon atoms at specific (conserved) sites along the protein matrix tunnel. The crystallographic results are in keeping with data from molecular dynamics simulations, where a dioxygen molecule is left free to diffuse within the protein matrix. Modulation of xenon binding over four main sites is related to the structural properties of the tunnel system in the three trHbs and may be connected to their functional roles. In a parallel crystallographic investigation on M. tuberculosis trHbN, we show that butyl isocyanide also binds within the apolar tunnel, in excellent agreement with concepts derived from the xenon binding experiments. These results, together with recent data on atypical CO rebinding kinetics to group I trHbs, underline the potential role of the tunnel system in supporting diffusion, but also accumulation in multiple copies, of low polarity ligands/molecules within group I trHbs.     

Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.     

Cavities and packing defects in the structural dynamics of myoglobin

Small globular proteins contain internal cavities and packing defects that reduce thermodynamic stability but seem to play a role in controlling function by defining pathways for the diffusion of the ligand/substrate to the active site. In the case of myoglobin (Mb), a prototype for structure–function relationship studies, the photosensitivity of the adduct of the reduced protein with CO, O₂ and NO allows events related to the migration of the ligand through the matrix to be followed. The crystal structures of intermediate states of wild-type (wt) and mutant Mbs show the photolysed CO to be located either in the distal heme pocket (primary docking site) or in one of two alternative cavities (secondary docking sites) corresponding to packing defects accessible to an atom of xenon. These results convey the general picture that pre-existing internal cavities are involved in controlling the dynamics and reactivity of the reactions of Mb with O₂ and other ligands, including NO.