Theoretical characterization of carbon monoxide vibrational spectrum in sperm whale myoglobin distal pocket

In this article we use the perturbed matrix method and an extended molecular dynamics sampling of the carbon monoxide (CO) in the myoglobin distal pocket to characterize the CO vibrational spectrum and hence to relate its spectroscopic features with the atomic-molecular behavior. Results show the accuracy of the method employed and confirm the assignment of the spectroscopic B1 and B2 states proposed by Lim et al.     

Evidence of quasi-linear gas transport through sperm whale myoglobin

The diffusion of molecular oxygen or its isosteric analogue, carbon monoxide, from the surface of myoglobin to its deeply imbedded haem appears to represent one of the simplest protein functions. Hence, it was chosen for the study of the possible role of a global controlling effect like an attractor. However, whereas the six statistical criteria of the classical non-linear dynamic analysis for the existence of an attractor in myoglobin were fulfilled and invariant to the Fourier transformation, the properties of this attractor were not as simple as anticipated. The parameters were tested and confirmed by alternative approaches, the interpoint distance method of Judd and Fourier transformation. If the diffusion were approximately linear, the order of the attractor would be expected to be near one. However, a clearly higher value, 1.46+/-0.03, was found, indicating the existence of additional steps. Later, the latter were identified as a 90 degrees rotation of CO followed by a translocation by 0.4 A to a transient pocket. These additional steps may explain the high number of regulatory factors found, 10+/-1. The autocorrelation function was damped with a correlation length of at least 20 residues. The Poincaré plot showed a dense domain compatible with the cross-section of a quasi-spherical attractor. The first Lyapunov exponent, lambda(1), was clearly positive. The Hurst fractal coefficient was 1.90+/-0.22, indicating a clear departure from simple linear diffusion.     

Functional effects of heme orientational disorder in sperm whale myoglobin

The optical absorption and ligand binding properties of newly reconstituted sperm whale myoglobin were examined systematically at pH 8, 20 degrees C. The conventional absorbance and magnetic circular dichroism spectra of freshly reconstituted samples were identical to those of the native protein. In contrast, reconstituted azide or CO myoglobin initially exhibited less circular dichroism in the Soret wavelength region than native myoglobin. These data support the theory proposed by La Mar and co-workers (La Mar, G. N., Davis, N. L., Parish, D. W., and Smith, R. M. (1983) J. Mol. Biol. 168, 887-896) that protoheme inserts into apomyoglobin in two distinct orientations. The equilibrium and kinetic parameters for O2 and CO binding to newly reconstituted myoglobin were observed to be identical to those of the native protein. Thus, the orientation of the heme group has no effect on the physiological properties of myoglobin. This result is in disagreement with the preliminary report of Livingston et al. (Livingston, D. J., Davis, N. L., La Mar, G. N., and Brown, W. D. (1984) J. Am. Chem. Soc. 106, 3025-3026) which suggested that the abnormal heme conformation exhibited a 10-fold greater affinity and association rate constant for O2 binding. Significant kinetic heterogeneity was observed only for long-chain isonitrile binding to newly reconstituted myoglobin, and even in these cases, the rate constants for the abnormal and normal heme conformations differed by less than a factor of 4.     

pH-dependent stability of sperm whale myoglobin in water-guanidine hydrochloride solutions

An experimental-theoretical approach for the elucidation of protein stability is proposed. The theoretical prediction of pH-dependent protein stability is based on the macroscopic electrostatic model for calculation of the pH-dependent electrostatic free energy of proteins. As a test of the method we have considered the pH-dependent stability of sperm whale metmyoglobin. Two theoretical methods for evaluation of the electrostatic free energy and p K values are applied: the finite-difference Poisson-Boltzmann method and the semiempirical approach based on the modified Tanford-Kirkwood theory. The theoretical results for electrostatic free energy of unfolding are compared with the experimental data for guanidine hydrochloride unfolding under equilibrium conditions over a wide pH range. Using the optical parameters of the Soret absorbance to monitor conformational equilibrium and Tanford's method to estimate the resulting data, it was found that the conformational free energy of unfolding of metmyoglobin is 16.3 kcal mol(-1) at neutral pH values. The total unfolding free energies were calculated on the basis of the theoretically predicted electrostatic unfolding free energies and the experimentally measured midpoints (pH(1/2)) of acidic and alkaline denaturation transitions. Experimental data for alkaline denaturation were used for the first time in theoretical analysis of the pH-dependent unfolding of myoglobin. The present results demonstrate that the simultaneous application of appropriate theoretical and experimental methods permits a more complete analysis of the pH-dependent and pH-independent properties and stability of globular proteins.     

