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: an infrared kinetic study of ligand migration in mutants YQR and YQRF

Recombination of carbon monoxide to myoglobin mutants YQR and YQRF was studied using transient infrared absorption spectroscopy and Fourier transform infrared-temperature derivative spectroscopy (FTIR-TDS). Photoproduct states B, C', C" and D associated with ligands residing in different protein cavities have been identified. After photolysis, ligands migrate to primary docking site B and subsequently rebind or escape to a secondary site (C) within the Xe4 cavity. For YQR, a global analysis of the isothermal rebinding kinetics below 160 K and the TDS data reveal a correlation between the enthalpy barriers governing the two processes. Above 120 K, a protein conformational change in both YQR and YQRF converts photoproduct C' into C" with markedly slowed kinetics. Above approximately 180 K, ligands migrate to the proximal Xe1 site (D) and also exit into the solvent, from where they rebind in a bimolecular reaction.     

Structural dynamics of myoglobin: FTIR-TDS study of NO migration and binding

By using Fourier transform infrared photolysis difference spectroscopy combined with temperature derivative spectroscopy at cryogenic temperatures, we have measured infrared spectra of the stretching absorption on nitric oxide (NO) in the heme-bound and photodissociated states of ferrous and ferric nitrosyl myoglobin (MbNO) and a few site-specific Mb mutants. The NO absorption was utilized as a sensitive local probe of ligand interactions with active-site residues and movements within the protein. By comparison with results obtained in previous spectroscopic and structural studies of carbonmonoxy myoglobin (MbCO), the MbNO data were interpreted in structural terms. In the NO-bound state, conformational heterogeneity was inferred from the appearance of multiple bands, arising from different electrostatic interactions with active site residues, most importantly, His-64. In ferrous MbNO, a primary photoproduct site similar to site B of MbCO was found, as indicated by a characteristic NO stretching spectrum. In ferric MbNO, the His-64 side chain appears to interfere with trapping of NO in this site; only a very weak photoproduct spectrum was observed in Mb variants in which His-64 was present. Upon extended illumination, the photoproduct spectrum changed in a characteristic way, indicating that NO readily migrates to a secondary docking site C, the Xe4 cavity, in which the ligand performs librational motions on the picosecond time scale. This docking site may play a role in the physiological NO scavenging reaction. Surprisingly, NO cannot be trapped at all in secondary docking site D, the Xe1 cavity.     

Coupling of protein relaxation to ligand binding and migration in myoglobin

Protein relaxation, ligand binding, and ligand migration into a hydrophobic cavity in myoglobin are unified by a bounded diffusion model which produces an accurate fit to complex ligand rebinding data over eight decades in time and a 160 K temperature range, in qualitative agreement with time-resolved x-ray crystallography. Protein relaxation operates in a cyclic manner to move the ligand away from the binding site.     

Structural dynamics of distal histidine replaced mutants of myoglobin accompanied with the photodissociation reaction of the ligand

Protein dynamics observed by the transient grating (TG) method are studied for some site-directed mutants at the distal histidine of myoglobin (H64L, H64Q, H64V). The time profiles of the TG signals are very sensitive to the amino acid residue of the 64 position. It was found that the sensitivity is mostly caused by the different rates of the ligand escape from the protein to solvent and the magnitude of the molecular volume change. Several molecular origins of the volume difference between MbCO and Mb, such as the electrostatic interaction in the distal pocket, movement of helices, and distal water, are proposed. Interestingly, the volume difference between the CO-trapped Mb inside the protein interior and Mb is similar to that of the partial molar volume of CO in organic solvent. The effect of mutation on the nature of the CO trapped site is discussed.     

Quantitative delineation of how breathing motions open ligand migration channels in myoglobin and its mutants

Ligand migration and binding are central to the biological functions of many proteins such as myoglobin (Mb) and it is widely thought that protein breathing motions open up ligand channels dynamically. However, how a protein exerts its control over the opening and closing of these channels through its intrinsic dynamics is not fully understood. Specifically, a quantitative delineation of the breathing motions that are needed to open ligand channels is lacking. In this work, we present and apply a novel normal mode‐based method to quantitatively delineate what and how breathing motions open ligand migration channels in Mb and its mutants. The motivation behind this work springs from the observation that normal mode motions are closely linked to the breathing motions that are thought to open ligand migration channels. In addition, the method provides a direct and detailed depiction of the motions of each and every residue that lines a channel and can identify key residues that play a dominating role in regulating the channel. The all‐atom model and the full force‐field employed in the method provide a realistic energetics on the work cost required to open a channel, and as a result, the method can be used to efficiently study the effects of mutations on ligand migration channels and on ligand entry rates. Our results on Mb and its mutants are in excellent agreement with MD simulation results and experimentally determined ligand entry rates. Proteins 2015; 83:757–770. © 2015 Wiley Periodicals, Inc.     

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 of myoglobin

