Electron transfer in chemically and genetically modified myoglobins

The temperature dependences of the reduction potentials (E^o') of wildtype human myoglobin (Mb) and three site-directed mutants have been measured by using thin-layer spectroelectrochemistry. Residue Val68, which is in van der Waals contact with the heme in Mb, has been replaced by Glu, Asp, and Asn. At pH 7.0, reduction of the heme iron (III) in the former two proteins is accompanied by uptake of a proton by the protein. The changes in E^o', and the standard entropy (ΔS^o') and enthalpy (ΔH^o') of reduction in the mutant proteins were determined relative to values for wild-type; the change in E^o' at 25°C was about -200 millivolts for the Glu and Asp mutants, and about -80 millivolts for the Asn mutant. Reduction of Fe(III) to Fe(II) in the Glu and Asp mutants is accompanied by uptake of a proton. These studies demonstrate that Mb can tolerate substitution of a buried hydrophobic group by potentially charged and polar residues, and that such amino acid replacements can lead to substantial changes in the redox thermodynamics of the protein.

Through analysis of the temperature dependence and shapes of NMR dispersion signals, it is determined that a water molecule is bound to the sixth coordination site of the ferric heme in the Val68Asp and in the Val68Asn recombinant proteins while the carboxyl group of the sidechain of Glu68 occupies this position in Val68Glu. The relative rhombic distortions in the ESR spectra of these mutant proteins combined with H₂¹⁷O and spin interconversion experiments performed on them confirm the conclusions of the NMRD study.

The rates of intramolecular electron transfer (ET) of (NH₃)₅Ru-His48 (Val68Asp, His81GIn, Cys110AIa)Mb and (NH₃)₅Ru-His48 (Val68GIu,His81GIn,Cys110Ala)Mb were measured to be .85(3)s^-1 and .30(2)s^-1, respectively. This data supports the hypothesis that entropy of 111 reduction and reorganization energy of ET are inversely related. The rates of forward and reverse ET for (NH₃)₅ Ru-His48 (Val68GIu, His81 GIn, Cys110AIa)ZnMb -7.2(5)•10⁴s^-1and 1.4(2)•10⁵s^-1, respectively- demonstrate that the placement of a highly polar residue nearby does not significantly change the reorganization energy of the photoactive Zn porphyrin.

The distal histidine imidazoles of (NH₃)₄isnRu-His48 SWMb and (NH₃)₅Ru-His48 SWMb were cyanated with BrCN. The intramolecular ET rates of these BrCN-modified Mb derivatives are 5.5(6)s^-1 and 3.2(5)s^-1, respectively. These respective rates are 20 and 10 times faster than those of their noncyanated counterparts after the differences in ET rate from driving force are scaled according to the Marcus equation. This increase in ET rate of the cyanated Mb derivatives is attributed to lower reorganization energy since the cyanated Mb heme is pentacoordinate in both oxidation states; whereas, the native Mb heme loses a water molecule upon reduction so that it changes from six to five coordinate. The reorganization energy from Fe-OH₂ dissociation is estimated to be .2eV. This conclusion is used to reconcile data from previous experiments in our lab. ET in photoactive porphyrin-substituted myoglobins proceed faster than predicted by Marcus Theory when it is assumed that the only difference in ET parameters between photoactive porphyrins and native heme systems is driving force. However, the data can be consistently fit to Marcus Theory if one corrects for the smaller reorganization in the photoactive porphyrin systems since they do not undergo a coordination change upon ET.

Finally, the intramolecular ET rate of (NH₃)₄isnRu-His48 SWMb was measured to be 3.0(4)s^-1. This rate is within experimental error of that for (NH₃)₄pyrRu-His48 SWMb even though the former has 80mV more driving force. One likely possibility for this observation is that the tetraamminepyridineruthenium group undergoes less reorganization upon ET than the tetraammineisonicotinamideruthenium group. Moreover, analysis of the (NH₃)₄isnRu-His48 SWMb experimental system gives a likely explanation of why ET was not observed previously in (NH₃)₄isnRu-Cytochrome C.

I. Electron spin relaxation studies of manganese(II) complexes in acetonitrile. II. Nuclear spin relaxation studies of bromine in tetrabromomanganese(II) in acetonitrile

Magnetic resonance techniques have given us a powerful means for investigating dynamical processes in gases, liquids and solids. Dynamical effects manifest themselves in both resonance line shifts and linewidths, and, accordingly, require detailed analyses to extract desired information. The success of a magnetic resonance experiment depends critically on relaxation mechanisms to maintain thermal equilibrium between spin states. Consequently, there must be an interaction between the excited spin states and their immediate molecular environment which promote changes in spin orientation while excess magnetic energy is coupled into other degrees of freedom by non-radiative processes. This is well known as spin-lattice relaxation. Certain dynamical processes cause fluctuations in the spin state energy levels leading to spin-spin relaxation and, here again, the environment at the molecular level plays a significant role in the magnitude of interaction. Relatively few electron spin relaxation studies of solutions have been conducted and the present work is addressed toward the extension of our knowledge in this area and the retrieval of dynamical information from line shape analyses on a time scale comparable to diffusion controlled phenomena.

Specifically, the electron spin relaxation of three Mn⁺²3d⁵ complexes, Mn(CH₃CN)₆⁺², MnCl₄^-2 in acetonitrile has been studied in considerable detail. The effective spin Hamiltonian constants were carefully evaluated under a wide range of experimental conditions. Resonance widths of these Mn⁺² complexes were studied in the presence of various excess ligand ions and as a function of concentration, viscosity, temperature and frequency (X-band, ~9.5 Ԍ Hz and K-band, ~35 Ԍ Hz).

A number of interesting conclusions were drawn from these studies. For the Et₄NCl-₄^-2 system several relaxation mechanisms leading to resonance broadening were observed. One source appears to arise through spin-orbit interactions caused by modulation of the ligand field resulting from transient distortions of the complex imparted by solvent fluctuations in the immediate surroundings of the paramagnetic ion. An additional spin relaxation was assigned to the formation of ion pairs [Et₄N⁺…MnCl₄^-2] and it was possible to estimate the dissociation constant for this specie in acetonitrile.

The Bu₄NBr-MnBr₄^-2 study was considerably more interesting. As in the former case, solvent fluctuations and ion-pairing of the paramagnetic complex [Bu₄N⁺…MnBr₄^-2] provide significant relaxation for the electronic spin system. Most interesting, without doubt, is the onset of a new relaxation mechanism leading to resonance broadening which is best interpreted as chemical exchange. Thus, assuming that resonance widths were simply governed by electron spin state lifetimes, we were able to extract dynamical information from an interaction in which the initial and final states are the same

MnBr₄^-2 + Br^- = MnBr₄^-2 + Br^-.

The bimolecular rate constants were obtained at six different temperatures and their magnitudes suggested that the exchange is probably diffusion controlled with essentially a zero energy of activation. The most important source of spin relaxation in this system stems directly from dipolar interactions between the manganese 3d⁵ electrons. Moreover, the dipolar broadening is strongly frequency dependent indicating a deviation between the transverse and longitudinal relaxation times. We are led to the conclusion that the 3d⁵ spin states of ion-paired MnBr₄^-2 are significantly correlated so that dynamical processes are also entering the picture. It was possible to estimate the correlation time, ^Td, characterizing this dynamical process.

In Part II we study nuclear magnetic relaxation of bromine ions in the MnBr4-2-Bu4NBr-acetonitrile system. Essentially we monitor the 79Br and 81Br linewidths in response to the [MnBr4-2]/[Br-] ratio with the express purpose of supporting our contention that exchange is occurring between "free" bromine ions in the solvent and bromine in the first coordination sphere of the paramagnetic anion. The complexity of the system elicited a two-part study: (1) the linewidth behavior of Bu4NBr in anhydrous CH3CN in the absence of MnBr4-2 and (2) in the presence of MnBr4-2. It was concluded in study (1) that dynamical association, Bu4NBr k1= Bu4N+ + Br-, was modulating field-gradient interactions at frequencies high enough to provide an estimation of the unimolecular rate constant, k1. A comparison of the two isotopic bromine linewidth-mole fraction results led to the conclusion that quadrupole interactions provided the dominant relaxation mechanism. In study (2) the "residual" bromine linewidths for both 79Br and 81Br are clearly controlled by quadrupole interactions which appear to be modulated by very rapid dynamical processes other than molecular reorientation. We conclude that the "residual" linewidth has its origin in chemical exchange and that bromine nuclei exchange rapidly between a "free" solvated ion and the paramagnetic complex, MnBr4-2.

I. Quantum mechanics of one-dimensional two-particle models. II. A quantum mechanical treatment of inelastic collisions

Part I

Solutions of Schrödinger’s equation for system of two particles bound in various stationary one-dimensional potential wells and repelling each other with a Coulomb force are obtained by the method of finite differences. The general properties of such systems are worked out in detail for the case of two electrons in an infinite square well. For small well widths (1-10 a.u.) the energy levels lie above those of the noninteresting particle model by as much as a factor of 4, although excitation energies are only half again as great. The analytical form of the solutions is obtained and it is shown that every eigenstate is doubly degenerate due to the “pathological” nature of the one-dimensional Coulomb potential. This degeneracy is verified numerically by the finite-difference method. The properties of the square-well system are compared with those of the free-electron and hard-sphere models; perturbation and variational treatments are also carried out using the hard-sphere Hamiltonian as a zeroth-order approximation. The lowest several finite-difference eigenvalues converge from below with decreasing mesh size to energies below those of the “best” linear variational function consisting of hard-sphere eigenfunctions. The finite-difference solutions in general yield expectation values and matrix elements as accurate as those obtained using the “best” variational function.

The system of two electrons in a parabolic well is also treated by finite differences. In this system it is possible to separate the center-of-mass motion and hence to effect a considerable numerical simplification. It is shown that the pathological one-dimensional Coulomb potential gives rise to doubly degenerate eigenstates for the parabolic well in exactly the same manner as for the infinite square well.

Part II

A general method of treating inelastic collisions quantum mechanically is developed and applied to several one-dimensional models. The formalism is first developed for nonreactive “vibrational” excitations of a bound system by an incident free particle. It is then extended to treat simple exchange reactions of the form A + BC →AB + C. The method consists essentially of finding a set of linearly independent solutions of the Schrödinger equation such that each solution of the set satisfies a distinct, yet arbitrary boundary condition specified in the asymptotic region. These linearly independent solutions are then combined to form a total scattering wavefunction having the correct asymptotic form. The method of finite differences is used to determine the linearly independent functions.

The theory is applied to the impulsive collision of a free particle with a particle bound in (1) an infinite square well and (2) a parabolic well. Calculated transition probabilities agree well with previously obtained values.

Several models for the exchange reaction involving three identical particles are also treated: (1) infinite-square-well potential surface, in which all three particles interact as hard spheres and each two-particle subsystem (i.e. BC and AB) is bound by an attractive infinite-square-well potential; (2) truncated parabolic potential surface, in which the two-particle subsystems are bound by a harmonic oscillator potential which becomes infinite for interparticle separations greater than a certain value; (3) parabolic (untruncated) surface. Although there are no published values with which to compare our reaction probabilities, several independent checks on internal consistency indicate that the results are reliable.

Electron microprobe investigations of metamorphic reactions and mineral growth histories, Kwoiek area, British Columbia

The Kwoiek Area of British Columbia contains a pendant or screen of metamorphosed sedimentary and volcanic rocks almost entirely surrounded by a portion of the Coast Range Batholith, and intruded by several dozen stocks. The major metamorphic effects were produced by the quartz diorite batholithic rocks, with minor and later effects by the quartz diorite stocks. The sequence of important metamorphic reactions in the metasedimentary and metavolcanic rocks, ranging in grade from chlorite to sillimanite, is:

1. chlorite + carbonate + muscovite → epidote + biotite

2. chlorite + carbonate → actinolite + epidote

3. chlorite + muscovite → garnet + biotite

4. chlorite + epidote → garnet + hornblende

5. chlorite + muscovite → garnet + staurolite + biotite

6. chlorite + muscovite → aluminum silicate + biotite

7. muscovite + staurolite → garnet + aluminum silicate + biotite

8. staurolite → garnet + aluminum silicate

Continuous reactions, occurring between reactions 5 and 7, are:

A. chlorite + (high Ti) biotite + Al₂O₃ (from plagioclase?)→ garnet + staurolite + (low Ti) biotite + O₂

B. muscovite (phengitic) → garnet + staurolite +muscovite (less phengitic) + O₂ (?)

Detailed electron microprobe work on garnet, staurolite, biotite, and chlorite shows that:

(1) The garnet porphyroblasts are zoned according to a depletion model, called the Rayleigh depletion model, which assumes equilibrium between the edge of a growing garnet and the minerals which are unzoned, notably biotite, chlorite, and muscovite, but which assumes disequilibrium within the garnet.

(2) The staurolite porphyroblasts are also zoned, and from their zoning patterns reactions A, B, and 5 are documented. Progressive reduction of iron with increasing grade of metamorphism is also inferred from the staurolite zoning patterns.

(3) During a late period of falling temperature garnet continued to grow and the biotite and chlorite reequilibrated. The biotite, chlorite, and garnet edge compositions can vary from point to point in a given thin section, indicating that the volume of equilibrium at the final stage of metamorphism was only a few cubic microns.

(4) The horizon within the garnet that grew at maximum temperature can be identified. The Mg/Fe ratio of this horizon, if the garnet composition is a limiting composition in the Al₂O₃ - K₂O - FeO - MgO tetrahedron, increases systematically with increasing metamorphic grade. Biotite and chlorite compositions also show a general increase in Mg/Fe ratio with increasing metamorphic grade, but staurolite appears to show the reverse effect.

(5) The Mg/Fe ratio at the maximum temperature horizon of the garnet porphyroblasts is a function of its Mn content as evidenced from the study of five garnet-bearing rocks, collected from one outcrop area, with the same assemblage but with differing proportions of minerals.

An important implication of zoned minerals is that the effective composition of a system in a phase lies on the join between the homogeneous minerals (if there are two) and not within three-or- four-phase fields when a zoned mineral, such as garnet or staurolite, is present in the assemblage.

Study of the three aluminum silicates found in the Kwoiek Area showed that a constant pressure change in polymorphs from andalusite to kyanite to sillimanite took place with increasing temperature. This transition series is best explained by the metastable formation of andalusite.

Photographic materials on pages 15, 121, 160, 162, and 164 are essential and will not reproduce clearly on Xerox copies. Photographic copies should be ordered.