pH-dependence of photo-induced electron transfer in zinc-substituted sperm whale myoglobin

The electron transfer between the excited triplet state of zinc-substituted sperm whale myoglobin and Cu2+ has been studied by following the decay rate of delayed fluorescence. The Cu2+ bound on the surface of the myoglobin molecule are efficient quenchers of the excited electron state of Zn-myoglobin. Two bimolecular rate constants of quenching (KQ) for every pH investigated have been calculated. The pH-dependence of KQ1 indicates that the protonation of one amino acid residue (His-GH1 (119] is important for the process. Our results support the idea of the common nature of the mechanism of quenching by Cu2+ and oxidation of oxymyoglobin by Cu2+.     

A kinetic description of ligand binding to sperm whale myoglobin

Nanosecond recombination time courses were measured by photolyzing O2, NO, CO, methyl, ethyl, n-propyl, n-butyl, and tert-butyl isocyanide complexes of sperm whale myoglobin with a 30-ns laser pulse at pH 7, 20 degrees C. Absorbance was measured both during and after the excitation pulse and as a function of laser light intensity. The results were analyzed quantitatively in terms of a three-step reaction scheme, MbX in equilibrium B in equilibrium C in equilibrium Mb + X, where Mb is myoglobin, B represents a geminate state in which the ligand is present in the distal pocket but not covalently bound to the iron atom, and C, a state in which the ligand is still embedded in the protein but further away from the heme group. The fitted rate parameters were required to be consistent with the observed overall quantum yield, Q, which had been measured independently using much longer (approximately 0.5 ms) xenon flash pulses. Three major conclusions were derived from these analyses. First, the overall quantum yield of the ligand complex is determined primarily by the competition between the rate of iron-ligand bond formation from the initial photoproduct, kB----MbX, and the rate of migration away from state B, kB----C. For example, kB----C approximately equal to 30-100 microseconds-1 for all three gaseous ligands, whereas both Q and kB----MbX vary over 3 orders of magnitude (i.e. NO, Q = 0.001, kB----MbX approximately equal to 16,000 microseconds-1; O2, Q = 0.1, kB----MbX approximately equal to 500 microseconds-1; CO, Q = 1.0, kB----MbX approximately equal to 2 microseconds-1). Second, for NO, O2, and the isonitriles, the rate-limiting step in the overall association reaction starting from ligand in solution is the formation of state B. The rate constant for this process varies from 2 X 10(7) M-1 s-1 for the gaseous ligands to 0.02-1.4 X 10(5) M-1 s-1 for the isonitriles. In contrast, the B to MbX transition is limiting for CO binding. Third, for all the ligands except CO, the overall rate of dissociation is limited significantly both by the rate of thermal bond disruption, kMbX----B, and the competition between geminate recombination and migration away from the distal pocket (i.e. kB----C/(kB----MbX + kB----C]. In the case of CO, the rate of bond disruption is equal to the observed dissociation rate constant.     

Dielectric properties of hemoglobin and myoglobin. II. Dipole moment of sperm whale myoglobin

This paper is concerned with the molecular origin of the dipole moment of sperm whale myoglobin as it can be calculated from the dielectric dispersion at 1 Mcps on the basis of a mechanism of orientational polarization. It was possible to compare the dielectric increment of native myoglobin and its change during the reaction with bromo acetate with dipole moments calculated according to the known coordinates of the charged groups of the molecule. The agreement between the two shows that in myoglobin only the permanent dipole moment due to these charged groups is important, and that contributions from other possible sources remain within the limits of experimental error.     

Alteration of sperm whale myoglobin heme axial ligation by site-directed mutagenesis

Three mutant proteins of sperm whale myoglobin (Mb) that exhibit altered axial ligations were constructed by site-directed mutagenesis of a synthetic gene for sperm whale myoglobin. Substitution of distal pocket residues, histidine E7 and valine E11, with tyrosine and glutamic acid generated His(E7)Tyr Mb and Val(E11)Glu Mb. The normal axial ligand residue, histidine F8, was also replaced with tyrosine, resulting in His(F8)Tyr Mb. These proteins are analogous in their substitutions to the naturally occurring hemoglobin M mutants (HbM). Tyrosine coordination to the ferric heme iron of His(E7)Tyr Mb and His(F8)Tyr Mb is suggested by optical absorption and EPR spectra and is verified by similarities to resonance Raman spectral bands assigned for iron-tyrosine proteins. His(E7)Tyr Mb is high-spin, six-coordinate with the ferric heme iron coordinated to the distal tyrosine and the proximal histidine, resembling Hb M Saskatoon [His(beta E7)Tyr], while the ferrous iron of this Mb mutant is high-spin, five-coordinate with ligation provided by the proximal histidine. His(F8)Tyr Mb is high-spin, five-coordinate in both the oxidized and reduced states, with the ferric heme iron liganded to the proximal tyrosine, resembling Hb M Iwate [His(alpha F8)Tyr] and Hb M Hyde Park [His(beta F8)Tyr]. Val(E11)Glu Mb is high-spin, six-coordinate with the ferric heme iron liganded to the F8 histidine. Glutamate coordination to the ferric iron of this mutant is strongly suggested by the optical and EPR spectral features, which are consistent with those observed for Hb M Milwaukee [Val(beta E11)Glu]. The ferrous iron of Val(E11)Glu Mb exhibits a five-coordinate structure with the F8 histidine-iron bond intact.(ABSTRACT TRUNCATED AT 250 WORDS)     

Characterization of the reaction of methyl acetimidate with sperm whale myoglobin.

The effects of pH, acetimidate concentration, temperature, and reaction time of methyl acetimidate with sperm whale myoglobulin have been assessed. Reaction at pH 9.8 and 15 degrees C for 30 min with a sixfold excess of methyl acetimidate relative to each amino group yielded six acetimidomyoglobin derivatives which were separated and purified. Reaction with tetrahydrophthalic anhydride revealed the number of amino groups that remained unreacted in each separated component and made possible further subractionation. Modification at the NH2 terminus was quantitated by automated stepwise Edman degradation. The acetimidyl and tetrahydrophthalyl groups, were readily removable. The potentiometric titration of three of the completely deprotected components showed identity with the parent untreated sperm whale myoglobin. The first of two major products was acetimidated at all 19 epsilon-amino groups but not at the NH2 terminus. The second major product bore a blocked NH2 terminus but retained one unmodified epsilon-amino group, identified after modification by trinitrobenzenesulfonate as lysine residue 77. Of the minor components, one was identified as completely acetimidated at all 20 amino groups. The other three minor components appeared to contain irreversible by-products.     

New transient species of sperm whale myoglobin in photodissociation of dioxygen from oxymyoglobin

We carried out the flash photolysis of oxy complexes of sperm whale myoglobin, cobalt-substituted sperm whale myoglobin, and Aplysia myoglobin. When the optical absorption spectral changes associated with the O2 rebinding were monitored on the nanosecond to millisecond time scale, we found that the transient spectra of the O2 photoproduct of sperm whale myoglobin were significantly different from the static spectra of deoxy form. This was sharply contrasted with the observations that the spectra of the CO photoproduct of sperm whale myoglobin and of the O2 photoproducts of cobalt-substituted sperm whale myoglobin and Aplysia myoglobin are identical to the corresponding spectra of their deoxy forms. These results led us to suggest the presence of a fairly stable transient species in the O2 photodissociation from the oxy complex of sperm whale myoglobin, which has a protein structure different from the deoxy form. We denoted the O2 photo-product to be Mb*. In the time-resolved resonance Raman measurements, the nu Fe-His mode of Mb* gave the same value as that of the deoxy form, indicating that the difference in the optical absorption spectra is possibly due to the structural difference at the heme distal side rather than those of the proximal side. The structure of Mb* is discussed in relation to the dynamic motion of myoglobin in the O2 entry to or exit from the heme pocket. Comparing the structural characteristics of several myoglobins employed, we suggested that the formation of Mb* relates to the following two factors: a hydrogen bonding of O2 with the distal histidine, and the movement of iron upon the ligation of O2.     

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.     

Binding of alkylisocyanides with soybean leghemoglobin. Comparisons with sperm whale myoglobin

The binding of various linear and branched chain alkylisocyanides to soybean leghemoglobin has been studied with respect to association and dissociation kinetics and the results compared with those obtained in parallel on sperm whale and horse heart myoglobins; the linear ligands used (methyl to n-heptyl) cover a greater distribution of chain lengths than hitherto used. The association rate constants are much higher for leghemoglobin than for myoglobin, while the dissociation rates are slower. For a given protein, the dissociation rate constants are not much different when different isocyanides are used (except for methyl), whereas the association rates show complex behavior in relation with the alkyl chain length; singular differences are observed between leghemoglobin and sperm whale myoglobin in this regard. For myoglobin, the binding rate constants decrease from methyl to n-propyl, but remain approximately the same when the ligand carries a still longer alkyl chain. In contrast, for leghemoglobin, although the rate constants decrease from methyl to n-propyl, they show a progressive and important rise with longer alkyl substituents: n-butyl and n-pentyl.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.