Infrared Study of Carbon Monoxide Migration among Internal Cavities of Myoglobin Mutant L29W

Myoglobin, a small globular heme protein that binds gaseous ligands such asO₂, CO and NO reversibly at the heme iron, provides an excellent modelsystem for studying structural and dynamic aspects of protein reactions. Flashphotolysis experiments, performed over wide ranges in time and temperature, reveal a complex ligand binding reaction with multiple kinetic intermediates, resulting from protein relaxation and movements of the ligand within the protein. Our recent studies of carbonmonoxy-myoglobin (MbCO) mutant L29W, using time-resolved infrared spectroscopy in combination with x-ray crystallography, have correlated kinetic intermediates with photoproduct structures that are characterized by the CO residing in different internal protein cavities, so-called xenon holes. Here we have used Fourier transform infrared temperature derivative spectroscopy (FTIR-TDS) to further examine the role of internal cavities in the dynamics. Different cavities can be accessed by the CO ligands at different temperatures, and characteristic infrared absorption spectra have been obtained for the different locations of the CO ligand within the protein, enabling us to monitor ligand migration through the protein as well as conformational changes of the protein.     

Protein Collective Motions Coupled to Ligand Migration in Myoglobin

Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.     

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.     

Unraveling the molecular basis for ligand binding in truncated hemoglobins: The trHbO Bacillus subtilis case

Truncated hemoglobins (trHbs) are heme proteins present in bacteria, unicellular eukaryotes, and higher plants. Their tertiary structure consists in a 2‐over‐2 helical sandwich, which display typically an inner tunnel/cavity system for ligand migration and/or storage. The microorganism Bacillus subtilis contains a peculiar trHb, which does not show an evident tunnel/cavity system connecting the protein active site with the solvent, and exhibits anyway a very high oxygen association rate. Moreover, resonant Raman results of CO bound protein, showed that a complex hydrogen bond network exists in the distal cavity, making it difficult to assign unambiguously the residues involved in the stabilization of the bound ligand. To understand these experimental results with atomistic detail, we performed classical molecular dynamics simulations of the oxy, carboxy, and deoxy proteins. The free energy profiles for ligand migration suggest that there is a key residue, GlnE11, that presents an alternate conformation, in which a wide ligand migration tunnel is formed, consistently with the kinetic data. This tunnel is topologically related to the one found in group I trHbs. On the other hand, the results for the CO and O₂ bound protein show that GlnE11 is directly involved in the stabilization of the cordinated ligand, playing a similar role as TyrB10 and TrpG8 in other trHbs. Our results not only reconcile the structural data with the kinetic information, but also provide additional insight into the general behaviour of trHbs. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

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.     

The effect of ligand dynamics on heme electronic transition band III in myoglobin.

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.     

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.     

Different roles of protein dynamics and ligand migration in non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana

The ligand rebinding kinetics after photolysis of the CO complexes of Arabidopsis thaliana hemoglobins AHb1 and AHb2 in solution show very different amplitudes in the geminate phase, reflecting different migration pathways of the photodissociated ligand in the system of internal cavities accessible from the heme. The dependence of the geminate phase on CO concentration, temperature, encapsulation in silica gels and presence of glycerol confirms a remarkable difference in the internal structure of the two proteins and a dramatically different role of protein dynamics in regulating the reactivity with CO. This finding strongly supports the idea that they have distinct physiological functions.     

Nitric Oxide Dynamics in Truncated Hemoglobin: Docking Sites, Migration Pathways, and Vibrational Spectroscopy from Molecular Dynamics Simulations

