Equilibrium O2 distribution in the Zn2+-protoporphyrin IX deoxymyoglobin mimic: application to oxygen migration pathway analysis

Nuclear magnetic resonance-based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine-free wild-type T4 lysozyme was investigated using O₂ gas-pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂ -induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C-terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid-search. Two O₂ -accessible sites were identified in internal cavities: One site was consistent with the xenon-binding site in the protein in crystal, and the other site was established to be a novel ligand-binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground-state ensemble transiently develop channels for O₂ association with the internal protein cavities.

Analysis of O2-binding Sites in Proteins Using Gas-Pressure NMR Spectroscopy: Outer Surface Protein A

Internal cavities in proteins produce conformational fluctuations and enable the binding of small ligands. Here, we report a NMR analysis of O₂-binding sites by O₂-induced paramagnetic relaxation enhancements (PREs) on amide groups of proteins in solution. Outer surface protein A contains a nonglobular single-layer β-sheet that connects the N- and C-terminal globular domains. Several cavities have been observed in both domains of the crystallized protein structure. The receptor-binding sites are occluded and line the largest cavity of the C-terminal domain. We observed significant O₂-induced PREs for amide protons located around the largest cavity and at the central β-sheet. We suggested three potential O₂-accessible sites in the protein based on the 1/r⁶ distance dependence of the PRE. Two sites were in or close to the largest cavity and the third site was in the surface crevice of the central β-sheet. These results provide, to our knowledge, the first evidence of ligand binding to the surface crevice and cavity of the protein in solution. Because O₂ generally binds more specifically to hydrophobic rather than hydrophilic cavities within a protein, the results also indicated that the receptor-binding sites lining the largest cavity were in the hydrophobic environment in the ground-state conformation. Molecular dynamics simulations permitted the visualization of the rotational and translational motions of O₂ within the largest cavity, egress of O₂ from the cavity, and ingress of O₂ in the surface crevice of the β-sheet. These molecular dynamics simulation results qualitatively explained the O₂-induced changes in NMR observations. Exploring cavities that are sufficiently dynamic to enable access by small molecules can be a useful strategy for the design of stable proteins and their ligands.

Detecting O2 binding sites in protein cavities

Internal cavities are important elements in protein structure, dynamics, stability and function. Here we use NMR spectroscopy to investigate the binding of molecular oxygen (O₂) to cavities in a well-studied model for ligand binding, the L99A mutant of T4 lysozyme. On increasing the O₂ concentration to 8.9 mM, changes in ¹H, ¹⁵N, and ¹³C chemical shifts and signal broadening were observed specifically for backbone amide and side chain methyl groups located around the two hydrophobic cavities of the protein. O₂-induced longitudinal relaxation enhancements for amide and methyl protons could be adequately accounted for by paramagnetic dipolar relaxation. These data provide the first experimental demonstration that O₂ binds specifically to the hydrophobic, and not the hydrophilic cavities, in a protein. Molecular dynamics simulations visualized the rotational and translational motions of O₂ in the cavities, as well as the binding and egress of O₂, suggesting that the channel consisting of helices D, E, G, H, and J could be the potential gateway for ligand binding to the protein. Due to strong paramagnetic relaxation effects, O₂ gas-pressure NMR measurements can detect hydrophobic cavities when populated to as little as 1%, and thereby provide a general and highly sensitive method for detecting oxygen binding in proteins.

Modulating distal cavities in the α and β subunits of human HbA reveals the primary ligand migration pathway