Conformational fluctuations have been invoked to explain the observation that the diffusion of small ligands through a protein is a global phenomenon, as suggested (for example) by the oxygen induced fluorescence quenching of buried tryptophans. In enzymes processing large substrates, a channel to the catalytic site is often seen in the crystal structure; on the other hand in small globular proteins, it is not known if the cavities identified in the interior space are important in controlling their function by defining specific pathways in the diffusion to the active site. This point is addressed in this paper, which reports some relevant results obtained on myoglobin, the hydrogen atom of molecular biology. Protein conformational relaxations have been extensively investigated with myoglobin because the photosensivity of the adduct with CO, O2 and NO allows us to follow events related to the migration of the ligand through the matrix. Results obtained by laser photolysis, molecular dynamics simulations, X-ray diffraction of intermediate states of wt type and mutant myoglobins are briefly summarized. Crystallographic data on the photochemical intermediate of a new triple mutant of sperm whale myoglobin (Mb-YQR) show, for the first time, the photolyzed CO* sitting in one of the Xe-binding cavities, removed from the heme group. These results support the viewpoint that pre-existing 'packing defects' in the protein interior play a major role in controlling the dynamics of ligand binding, including oxygen, and thereby acquire a survival value.     

Dynamics of ligand binding to myoglobin.

Myoglobin rebinding of carbon monoxide and dioxygen after photodissociation has been observed in the temperature range between 40 and 350 K. A system was constructed that records the change in optical absorption at 436 nm smoothly and without break between 2 musec and 1 ksec. Four different rebinding processes have been found. Between 40 and 160 K, a single process is observed. It is not exponential in time, but approximately given by N(t) = (1 + t/to)-n, where to and n are temperature-dependent, ligand-concentration independent, parameters. At about 170 K, a second and at 200 K, a third concentration-independent process emerge. At 210 K, a concentration-dependent process sets in. If myoglobin is embedded in a solid, only the first three can be seen, and they are all nonexponential. In a liquid glycerol-water solvent, rebinding is exponential. To interpret the data, a model is proposed in which the ligand molecule, on its way from the solvent to the binding site at the ferrous heme iron, encounters four barriers in succession. The barriers are tentatively identified with known features of myoglobin. By computer-solving the differential equation for the motion of a ligand molecule over four barriers, the rates for all important steps are obtained. The temperature dependences of the rates yield enthalpy, entropy, and free-energy changes at all barriers. The free-energy barriers at 310 K indicate how myoglobin achieves specificity and order. For carbon monoxide, the heights of these barriers increase toward the inside; carbon monoxide consequently is partially rejected at each of the four barriers. Dioxygen, in contrast, sees barriers of about equal height and moves smoothly toward the binding site. The entropy increases over the first two barriers, indicating a rupturing of bonds or displacement of residues, and then smoothly decreases, reaching a minimum at the binding site. The magnitude of the decrease over the innermost barrier implies participation of heme and/or protein. The nonexponential rebinding observed at low temperatures and in solid samples implies that the innermost barrier has a spectrum of activation energies. The shape of the spectrum has been determined; its existence can be explained by assuming the presence of many conformational states for myoglobin. In a liquid at temperatures above about 230 K, relaxation among conformational states occurs and rebinding becomes exponential.     

Trapping intermediates in the crystal: ligand binding to myoglobin

Crystal structures of the reactive short-lived species that occur in chemical or binding reactions can be determined using X-ray crystallography via time-resolved or kinetic trapping approaches. Recently, various kinetic trapping methods have been used to determine the structure of intermediates in ligand binding to myoglobin.     

Structure and ligand binding properties of leucine 29(B10) mutants of human myoglobin.

