Assignment of the ferriheme resonances of the high-spin forms of nitrophorins 1 and 4 by 1H NMR spectroscopy: comparison to structural data obtained from X-ray crystallography

Dimerization of nitrophorin 4 at low pH and comparison to the K1A mutant of nitrophorin 1

Nitrophorin 4, one of the four NO-carrying heme proteins from the salivary glands of Rhodnius prolixus, forms a homodimer at pH 5.0 with a Kd of ∼8 μM. This dimer begins to dissociate at pH 5.5 and is completely dissociated to monomer at pH 7.3, even at 3.7 mM. The dimer is significantly stabilized by binding NO to the heme and at pH 7.3 would require dilution to well below 0.2 mM to completely dissociate the NP4-NO homodimer. The primary techniques used for investigating the homodimer and the monomer-dimer equilibrium were size-exclusion fast protein liquid chromatography at pH 5.0 and (1)H{(15)N} heteronuclear single-quantum coherence spectroscopy as a function of pH and concentration. Preparation of site-directed mutants of NP4 (A1K, D30A, D30N, V36A/D129A/L130A, K38A, R39A, K125A, K125E, D132A, L133V, and K38Q/R39Q/K125Q) showed that the N-terminus, D30, D129, D132, at least one heme propionate, and, by association, likely also E32 and D35 are involved in the dimerization. The "closed loop" form of the A-B and G-H flexible loops of monomeric NP4, which predominates in crystal structures of the monomeric protein reported at pH 5.6 but not at pH 7.5 and which involves all of the residues listed above except D132, is required for dimer formation. Wild-type NP1 does not form a homodimer, but NP1(K1A) and native N-terminal NP1 form dimers in the presence of NO. The homodimer of NP1, however, is considerably less stable than that of NP4 in the absence of NO. This suggests that additional aspartate or glutamate residues present in the C-terminal region of NP4, but not NP1, are also involved in stabilizing the dimer.

NMR investigations of nitrophorin 2 belt side chain effects on heme orientation and seating of native N-terminus NP2 and NP2(D1A)

Nitrophorin 2 (NP2), one of the four NO-storing and NO-releasing proteins found in the saliva of the blood-sucking bug Rhodnius prolixus, has a more ruffled heme and a high preference for a particular heme orientation (B) compared with nitrophorin 1 and nitrophorin 4, which show not a preference (A to B ratio of approximately 1:1), suggesting that it fits more tightly in the β-barrel protein. In this work we have prepared a series of "belt" mutants of NP2(D1A) and (ΔM0)NP2 aimed at reducing the size of aromatic or other residues that surround the heme, and investigated them as the high-spin aqua and low-spin N-methylimidazole complexes. The belt mutants included Y38A, Y38F, F42A, F66A, Y85A, Y85F, Y104A, I120T, and a triple mutant of NP2(D1A), the F42L, L106F, I120T mutant. Although I120 has been mainly considered to be a distal pocket residue, CδH3 of I120 lies directly above the heme 3-methyl, at 2.67 Å, of heme orientation B, or the 2-vinyl of A, and it thus plays a role as a belt mutant, a role that turns out to be extremely important in creating the strong favoring of the B heme orientation [A to B ratio of 1:14 for NP2(D1A) or 1:12 for (ΔM0)NP2]. The results show that the 1D (1)H NMR spectra of the high-spin forms are quite sensitive to changes in the shape of the heme binding cavity. The single mutation I120T eliminates the favorability of the B heme orientation by producing a heme A to B orientation ratio of 1:1, whereas the single mutation F42A reverses the heme orientation from an A to B ratio of 1:14 seen for NP2(D1A) to 10:1 for NP2(D1A,F42A). The most extreme ratio was found for the triple mutant of NP2(D1A), NP2(D1A,F42L,L105F,I120T), in which the A to B ratio is approximately 25:1, a ΔG change of about -3.5 kcal/mol or -14.1 kJ/mol with respect to NP2(D1A). The seating of the heme is modified as well in that mutant and in several others, by rotations of the heme by up to 4° from the seating observed in NP2(D1A), in order to relieve steric interactions between a vinyl β-carbon and a protein side chain, or to fill a cavity created by replacing a large protein side chain by a much smaller one; the latter was observed for all tyrosine to alanine mutants. These relatively small changes in seating have a measurable effect on the NMR spectra of the mutants, but are indeed minor in terms of overall seating and reactivity of the NP2(D1A) protein. The (1)H NMR resonances of the hemin substituents of the low-spin N-methylimidazole complexes of NP2(D1A,F42L,L105F,I120T) as well as NP2(D1A,I120T), NP2(D1A,Y104A), and NP2(D1A,F42A) have been assigned using natural abundance (1)H{(13)C} heteronuclear multiple quantum correlation and (1)H-(1)H nuclear Overhauser effect spectroscopy spectra.