The free volume in the active site of human HbA plays a crucial role in governing the bimolecular rates of O(2), CO, and NO binding, the fraction of geminate ligand recombination, and the rate of NO dioxygenation by the oxygenated complex. We have decreased the size of the distal pocket by mutating Leu(B10), Val(E11), and Leu(G8) to Phe and Trp and that of other more internal cavities by filling them with Xe at high gas pressures. Increasing the size of the B10 side chain reduces bimolecular rates of ligand binding nearly 5000-fold and inhibits CO geminate recombination due to both reduction of the capture volume in the distal pocket and direct steric hindrance of Fe-ligand bond formation. Phe and Trp(E11) mutations also cause a decrease in distal pocket volume but, at the same time, increase access to the Fe atom because of the loss of the γ2 CH(3) group of the native Val(E11) side chain. The net result of these E11 substitutions is a dramatic increase in the rate of geminate recombination because dissociated CO is sequestered close to the Fe atom and can rapidly rebind without steric resistance. However, the bimolecular rate constants for binding of ligand to the Phe and Trp(E11) mutants are decreased 5-30-fold, because of a smaller capture volume. Geminate and bimolecular kinetic parameters for Phe and Trp(G8) mutants are similar to those for the native HbA subunits because the aromatic rings at this position cause little change in distal pocket volume and because ligands do not move past this position into the globin interior of wild-type HbA subunits. The latter conclusion is verified by the observation that Xe binding to the α and β Hb subunits has little effect on either geminate or bimolecular ligand rebinding. All of these experimental results argue strongly against alternative ligand migration pathways that involve movements through the protein interior in HbA. Instead, ligands appear to enter through the His(E7) gate and are captured directly in the distal cavity.

The stretching frequencies of bound alkyl isocyanides indicate two distinct ligand orientations within the distal pocket of myoglobin

The FTIR spectra for alkyl isocyanides (CNRs) change from a single nu(CN) band centered at approximately 2175 cm(-1) to two peaks at approximately 2075 and approximately 2125 cm(-1) upon binding to sperm whale myoglobin (Mb). The low- and high-frequency peaks have been assigned to in and out conformations, respectively. In the in conformation, the ligand is pointing toward the protein interior, and the distal His64(E7) is in a closed position, donates a H-bond to the bound isocyano group, enhances back-bonding, and lowers the C-N bond order. In the out conformation, the ligand side chain points toward solvent through a channel opened by outward rotation of His64. Loss of positive polarity near the binding site causes an increase in C-N bond order. Support for this interpretation is threefold: (1) similar shifts to lower frequency occur for MbCO complexes when H-bond donation from His64(E7) occurs; (2) only one peak at approximately 2125 cm(-1), indicative of an apolar environment, is observed for CNRs bound to H64A or H64L Mb mutants or to chelated protoheme in soap micelles; and (3) the fraction of in conformation based on FTIR spectra correlates strongly with the fraction of geminate recombination after nanosecond laser photolysis. The in alkyl side chain conformation causes the photodissociated ligand to be "stuck" in the distal pocket, promoting internal rebinding, whereas the out conformation inhibits geminate recombination because part of the ligand is already in an open E7 channel, poised for rapid escape.

Optical detection of disordered water within a protein cavity

Internal water molecules are important to protein structure and function, but positional disorder and low occupancies can obscure their detection by X-ray crystallography. Here, we show that water can be detected within the distal cavities of myoglobin mutants by subtle changes in the absorbance spectrum of pentacoordinate heme, even when the presence of solvent is not readily observed in the corresponding crystal structures. A well-defined, noncoordinated water molecule hydrogen bonded to the distal histidine (His64) is seen within the distal heme pocket in the crystal structure of wild type (wt) deoxymyoglobin. Displacement of this water decreases the rate of ligand entry into wt Mb, and we have shown previously that the entry of this water is readily detected optically after laser photolysis of MbCO complexes. However, for L29F and V68L Mb no discrete positions for solvent molecules are seen in the electron density maps of the crystal structures even though His64 is still present and slow rates of ligand binding indicative of internal water are observed. In contrast, time-resolved perturbations of the visible absorption bands of L29F and V68L deoxyMb generated after laser photolysis detect the entry and significant occupancy of water within the distal pockets of these variants. Thus, the spectral perturbation of pentacoordinate heme offers a potentially robust system for measuring nonspecific hydration of the active sites of heme proteins.

Spin relaxation measurements of electrostatic bias in intermolecular exploration