Site-specific mutants of human myoglobin (Mb) have been prepared, in which Leu29 (B10) is replaced by Ala(L29A) or Ile(L29I), in order to examine the influence of this highly conserved residue in the hydrophobic clusters of the heme distal site on the heme environmental structure and ligand binding properties of Mb. Structural characterizations of these recombinant Mbs are studied by electronic absorption, infrared (IR), one- and two-dimensional proton nuclear magnetic resonance spectroscopies, and ligand-binding kinetics by laser photolysis measurements under ambient and high pressures (up to 2000 bar). Multiple split carbon monoxide (CO) stretch bands in the IR spectra of mutant Mbs exhibit a relative decrease of the 1945 cm-1 band (approximately 50%) which is associated with an upright binding geometry of CO, accompanied by an increase of the tilted CO conformer at 1932 cm-1. On the basis of these results, replacement of Leu29(B10) by Ala or Ile appears to allow bound CO to rotate from a conformation pointing toward the beta meso carbon of the heme group to the one pointing toward the alpha meso carbon atom, presumably filling the space left by removal of the delta 2 carbon atom of Leu29(B10). These substitutions cause the rate constants for CO and O2 association to decrease almost 3-5-fold. Present results show that CO and O2 bindings to the heme iron of Mb are controlled by Leu29(B10) by influencing the structure of close vicinity of the heme and the geometry of iron-bound ligand. Further, mutant Mbs (Leu72(E15)----Ala and Leu104 (G5)----Ala) which have altered residues in another hydrophobic clusters around proximal and distal site are also examined.     

Structural dynamics of liganded myoglobin.

X-ray crystallography can reveal the magnitudes and principal directions of the mean-square displacements of every atom in a protein. This structural information is complementary to the temporal information obtainable by spectroscopic techniques such as nuclear magnetic resonance. Determination of the temperature dependence of the mean-square displacements makes it possible to separate large conformational motions from simple thermal vibrations. The contribution of crystal lattice disorder to the overall apparent displacement can be estimated by Mössbauer spectroscopy. This technique has been applied to high resolution x-ray diffraction data from sperm whale myoglobin in its Met iron and oxy cobalt forms. Both crystal structures display regions of large conformational motions, particularly at the chain termini and in the region of the proximal histidine. Overall, the mean-square displacement increases with increasing distance from the center of gravity of the molecule. Some regions of the heme pocket in oxy cobalt myoglobin are more rigid than the corresponding regions in Met myoglobin.     

Controlling ligand binding in myoglobin by mutagenesis.

A quadruple mutant of sperm whale myoglobin was constructed to mimic the structure found in Ascaris suum hemoglobin. The replacements include His(E7)-->Gln, Leu(B10)-->Tyr, Thr(E10)--> Arg, and Ile(G8)-->Phe. Single, double, and triple mutants were characterized to dissect out the effects of the individual substitutions. The crystal structures of the deoxy and oxy forms of the quadruple mutant were determined and compared with that of native Ascaris hemoglobin. Tyr(B10) myoglobin displays low O(2) affinity, high dissociation rate constants, and heterogeneous kinetic behavior, suggesting unfavorable steric interactions between the B10 phenol side chain and His(E7). In contrast, all mutants containing the Tyr(B10)/Gln(E7) pair show high O(2) affinity, low dissociation rate constants, and simple, monophasic kinetic behavior. Replacement of Ile(107) with Phe enhances nanosecond geminate recombination singly and in combination with the Tyr(B10)/Gln(E7)/Arg(E10) mutation by limiting access to the Xe4 site. These kinetic results and comparisons with native Ascaris hemoglobin demonstrate the importance of distal pocket cavities in governing the kinetics of ligand binding. The approximately 150-fold higher O(2) affinity of Ascaris hemoglobin compared with that for Tyr(B10)/Gln(E7)-containing myoglobin mutants appears to be the result of favorable proximal effects in the Ascaris protein, due to a staggered orientation of His(F8), the lack of a hydrogen bonding lattice between the F4, F7, and F8 residues, and the presence of a large polar Trp(G5) residue in the interior portion of the proximal heme pocket.     

On the mechanism of ligand binding to myoglobin

The association reaction of CO and O₂ with heme is expected to reflect the differences in the electronic structures of the two ligands. CO binding should be controlled by a high spin/low spin transition while oxygen binding is spin-allowed. Dioxygen should thus bind substantially faster than CO. The experimental association rates of the two ligands are, however, almost identical. We propose that the reaction is triggered in both cases by a fast structural intermediate which allows the CO molecule to bind adiabatically. A suitable structural transition has been identified recently by inelastic neutron scattering.     