Atomistic simulations of nitric oxide (NO) dynamics and migration in the trHbN of Mycobacterium tuberculosis are reported. From extensive molecular dynamics simulations (48 ns in total), the structural and energetic properties of the ligand docking sites in the protein have been characterized and a connectivity network between the ligand docking sites has been built. Several novel migration and exit pathways are found and are analyzed in detail. The interplay between a hydrogen-bonding network involving residues Tyr³³ and Gln⁵⁸ and the bound O₂ ligand is discussed and the role of Phe⁶² residue in ligand migration is examined. It is found that Phe⁶² is directly involved in controlling ligand migration. This is reminiscent of His⁶⁴ in myoglobin, which also plays a central role in CO migration pathways. Finally, infrared spectra of the NO molecule in different ligand docking sites of the protein are calculated. The pocket-specific spectra are typically blue-shifted by 5-10 cm⁻¹, which should be detectable in future spectroscopic experiments.     

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.     

Photothermal studies of CO photodissociation from mixed valence Escherichia coli cytochrome bo3

Here, we report the volume and enthalpy changes accompanying CO photodissociation from the mixed valence form of cytochrome bo3 oxidase from Escherichia coli. The results of photoacoustic calorimetry indicate two kinetic phases with distinct volume and enthalpy changes accompanying CO photodissociation from heme o3 and its transfer to CuB. The first phase occurring on a timescale of <50 ns is characterized by a volume decrease of -1.3+/-0.3 mL mol-1 and enthalpy change of 32+/-1.6 kcal mol-1. Subsequently, a volume increase of 2.9 mL mol-1 with an enthalpy change of -5.3+/-2.5 kcal mol-1 is observed with the lifetime of approximately 250 ns (this phase has not been detected in previous optical studies). These volume and enthalpy changes differ from the volume and enthalpy changes observed for CO dissociation from fully reduced cytochrome bo3 oxidase indicating that the heme o3/CuB active site dynamics are affected by the redox state of heme b.     

Low-temperature flash photolysis studies of cytochrome oxidase and its environment.

The CO-binding kinetics of cytochrome a3, in isolated, detergent-solubilized cytochrome oxidase have been studied by flash photolysis over wide ranges of CO concentration and temperature. The results strongly suggest that CO has an intermediate bound state in its path to the final bound state at the heme iron. In the temperture range 230-273 K in frozen aqueous solutions, the recombination rates depend upon CO concentration; at low CO concentrations the kinetics are biphasic. The rate of the faster process depends upon the detergent concentration, that of the slower process upon the salt concentration. In addition, the faster process depends upon the amount of CO photodissociated. It is concluded that the cytochrome oxidase molecules are aggregated in regions that contain detergent and possibly some lipids. The regions retain considerable fluid character well below the macroscopic freezing point of the solution. The faster phase of the recombination is interpreted as the rebinding of CO molecules that remain in the fluid region after photodissociation. The slower phase would then be due to the migration of some dissociated CO out into surrounding frozen solvent. The non-Arrhenius behavior of both phase probably represents partial melting of the medium; preliminary NMR measurements of mobile protons support this hypothesis. Many of the kinetic features described here are also seen in mitochondria; thus the detergent-solubilized cytochrome oxidase may be a useful model system for the ligand-binding behavior of the enzyme in the mitochondrial membrane.     

Quantitative analysis of ligand migration from transition networks

In this work we use transition network analysis for the first time to investigate ligand migration in truncated hemoglobin (trHbN) and obtain kinetic information about the docking-site dynamics in the protein. A comparison with explicit water molecular dynamics simulations (100 ns in total) shows that the rate constants derived from the network analysis are realistic. The transition network analysis provides 1) The time-resolved connectivity network in the protein; 2) The half-lives of the docking sites; 3) The transition timescales between two given docking sites; and 4) The extent of population transfer among different docking sites of the protein as a function of lag time. We investigate the role of the Tyr³³ and Gln⁵⁸ residues in ligand migration by studying ligand migration in four mutants of trHbN. The mutation study suggests that residues Tyr³³ and Gln⁵⁸ stabilize the NO ligand in the Xe2 docking site of trHbN, thus facilitating the efficiency of the NO detoxification reaction.     