Expression, purification, and solid-state NMR characterization of the membrane binding heme protein nitrophorin 7 in two electronic spin states

The nitrophorins (NPs) comprise a group of NO transporting ferriheme b proteins found in the saliva of the blood sucking insect Rhodnius prolixus . In contrast to other nitrophorins (NP1-4), the recently identified membrane binding isoform NP7 tends to form oligomers and precipitates at higher concentrations in solution. Hence, solid-state NMR (ssNMR) was employed as an alternative method to gain structural insights on the precipitated protein. We report the expression and purification of (13)C,(15)N isotopically labeled protein together with the first ssNMR characterization of NP7. Because the size of NP7 (21 kDa) still provides a challenge for ssNMR, the samples were reverse labeled with Lys and Val to reduce the number of crosspeaks in two-dimensional spectra. The two electronic spin states with S = 1/2 and S = 0 at the ferriheme iron were generated by the complexation with imidazole and NO, respectively. ssNMR spectra of both forms are well resolved, which allows for sequential resonance assignments of 22 residues. Importantly, the ssNMR spectra demonstrate that aggregation does not affect the protein fold. Comparison of the spectra of the two electronic spin states allows the determination of paramagnetically shifted cross peaks due to pseudocontact shifts, which assists the assignment of residues close to the heme center.

Electron spin density on the axial His ligand of high-spin and low-spin nitrophorin 2 probed by heteronuclear NMR spectroscopy

The electronic structure of heme proteins is exquisitely tuned by the interaction of the iron center with the axial ligands. NMR studies of paramagnetic heme systems have been focused on the heme signals, but signals from the axial ligands have been rather difficult to detect and assign. We report an extensive assignment of the (1)H, (13)C and (15)N resonances of the axial His ligand in the NO-carrying protein nitrophorin 2 (NP2) in the paramagnetic high-spin and low-spin forms, as well as in the diamagnetic NO complex. We find that the high-spin protein has σ spin delocalization to all atoms in the axial His57, which decreases in size as the number of bonds between Fe(III) and the atom in question increases, except that within the His57 imidazole ring the contact shifts are a balance between positive σ and negative π contributions. In contrast, the low-spin protein has π spin delocalization to all atoms of the imidazole ring. Our strategy, adequately combined with a selective residue labeling scheme, represents a straightforward characterization of the electron spin density in heme axial ligands.

Native N-terminus nitrophorin 2 from the kissing bug: similarities to and differences from NP2(D1A)

The first amino acid of mature native nitrophorin 2 is aspartic acid, and when expressed in E. coli, the wild-type gene of the mature protein retains the methionine-0, which is produced by translation of the start codon. This form of NP2, (M0)NP2, has been found to have different properties from its D1A mutant, for which the Met0 is cleaved by the methionine aminopeptidase of E. coli (R. E. Berry, T. K. Shokhireva, I. Filippov, M. N. Shokhirev, H. Zhang, F. A. Walker, Biochemistry 2007, 46, 6830). Native N-terminus nitrophorin 2 ((ΔM0)NP2) has been prepared by employing periplasmic expression of NP2 in E. coli using the pelB leader sequence from Erwinia carotovora, which is present in the pET-26b expression plasmid (Novagen). This paper details the similarities and differences between the three different N-terminal forms of nitrophorin 2, (M0)NP2, NP2(D1A), and (ΔM0)NP2. It is found that the NMR spectra of high- and low-spin (ΔM0)NP2 are essentially identical to those of NP2(D1A), but the rate and equilibrium constants for histamine and NO dissociation/association of the two are different.

