Solution and Crystal Structures of a Sperm Whale Myoglobin Triple Mutant That Mimics the Sulfide-binding Hemoglobin from Lucina pectinata*

Bao D. Nguyen, Xuefeng Zhao, Krishnamurthi Vyas, and Gerd N. La MarDagger

From the Department of Chemistry, University of California, Davis, California 95616

R. Ashley Lile, Eric Allen Brucker§, George N. Phillips Jr., and John S. Olson

From the Department of Biochemistry and Cell Biology and the W. M. Keck Center for Computational Biology, Rice University, Houston, Texas 77005-1892

Jonathan B. Wittenberg

From the Department of Physiology and Biophysics, Albert Einstein College of Medicine, Bronx, New York 10461

    ABSTRACT
Top
Abstract
Introduction
Procedures
Results
Discussion
References

The bivalve mollusc Lucina pectinata harbors sulfide-oxidizing chemoautotrophic bacteria and expresses a monomeric hemoglobin I, HbI, with normal O2, but extraordinarily high sulfide affinity. The crystal structure of aquomet Lucina HbI has revealed an active site with three residues not commonly found in vertebrate globins: Phe(B10), Gln(E7), and Phe(E11) (Rizzi, M., Wittenberg, J. B., Coda, A., Fasano, M., Ascenzi, P., and Bolognesi, M. (1994) J. Mol. Biol. 244, 86-89). Engineering these three residues into sperm whale myoglobin results in a triple mutant with ~700-fold higher sulfide affinity than for wild-type. The single crystal x-ray structure of the aquomet derivative of the myoglobin triple mutant and the solution 1H NMR active site structures of the cyanomet derivatives of both the myoglobin mutant and Lucina HbI have been determined to examine further the structural origin of their unusually high sulfide affinities. The major differences in the distal pocket is that in the aquomet form the carbonyl of Gln64(E7) serves as a H-bond acceptor, whereas in the cyanomet form the amido group acts as H-bond donor to the bound ligand. Phe68(E11) is rotated ~90° about chi 2 and located ~1-2 Å closer to the iron atom in the myoglobin triple mutant relative to its conformation in Lucina HbI. The change in orientation potentially eliminates the stabilizing interaction with sulfide and, together with the decrease in size of the distal pocket, accounts for the 7-fold lower sulfide affinity of the myoglobin mutant compared with that of Lucina HbI.

    INTRODUCTION
Top
Abstract
Introduction
Procedures
Results
Discussion
References

Despite substantial amino acid sequence variability, the reversible binding of molecular oxygen in myoglobin (Mb)1 and hemoglobin (Hb) is achieved by a surprisingly invariant protein folding topology: a heme group imbedded within 7-8 packed alpha -helices (1, 2). Examination of the more than 300 known sequences of vertebrate Hbs and Mbs and the >130 nonvertebrate globin-like sequences demonstrates that only two residues appear to be conserved, the proximal His at helical position F-8 and the Phe parallel to the heme surface at an interhelical loop, position CD1 (1, 3). The residue considered to be pivotal in the stabilization of the bound O2 is the one whose side chain can provide a neutral donor for a hydrogen bond with the ligand, such as His or Gln at the distal position E-7. Among vertebrate Mbs and Hbs, His predominates, with Gln occurring only in elephant Mb (4), shark Mb (5), hagfish Hb (6), and the alpha  chain of opossum Hb (7). Although His is also found in the majority of nonvertebrate Hbs and Mbs, Gln is not uncommon and, in addition, other residues may occupy the E-7 position, e.g. Val in Aplysia (8) and several other mollusc Mbs (9), Tyr in the Mb from Paramphistomun epiclitum (10, 11), and Leu in the monomeric components of Glycera Hb (12). Although capable of O2 binding, a number of the variant globins exhibit some anomalous functional properties. Parallel functional and structural studies on genetically engineered mutants of mammalian Mb have shown that changes in only a very limited number of residues in the distal pocket often transfers the essence of the unusual behavior of the natural genetic variant to the appropriate reference protein mutant (13-17).

Cytoplasmic Hbs occur in many symbiotic associations between bacteria and invertebrates or plants and may have a physiological role in the symbiosis (18-20). The bivalve mollusc Lucina pectinata, found in sulfide-rich coastal sediments, harbors sulfide-oxidizing chemoautotrophic bacteria (21, 22). Its abundant cytoplasmic hemoglobins consist of three single-chain components, each having a moderately high affinity for oxygen (P50 = 0.1-0.2 torr), which is achieved, however, by a very different balance of combination and dissociation rates. Reactions of the monomeric HbI with oxygen are rapid, whereas those of HbII and HbIII are extraordinarily slow (23).

Recent crystal structures of Lucina aquometHbI and the sulfide bound form (24, 25) have revealed a distal Gln64(E7) and Phe29(B10), neither of which is rare among invertebrate Mbs, Hbs. However, the E-11 position is uniquely occupied by Phe rather than the aliphatic amino acids Val, Ile, or Leu. To examine the structural origin of strong sulfide binding in Lucina HbI, we have systematically replaced His64(E7), Leu29(B10), and Val64(E11) with Gln, Phe, and Phe, respectively, in sperm whale myoglobin and then determined the structures of the resulting multiple mutants by NMR and x-ray crystallography. This approach is similar to that used to study the unusual functional properties of Aplysia Mb, Ascaris Hb, and elephant Mb, and provides comparison between the native protein and the synthetic mimic in two ligation/spin states and between protein structures in crystal and solution (13-17, 26, 27). The NMR studies were pursued to define more carefully the positions of the labile protons of Gln64(E7) and their role as hydrogen bond donors to the bound ligand. Both the paramagnetic relaxation and induced dipolar shifts in metMbCN and metHbICN can be used to model the position of distal residues (14-17, 26, 29, 30). The combined NMR and crystallographic results reveal systematic differences between the orientation of Gln64(E7) in the aquomet and cyanomet complexes of Lucina HbI and L29F/H64Q/V68F-Mb, which demonstrate that Gln64(E7) may serve as either a H-bond acceptor or donor to bound ligands and that a difference in the orientation of the Phe68(E11) ring may account for the difference in sulfide affinity between Lucina HbI and the sperm whale triple mutant Mb mimic.

    EXPERIMENTAL PROCEDURES
Top
Abstract
Introduction
Procedures
Results
Discussion
References

Protein Preparation-- The pUC19 plasmid containing V68F sperm whale myoglobin was made by Egeberg et al. (31). The plasmid containing L29F/H64Q/V68F sperm whale Mb was constructed from pEMBL19 containing L29F Mb using cassette mutagenesis as described in Springer et al. (32). Vectors harboring the mutated gene were transformed into Escherichia coli strain TB1 and expressed constitutively using a 100-liter fermentor. The resulting soluble holomyoglobin protein was purified as described previously (33, 34).