Rapid kinetic studies of acetyl-CoA synthesis: evidence supporting the catalytic intermediacy of a paramagnetic NiFeC species in the autotrophic Wood-Ljungdahl pathway

CO dehydrogenase/acetyl-CoA synthase (CODH/ACS), a key enzyme in the Wood-Ljungdahl pathway of anaerobic CO(2) fixation, is a bifunctional enzyme containing CODH, which catalyzes the reversible two-electron oxidation of CO to CO(2), and ACS, which catalyzes acetyl-CoA synthesis from CoA, CO, and a methylated corrinoid iron-sulfur protein (CFeSP). ACS contains an active site nickel iron-sulfur cluster that forms a paramagnetic adduct with CO, called the nickel iron carbon (NiFeC) species, which we have hypothesized to be a key intermediate in acetyl-CoA synthesis. This hypothesis has been controversial. Here we report the results of steady-state kinetic experiments; stopped-flow and rapid freeze-quench transient kinetic studies; and kinetic simulations that directly test this hypothesis. Our results show that formation of the NiFeC intermediate occurs at approximately the same rate as, and its decay occurs 6-fold faster than, the rate of acetyl-CoA synthesis. Kinetic simulations of the steady-state and transient kinetic results accommodate the NiFeC species in the mechanism and define the rate constants for the elementary steps in acetyl-CoA synthesis. The combined results strongly support the kinetic competence of the NiFeC species in the Wood-Ljungdahl pathway. The results also imply that the methylation of ACS occurs by attack of the Ni(1+) site in the NiFeC intermediate on the methyl group of the methylated CFeSP. Our results indicate that CO inhibits acetyl-CoA synthesis by inhibiting this methyl transfer reaction. Under noninhibitory CO concentrations (below 100 microM), formation of the NiFeC species is rate-limiting, while at higher inhibitory CO concentrations, methyl transfer to ACS becomes rate-limiting.     

Ligand migration and protein fluctuations in myoglobin mutant L29W

We have determined eight X-ray structures of myoglobin mutant L29W at various experimental conditions. In addition, infrared spectroscopic experiments are presented, which are discussed in the light of the X-ray structures. Two distinct conformations of the CO-ligated protein were identified, giving rise to two stretching bands of heme-bound CO. If L29W MbCO crystals are illuminated around 180 K, a deoxy species is formed. The CO molecules migrate to the proximal side of the heme and remain trapped in the so-called Xe1 cavity upon temperature decrease to 105 K. The structure of this photoproduct is almost identical to the equilibrium high-temperature deoxy Mb structure. If the temperature is cycled to increasingly higher values, CO recombination is observed. Three intermediate structures have been determined during the rebinding process. Efficient recombination occurs only above 180 K, the characteristic temperature for the onset of protein dynamics. Rebinding is remarkably slow because bulky residues His64 and Trp29 block important migration pathways of the CO molecule.     

Quantitative studies of bacterial chemotaxis and microbial population dynamics

Although there is a long history of conjecture regarding the role and significance of bacterial chemotaxis in microbial ecology, only recently has a significant body of work appeared attempting to address this issue. The purpose of this paper is to provide a concise overview of this work, which combined mathematical modeling of bacterial population migration and experimental measurement of the model parameters with modeling of competitive microbial population dynamics in a nonmixed environment. Predictions from the population dynamics models, based on experimental estimates of the various motility and growth parameter values, are related to the small number of experimental observations available to date dealing with the effects of bacterial motility on competition in a nonmixed environment. Current results indicate that cell motility and chemotaxis properties can be as important to population dynamics as cell growth kinetic properties, so that greater attention to this aspect of microbial behavior is warranted in future studies of microbial ecology.     

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 MbO₂ 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 MbO₂ 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 O₂ and CO escape appears to take place exclusively from the primary site. In contrast, at T ≈ 250 K, roughly 50% of the CO molecules that have escaped from the protein originate from the Xe1 secondary site.