We utilize the paramagnetic contribution to proton spin-lattice relaxation rate constants induced by freely diffusing charged paramagnetic centers to investigate the effect of charge on the intermolecular exploration of a protein by the small molecule. The proton NMR spectrum provided 255 resolved resonances that report how the explorer molecule local concentration varies with position on the surface. The measurements integrate over local dielectric constant variations, and, in principle, provide an experimental characterization of the surface free energy sampling biases introduced by the charge distribution on the protein. The experimental results for ribonuclease A obtained using positive, neutral, and negatively charged small nitroxide radicals are qualitatively similar to those expected from electrostatic calculations. However, while systematic electrostatic trends are apparent, the three different combinations of the data sets do not yield internally consistent values for the electrostatic contribution to the intermolecular free energy. We attribute this failure to the weakness of the electrostatic sampling bias for charged nitroxides in water and local variations in effective translational diffusion constant at the water-protein interface, which enters the nuclear spin relaxation equations for the nitroxide-proton dipolar coupling.

Water and ligand entry in myoglobin: assessing the speed and extent of heme pocket hydration after CO photodissociation

A previously undescribed spectrokinetic assay for the entry of water into the distal heme pocket of wild-type and mutant myoglobins is presented. Nanosecond photolysis difference spectra were measured in the visible bands of sperm whale myoglobin as a function of distal pocket mutation and temperature. A small blue shift in the 560-nm deoxy absorption peak marked water entry several hundred nanoseconds after CO photodissociation. The observed rate suggests that water entry is rate-limited by the escape of internal dissociated CO. The heme pocket hydration and geminate recombination yields were found to be the primary factors controlling the overall bimolecular association rate constants for CO binding to the mutants studied. The kinetic analysis provides estimates of 84%, 60%, 40%, 0%, and 99% for the steady-state hydrations of wild-type, H64Q, H64A, H64L, and V68F deoxymyoglobin, respectively. The second-order rate constants for CO and H(2)O entry into the empty distal pocket of myoglobin are markedly different, 8 x 10(7) and 2 x 10(5) M(-1).s(-1), respectively, suggesting that hydrophobic partitioning of the apolar gas from the aqueous phase into the relatively apolar protein interior lowers the free energy barrier for CO entry.

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.

Ligand migration pathway and protein dynamics in myoglobin: a time-resolved crystallographic study on L29W MbCO

By using time-resolved x-ray crystallography at room temperature, structural relaxations and ligand migration were examined in myoglobin (Mb) mutant L29W from nanoseconds to seconds after photodissociation of carbon monoxide (CO) from the heme iron by nanosecond laser pulses. The data were analyzed in terms of transient kinetics by fitting trial functions to integrated difference electron density values obtained from select structural moieties, thus allowing a quantitative description of the processes involved. The observed relaxations are linked to other investigations on protein dynamics. At the earliest times, the heme has already completely relaxed into its domed deoxy structure, and there is no photo-dissociated CO visible at the primary docking site. Initial relaxations of larger globin moieties are completed within several hundred nanoseconds. They influence the concomitant migration of photo-dissociated CO to the Xe1 site, where it appears at approximately 300 ns and leaves again at approximately 1.5 ms. The extremely long residence time in Xe1 as compared with wild-type MbCO implies that, in the latter protein, the CO exits the protein from Xe1 predominantly via the distal pocket. A well-defined deligated state is populated between approximately 2 micros and approximately 1 ms; its structure is very similar to the equilibrium deoxy structure. Between 1.5 and 20 ms, no CO is visible in the protein interior; it is either distributed among many sites within the protein or has escaped to the solvent. Finally, recombination at the heme iron occurs after >20 ms.

Increased peptide deformylase activity for N-formylmethionine processing of proteins overexpressed in Escherichia coli: application to homogeneous rubredoxin production