Protein Crystallization-- Crystals of the recombinant L29F/H64Q/V68F-Mb were grown by the batch method (35, 36) in a temperature-controlled environment (17 °C). Large single crystals developed out of solutions ranging from 2.4 to 2.7 M ammonium sulfate buffered at pH 9 (20 mM Tris-HCl, 1 mM EDTA).

X-ray Data Collection and Structure Determination-- Crystals were hexagonal, space group P6, with one molecule per asymmetric unit as originally solved by Phillips et al. (35). They were mounted in sealed quartz capillary tubes prior to data collection. A single data set was collected from one crystal at room temperature using a Rigaku R-axis IIc imaging plate system and copper Kalpha radiation from a Siemens rotating anode operated at 50 mV and 90 mA. The unit cell dimensions of the L29F/H64Q/V68F-metMbH2O were a = 91.39, c = 45.75. A total of 74,100 measured reflections from 60 images (Delta omega  = 1.5°) with an Rmerge of 5.7% for all data were reduced to 15,526 unique intensities (5.0 to 1.85 Å, 87.3% complete) using the program XDS.

Crystal Structure Refinement-- Starting coordinates were generated using the L29F/H64Q-MbCO (Protein Data Bank entry 1MCY, Brookhaven National Laboratory) (17) to calculate initial phases. The coordinates for the additional mutated residue (V68F) were built from an electron density omit map. Solvent atoms were introduced in positive density peaks over four sigma . Cycles of conventional positional refinement were carried out with XPLOR (37-39) alternated with manual fitting using the CHAIN software package (40). The crystal structures were all refined with Engh and Huber (41) topology and parameters with no restraints on the iron atom to remove bias in metal position. The crystallographic refinements converged to an R-factor of 16.8% for the triple mutant metmyoglobin structure with root mean square bond deviations of 0.017 Å. Coordinates have been deposited in the Protein Data Bank (entry 1OBM).

1H NMR Measurements-- All the 1H NMR spectra were collected on the GE Omega 500 MHz spectrometer. The strongly relaxed signals were optimally detected in water-eliminated Fourier transform spectra (42). Nonselective T1s for the resolved strongly relaxed protons were measured via inversion-recovery experiment. Steady state NOEs were recorded as described previously (43). The phase-sensitive TOCSY (44, 45), NOESY (46), and conventional magnitude COSY (MCOSY) (47) employed the method described by States et al. (48) to provide quadrature detection in the t1 dimension. Solvent suppression, when required, was achieved by direct saturation in the relaxation delay period. 512 blocks were collected with 25.0 kHz spectral widths to include all resonances, and 10 kHz to improve resolution for the diamagnetic envelope. 128 to 256 scans were accumulated with repetition rate of 0.7 s-1 or 1.2 s-1 for each block with free induction decays of 2048 complex points. The data were processed as described previously (49); details are given in the figure captions. All two-dimensional data were processed on Silicon Graphics workstation using the software package Felix from Biosym/MSI (San Diego).

Magnetic Axes Determination-- The magnetic axes were determined as described previously (15, 16, 29, 30, 50). Experimental dipolar shifts for the structurally conserved proximal side of the heme were used as input to search for the Euler rotation angles, Gamma (alpha ,beta ,gamma ), that transforms the molecular pseudo-symmetry coordinates (x', y', z', or R, theta ', Omega ' (Fig. 1)) readily obtained from crystal coordinates (17, 24, 25, 51) into magnetic axes, x, y, and z, by minimizing the global error function,
<FR><NU><UP>F</UP></NU><DE>n</DE></FR>=<LIM><OP>∑</OP><UL>n</UL></LIM>‖&dgr;<SUB><UP>dip</UP></SUB>(<UP>obs</UP>)−&dgr;<SUB><UP>dip</UP></SUB>(<UP>calc</UP>)&Ggr;(&agr;, &bgr;, &ggr;)‖<SUP>2</SUP> (Eq. 1)
where
&dgr;<SUB><UP>dip</UP></SUB>(<UP>calc</UP>)=<FR><NU>1</NU><DE>3<UP>N</UP></DE></FR><FENCE>&Dgr;&khgr;<SUB><UP>ax</UP></SUB>(3<UP>cos</UP><SUP>2</SUP>&thgr;′−1)<UP>R</UP><SUP><UP>−</UP>3</SUP>+<FR><NU>3</NU><DE>2</DE></FR>(&Dgr;&khgr;<SUB><UP>rh</UP></SUB> <UP>sin</UP><SUP>2</SUP>&thgr;′<UP>cos</UP>2&OHgr;′)<UP>R</UP><SUP><UP>−</UP>3</SUP></FENCE> (Eq. 2)
and
&dgr;<SUB><UP>dip</UP></SUB>(<UP>obs</UP>)=&dgr;<SUB><UP>DSS</UP></SUB>(<UP>obs</UP>)−&dgr;<SUB><UP>dia</UP></SUB> (Eq. 3)
Delta chi ax and Delta chi rh are axial and rhombic anisotropies, and delta DSS(obs) is the observed chemical shift referenced to DSS. delta dia is the shift in the isostructural diamagnetic MbCO complex (52, 53), or calculated for protons whose delta dia are not available as described in detail by Qin et al. (15, 16). Minimizing the error function F/n in Equation 1 was performed over three parameters, alpha , beta , and gamma , using available Delta chi ax and Delta chi rh, or extended to all five parameters to yield both the Euler angles and anisotropies as described in detail previously (29).


View larger version (32K):
[in this window]
[in a new window]
 
Fig. 1.   Schematic structure of the heme cavity in sperm whale mutant L29F/H64Q/V68F-Mb and Lucina HbI with the proximal (squares) and distal (circles) residues in contact with heme, axial His, and ligand identified by the position on helices. In cases where the residues differ for the two proteins, the residue is labeled first for L29F/H64Q/V68F-Mb and second for Lucina HbI. The solid lines reflect dipolar (NOESY) contacts observed for both proteins; dotted and dashed lines represent such dipolar contacts observed only in the triple mutant and Lucina HbI, respectively. For all but Gln64(E7) in Lucina HbI, the observed dipolar contacts are completely consistent with the crystal structure of their respective aquo-met forms. The x',y',z' coordinate system represents an iron-centered coordinate system derived from crystal coordinates of MbCO. The magnetic axes, x, y, z, in which the paramagnetic susceptibility tensor is diagonal, are related by the Euler rotation Gamma (alpha , beta , gamma ) via (x, y, z) = (x', y', z')Gamma (alpha , beta , gamma ), where beta  represents the tilt of major magnetic axes from the heme normal, alpha  is the direction of this tilt, defined as the angle between the projection of -z on the x',y' plane and x' axis, and kappa  ~ alpha  + gamma  locates the projection of the rhombic axes onto the heme plane.

