Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

X-ray structure determination of a metastable state of carbonmonoxy myoglobin after photodissociation.

The x-ray structure of carbon monoxide (CO)-ligated myoglobin illuminated during data collection by a laser diode at the wavelength lambda = 690 nm has been determined to a resolution of 1.7 A at T = 36 K. For comparison, we also measured data sets of deoxymyoglobin and CO-ligated myoglobin. In the photon-induced structure the electron density associated with the CO ligand can be described by a tube extending from the iron into the heme pocket over more than 4 A. This density can be interpreted by two discrete positions of the CO molecule. One is close to the heme iron and can be identified to be bound CO. In the second, the CO is dissociated from the heme iron and lies on top of pyrrole ring C. At our experimental conditions the overall structure of myoglobin in the metastable state is close to the structure of a CO-ligated molecule. However, the iron has essentially relaxed into the position of deoxymyoglobin. We compare our results with those of Schlichting el al. [Schlichting, I., Berendzen, J., Phillips, G. N., Jr., & Sweet, R. M. (1994) Nature 317, 808-812], who worked with the myoglobin mutant (D122N) that crystallizes in the space group P6 and Teng et al. [Teng, T. Y., Srajer, V. & Moffat, K. (1994) Nat. Struct. Biol. 1, 701-705], who used native myoglobin crystals of the space group P2(1). Possible reasons for the structural differences are discussed.

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Mapping the Intramolecular Vibrational Energy Flow in Proteins Reveals Functionally Important Residues