Assignment of the 1H NMR resonances of protein residues in close proximity to the heme of the nitrophorins: similarities and differences among the four proteins from the saliva of the adult blood-sucking insect Rhodnius prolixus

The nuclear Overhauser effects (NOEs) observed between heme substituent protons and a small number of nearby protein side chain protons in the water-elimination Fourier transform NOE spectroscopy (WEFT-NOESY) spectra of high- and low-spin wild-type nitrophorin (NP) 2 and its ligand complexes have been analyzed and compared with those observed for the same complexes of wild-type NP3. These assignments were made on naturally abundant isotope samples, with the most useful protein side chains being those of Ile120, Leu122, and Leu132 for NP2 and NP3, and Thr121, Leu123, and Leu133 for NP1 and NP4. It is found that the NOEs observed are identical, with extremely similar protein side chain proton chemical shifts. This is strong evidence that the structure of NP3, for which no X-ray crystal structures are available, is essentially identical to that of NP2, at least near the heme binding pocket. Similarly, the NOEs observed between heme substituents and protein side chains for NP1 and NP4 also indicate that the structures of the protein having both A and B heme orientations are very similar to each other, as well as to the proteins with major B heme orientation of NP2 and NP3. These A and B connectivities can be seen, even though the two heme orientations have similar populations in NP1 and NP4, which complicates the analysis of the NOESY spectra. The histamine complex of wild-type NP2 shows significant shifts of the Leu132 side chain protons relative to all other ligand complexes of NP1-NP4 because of the perturbation of the structure near Leu132 caused by the histamine's side chain ammonium hydrogen bond to the Asp29 side chain carboxylate.

NMR studies of nitrophorin distal pocket side chain effects on the heme orientation and seating of NP2 as compared to NP1

The nitrophorins (NP) of the adult blood-sucking insect Rhodnius prolixus fall into two pairs based on sequence identity (NP1,4 (90%) and NP2,3 (79%)), which differ significantly in the size of side chains of residues which contact the heme. These residues include those in the distal pocket of NP2 (I120) and NP1 (T121) and the "belt" that surrounds the heme of NP2 (S40, F42), and NP1(A42, L44). To determine the importance of these residues and others conserved or very similar for the two pairs, including L122(123), L132(133), appropriate mutants of NP2 and NP1 have been prepared and studied by (1)H NMR spectroscopy. Wild-type NP2 has heme orientation ratio (A:B) of 1:8 at equilibrium, while wild-type NP1 has A:B ~1:1 at equilibrium. Another difference between NP2 and NP1 is in the heme seating with regard to His57(59). It is found that among the distal pocket residues investigated, the residue most responsible for heme orientation and seating is I120(T121). F42(L44) and L106(F107) may also be important, but must be investigated in greater detail.

Ligand orientation of human neuroglobin obtained from solution NMR and molecular dynamics simulation as compared with X-ray crystallography

Neuroglobin, a new member of hemoprotein family, can reversibly bind oxygen and take part in many biological processes such as enzymatic reaction, signal transduction and the mitochondria function. Different from myoglobin and hemoglobin, it has a hexacoordinated heme environment, with histidyl imidazole of proximal His(96)(F8) and distal His(64)(E7) directly bound to the metal ion. In the present work, solution (1)H NMR spectroscopy was employed to investigate the electronic structure of heme center of wild-type met-human neuroglobin. The resonances of heme protons and key residues in the heme pocket were assigned. Two heme orientations resulting from a 180 degrees rotation about the alpha-gamma-meso axis with a population ratio about 2:1 were observed. Then the (1)H NMR chemical shifts of the ferriheme methyl groups were used to predict orientations of the axial ligand. The obtained axial ligand plane angle phi is consistent with that from the molecular dynamics simulation but not with those from the crystal data. Compared with mouse neuroglobin, the obtained average ligand orientation of human neuroglobin reflects the changeability of heme environment for the Ngb family.