Dipolar Shift Simulations-- The position of a substituted or perturbed residue can be determined by minimizing a local error function. This local error function, designated F*(residue)/n' to distinguish it from that global error function in Equation 1, is given by Ref. 30,
F*(<UP>residue</UP>)/n′=<LIM><OP>∑</OP><UL>n</UL></LIM>‖&dgr;<SUB><UP>dip</UP></SUB>(<UP>obs</UP>)−&dgr;<SUB><UP>dip</UP></SUB>(&phgr;, &OHgr;, <UP>r</UP>)‖<SUP>2</SUP> (Eq. 4)
where delta dip(phi , Omega , r) represents the delta dip(calc) as a function of a bond rotation phi  or a translation of a residue by a distance r in a direction defined by the angle Omega , using the magnetic axes derived from conserved structural elements; n' is the number of residue dipolar shifts. The bond angle, phi , or movement defined by Omega  angle and distance r, that mimize the residual error function F*(residue)/n', defines structural changes as described in detail previously (15-17, 29). When available, the influence of paramagnetic relaxation, T1-1 proportional to  RFe-6, allowed estimate of RFe using the relation,
T<SUB><UP>1i</UP></SUB>/T<SUB><UP>1j</UP></SUB>=R<SUP>6</SUP><SUB><UP>Fe−i</UP></SUB>/<UP>R</UP><SUP>6</SUP><SUB><UP>Fe−j</UP></SUB> (Eq. 5)
where the His(F8) ring NH (RFe = 5.0 Å) and heme methyl (RFe = 6.1 Å) provide upper and lower limits to a proton of interest (54). The molecular modeling was carried out on a Silicon Graphics INDIGO from available crystal coordinates (17, 24, 25, 51) using the INSIGHT II (Biosym/MSI, San Diego) and MIDAS (UCSF) programs.

    RESULTS
Top
Abstract
Introduction
Procedures
Results
Discussion
References

Hydrogen Sulfide and O2 Binding to Lucina HbI and Sperm Whale Mb

A brief summary of the affinities of the myoglobin mutants and Lucina HbI for O2 and hydrogen sulfide2 is shown in Table I (17, 23, 28). Lucina HbI shows a 5,000-fold higher affinity for sulfide than recombinant wild-type sperm whale myoglobin, whereas the two proteins have similar affinities for oxygen. There are progressive increases in KH2S as the active site of the whale myoglobin is altered to resemble that of the mollusc hemoglobin. All three single mutations, L29F, H64Q, and V68F produce large, 6-20-fold increases in sulfide affinity, and the effects on KH2S are roughly additive in the double and triple mutants. In contrast, these mutations produce markedly different effects on O2 affinity. For example, the single H64Q and L29F mutations produces a 5-fold decrease and 15-fold increase, respectively, in O2 affinity. As a result, KO2 for the triple mutant is only 3-4-fold higher than wild-type myoglobin and very similar to that of Lucina HbI (Table I).

                              
View this table:
[in this window]
[in a new window]
 
Table I
Equilibrium dissociation constants for O2 and H2S binding to Lucina HbI and mutants of recombinant sperm whale myoglobin

High Resolution Crystal Structure of Sperm Whale L29F/H64Q/V68F-metMbH2O

The three mutated side chains in the recombinant myoglobin do not alter the overall tertiary structure of myoglobin nor do they appear to sterically hinder the bound ligand. Two of these mutated residues, Gln64(E7) and Phe68(E11), occupy positions equivalent to those observed in the structures of corresponding H64Q-Mb and V68F-Mb single mutants (36, 55). In contrast, the aromatic ring at position 29 in the triple mutant is perpendicular to the orientation found in sperm whale L29F- and L29F/H64Q-MbCO (Fig. 2) (17, 56). A detailed comparison of the distal pockets of the L29F/H64Q/V68F-metMbH2O and Lucina metHbIH2O (24) is presented in Table II and Fig. 2. The overall positions of the three mutated residues in the recombinant sperm whale metmyoglobin are close to those of the corresponding residues in the clam protein, but the exact orientations of the phenyl side chains are significantly different. The distal glutamines in both proteins form hydrogen bonds with the coordinated water molecule.


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 2.   Stereo diagrams of the orientations of the distal Phe(B10), Gln(E7), and Phe(E11) residues. A, L29F/H64Q/V68F-metMbCN in solution (dark lines), and L29F/H64Q/V68F-metMbH2O in a crystal (light lines); B, Lucina metHbICN in solution (dark lines) and Lucina metHbIH2O (24) in a crystal (light lines); C, L29F/H64Q/V68F-metMbH2O in a crystal (light lines) and Lucina metHbIH2O (24) in a crystal (dark lines); the conserved Phe43(CD1) is also included in this comparison. The side chain oxygen atom of the distal glutamine is shown as a ball in each case.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Conformational parameters in the distal pockets for the aquo-met (crystal) and cyanomet (solution) structures of sperm whale L29F/H64Q/V68F-Mb and Lucina pectinata HbI

Heme Pocket Structure of L29F/H64Q/V68F-metMbCN and Lucina metHbICN

Resonance Assignments-- The 500 MHz 1H NMR spectra for L29F/H64Q/V28/F-metMbCN, V68F-metMbCN, and Lucina-metHbICN in 1H2O at 25 °C are shown in Fig. 3, A-C, respectively. The signals for the heme for each protein were located and assigned by the characteristic pattern of dipolar contacts about the heme periphery among the TOCSY identified vinyl, propionate groups, and the pyrrole methyl, meso-Hs in a fashion standard for both diamagnetic and paramagnetic heme proteins. The pattern of heme chemical shifts is similar to that observed in wild-type sperm whale metMbCN (not shown). Resonances for residues near the heme are assigned to the degree possible, as limited by spectral congestion and paramagnetic relaxation, and by standard backbone dipolar connectivities together with characteristic TOCSY and/or COSY connectivities for several side chains (15, 16, 49). Other residues not addressable by these approaches are assigned on the basis of dipolar contacts to the heme and/or to other assigned residues. Inasmuch as this assignment strategy has been reported in detail for both wild-type and mutant metMbCN (17, 29, 30, 43), two-dimensional data are shown only as relevant to the placement of perturbed distal pocket residues in both protein complexes.


View larger version (25K):
[in this window]
[in a new window]
 
Fig. 3.   500 MHz 1H NMR spectra of (A) L29/H64Q/V68F-metMbCN, (B) V68F-metMbCN, and (C) Lucina metHbICN in 1H2O at 25 °C, pH 7.6. Heme resonances are labeled by the Fischer notation and protein signals by the standard one-letter amino acid code. Steady-state NOE difference spectra upon saturating the low-field Gln64(E7) Nepsilon H peaks in: D, L29F/H64Q/V68F-metMbCN, and E, Lucina metHbICN.

The fingerprint portion of the NOESY and COSY maps (not shown) for each protein locates two extended helical fragments each, the shorter one corresponding to F helix residues F4-F9 that includes the strongly hyperfine shifted axial His(F8). The strongly relaxed and hyperfine shifted proximal FG corner residues (His97(FG3) and Ile99(FG5) in Mb mutants, Arg99(FG1) and Ile101(FG4) in Lucina HbI) are identified by their characteristic NOESY cross-peaks to the 5-CH3. For each protein, the pattern of heme residue and inter residue dipolar contacts and paramagnetic relaxation reveal a proximal side heme pocket structure that is indistinguishable from that in the respective reference crystal structures (Fig. 1). The chemical shifts for relevant assigned resonances for the three protein complexes are listed in Table III.

                              
View this table:
[in this window]
[in a new window]
 