Ligand and proton exchange dynamics in recombinant human myoglobin mutants

Site-specific mutants of human myoglobin have been prepared in which lysine 45 is replaced by arginine (K45R) and aspartate 60 by glutamate (D60E), in order to examine the influence of these residues and their interaction on the dynamics of the protein. These proteins were studied by a variety of methods, including one and two-dimensional proton nuclear magnetic resonance spectroscopy, exchange kinetics for the distal and proximal histidine NH protons as a function of pH in the met cyano forms, flash photolysis of the CO forms, and ligand replacement kinetics. The electronic absorption and proton nuclear magnetic resonance spectra of the CO forms of these proteins are virtually identical, indicating that the structure of the heme pocket is unaltered by these mutations. There are, however, substantial changes in the dynamics of both CO binding and proton exchange for the mutant K45R, whereas the mutant D60E exhibits behavior indistinguishable from the reference human myoglobin. K45R has a faster CO bimolecular recombination rate and slower CO off-rate relative to the reference. The kinetics for CO binding are independent of pH (6.5 to 10) as well as ionic strength (0 to 1 M-NaCl). The exchange rate for the distal histidine NH is substantially lower for K45R than the reference, whereas the proximal histidine NH exchange rate is unaltered. The exchange behavior of the human proteins is similar to that reported for a comparison of the exchange rates for myoglobins having lysine at position 45 with sperm whale myoglobin, which has arginine at this position. This indicates that the differences in exchange rates reflects largely the Lys----Arg substitution. The lack of a simple correlation for the CO kinetics with this substitution means that these are sensitive to other factors as well. Specific kinetic models, whereby substitution of arginine for lysine at position 45 can affect ligand binding dynamics, are outlined. These experiments demonstrate that a relatively conservative change of a surface residue can substantially perturb ligand and proton exchange dynamics in a manner that is not readily predicted from the static structures.     

Distal pocket polarity in ligand binding to myoglobin: structural and functional characterization of a threonine68(E11) mutant

Site-directed mutagenesis studies have confirmed that the distal histidine in myoglobin stabilizes bound O2 by hydrogen bonding and have suggested that it is the polar character of the imidazole side chain rather than its size that limits the rate of ligand entry into the protein. We constructed an isosteric Val68 to Thr replacement in pig myoglobin (i) to investigate whether the O2 affinity could be increased by the introduction of a second hydrogen-bonding group into the distal heme pocket and (ii) to examine the influence of polarity on the ligand binding rates more rigorously. The 1.9-A crystal structure of Thr68 aquometmyoglobin confirms that the mutant and wild-type proteins are essentially isostructural and reveals that the beta-OH group of Thr68 is in a position to form hydrogen-bonding interactions both with the coordinated water molecule and with the main chain greater than C=O of residue 64. The rate of azide binding to the ferric form of the Thr68 mutant was 60-fold lower than that for the wild-type protein, consistent with the proposed stabilization of the coordinated water molecule. However, bound O2 is destabilized in the ferrous form of the mutant protein. The observed 17-fold lowering of the O2 affinity may be a consequence of the hydrogen-bonding interaction made between the Thr68 beta-OH group and the carbonyl oxygen of residue 64. Overall association rate constants for O2, NO, and alkyl isocyanide binding to ferrous pig myoglobin were 3-10-fold lower for the mutant compared to the wild-type protein, whereas that for CO binding was little affected.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Quantifying the importance of protein conformation on ligand migration in myoglobin

Myoglobin (Mb) is a model system for ligand binding and migration. The energy barriers (ΔG) for ligand migration in Mb have been studied in the past by experiment and theory and significant differences between different approaches were found. From experiment, it is known that Mb can assume a large number of conformational substates. In this work, these substates are investigated as a possible source of the differences in migration barriers. We show that the initial structure significantly affects the calculated ΔG for a particular transition and that fluctuations in barrier heights δΔG are of similar magnitude as the free energy barriers themselves. The sensitivity of ΔG to the initial structure is compared to other sources of errors. Different protein structures can affect the calculated ΔG by up to 4 kcal/mol, whereas differences between simple point charge models and more elaborate multipolar charge models for the CO-ligand are smaller by a factor of two. Analysis of the structural changes underlying the large effect of the conformational substate reveals the importance of coupling between protein and ligand motion for migration.