1H and 13C NMR spectroscopic studies of the ferriheme resonances of three low-spin complexes of wild-type nitrophorin 2 and nitrophorin 2(V24E) as a function of pH

The ferriheme resonances of the low-spin (S = 1/2) complexes of wild-type (wt) nitrophorin 2 (NP2) and its heme pocket mutant NP2(V24E) with imidazole (ImH), histamine (Hm), and cyanide (CN⁻) as the sixth ligand have been investigated by NMR spectroscopy as a function of pH (4.0–7.5). For the three wt NP2 complexes, the ratio of the two possible heme orientational isomers, A and B, remains almost unchanged (ratio of A:B approximately 1:6 to 1:5) over this wide pH range. However, strong chemical exchange cross peaks appear in the nuclear Overhauser effect spectroscopy/exchange spectroscopy (NOESY/EXSY) spectra for the heme methyl resonances at low pH (pH* 4.0–5.5), which indicate chemical exchange between two species. We have shown these to be two different exogenous ImH or Hm orientations that are denoted B and B′, with the ImH plane nearly parallel and perpendicular to the ImH plane of the protein-provided His57, respectively. The wt NP2–CN complex also shows EXSY cross peaks due to chemical exchange, which is shown to be a result of interchange between two ruffling distortions of the heme. The same ruffling distortion interchange is also responsible for the ImH and Hm chemical exchange. For the three NP2(V24E) ligand complexes, no EXSY cross peaks are observed, but the A:B ratios change dramatically with pH. The fact that heme favors the A orientation highly for NP2(V24E) at low pH as compared with wt NP2 is believed to be due to the steric effect of the V24E mutation. The existence of the B′ species at lower pH for wt NP2 complexes and the increase in A heme orientation at lower pH for NP2(V24E) are believed to be a result of a change in structure near Glu53 when it is protonated at low pH. ¹H{¹³C} heteronuclear multiple quantum coherence (HMQC) spectra are very helpful for the assignment of heme and nearby protein side chain resonances.

Effect of mutation of carboxyl side-chain amino acids near the heme on the midpoint potentials and ligand binding constants of nitrophorin 2 and its NO, histamine, and imidazole complexes

Nitrophorins (NPs) are a group of NO-carrying heme proteins found in the saliva of a blood-sucking insect from tropical Central and South America, Rhodnius prolixus, the "kissing bug". NO is kept stable for long periods of time by binding it as an axial ligand to a ferriheme center. The fact that the nitrophorins are stabilized as Fe(III)-NO proteins is a unique property because most heme proteins are readily autoreduced by excess NO and bind NO to the Fe(II) heme irreversibly (K(d)s in the picomolar range). In contrast, the nitrophorins, as Fe(III) heme centers, have K(d)s in the micromolar to nanomolar range and thus allow NO to dissociate upon dilution following injection into the tissues of the victim. This NO can cause vasodilation and thereby allow more blood to be transported to the site of the wound. We prepared 13 site-directed mutants of three major nitrophorins, NP2, NP1, and NP4, to investigate the stabilization of the ferric-NO heme center and preservation of reversible binding that facilitates these proteins' NO storage, transport, and release functions. Of the mutations in which Glu and/or Asp were replaced by Ala, most of these carboxyls show a significant role stabilizing Fe(III)-NO over Fe(II)-NO, with buried E53 of NP2 or E55 of NP1 and NP4 being the most important and partially buried D29 of NP2 or D30 of NP4 being second in importance. The pK(a)s of the carboxyl groups studied vary significantly but all are largely deprotonated at pH 7.5 except E124.

Assignment of the ferriheme resonances of high- and low-spin forms of the symmetrical hemin-reconstituted nitrophorins 1-4 by 1H and 13C NMR spectroscopy: the dynamics of heme ruffling deformations