Table III
1H NMR chemical shifts of the amino acid residues in sperm whale L29F/H64Q/V68F-metMbCN and V68F-metMbCN and in Lucina metHbICN
Chemical shifts in ppm at 25 °C (L29F/H64Q/V68F-metMbCN and Lucina metHbICN) and at 30 °C (V68F-metMbCN), pH 7.2, reference to DSS.

The longer helical fragment represents the E helix residues E-7-E-14 in each protein, with residues E-10, E-11 and E-14 exhibiting the expected NOESY cross-peaks in the heme, including that of the Phe68(E11) ring to the 2-vinyl group. The TOCSY/MCOSY detected fragment NHCalpha H-Cbeta H for Gln64(E7) in each protein shows strong NOESY cross-peaks to a moderately relaxed and hyperfine shifted CH-CH2 fragment which is assigned to the remainder of Gln64(E7); the Cbeta H-Cbeta 'H TOCSY/COSY connectivities are not observed because the cross-peak is too close to the diagonal in each protein (Figs. 4 and 5). A strongly relaxed (T1 ~ 12 ms) and strongly low-field hyperfine shifted labile proton, when saturated, exhibits an intense NOE to another upfield shifted and moderately relaxed (T1 ~ 50 ms) labile proton in Lucina metHbICN (Fig. 3E), but not detected because of overlap with the solvent signals in the triple mutant, as well as to the Gln64(E7) Cgamma Hs (and for the triple mutant, to Thr67(E10); Fig. 3D) and confirms the two labile protons as arising from the Nepsilon Hs of Gln64(E7). The T1s via Equation 5 reveal, RFe(Nepsilon H) = 4.2 ± 0.2 Å for the low-field NH for both proteins, and for Lucina metHbICN, the upfield resolved RFe(Nepsilon H') = 4.7 ± 0.2 Å. The resulting RFe are listed in Table II.


View larger version (45K):
[in this window]
[in a new window]
 
Fig. 4.   Portions of the MCOSY (A, A') and NOESY (tau m = 50 ms) (B, B') spectra of Lucina metHbICN in 2H2O at 25 °C, pH 7.0, showing the identification of the Phe29(B10) and Gln64(E7) spin system and the relevant dipolar contact between the two residues. The MCOSY data were processed by applying 512 t1 × 512 t2 points prior to zero filling to 2048 × 2048 data points and Fourier transformation. The NOESY data were processed by applying 30°-shifted sine-bell squared window over 1024 t1 × 256 t2 prior to zero-filling to 2048 × 2048 data points and Fourier transformation.


View larger version (44K):
[in this window]
[in a new window]
 
Fig. 5.   Portions of the MCOSY (A) and NOESY (tau m = 50 ms) (B) spectra of L29F/H64Q/V68F-metMbCN in 1H2O at 25 °C, pH 7.0, showing the identification of the Phe29(B10) and Gln64(E7) spin system and the relevant dipolar contact between the two residues. MCOSY and NOESY data were processed as described in detail in the legend to Fig. 4.

One strongly relaxed (for Czeta H) and hyperfine shifted aromatic ring with NOESY contact to 5-CH3 and Gln64(E7) (not shown) which uniquely identifies Phe43(CD1). Another strongly hyperfine shifted and relaxed (Czeta H 19 ppm, T1 ~ 33 ms) aromatic ring exhibits NOESY cross-peaks to Phe68(E11), Gln64(E7), Gly/Ala65(E8), and Phe43(CD1) and definitively assigns Phe29(B10) for each protein (Figs. 4 and 5). The chemical shifts of relevant assigned distal residues for the three proteins are included in Table III. The pattern of dipolar contacts between distal residues and the heme, and among distal residues, are predicted by the crystal structure, as shown schematically in Fig. 1. The difference that are noted, as described below, involve primarily the residues Gln64(E7).

Magnetic Axes Determination-- All dipolar shifted residues except Phe68(E11) exhibit a good correlation between delta dip(obs) and the slope of the chemical shift in a Curie (chemical shift versus T-1) plot that is indicative of well defined orientations to the heme (not shown) (49). A variety of proximal proton dipolar shifts which correlate with their Curie slopes for each protein were used.3 To determine the magnetic axes, both as three-parameter searches for Gamma (alpha ,beta ,gamma ) using the wild-type metMbCN anisotropies, Delta chi ax = 2.04 × 10-9 m3/mol, Delta chi rh = -0.48 × 10-9 m3/mol (29), and as five-parameter searches for both Gamma (alpha ,beta ,gamma ) and the anisotropies. Both mutants exhibit excellent correlation between delta dip(obs) and delta dip(calc) with very low residual F/n in Equation 1 (not shown) typical of previous NMR studies of mutant metMbCNs (15, 16, 29, 30, 50). The chemical shifts for numerous distal residues are also well predicted. The resulting Euler angles using different input data sets were highly clustered for each protein, although distinct for each of the proteins. For L29F/H64Q/V68F-metMbCN, the optimized angles (and ranges for the various fits) are: a tilt of the major axes from the heme normal, beta  ~ 6.3° (6.1-6.6°), direction of tilt, alpha  ~ -40° (-40° to -50°) and rhombic axes, kappa  ~ alpha  + gamma  = 30° (20° to 40°) using MbCO crystal coordinates. The magnetic axes for V68F-metMbCN are alpha  = 5° (0-10°), beta  = 7.5° (7-8°), and kappa  = 40° (30-50°). The optimized anisotropies for each five-parameter fit differ inconsequentally for those of wild-type metMbCN and yield Euler angles within the ranges obtained by the 3-parameter fit (not shown). When the coordinates for L29F/H64Q/V68F-metMbCN are taken from the L29F/H64Q/V68F-metMbH2O crystal structure presented herein, we obtain values for the angles and anisotropies that are well within the ranges obtained using the wild-type coordinates. The magnetic axes (and ranges in their values), obtained similarly for Lucina metMbICN, using either the Lucina metHbIH2O (24) or metHbIH2S (25) crystal coordinates, are beta  = 7° (6-8°), alpha  = 155 (150-160), kappa  ~ alpha  + gamma  = 255° (250-260°), and optimized anisotropies are unchanged from those of wild-type metMbCN (not shown).