Deformylation of the initiator N-formylmethionine does not always proceed to completion for proteins overexpressed in Escherichia coli. To overcome this limitation, the def gene encoding the Escherichia coli peptide deformylase was cloned into the plysS plasmid under the tetracycline (Tc) promoter control. The efficiency of this constitutive level of peptide deformylase expression was demonstrated for the case of the rubredoxins from both mesophilic and hyperthermophilic organisms which normally retain a majority of their N-formyl terminal form. Indicating the potential structural/functional significance of residual formylation, the presence of a highly solvent exposed N-formyl group in rubredoxin is discernable in the amide NMR chemical shifts for the active site metal-coordinating cysteines more than 21A away.

An NMR study of the origin of dioxygen-induced spin-lattice relaxation enhancement and chemical shift perturbation

Due to its depth-dependent solubility, oxygen exerts paramagnetic effects which become progressively greater toward the hydrophobic interior of micelles, and lipid bilayer membranes. This paramagnetic gradient, which is manifested as contact shift perturbations (19F and 13C NMR) and spin-lattice relaxation enhancement (19F and 1H NMR), has been shown to be useful for precisely determining immersion depth, membrane protein secondary structure, and overall topology of membrane proteins. We have investigated the influence of oxygen on 19F and 13C NMR spectra and spin-lattice relaxation rates of a semiperfluorinated detergent, (8,8,8)-trifluoro (3,3,4,4,5,5,6,6,7,7)-difluoro octylmaltoside (TFOM) in a model membrane system, to determine the dominant paramagnetic spin-lattice relaxation and shift-perturbation mechanism. Based on the ratio of paramagnetic spin-lattice relaxation rates of 19F and directly bonded 13C nuclei, we conclude that the dominant relaxation mechanism must be dipolar. Furthermore, the temperature dependence of oxygen-induced chemical shift perturbations in 9F NMR spectra suggests a contact interaction is the dominant shift mechanism. The respective hyperfine coupling constants for 19F and 13C nuclei can then be estimated from the contact shifts <(deltav/v0)19F> and <(deltav/v0)13C>, allowing us to estimate the relative contribution of scalar and dipolar relaxation to 19F and 13C nuclei. We conclude that the contribution to spin-lattice relaxation from the oxygen induced paramagnetic scalar mechanism is negligible.

Mapping oxygen accessibility to ribonuclease a using high-resolution NMR relaxation spectroscopy

Paramagnetic contributions to nuclear magnetic spin-lattice relaxation rate constant induced by freely diffusing molecular oxygen measured at hundreds of different protein proton sites provide a direct means for characterizing the exploration of the protein by oxygen. This report focuses on regions of ribonuclease A where the rate constant enhancements are either quite large or quite small. We find that there are several regions of enhanced oxygen affinity for the protein both on the surface and in interior pockets where sufficient free volume permits. Oxygen has weak associative interactions with a number of surface crevices that are generally between secondary structural elements of the protein fold. Several regions near the surface have higher than expected accessibility to oxygen indicating that structural fluctuations in the protein provide intermolecular access. Oxygen penetrates part of the hydrophobic interior, but affinity does not correlate simply with hydrophobicity indices. Oxygen is excluded from regions of high interior packing density and a few surface sites where x-ray diffraction data have indicated the presence of specific hydration with high occupancy.

Competition with xenon elicits ligand migration and escape pathways in myoglobin

Evidence for ligand migration toward the xenon-binding cavities in myoglobin comes from a number of laser photolysis studies of MbO2 including mutants and from cryo- and time-resolved crystallography of MbCO. To explore ligand migration in greater detail, we investigated the rebinding kinetics of both MbO2 and MbCO under a xenon partial pressure ranging from 1 to 16 atm over the temperature range (293-77 K). Below 180 K xenon affects to a significant, but minor, extent the thermodynamic parameters for rebinding from the primary docking site in each Mb taxonomic substate. Above 200 K the ligand migrates to the proximal Xe1 site but when the latter is occupied by xenon a new kinetic process appears. It is attributed to rebinding from transient docking sites located on the path between the primary and the secondary docking site of both ligands. Ligand escape exhibits a more complicated pattern than expected. At room temperature O2 and CO escape appears to take place exclusively from the primary site. In contrast, at T approximately 250 K, roughly 50% of the CO molecules that have escaped from the protein originate from the Xe1 secondary site.