The four major nitrophorins (NPs) of the adult blood-sucking insect Rhodnius prolixus have been reconstituted with the "symmetrical hemin" 2,4-dimethyldeuterohemin, and their NMR spectra have been investigated as the high-spin (S=5/2) aqua and low-spin (S=1/2) N-methylimidazole (NMeIm) and cyanide complexes. The NMeIm complexes allow assignment of the high-spin hemin resonances by saturation transfer difference spectroscopy. The cyanide complexes were investigated as paramagnetic analogues of the NO complexes. It is shown that the hemin ring is highly distorted from planarity, much more so for NP2 than for NP1 and NP4 (with ruffling being the major distortion mode), for both high- and low-spin forms. For the cyanide complexes, the conformation of the distorted ring changes on the NMR timescale to yield chemical exchange (exchange spectroscopy, EXSY) cross peaks for NP1sym(CN), NP3sym(CN) and NP4sym(CN) but not for NP2sym(CN). These changes in nonplanar conformation are visualized as a "rolling" of the ruffled macrocycle ridges through some number of degrees, the lowest-energy ruffling mode. This probably occurs in response to slow protein dynamics that cause the I120 and L132 side chains in the distal heme pocket to move in opposite directions (up and away vs. down and toward the hemin ring). This in turn changes the out-of-plane displacements of the 2M and 3M of the symmetrical hemin on the NMR timescale. Two other types of dynamics, i.e., changes in heme seating and NMeIm rotation, are also observed. The highly distorted heme and the dynamics it causes are unique to the NPs and a few other heme proteins with highly distorted macrocycles.

Assignment of Ferriheme Resonances for High- and Low-Spin Forms of Nitrophorin 3 by H and C NMR Spectroscopy and Comparison to Nitrophorin 2: Heme Pocket Structural Similarities and Differences

Nitrophorin 3 (NP3) is the only one of the four major NO-binding heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus (also called the Kissing Bug) for which it has not been possible to obtain crystals of diffraction quality for structure determination by X-ray crystallography. Thus we have used NMR spectroscopy, mainly of the hyperfine-shifted ferriheme substituent resonances, to learn about the similarities and differences in the heme pocket and the iron active site of NP3 as compared to NP2, which has previously been well-characterized by both X-ray crystallography and NMR spectroscopy. Only one residue in the heme pocket differs between the two, F27 of NP2 is Y27 for NP3; in both cases this residue is expected to interact strongly with the 2-vinyl side chain of the B heme rotational isomer or the 4-vinyl of the A heme rotational isomer. Both the high-spin (S = 5/2) aquo complex, NP3-H(2)O, and the low-spin (S = 1/2) N-methylimidazole (NMeIm) complex of NP3 have been studied. It is found that the chemical shifts of the protons of both forms are similar to those of the corresponding NP2 complexes, but with minor differences that indicate a slightly different angle for the proximal histidine (H57) ligand plane. The B heme rotational isomer is preferred by both NP3 and NP2 in both spin states, but to a greater extent when phenylalanine is present at position 27 (A:B = 1:8 for NP2, 1:6 for NP3-Y27F, 1:4 for NP3, and 1:3 for NP2-F27Y). Careful analysis of the 5Me and 8Me shifts of the A and B isomers of the two high-spin nitrophorins leads to the conclusion that the heme environment for the two isomers differs in some way that cannot be explained at the present time. The kinetics of deprotonation of the high-spin complexes of NP2 and NP3 are very different, with NP2 giving well-resolved high-spin aquo and "low-spin" hydroxo proton NMR spectra until close to the end of the titration, while NP3 exhibits broadened (1)H NMR spectra indicative of an intermediate rate of exchange on the NMR timescale between the two forms throughout the titration. The heme methyl shifts of NP2-OH are similar in magnitude and spread to those of NP2-CN, while those of metmyoglobin-hydroxo complexes are much larger in magnitude but not spread. It is concluded that the hydroxo complex of NP2 is likely S = 1/2 with a mixed (d(XY))(2)(d(XZ),d(YZ))(3)/(d(xy))(1)(d(xz),d(yz))(4) electron configuration, while those of met-Mb-OH are likely S = 1/2,3/2 mixed spin systems.

Effect of the N-terminus on heme cavity structure, ligand equilibrium, rate constants, and reduction potentials of nitrophorin 2 from Rhodnius prolixus