Orientation of Distal Residues in L29F/H64Q/V68F-metMbCN-- Using either the wild-type MbCO crystal coordinates with the mutated Phe29(B10) inserted with an orientation as found in the crystal structure of L29F/H64Q-MbCO (half-open squares in Fig. 6A), or directly the crystal coordinates of the presently characterized L29F/H64Q/V68F-metHbH2O (open squares in Fig. 6A), a qualitative, but not quantitative, correlation between delta dip(calc) and delta dip(obs) are observed, and the RFe (Czeta H = 4.7 Å) is shorter than indicated by T1 = 33 ms (RFe = 5.1 ± 0.3 Å). However, altering chi 1 by 5° from that in the L29F/H64Q-metMbCO and a 50° rotation of chi 2 leads to a Phe29(B10) orientation whose delta dip(calc) correlate very well with delta dip(obs) (closed squares in Fig. 6A) and for which the RFe (Czeta H) = 5.1 Å is consistent with the T1 value.


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 6.   Plot of delta dip(obs) versus delta dip(calc) for key distal residues. A, L29F/H64Q/V68F-metMbCN for Phe29(B10) (squares) and Gln64(E7) (circles) based on the L29F/H64Q-MbCO crystal structure (half-filled marker), L29F/H64Q/V68F-metMbH2O (open marker) crystal structures, and the latter crystal structures with Gln64(E7) chi 3 rotated by 180° (open marker with asterisk) as determined by NMR (closed markers); B, Lucina metHbICN for Phe29(B10) (squares), Gln64(E7) (circles) and Phe43(CD1) (triangles) based on the crystal structures of Lucina metHbIH2O (24) (open markers), the crystal structure of Lucina metHbIH2S (25) (half-closed markers) and with Gln64(E7) chi 3 rotated by 180° (half-closed markers with asterisk), and as determined by NMR (closed markers).

Introduction of Gln64(E7) into the wild-type MbCO crystal structure based on its orientation in the H64Q-MbCO (or L29F/H64Q-MbCO) crystal structure results in reasonable agreement with both T1 data and delta dip(calc) (half-closed circles in Fig. 6A). However, when directly using the crystal coordinates for L29F/H64Q/V68F-metMbH2O, neither the delta dip(calc) (open circles in Fig. 6A) nor T1 value for Nepsilon H (RFe = 5.9 Å) are reasonably predicted (T1 = 14 ms, RFe = 4.2 ± 0.2). Small rotation of chi 1 (~5°) and chi 2 (~8°), followed by a 180° rotation of chi 3 for Gln64(E7) in the L29F/H64Q/V68F-metMbH2O crystal structure leads to an excellent fit for delta dip(obs) versus delta dip(calc) (closed circles in Fig. 6A) and yields RFe(Nepsilon H) = 4.2 Å, in good agreement with the 4.2-Å value obtained from the T1 = 14 ms via Equation 5. The other Nepsilon H is correctly predicted to resonate close to the solvent signal by the optimized delta dip(calc). Similarly small variations in chi 1, chi 2 starting with the Gln64(E7) orientation in the single mutant lead to a Gln orientation essentially the same as that obtained starting with the coordinate of the triple mutant (not shown).

The anomalous variable temperature slopes (see above), small delta dip(calc) which are relatively insensitive to ring orientation, and the presence of comparable magnitude but opposed sign delta dip(calc) and ring current shift rendered the optimization of the orientation of Phe68(E11) ring on the basis of delta dip(obs) completely impractical (50). The pattern of small delta dip for Phe68(E11), nevertheless, are consistent with observed dipolar shifts and NOESY contacts. The orientation for these three mutated residues in L29F/H64Q/V68F-metMbCN are shown in Fig. 2A, and the relevant side chain bond angles and distances to the iron are listed in Table II.

Orientation of Distal Residues in Lucina MetHbICN-- The orientations for Gln64(E7) in the Lucina metMbH2O (24) or metHbIH2S (25) crystal structures yielded completely unacceptable fits for delta dip(obs) versus delta dip(calc) for Lucina-metMbICN in solution (as shown by open circles and half-closed circles, respectively, in Fig. 6B); the correlation is particularly poor for the labile Nepsilon Hs where both models predict even the wrong sign for delta dip(calc). Moreover, the low-field Nepsilon H closest to the iron has RFe ~3.5 Å and 5.2 Å in the two crystal structures while the observed T1 = 12 ms predicts RFe ~ 4.2 ± 0.2 Å via Equation 5. Since the Gln64(E7) orientation in Lucina metHbIH2O had the Nepsilon H oriented toward the ligated water, which is expected to be a H-bond donor, the Gln orientation with chi 3 rotated by 180° to interchange the carbonyl and amide groups was considered. This altered orientation leads to different Nepsilon H2 dipolar shifts (not shown) that fit delta dip(obs) better, but still unacceptably. Attempts to obtain a fit for delta dip(obs) by sequential chi 1, chi 2, chi 3 rotation failed since delta dip(calc) for Cbeta Hs is quite insensitive to the chi 1. Instead, the combination of the large low-field delta dip(obs) and RFe ~ 4.2 Å (from the T1 = 12 ms) was used to uniquely locate the low-field Nepsilon H in the crystal coordinates, and a search made for the position in space of the other Nepsilon H (with intra Nepsilon H distance 1.88 Å) to satisfy both the upfield delta dip(obs) and RFe ~ 4.6 ± 0.2 Å obtained from its T1 = 22 ms. Upon obtaining a reasonable fit for the Nepsilon Hs that satisfy both relaxation and delta dip(obs) constraints, a search was pursued for the range of chi 1, chi 2, and chi 3 allowed by the fixed Calpha H and NH2 positions. The good correlation between delta dip(calc) and delta dip(obs) for this optimized Gln orientation is shown in Fig. 6B (closed circles). The resulting tortional angles for Gln64(E7) are listed in Table II where they are compared with values in the crystal structure of Lucina metHbIH2O.