Unveiling the mechanisms of energy relaxation in biomolecules is key to our understanding of protein stability, allostery, intramolecular signaling, and long-lasting quantum coherence phenomena at ambient temperatures. Yet, the relationship between the pathways of energy transfer and the functional role of the residues involved remains largely unknown. Here, we develop a simulation method of mapping out residues that are highly efficient in relaxing an initially localized excess vibrational energy and perform site-directed mutagenesis functional assays to assess the relevance of these residues to protein function. We use the ligand binding domains of thyroid hormone receptor (TR) subtypes as a test case and find that conserved arginines, which are critical to TR transactivation function, are the most effective heat diffusers across the protein structure. These results suggest a hitherto unsuspected connection between a residue`s ability to mediate intramolecular vibrational energy redistribution and its functional relevance.

Carbohydrate supplementation and alterations in neutrophils, and plasma cortisol and myoglobin concentration after intense exercise

The present study examined the effect of carbohydrate supplementation on changes in neutrophil counts, and the plasma concentrations of cortisol and myoglobin after intense exercise. Eight well-trained male runners ran on a treadmill for 1 h at 85% maximal oxygen uptake on two separate occasions. In a double-blind cross-over design, subjects consumed either 750 ml of a 10% carbohydrate (CHO) drink or a placebo drink on each occasion. The order of the trials was counter-balanced. Blood was drawn immediately before and after exercise, and 1 h after exercise. Immediately after exercise, neutrophil counts (CHO, 49%; placebo, 65%; P<0.05), plasma concentrations of glucose (CHO, 43%; P<0.05), lactate (CHO, 130%; placebo, 130%; P<0.01), cortisol (CHO, 100%; placebo, 161%; P<0.01), myoglobin (CHO, 194%; placebo, 342%; P<0.01) all increased significantly. One hour post-exercise, plasma myoglobin concentration (CHO, 331%; placebo, 482%; P<0.01) and neutrophil count (CHO, 151%; placebo, 230% P<0.01) both increased further above baseline. CHO significantly attenuated plasma myoglobin concentration and the neutrophil count after exercise (P<0.01), but did not affect plasma cortisol concentration. The effects of CHO on plasma myoglobin concentration may be due to alterations in cytokine synthesis, insulin responses or myoglobin clearance rates from the bloodstream during exercise. Plasma cortisol responses to CHO during exercise may depend on the intensity of exercise, or the amount of CHO consumed. Lastly, cortisol appears to play a minor role in the mobilisation of neutrophils after intense exercise.

Action of N-bromosuccinimide on sperm whale myoglobin

Previously, it was reported from this laboratory that the heme groups of hemoglobin are “buried” within globin at pH 4.0 and not dissociated, on the basis of the obiligatory requirement of urea for the reaction of N-bromosuccinimide with the heme groups of hemoglobin at pH4.0, and also on the basis of the “normalization” of the spectrum of hemoglobin at this pH in the presence of urea or sucrose. In the present study, it has been shown that the behaviour of sperm whale myoglobin with respect to its reaction with N-bromosuccinimide and with respect to spectral “normalization” in urea or sucrose are essentially similar to that of hemoglobin. It has also been demonstrated that the spectral “normalization” obtained with crystalline hemin is not identical with that obtained with either hemoglobin or myoglobin. The bearing of the results of the present study on the earlier work on hemoglobin is indicated.

Employing denaturation for rapid electrochemical detection of myoglobin using TiO2 nanotubes

An alternative antibody-free strategy for the rapid electrochemical detection of cardiac myoglobin has been demonstrated here using hydrothermally synthesized TiO2 nanotubes (Ti-NT). The denaturant induced unfolding of myoglobin led to easy access of the deeply buried electroactive heme center and thus the efficient reversible electron transfer from protein to electrode surface. The sensing performance of the Ti-NT modified electrodes were compared vis a vis commercially available titania and GCEs. The tubular morphology of the Ti-NT led to facile transfer of electrons to the electrode surface, which eventually provided a linear current response (obtained from cyclic voltammetry) over a wide range of Mb concentration. The sensitivity of the Ti-NT based sensor was remarkable and was equal to 18 mu A mg(-1) ml (detection limit = 50 nM). This coupled with the rapid analysis time of a few tens of minutes (compared to a few days for ELISA) demonstrates its potential usefulness for the early detection of acute myocardial infarction (AMI).

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.

Direct electrochemistry and surface plasmon resonance characterization of alternate layer-by-layer self-assembled DNA-myoglobin thin films on chemically modified gold surfaces

Alternate layer-by-layer (L-by-L) polyion adsorption onto gold electrodes coated with chemisorbed cysteamine gave stable, electroactive multilayer films containing calf thymus double stranded DNA (CT ds-DNA) and myoglobin (Mb). Direct, quasi-reversible electron exchange between gold electrodes and proteins involved the Mb heme Fe2+/Fe3+ redox couple. The formation of L-by-L (DNA/Mb), films was characterized by both in situ surface plasmon resonance (SPR) monitoring and cyclic voltammetry (CV). The effective thickness of DNA and Mb monolayers in the (DNA/Mb)l bilayer were 1.0 +/- 0.1 and 2.5 +/- 0.1 mn, corresponding to the surface coverage of similar to65% and similar to89% of its full packed monolayer, respectively. A linear increase of film thickness with increasing number of layers was confirmed by SPR characterizations. At pH 5.5, the electroactive Mb in films are those closest to the electrode surface; additional protein layers did not communicate with the electrode. CV studies showed that electrical communication might occur through hopping conduction via the electrode/base pair/Mb channel, thanks to the DNA-Mb interaction. After the uptake of Zn2+, a special electrochemical behavior, where MbFe(2+) acts as a DNA-binding reduction catalyst in the Zn2+-DNA/Mb assembly, takes place.

Application of resonance raman spectroscopy to study the effect of rare-earth ion Er3+ on myoglobin

The effect of rare-earth ion Er3+ On myoglobin(Mb) was studied by using Resonance Raman spectroscopy. The results show that with the variation of Er3+ concentrations, both the oxidation state and spin state of Mb are sensitive to the perturbation of Er3+. Er3+ added to Mb affects the oxidation and spin state synchronously. The structure-sensitive groups of Mb are more accessible to the Er3+ than other groups. According to the fluorometry and CD spectra studied and our results as mentioned above, we considered that Er3+ does not interact with heme directly, and Er3+ probably leads to the conformational changes of Mb due to the change of oxidation and spin state of Heme.

Synchronous fluorescence spectra of myoglobin

The synchronous fluorescence spectra of myoglobin were studies for the first time. The fluorescence peals observed in the spectra were assigned, When the wavelength interval (Delta lambda) is 80 nm, the main peak at 335 nm is originated from the tryptophan residues in the myoglobin molecule. When Delta lambda is 20 nn, the peak at 308 nm is mainly due to the tyrosine residues in the myoglobin molecule and in a small part due to the tryptophan residues. Two peaks at 322 and 596 nm were observed in the spectrum of myoglobin for Delta lambda = 40 nm. The peak at 322 nm is due to both tyrosine and tryptophan residues. The peak at 596 nm is attributed to the heme group in the myoglobin molecule.

The direct electrochemistry of cryo-hydrogel immobilized myoglobin at a glassy carbon electrode

A cryo-hydrogel membrane (CHM) immobilized at a glassy carbon (GC) electrode is reported for the direct electron transfer of redox proteins. The most attractive characteristics of this CHM were its hydrophilic micro-environment for incorporated proteins to retain their activities, its high ability for protection against interference of denatured and adsorbed proteins at the electrode, its potential applications for various proteins or enzymes, as well as its high mechanical strength and thermal stability. A clear well developed and stable redox wave was obtained for commercially available horse heart myoglobin without further purification, giving a peak to peak separation Delta E(p) = 93 mV at 5 mV s(-1) and the formal electrode potential E(0)' = -0.158 V (vs. Ag/AgCl). The formal heterogeneous electron transfer rate constant was calculated as k(0)' = 5.7 X 10(-4) cm s(-1) at pH 6.5, showing rapid electron transfer was achieved. The pH controlled conformational equilibria, acid state --> natural state --> basic I state --> basic II state, of myoglobin at the CHM GC electrode in the pH range 0-13.8 were also observed and are discussed in detail.

A cryo-hydrogel immobilized protein electrode for direct electron transfer of myoglobin

The direct electron transfer process of horse heart myoglobin, which was immobilized into a new type of cryo-hydrogel membrane on a glassy carbon electrode surface, was studied and the characteristics of this cryo-hydrogel immobilized protein electrode were discussed.