The D1A mutant of recombinant NP2 has been prepared and shown to have the expression-initiation methionine-0 cleaved during expression in E. coli, as is the case for recombinant NP4, where Ala is the first amino acid for the recombinant protein as well as for the mature native protein. The heme substituent 1H NMR chemical shifts of NP2-D1A and those of its imidazole, N-methylimidazole, and cyanide complexes are rather different from those of NP2-M0D1. This difference is likely due to the much smaller size of the N-terminal amino acid (A) of NP2-D1A, which allows the formation of the closed loop form of this protein, as it does for NP4 (Weichsel, A., Andersen, J. F., Roberts, S. A., and Montfort, W. R. (2000) Nature Struct. Biol. 7, 551-554). The ratio of the two hemin rotational isomers A and B is different for the two proteins, and the rate at which the A:B ratio reaches equilibrium is strikingly different (NP2-M0D1 t1/2 for heme rotation approximately 2 h, NP2-D1A t1/2 approximately 43 h). This difference is consistent with the high stability of the closed loop form of the NP2-D1A protein and infrequent opening of the loops that could allow heme to at least partially exit the binding pocket in order to rotate about its alpha,gamma-meso axis. Consistent with this, the rates of histamine binding and release to/from NP2-D1A are significantly slower than those for NP2-M0D1 at pH 7.5. This work suggests that care must be taken in interpreting data obtained from proteins that carry the expression-initiation M0.

Assignment of the ferriheme resonances of the low-spin complexes of nitrophorins 1 and 4 by (1)H and (13)C NMR spectroscopy: comparison to structural data obtained from X-ray crystallography

In this work we report the assignment of the majority of the ferriheme resonances of low-spin nitrophorins (NP) 1 and 4 and compare them to those of NP2, published previously. It is found that the structure of the ferriheme complexes of NP1 and NP4, in terms of the orientation of the ligand(s), can be determined with good accuracy by NMR techniques in the low-spin forms and that angle plots proposed previously (Shokhirev, N. V.; Walker, F. A. J. Biol. Inorg. Chem. 1998, 3, 581-594) describe the angle of the effective nodal plane of the axial ligands in solution. The effective nodal plane of low-spin NP1, NP4, and NP2 complexes is in all cases of imidazole and histamine complexes quite similar to the average of the His-59 or -57 and the exogenous ligand angles seen in the X-ray crystal structures. For the cyanide complexes of the nitrophorins, however, the effective nodal plane of the axial ligand does not coincide with the actual histidine-imidazole plane orientation. This appears to be a result of the contribution of an additional source of asymmetry, the orientation of one of the zero-ruffling lines of the heme. Probably this effect exists for the imidazole and histamine complexes as well, but because the effect of asymmetry that occurs from planar exogenous axial ligands is much larger than the effect of heme ruffling the effect of the zero-ruffling line can only be detected for the cyanide complexes, where the only ligand plane is that of the proximal histidine. The three-dimensional structures of the three NP-CN complexes, including that of NP2-CN reported herein, confirm the high degree of ruffling of these complexes. There is an equilibrium between the two heme orientations (A and B) that depends on the heme cavity shape and changes somewhat with exogenous axial ligand. The A:B ratio can be much more accurately measured by NMR spectroscopy than by X-ray crystallography.

Assignment of the ferriheme resonances of the high-spin forms of nitrophorins 1 and 4 by 1H NMR spectroscopy: comparison to structural data obtained from X-ray crystallography

In this work, we report the assignment of the majority of the ferriheme resonances of high-spin nitrophorins (NPs) 1 and 4 and compare them to those of NP2, published previously. It is found that the structures of the ferriheme complexes of NP1 and NP4, in terms of the orientation of the histidine imidazole ligand, can be described with good accuracy by NMR techniques and that the angle plot proposed previously for the high-spin form of the NPs (Shokhireva, T. Kh.; Shokhirev, N. V.; Walker, F. A. Biochemistry 2003, 42, 679-693) describes the angle of the effective nodal plane of the axial histidine imidazole in solution. There is an equilibrium between the two heme orientations (A and B), which depends on the heme cavity shape, which can be altered by mutation of amino acids with side chains (phenyl vs tyrosyl) near the potential position where a heme vinyl group would be in one of the isomers. The A:B ratio can be much more accurately measured by NMR spectroscopy than by X-ray crystallography.