The orientation of Phe29(B10) in the crystal structures of either Lucina metHbIH2O (24) or metHbIH2S (25) yield very similar ring proton delta dip(calc) that correlate equally well with delta dip(obs) (open and closed squares in Fig. 6B, respectively). The crystallographic RFe(Czeta H) (5.6 Å), however, is somewhat larger than the value reflected in its T1 ~ 40 ms (RFe = 5.2 ± 0.2 Å via Equation 5. A ~5° change in chi 1 alters the delta dip(calc) inconsequentially but results in a reasonable RFe (Czeta H) ~5.3 Å. Phe68(E11) exhibits only minor delta dip(obs) and, like in the triple mutant Mb, yields anomalous Curie plots that reflect variable population of alternate orientations (50). The crystallographic orientation, however, is consistent both with delta dip(calc) (not shown) and NOESY cross-peak pattern (Fig. 1). It is noted that for Lucina metHbICN, in contrast to L29F/H64Q/V68F-metMbCN, the crystallographic distal Phe43(CD1) placement in either crystal structure yields similarly poor fits for delta dip(obs) versus delta dip(calc) (open triangles in Fig. 6B). Moreover, the Gln64(E7) orientation deduced above results in a ~0.3 Å van der Waals overlap with the crystallographically defined Phe43(CD1) ring. Moving the Phe43(CD1) ring 0.3 Å parallel to the heme so as to abolish this steric interaction leads to a much better fit between delta dip(obs) and delta dip(calc) for Phe43(CD1) (closed triangles, Fig. 6B). The positions of the distal residues for metHbIH2O in the crystal and metHbICN in solution are compared in Fig. 2B, and reveal significant differences primarily for the Gln64(E7) orientation.

    DISCUSSION
Top
Abstract
Introduction
Procedures
Results
Discussion
References

Solution NMR Cavity Structure-- The conventional two-dimensional NMR approach allowed unambiguous assignment for all relevant residues in the heme cavity. The intrinsic paramagnetism of cyano-metMb, instead of being an impediment to assignment, was actually an advantage since the dipolar shift with minimal relaxation provide significantly enhanced resolution for the E and F helical segments. For the proximal side of the heme, the NOESY cross-peak patterns and paramagnetic relaxation for the two Mb mutants were indistinguishable from those observed for wild-type metMbCN and predicted by the wild-type MbCO structure. A similarly strongly conserved proximal structure is confirmed for Lucina metHbICN. These conserved proximal structures allow the determination of meaningful magnetic axes for each protein. The robust nature of the magnetic axes for L29F/H64Q/V68F-metMbCN is confirmed by the indistinguishable magnetic axes determined using either the wild-type MbCO or triple mutant crystal structures.

For the distal residues, whose relaxation properties, dipolar shifts and/or NOESY cross-peak pattern differed from that expected for the reference crystal structure, the combined restraints for dipolar shifts and paramagnetic relaxation allow the determination of the their orientations in a more definitive manner than allowed by the NOESY cross-peak pattern themselves. The accuracy of the distal orientations is limited only by the basic assumption that the helical backbone positions are inconsequentially perturbed in the mutant and/or ligation variant relative to the available crystal structure. It is noted that the labile protons of the distal Gln64(E7) are readily assigned. The combined effect of the dipolar shift and paramagnetic relaxation allow the placement of the labile protons in the distal pocket so as to more accurately describe hydrogen bonding between Gln64(E7) and the bound ligand than is possible in some of the crystal structures.

The aromatic ring for Phe68(E11) is oriented away from the expected binding site with the ring inserting into a largely pre-existing cavity in the distal pocket (55). The absence of a methyl group at the gamma 2 position results in less direct hindrance of the bound ligand than is observed for the naturally occurring Val side chain. Thus, the Phe68(E11) mutant behaves like V68A-metMbCN; both exhibit less tilt, beta  ~ 7-9°, than that in the wild-type protein (beta  = 16°) (30, 57). Last, the labile Gln64(E7) protons are in a position to form hydrogen bonds with bound cyanide in the distal pockets of both L29F/H64Q/V68F metMbCN and Lucina metHbICN. This orientation is probably similar to that in the oxy complexes.

Comparison of Distal Pockets in Solution and Crystal Structures-- For all but select distal residues, the heme pocket structure for each protein is essentially the same in cyanomet and aquo-met derivatives. In the case of the sperm whale triple mutant, the major difference in the Gln64(E7) orientation between the cyanomet and aquomet complexes is ~180° rotation about chi 3 that interchanges the carboxy and amide termini (Table II, Fig. 2A). The 180° rotation of the Gln side chain terminus is completely consistent with the stabilizing roles of the side chain for the different ligands. For the cyano complex, the Nepsilon H2 is oriented to serve as a H-bond donor to bound cyanide, whereas in the aquomet complex, the carbonyl oxygen is oriented to accept a H-bond from coordinated water. The small changes in chi 1, chi 2 for Gln64(E7) (Table II) help to accommodate the larger cyanide ligand.

As shown in Fig. 2B, the solution orientation of Gln64(E7) in Lucina metHbCN I differs substantially from that reported for the corresponding aquomet crystal. The altered orientation is firmly and independently established by the NOESY pattern from Cgamma Hs to Phe29(B10), the dipolar shift simulation, and the paramagnetic relaxation. However, both structures place the amino group in a position to interact with the bound ligand. The electron density in the crystal structure does not allow differentiation between the terminal carbonyl O atom and the Nepsilon atom.4 However, the carbonyl O atom rather than the reported NH2 must be pointing toward the coordinated water to serve as H-bond acceptor. Donation of proton to the bound water would give the ligand a partial positive charge (i.e. H3O+ character) causing substantial electrostatic repulsion with the net +1 charge on the hemin iron atom. Thus, the solution NMR data, through the dipolar shift and paramagnetic relaxation, provide the most definitive location of the NH2 group. The small movements of Phe29(B10) and Phe43(CD1) in Lucina HbI appear to reflect minor accommodation of the larger cyanide ligand. However, the orientations of Gln64(E7) do differ in solution and crystal, as shown in Fig. 2B.

The orientation of the Gln64(E7) side chain in the sulfide complex of Lucina metHbI (not shown) is intermediate between that for the aquomet crystal and that for metHbICN in solution (Fig. 2B). The CO versus NH2 terminus orientation in the sulfide complex of Lucina metHbI also differs from that in metHbICN in solution by a 180° rotation at chi 3. Again, this difference is consistent hydrogen bond donation to bound cyanide and acceptance from bound H2S.

Implications for Sulfide Binding-- All three single distal pocket substitutions in sperm whale Mb cause an enhancement of sulfide binding but have variable effects on oxygen affinity. The key mutation is H64Q, which causes a 25-fold increase in sulfide affinity by itself. The triple mutant, L29F/H64Q/V68F-metMb, has a sulfide affinity which is ~700 times greater than that of wild-type metMb and only ~7-fold less than that of Lucina HbI. Thus the triple mutant exhibits an increased sulfide binding free energy of ~3.8 kcal/mol or about ~75% of the stabilization of the bound sulfide (~5.1 kcal/mol) by native Lucina HbI relative to that by wild-type sperm whale Mb. Hence, a major portion, but not all, of the remarkable sulfide affinity of Lucina HbI can be transferred to a mammalian Mb by replacing the three residues in contact with bound sulfide with those found in Lucina HbI. The high affinity of Lucina HbI for sulfide was attributed (24, 25) to a combination of the hydrogen-bond acceptance by the Gln64(E7) side chain carbonyl, the hydrophobic pocket provided by the "cage" generated by the three Phe side chains, B10, CD1, and E11, and highly favorable electrostatic interactions between the sulfur and the edges of the three aromatic rings. While the distal pockets for the triple mutant Mb and Lucina HbI are similar (Fig. 2), there are some significant differences that relate to their differences in affinity for sulfide.

The side chains of Gln64(E7) and Phe43(CD1) have similar orientations in the sperm whale triple mutant and the clam protein (Fig. 2C). However, both Phe29(B10) and Phe68(E11) in the sperm whale mutant are rotated ~90° about chi 2 relative to their positions in Lucina HbI and their edges are closer to the iron atom in the mollusc Hb (Fig. 2C). As a result of the differences, Lucina hemoglobin has a significantly larger ligand-binding site than that found in the myoglobin mutants, which should facilitate the binding of the large H2S ligand. The change in orientation of Phe29(B10) should have little affect on stabilization of bound sulfide since the positive edge of the ring multipole is pointing toward the bound ligand in both proteins. In Lucina HbI, the positive edge of the Phe68(E11) ring also points toward the bound ligand (Fig. 2C). This orientation of the Phe(E11) side chain is caused by direct steric interactions Phe28(B9) in Lucina HbI (25). Stabilization of bound sulfide by Phe68(E11) ring does not appear to occur in the myoglobin triple mutant since the E-11 side chain is allowed to take a position more perpendicular to the plane of the heme since the B9 position is occupied by a smaller Ile residue. A more complete and quantitative interpretation of the sulfide binding data will require detailed volume and electrostatic calculations for the entire set of mutants listed in Table I; such studies are now in progress.

    ACKNOWLEDGEMENTS

We thank Eileen Singleton for the expression and purification of the triple mutant myoglobin, Mike Berry for help in the protein's crystallization and subsequent structural refinement, and Jun Qin for assistance with NMR experiments.

    FOOTNOTES

* This work was supported by United States National Institutes of Health Grants HL16087 (to G. N. L.), GM35649 and HL47020 (to J. S. O.), AR40252 (to G. N. P.), and Postdoctoral Fellowship AR08355 (to E. A. B.), the States of Texas Advanced Technology Program Grant 003604-025 (to G. N. P. and J. S. O.), Robert A. Welch Foundation Grants C-612 (to J. S. O.) and C-1142 (to G. N. P.), and the W. M. Keck Center for Computational Biology.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Dagger To whom correspondence should be addressed.

§ Current address: Somatogen, Inc., Boulder, CO 80301.

1 The abbreviations used are: Mb, myoglobin; Hb, hemoglobin; DSS, 2,2-dimethyl-2-silapentane-5-sulfonate; MbCO, carbonmonoxymyoglobin; MCOSY, magnitude COSY; metMbCN, cyanomet myoglobin; NMR, nuclear magnetic resonance; NOE, nuclear Overhauser effect; NOESY, two-dimensional nuclear Overhauser spectroscopy; ppm, parts per million; TOCSY, two-dimensional total correlation spectroscopy; V68F-Mb, Val68(E10)right-arrowPhe-Mb; L29F/H64Q/V68F-Mb, Leu29(B10) right-arrow Phe, His64(E7) right-arrow Gln, Val68(E11) right-arrow Phe-Mb; HbI, monomeric hemoglobin component I; metHbICN, cyanomet-hemoglobin I.

2  E. A. Brucker, G. N. Phillips, R. Lile, J. S. Olson, R. F. Eich, and J. B. Wittenberg, unpublished results.

3 The delta dip(obs) used for the magnetic axes for L29F/H64Q/V68F-metMbCN and V68F-metMbCN are: Leu89(F4) Calpha H; Ala90(F5) Calpha H, Cbeta H3; Phe138(H15) Cdelta Hs, Cepsilon Hs, Czeta H; Ile99(FG5) Calpha H, Cbeta H, Cgamma H, Cgamma H', Cgamma H3, Cdelta H3; His97(FG3) Cepsilon H, Cdelta H; and those used for Lucina metHbICN are: Ala62(E5)-Val72(E15), Phe92(F4)-Ala98(F10), Arg101(FG1), Gly102(FG2), Ile104(FG4), Ala108(G7)-Phe110(G9), Val130(H11), and Ala131(H12).

4 M. Bolognesi, personal communication.

    REFERENCES
Top
Abstract
Introduction
Procedures
Results
Discussion
References

  1. Lesk, A. M., and Chothia, C. (1980) J. Mol. Biol. 136, 225-270[Medline] [Order article via Infotrieve]
  2. Bashford, D., Chothia, C., and Lesk, A. M. (1987) J. Mol. Biol. 196, 199-216[Medline] [Order article via Infotrieve]
  3. Vinogradov, S. N., Walz, D. A., Pohajdak, B., Moens, L., Kapp, O. H., Susuzki, T., and Trotman, C. N. A. (1993) Comp. Biochem. Physiol. 106, 1-26[CrossRef]
  4. Romero-Herrera, A. E., Goodman, M., Dene, H., Bartiniki, D. E., and Mizukami, H. (1981) J. Mol. Evol. 17, 140-147[Medline] [Order article via Infotrieve]
  5. Suzuki, T., Muramatsu, R., Kisamori, T., and Furukohri, T. (1988) Zool. Sci. 5, 69-76
  6. Liljeqvist, G., Paléus, S., and Braunitzer, G. (1982) J. Mol. Evol. 18, 102-108[Medline] [Order article via Infotrieve]
  7. Stenzel, P., Brimhall, B., Jones, R. T., Black, J. A., McLachlan, A., and Gibson, D. (1979) J. Biol. Chem. 254, 2071-2076[Medline] [Order article via Infotrieve]
  8. Bolognesi, M., Onesti, S., Gatti, G., Coda, A., Ascenzi, P., Giacometti, A., and Brunori, M. (1989) J. Mol. Biol. 205, 529-544[CrossRef][Medline] [Order article via Infotrieve]
  9. Suzuki, T., Furukohri, T., and Okamoto, S. (1993) J. Protein Chem. 12, 45-50[Medline] [Order article via Infotrieve]
  10. Rashid, K. A., Van Hauwaert, M.-L., Haque, M., Siddiqi, A. H., Lasters, I., De Mayer, M., Griffon, N., Marden, M. C., Dewilde, S., Clauwaert, J., Vinogradov, N., and Moens, L. (1997) J. Biol. Chem. 272, 2992-2999[Abstract/Free Full Text]
  11. Zhang, W., Rashid, K. A., Haque, M., Siddiqi, A. H., Vinogradov, S. N., Moens, L., and La Mar, G. N. (1997) J. Biol. Chem. 272, 3000-3006[Abstract/Free Full Text]
  12. Arents, G., and Love, W. E. (1989) J. Mol. Biol. 210, 149-161[Medline] [Order article via Infotrieve]
  13. Travaglini Allocatelli, C., Cutruzzolá, F., Brancaccio, A., Brunori, M., Qin, J., and La Mar, G. N. (1993) Biochemistry 32, 6041-6049[Medline] [Order article via Infotrieve]
  14. Vyas, K., Rajarathnam, K., Yu, L. P., Emerson, S. D., La Mar, G. N., Krishnamoorthi, R., and Mizukami, H. (1993) J. Biol. Chem. 268, 14826-14835[Abstract/Free Full Text]
  15. Qin, J., La Mar, G. N., Ascoli, F., and Brunori, M. (1993) J. Mol. Biol. 231, 2009-1023
  16. Qin, J., La Mar, G. N., Cutruzzolá, F., Travaglini Allocatelli, C., Brancaccio, A., and Brunori, M. (1993) Biophys. J. 65, 2178-2190[Abstract]
  17. Zhao, X., Vyas, K., Nguyen, B. D., Rajarathnam, K., La Mar, G. N., Li, T., Phillips, G. N., Jr., Eich, R. F., Olson, J. S., Ling, J., and Bocian, D. F. (1995) J. Biol. Chem. 270, 20763-20774[Abstract/Free Full Text]
  18. Wittenberg, J. (1985) Bull. Biol. Soc. Wash. 6, 301-310
  19. Wittenberg, J. B., and Wittenberg, B. A. (1990) Annu. Rev. Biochem. Biophys. Commun. 19, 217-241
  20. Wittenberg, J. B., and Kraus, D. W. (1991) in Structure and Function of Invertebrate Oxygen Carriers (Vinogradov, S. N., and Kapp, O. H., eds), pp. 323-330, Springer Verlag, New York
  21. Fisher, C. R. (1990) Revs. Aquat. Sci. 2, 399-446
  22. Childress, J. J., and Fisher, C. R. (1992) Oceanogr. Marine Biol. Annu. Rev. 30, 337-441
  23. Kraus, D. W., and Wittenberg, J. B. (1990) J. Biol. Chem. 265, 16043-16053[Abstract/Free Full Text]
  24. Rizzi, M., Wittenberg, J. B., Coda, A., Fasano, M., Ascenzi, P., and Bolognesi, M. (1994) J. Mol. Biol. 244, 86-89[CrossRef][Medline] [Order article via Infotrieve]
  25. Rizzi, M., Wittenberg, J. B., Coda, A., Ascenzi, P., and Bolognesi, M. (1996) J. Mol. Biol. 258, 1-5[CrossRef][Medline] [Order article via Infotrieve]
  26. Zhang, W., Cutruzzolá, F., Travaglini Allocatelli, C., Brunori, M., and La Mar, G. N. (1997) Biophys. J. 73, 1019-1030[Abstract]
  27. Bisig, D. A., Di Iorio, E. E., Diederichs, K., Winterhalter, K. H., and Pionteck, K. (1995) J. Biol. Chem. 270, 20754-20762[Abstract/Free Full Text]
  28. Springer, B. A., Sligar, S. G., Olson, J. S., and Phillips, G. N., Jr. (1994) Chem. Rev. 94, 699-714
  29. Rajarathnam, K., La Mar, G. N., Chiu, M., and Sligar, S. G. (1992) J. Am. Chem. Soc. 114, 9048-9058
  30. Rajarathnam, K., Qin, J., La Mar, G. N., Chiu, M. L., and Sligar, S. G. (1993) Biochemistry 32, 5670-5680[Medline] [Order article via Infotrieve]
  31. Egeberg, K. D., Springer, B. A., Sligar, S. G., Carver, T. E., Rohlfs, R. J., and Olson, J. S. (1990) J. Biol. Chem. 265, 11788-11795[Abstract/Free Full Text]
  32. Springer, B. A., Egeberg, K. D., Sligar, S. G., Rohlfs, R. J., Mathews, A. J., and Olson, J. S. (1989) J. Biol. Chem. 264, 3057-3060[Abstract/Free Full Text]
  33. Carver, T. E., Brantley, R. E., Jr., Singleton, E. W., Arduini, R. M., Quillin, M. L., Phillips, G. N., Jr., and Olson, J. S. (1992) J. Biol. Chem. 267, 14443-14450[Abstract/Free Full Text]
  34. Springer, B. A., and Sligar, S. G. (1987) Proc. Natl. Acad. Sci. U. S. A. 84, 8961-8965[Abstract]
  35. Phillips, G. N., Jr., Arduini, R. M., Springer, B. A., and Sligar, S. G. (1990) Proteins Struct. Funct. Genet. 7, 358-365[Medline] [Order article via Infotrieve]
  36. Quillin, M. L., Arduini, R. M., Olson, J. S., and Phillips, G. N., Jr. (1993) J. Mol. Biol. 234, 140-155[CrossRef][Medline] [Order article via Infotrieve]
  37. Brünger, A. T., Kuriyan, J., and Karplus, M. (1987) Science 235, 458-460
  38. Brünger, A. T., Karplus, M., and Petsko, G. A. (1989) Acta Crystallogr. Sec. A 45, 50-61[CrossRef]
  39. Brünger, A. T. (1992) X-PLOR Version 3.1 Manual, Yale University
  40. Sack, J. S. (1988) J. Mol. Graph. 6, 244-245
  41. Engh, R. A., and Huber, R. (1991) Acta Crystallogr. Sec. A 47, 392-400[CrossRef]
  42. Gupta, R. K. (1976) J. Magn. Res. 24, 461-465
  43. Emerson, S. D., and La Mar, G. N. (1990) Biochemistry 29, 1545-1556[Medline] [Order article via Infotrieve]
  44. Braunschweiler, L., and Ernst, R. R. (1983) J. Magn. Reson. 53, 521-528
  45. Davis, D. G., and Bax, A. (1985) J. Am. Chem. Soc. 107, 2820-2821
  46. Jeener, J., Meier, B. H., Bachmann, P., and Ernst, R. R. (1979) J. Chem. Phys. 71, 4546-4553[CrossRef]
  47. Bax, A. (1982) Two-dimensional Nuclear Magnetic Resonance in Liquids, Delft University Press, Dordrecht, Holland
  48. States, D. J., Haberkorn, R. A., and Reuben, D. J. (1982) J. Magn. Reson. 48, 286-292
  49. Qin, J., and La Mar, G. N. (1992) J. Biomol. NMR 2, 597-618[Medline] [Order article via Infotrieve]
  50. Emerson, S. D., and La Mar, G. N. (1990) Biochemistry 29, 1556-1566[Medline] [Order article via Infotrieve]
  51. Kuriyan, J., Witz, S., Karplus, M., and Petsko, G. A. (1986) J. Mol. Biol. 192, 133-154[Medline] [Order article via Infotrieve]
  52. Dalvit, C. G., and Wright, P. E. (1987) J. Mol. Biol. 194, 313-327[Medline] [Order article via Infotrieve]
  53. Chiu, M. L. (1992) Structure and Dynamics of Myoglobin and its MutantsPh.D. thesis, University of Illinois, Urbana, IL
  54. Zhang, W., La Mar, G. N., and Gersonde, K. (1996) Eur. J. Biochem. 237, 841-853[Abstract]
  55. Quillin, M. L., Li, T., Olson, J. S., Phillips, G. N., Jr., Dou, Y., Ikeda-Saito, M., Regan, R., and Gibson, Q. H. (1995) J. Mol. Biol. 245, 416-436[CrossRef][Medline] [Order article via Infotrieve]
  56. Carver, T. E., Rohlfs, R. J., Olson, J. S., Gibson, Q. H., Blackmore, R. S., Springer, B. A., and Sligar, S. G. (1990) J. Biol. Chem. 265, 20007-20020[Abstract/Free Full Text]
  57. Rajarathnam, K., Qin, J., La Mar, G. N., Chiu, M. L., and Sligar, S. G. (1994) Biochemistry 33, 5493-5501[Medline] [Order article via Infotrieve]


Copyright © 1998 by The American Society for Biochemistry and Molecular Biology, Inc.