Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Development of Nonheme {FeNO}7 Complexes Based on the Pyrococcus furiosus Rubredoxin for Red-Light-Controllable Nitric Oxide Release

Nitric oxide (NO) is an essential biological messenger, contributing a significant role in a diverse range of physiological processes. The light-controllable NO releasers are of great interest because of their potential as agents for NO-related research and therapeutics. Herein, we developed a pair of red-light-controllable NO releasers, pfRd-C9A-{FeNO}⁷ and pfRd-C42A-{FeNO}⁷ (pfRd = Pyrococcus furiosus rubredoxin), by constructing a nonheme {FeNO}⁷ center within the redesigned iron-sulfur protein scaffolds. While shown to be both air and thermally stable, these complexes are highly sensitive to red-light irradiation with temporal precision, which was confirmed by electron paramagnetic resonance spin trapping and Griess assay. The temporally controlled NO release from these complexes was also demonstrated in DNA cleavage assay. Overall, this study demonstrates that such a protein-based nonheme iron nitrosyl system could be a viable chemical tool for precise NO administration.

Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin

The recently identified nonsymbiotic hemoglobin gene MtGlb1-2 of the legume Medicago truncatula possesses unique properties as it generates four alternative splice forms encoding proteins with one or two heme domains. Here we investigate the ligand binding kinetics of MtGlb1-2.1 and MtGlb1-2.4, bearing two hemes and one heme, respectively. Unexpectedly, the overall time-course of ligand rebinding was unusually fast. Thus, we complemented nanosecond laser flash photolysis kinetics with data collected with a hybrid femtosecond–nanosecond pump–probe setup. Most photodissociated ligands are rebound geminately within a few nanoseconds, which leads to rates of the bimolecular rebinding to pentacoordinate species in the 10⁸ M⁻¹s⁻¹ range. Binding of the distal histidine to the heme competes with CO rebinding with extremely high rates (k_h ~ 10⁵ s⁻¹). Histidine dissociation from the heme occurs with comparable rates, thus resulting in moderate equilibrium binding constants (K_H ~ 1). The rate constants for ligation and deligation of distal histidine to the heme are the highest reported for any plant or vertebrate globin. The combination of microscopic rates results in unusually high overall ligand binding rate constants, a fact that contributes to explaining at the mechanistic level the extremely high reactivity of these proteins toward the physiological ligands oxygen, nitric oxide and nitrite.

Iron transitions during activation of allosteric heme proteins in cell signaling

Allosteric heme proteins can fulfill a very large number of different functions thanks to the remarkable chemical versatility of heme through the entire living kingdom. Their efficacy resides in the ability of heme to transmit both iron coordination changes and iron redox state changes to the protein structure. Besides the properties of iron, proteins may impose a particular heme geometry leading to distortion, which allows selection or modulation of the electronic properties of heme. This review focusses on the mechanisms of allosteric protein activation triggered by heme coordination changes following diatomic binding to proteins as diverse as the human NO-receptor, cytochromes, NO-transporters and sensors, and a heme-activated potassium channel. It describes at the molecular level the chemical capabilities of heme to achieve very different tasks and emphasizes how the properties of heme are determined by the protein structure. Particularly, this reviews aims at giving an overview of the exquisite adaptability of heme, from bacteria to mammals.

Adiabatic Ligand Binding in Heme Proteins: Ultrafast Kinetics of Methionine Rebinding in Ferrous Cytochrome c

The dynamics of methionine geminate recombination following photodissociation in ferrous cytochrome c is investigated over a broad temperature range. The kinetic response, above the solvent glass transition ( T_g), is nearly monoexponential and displays a weak temperature dependence. Below T_g, the rebinding kinetics are nonexponential and can be explained using a quenched distribution of enthalpic rebinding barriers, arising from a relatively narrow distribution of heme out-of-plane displacements. The Arrhenius prefactor of this (Δ S = 2) reaction is ∼10¹¹ s^-1, which is similar to what has been found for the (Δ S = 1) NO binding reaction in heme proteins. This observation, along with other examples of ultrafast CO binding, provides strong evidence that ligand binding to heme is an adiabatic reaction with a spin-independent prefactor. In order to simultaneously account for the adiabatic nature of the reaction as well as the temperature dependence of both ultrafast CO and methionine geminate rebinding, it is proposed that a spin triplet state intersects and strongly couples to the reactant ( S = 2) and product ( S = 0) state surfaces in the transition state region along the reaction coordinate. It is also suggested that the nature of the intersecting triplet state and the reaction path may depend upon the proximity of the photolyzed ligand relative to the iron atom. At temperatures below ∼60 K, the kinetic data suggest that there is either an unexpected retardation of the heme photoproduct relaxation or that heavy atom quantum mechanical tunneling becomes significant.

Cold-Adaptation Signatures in the Ligand Rebinding Kinetics to the Truncated Hemoglobin of the Antarctic Bacterium Pseudoalteromonas haloplanktis TAC125

Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologues, resulting in enhanced reaction rates at low temperatures. In this context, protein-bound water molecules were suggested to play a major role, and their weaker interactions at protein active sites have been associated with cold adaptation. In this work, we tested this hypothesis on truncated hemoglobins (a family of microbial heme-proteins of yet-unclear function) applying molecular dynamics simulations and ligand-rebinding kinetics on a protein from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 in comparison with its thermophilic Thermobifida fusca homologue. The CO rebinding kinetics of the former highlight several geminate phases, with an unusually long-lived geminate intermediate. An articulated tunnel with at least two distinct docking sites was identified by analysis of molecular dynamics simulations and was suggested to be at the origin of the unusual geminate rebinding phase. Water molecules are present in the distal pocket, but their stabilization by TrpG8, TyrB10, and HisCD1 is much weaker than in thermophilic Thermobifida fusca truncated hemoglobin, resulting in a faster geminate rebinding. Our results support the hypothesis that weaker water-molecule interactions at the reaction site are associated with cold adaptation.

Electrostatic Tuning of the Ligand Binding Mechanism by Glu27 in Nitrophorin 7

Nitrophorins (NP) 1-7 are NO-carrying heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus. The isoform NP7 displays peculiar properties, such as an abnormally high isoelectric point, the ability to bind negatively charged membranes, and a strong pH sensitivity of NO affinity. A unique trait of NP7 is the presence of Glu in position 27, which is occupied by Val in other NPs. Glu27 appears to be important for tuning the heme properties, but its influence on the pH-dependent NO release mechanism, which is assisted by a conformational change in the AB loop, remains unexplored. Here, in order to gain insight into the functional role of Glu27, we examine the effect of Glu27 → Val and Glu27 → Gln mutations on the ligand binding kinetics using CO as a model. The results reveal that annihilation of the negative charge of Glu27 upon mutation reduces the pH sensitivity of the ligand binding rate, a process that in turn depends on the ionization of Asp32. We propose that Glu27 exerts a through-space electrostatic action on Asp32, which shifts the pKa of the latter amino acid towards more acidic values thus reducing the pH sensitivity of the transition between open and closed states.

Ultrafast CO Kinetics in Heme Proteins: Adiabatic Ligand Binding and Heavy Atom Tunneling

The ultrafast kinetics of CO rebinding to carbon monoxide oxidation activator protein (ChCooA) are measured over a wide temperature range and compared with the kinetics of CO binding in other heme systems such as myoglobin (Mb) and hemoglobin (Hb). The Arrhenius prefactor for CO binding to ChCooA and protoheme (∼10¹¹ s^-1) is similar to what is found for spin-allowed NO binding to heme proteins and is several orders of magnitude larger than the prefactor of Mb and Hb (∼10⁹ s^-1). This indicates that the CO binding reaction is adiabatic, in contrast to the commonly held view that it is nonadiabatic due to spin-forbidden (ΔS = 2) selection rules. Under the adiabatic condition, entropic factors, rather than spin-selection rules, are the source of the reduced Arrhenius prefactors associated with CO binding in Mb and Hb. The kinetic response of ChCooA-CO is nonexponential at all temperatures, including 298 K, and is described quantitatively using a distribution of enthalpic rebinding barriers associated with heterogeneity in the heme doming conformation. Above the solvent glass transition (T_g ∼ 180 K), the rebinding progress slows as temperature increases, and this is ascribed to an evolution of the distribution toward increased heme doming and larger enthalpic barriers. Between T_g and ∼60 K, the nonexponential rebinding slows down as the temperature is lowered and the survival fraction follows the predictions expected for a quenched barrier distribution. Below ∼60 K the rebinding kinetics do not follow these predictions unless quantum mechanical tunneling along the heme doming coordinate is also included as an active channel for CO binding.

Structure and Activation of Soluble Guanylyl Cyclase, the Nitric Oxide Sensor

SIGNIFICANCE

Soluble guanylyl/guanylate cyclase (sGC) is the primary receptor for nitric oxide (NO) and is central to the physiology of blood pressure regulation, wound healing, memory formation, and other key physiological activities. sGC is increasingly implicated in disease and is targeted by novel therapeutic compounds. The protein displays a rich evolutionary history and a fascinating signal transduction mechanism, with NO binding to an N-terminal heme-containing domain, which activates the C-terminal cyclase domains. Recent Advances: Crystal structures of individual sGC domains or their bacterial homologues coupled with small-angle x-ray scattering, electron microscopy, chemical cross-linking, and Förster resonance energy transfer measurements are yielding insight into the overall structure for sGC, which is elongated and likely quite dynamic. Transient kinetic measurements reveal a role for individual domains in lowering NO affinity for heme. New sGC stimulatory drugs are now in the clinic and appear to function through binding near or directly to the sGC heme domain, relieving inhibitory contacts with other domains. New sGC-activating drugs show promise for recovering oxidized sGC in diseases with high inflammation by replacing lost heme.

CRITICAL ISSUES

Despite the many recent advances, sGC regulation, NO activation, and mechanisms of drug binding remain unclear. Here, we describe the molecular evolution of sGC, new molecular models, and the linked equilibria between sGC NO binding, drug binding, and catalytic activity.

FUTURE DIRECTIONS

Recent results and ongoing studies lay the foundation for a complete understanding of structure and mechanism, and they open the door for new drug discovery targeting sGC. Antioxid. Redox Signal. 26, 107-121.

Structural changes and picosecond to second dynamics of cytochrome c in interaction with nitric oxide in ferrous and ferric redox states

After dissociation NO rebinds to Cyt c in 10 ps whereas Met80 rebinds in 5 μs after NO release from Cyt c. A complete view of heme – NO dynamics within 12 orders of magnitude of time in Cyt c is presented.

Kinetic Control of O2 Reactivity in H-NOX Domains

Transient absorption, resonance Raman, and vibrational coherence spectroscopies are used to investigate the mechanisms of NO and O2 binding to WT Tt H-NOX and its P115A mutant. Vibrational coherence spectra of the oxy complexes provide clear evidence for the enhancement of an iron-histidine mode near 217 cm(-1) following photoexcitation, which indicates that O2 can be dissociated in these proteins. However, the quantum yield of O2 photolysis is low, particularly in the wild type (≲3%). Geminate recombination of O2 and NO in both of these proteins is very fast (∼1.4 × 10(11) s(-1)) and highly efficient. We show that the distal heme pocket of the H-NOX system forms an efficient trap that limits the O2 off-rate and determines the overall affinity. The distal pocket hydrogen bond, which appears to be stronger in the P115A mutant, may help retard the O2 ligand from escaping into the solvent following either photoinduced or thermal dissociation. This, along with a strengthening of the Fe-O2 bond that is correlated with the significant heme ruffing and saddling distortions, explains the unusually high O2 affinity of WT Tt H-NOX and the even higher affinity found in the P115A mutant.

Ultrafast Spectroscopy Evidence for Picosecond Ligand Exchange at the Binding Site of a Heme Protein: Heme-Based Sensor YddV

An important question for the functioning of heme proteins is whether different ligands present within the protein moiety can readily exchange with heme-bound ligands. Studying the dynamics of the heme domain of the Escherichia coli sensor protein YddV upon dissociation of NO from the ferric heme by ultrafast spectroscopy, we demonstrate that when the hydrophobic leucine residue in the distal heme pocket is mutated to glycine, in a substantial fraction of the protein water replaces NO as an internal ligand in as fast as ∼4 ps. This process, which is near-barrierless and occurs orders of magnitude faster than the corresponding process in myoglobin, corresponds to a ligand swap of NO with a water molecule present in the heme pocket, as corroborated by molecular dynamics simulations. Our findings provide important new insight into ligand exchange in heme proteins that functionally interact with different external ligands.

Motion of proximal histidine and structural allosteric transition in soluble guanylate cyclase

We investigated the changes of heme coordination in purified soluble guanylate cyclase (sGC) by time-resolved spectroscopy in a time range encompassing 11 orders of magnitude (from 1 ps to 0.2 s). After dissociation, NO either recombines geminately to the 4-coordinate (4c) heme (τG1 = 7.5 ps; 97 ± 1% of the population) or exits the heme pocket (3 ± 1%). The proximal His rebinds to the 4c heme with a 70-ps time constant. Then, NO is distributed in two approximately equal populations (1.5%). One geminately rebinds to the 5c heme (τG2 = 6.5 ns), whereas the other diffuses out to the solution, from where it rebinds bimolecularly (τ = 50 μs with [NO] = 200 μM) forming a 6c heme with a diffusion-limited rate constant of 2 × 10(8) M(-1)⋅s(-1). In both cases, the rebinding of NO induces the cleavage of the Fe-His bond that can be observed as an individual reaction step. Saliently, the time constant of bond cleavage differs depending on whether NO binds geminately or from solution (τ5C1 = 0.66 μs and τ5C2 = 10 ms, respectively). Because the same event occurs with rates separated by four orders of magnitude, this measurement implies that sGC is in different structural states in both cases, having different strain exerted on the Fe-His bond. We show here that this structural allosteric transition takes place in the range 1-50 μs. In this context, the detection of NO binding to the proximal side of sGC heme is discussed.

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

NMR studies of the dynamics of high-spin nitrophorins: comparative studies of NP4 and NP2 at close to physiological pH

The β-barrel nitrophorin (NP) heme proteins are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands. NO is bound to iron of the NPs and is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1–NP4, which have sequence similarities in two pairs, NP1 and NP4 (90% identical) and NP2 and NP3 (80% identical). The available crystal structures of NP4 have been used to propose that pH-dependent changes in the conformation of two loops between adjacent β-strands at the front opening of the protein, the A–B and G–H loops, determine the rate of NO release. At pH 7.3, NP4 releases NO 17 times faster than NP2 does. In this work, the aqua complexes of NP4 and NP2 have been investigated by nuclear magnetic resonance (NMR) relaxation measurements to probe the pico- to nanosecond and micro- to millisecond time scale motions at two pH values, 6.5 and 7.3. It is found that NP4-OH₂ is fairly rigid and only residues in the loop regions show dynamics at pH 6.5; at pH 7.3, much more dynamics of the loops and most of the β-strands are observed while the α-helices remain fairly rigid. In comparison, NP2-OH₂ shows much less dynamics, albeit somewhat more than that of the previously reported NP2-NO complex [Muthu, D., Berry, R. E., Zhang, H., and Walker, F. A. (2013) Biochemistry 52, 7910–7925]. The reasons for this major difference between NP4 and NP2 are discussed.

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.

Underlying thermodynamics of pH-dependent allostery

Understanding the effects of coupling protein protonation and conformational states is critical to the development of drugs targeting pH sensors and to the rational engineering of pH switches. In this work, we address this issue by performing a comprehensive study of the pH-regulated switch from the closed to the open conformation in nitrophorin 4 (NP4) that determines its pH-dependent activity. Our calculations show that D30 is the only amino acid that has two significantly different pKas in the open and closed conformations, confirming its critical role in regulating pH-dependent behavior. In addition, we describe the free-energy landscape of the conformational change as a function of pH, obtaining accurate estimations of free-energy barriers and equilibrium constants using different methods. The underlying thermodynamic model of the switch workings suggests the possibility of tuning the observed pKa only through the conformational equilibria, keeping the same conformation-specific pKas, as evidenced by the proposed K125L mutant. Moreover, coupling between the protonation and conformational equilibria results in efficient regulation and pH-sensing around physiological pH values only for some combinations of protonation and conformational equilibrium constants, placing constraints on their possible values and leaving a narrow space for protein molecular evolution. The calculations and analysis presented here are of general applicability and provide a guide as to how more complex systems can be studied, offering insight into how pH-regulated allostery works of great value for designing drugs that target pH sensors and for rational engineering of pH switches beyond the common histidine trigger.

Investigations of heme distortion, low-frequency vibrational excitations, and electron transfer in cytochrome c

Cytochrome (cyt) c is an important electron transfer protein. The ruffling deformation of its heme cofactor has been suggested to relate to its electron transfer rate. However, there is no direct experimental evidence demonstrating this correlation. In this work, we studied Pseudomonas aeruginosa cytochrome c551 and its F7A mutant. These two proteins, although similar in their X-ray crystal structure, display a significant difference in their heme out-of-plane deformations, mainly along the ruffling coordinate. Resonance Raman and vibrational coherence measurements also indicate significant differences in ruffling-sensitive modes, particularly the low-frequency γa mode found between ∼50-60 cm(-1). This supports previous assignments of γa as having a large ruffling content. Measurement of the photoreduction kinetics finds an order of magnitude decrease of the photoreduction cross-section in the F7A mutant, which has nearly twice the ruffling deformation as the WT. Additional measurements on cytochrome c demonstrate that heme ruffling is correlated exponentially with the electron transfer rates and suggest that ruffling could play an important role in redox control. A major relaxation of heme ruffling in cytochrome c, upon binding to the mitochondrial membrane, is discussed in this context.

YC-1 binding to the β subunit of soluble guanylyl cyclase overcomes allosteric inhibition by the α subunit

Soluble guanylate cyclase (sGC) is a heterodimeric heme protein and the primary nitric oxide receptor. NO binding stimulates cyclase activity, leading to regulation of cardiovascular physiology and making sGC an attractive target for drug discovery. YC-1 and related compounds stimulate sGC both independently and synergistically with NO and CO binding; however, where the compounds bind and how they work remain unknown. Using linked equilibrium binding measurements, surface plasmon resonance, and domain truncations in Manduca sexta and bovine sGC, we demonstrate that YC-1 binds near or directly to the heme-containing domain of the β subunit. In the absence of CO, YC-1 binds with a Kd of 9-21 μM, depending on the construct. In the presence of CO, these values decrease to 0.6-1.1 μM. Pfizer compound 25 bound ∼10-fold weaker than YC-1 in the absence of CO, whereas compound BAY 41-2272 bound particularly tightly in the presence of CO (Kd = 30-90 nM). Additionally, we found that CO binds much more weakly to heterodimeric sGC proteins (Kd = 50-100 μM) than to the isolated heme domain (Kd = 0.2 μM for Manduca β H-NOX/PAS). YC-1 greatly enhanced binding of CO to heterodimeric sGC, as expected (Kd ∼ 1 μM). These data indicate the α subunit induces a heme pocket conformation with a lower affinity for CO and NO. YC-1 family compounds bind near the heme domain, overcoming the α subunit effect and inducing a heme pocket conformation with high affinity. We propose this high-affinity conformation is required for the full-length protein to achieve high catalytic activity.

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.

NMR studies of the dynamics of nitrophorin 2 bound to nitric oxide

The Rhodnius nitrophorins are β-barrel proteins of the lipocalin fold with a heme protruding from the open end of the barrel. They are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands, where NO is bound to iron. NO is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1-NP4. At pH 7.3, NP4 releases NO 17 times faster than NP2 does, as measured by stopped-flow kinetics. A number of crystal structures of the least abundant protein, NP4, are available. These structures have been used to propose that two loops between adjacent β-strands at the front opening of the protein, the A-B and G-H loops, determine the rate of NO release. To learn how the protein loops contribute to the release of NO for each of the nitrophorins, the dynamics of these proteins are being studied in our laboratory. In this work, the NP2-NO complex has been investigated by nuclear magnetic resonance relaxation measurements to probe the picosecond-to-nanosecond and microsecond-to-millisecond time scale motions at three pH values, 5.0, 6.5, and 7.3. It is found that at pH 5.0 and 6.5, the NP2-NO complex is rigid and only a few residues in the loop regions show dynamics, while at pH 7.3, somewhat more dynamics, particularly of the A-B loop, are observed. Comparison to other lipocalins shows that all are relatively rigid, and that the dynamics of lipocalins in general are much more subtle than those of mainly α-helical proteins.

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.

pH-dependent picosecond structural dynamics in the distal pocket of nitrophorin 4 investigated by 2D IR spectroscopy

Nitrophorin 4 (NP4) belongs to a family of pH-sensitive, nitric oxide (NO) transporter proteins that undergo a large structural change from a closed to an open conformation at high pH to allow for NO delivery. Measuring the pH-dependent structural dynamics in NP4-NO around the ligand binding site is crucial for developing a mechanistic understanding of NO binding and release. In this study, we use coherent two-dimensional infrared (2D IR) spectroscopy to measure picosecond structural dynamics sampled by the nitrosyl stretch in NP4-NO as a function of pH at room temperature. Our results show that both the closed and open conformers of the protein are present at low (pD 5.1) and high (pD 7.9) pH conditions. The closed and open conformers are characterized by two frequencies of the nitrosyl stretching vibration labeled A0 and A1, respectively. Analysis of the 2D IR line shapes reveals that at pD 5.1, the closed conformer experiences structural fluctuations arising from solvation dynamics on a ∼3 ps time scale. At pD 7.9, both the open and closed conformers exhibit fluctuations on a ∼1 ps time scale. At both pD conditions, the closed conformers maintain a static distribution of structures within the experimental time window of 100 ps. This is in contrast to the open conformer, which is able to interconvert among its substates on a ∼100 ps time scale. Our results directly measure the time scales of solvation dynamics in the distal pocket, the flexibility of the open conformation at high pH, and the rigidity of the closed conformers at both pH conditions. We discuss how the pH-dependent equilibrium structural fluctuations of the nitrosyl ligand measured in this study are related to the uptake and delivery of nitric oxide in NP4.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Kinetics and computational studies of ligand migration in nitrophorin 7 and its Δ1-3 mutant

Nitrophorins (NPs) are nitric oxide (NO)-carrying heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus. Though NP7 exhibits a large sequence resemblance with other NPs, two major differential features are the ability to interact with negatively charged cell surfaces and the presence of a specific N-terminus composed of three extra residues (Leu1-Pro2-Gly3). The aim of this study is to examine the influence of the N-terminus on the ligand binding, and the topological features of inner cavities in closed and open states of NP7, which can be associated to the protein structure at low and high pH, respectively. Laser flash photolysis measurements of the CO rebinding kinetics to NP7 and its variant NP7(Δ1-3), which lacks the three extra residues at the N-terminus, exhibit a similar pattern and support the existence of a common kinetic mechanism for ligand migration and binding. This is supported by the existence of a common topology of inner cavities, which consists of two docking sites in the heme pocket and a secondary site at the back of the protein. The ligand exchange between these cavities is facilitated by an additional site, which can be transiently occupied by the ligand in NP7, although it is absent in NP4. These features provide a basis to explain the enhanced internal gas hosting capacity found experimentally in NP7 and the absence of ligand rebinding from secondary sites in NP4. The current data allow us to speculate that the processes of docking to cell surfaces and NO release may be interconnected in NP7, thereby efficiently releasing NO into a target cell. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

pH-Dependent conformational changes in proteins and their effect on experimental pK(a)s: the case of Nitrophorin 4

The acid-base behavior of amino acids is an important subject of study due to their prominent role in enzyme catalysis, substrate binding and protein structure. Due to interactions with the protein environment, their pK_as can be shifted from their solution values and, if a protein has two stable conformations, it is possible for a residue to have different “microscopic”, conformation-dependent pK_a values. In those cases, interpretation of experimental measurements of the pK_a is complicated by the coupling between pH, protonation state and protein conformation. We explored these issues using Nitrophorin 4 (NP4), a protein that releases NO in a pH sensitive manner. At pH 5.5 NP4 is in a closed conformation where NO is tightly bound, while at pH 7.5 Asp30 becomes deprotonated, causing the conformation to change to an open state from which NO can easily escape. Using constant pH molecular dynamics we found two distinct microscopic Asp30 pK_as: 8.5 in the closed structure and 4.3 in the open structure. Using a four-state model, we then related the obtained microscopic values to the experimentally observed “apparent” pK_a, obtaining a value of 6.5, in excellent agreement with experimental data. This value must be interpreted as the pH at which the closed to open population transition takes place. More generally, our results show that it is possible to relate microscopic structure dependent pKa values to experimentally observed ensemble dependent apparent pK_as and that the insight gained in the relatively simple case of NP4 can be useful in several more complex cases involving a pH dependent transition, of great biochemical interest.

The interaction of an amino acid with its protein environment can result in an acid-base behavior that is very different from what would be observed in solution. This environment can be greatly altered when the protein changes conformation. As a result, the amino acid will have two different “microscopic” pK_a values. Nitrophorin 4 is a good case study to explore this behavior, because it undergoes a pH-dependent conformational change that is well characterized experimentally. Using computer simulation tools, we found that the key titratable Aspartic acid 30, has two very different microscopic pK_as: 4.3 and 8.5, which are significantly different to the observed transition pK_a in solution. However, using a simple model, we were able to understand how this causes the conformational change to take place at pH∼6.5, as measured experimentally. The insight gained in this relatively simple case can be useful in other more complex cases where the apparent pK_a is also a result of the interplay of different conformations where some amino acids experience very different environments.

Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca

CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ∼100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ∼2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ∼100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ∼3 to ∼100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb.

Heterogeneous kinetics of the carbon monoxide association and dissociation reaction to nitrophorin 4 and 7 coincide with structural heterogeneity of the gate-loop

NO is an important signaling molecule in human tissue. However, the mechanisms by which this molecule is controlled and directed are currently little understood. Nitrophorins (NPs) comprise a group of ferriheme proteins originating from blood-sucking insects that are tailored to protect and deliver NO via coordination to and release from the heme iron. Therefore, the kinetics of the association and dissociation reactions were studied in this work using the ferroheme-CO complexes of NP4, NP4(D30N), and NP7 as isoelectronic models for the ferriheme-NO complexes. The kinetic measurements performed by nanosecond laser-flash-photolysis and stopped-flow are accompanied by resonance Raman and FT-IR spectroscopy to characterize the carbonyl species. Careful analysis of the CO rebinding kinetics reveals that in NP4 and, to a larger extent, NP7 internal gas binding cavities are located, which temporarily trap photodissociated ligands. Moreover, changes in the free energy barriers throughout the rebinding and release pathway upon increase of the pH are surprisingly small in case of NP4. Also in case of NP4, a heterogeneous kinetic trace is obtained at pH 7.5, which corresponds to the presence of two carbonyl species in the heme cavity that are seen in vibrational spectroscopy and that are due to the change of the distal heme pocket polarity. Quantification of the two species from FT-IR spectra allowed the fitting of the kinetic traces as two processes, corresponding to the previously reported open and closed conformation of the A-B and G-H loops. With the use of the A-B loop mutant NP4(D30N), it was confirmed that the kinetic heterogeneity is controlled by pH through the disruption of the H-bond between the Asp30 side chain and the Leu130 backbone carbonyl. Overall, this first study on the slow phase of the dynamics of diatomic gas molecule interaction with NPs comprises an important experimental contribution for the understanding of the dynamics involved in the binding/release processes of NO/CO in NPs.

Effect of DNA binding on geminate CO recombination kinetics in CO-sensing transcription factor CooA

Carbon monoxide oxidation activator (CooA) proteins are heme-based CO-sensing transcription factors. Here we study the ultrafast dynamics of geminate CO rebinding in two CooA homologues, Rhodospirillum rubrum (RrCooA) and Carboxydothermus hydrogenoformans (ChCooA). The effects of DNA binding and the truncation of the DNA-binding domain on the CO geminate recombination kinetics were specifically investigated. The CO rebinding kinetics in these CooA complexes take place on ultrafast time scales but remain non-exponential over many decades in time. We show that this non-exponential kinetic response is due to a quenched enthalpic barrier distribution resulting from a distribution of heme geometries that is frozen or slowly evolving on the time scale of CO rebinding. We also show that, upon CO binding, the distal pocket of the heme in the CooA proteins relaxes to form a very efficient hydrophobic trap for CO. DNA binding further tightens the narrow distal pocket and slightly weakens the iron-proximal histidine bond. Comparison of the CO rebinding kinetics of RrCooA, truncated RrCooA, and DNA-bound RrCooA proteins reveals that the uncomplexed and inherently flexible DNA-binding domain adds additional structural heterogeneity to the heme doming coordinate. When CooA forms a complex with DNA, the flexibility of the DNA-binding domain decreases, and the distribution of the conformations available in the heme domain becomes restricted. The kinetic studies also offer insights into how the architecture of the heme environment can tune entropic barriers in order to control the geminate recombination of CO in heme proteins, whereas spin selection rules play a minor or non-existent role.

Spectroscopic analysis of protein Fe–NO complexes

The toxic free radical NO (nitric oxide) has diverse biological roles in eukaryotes and bacteria, being involved in signalling, vasodilation, blood clotting and immunity, and as an intermediate in microbial denitrification. The predominant biological mechanism of detecting NO is through the formation of iron nitrosyl complexes, although this is a deleterious process for other iron-containing enzymes. We have previously applied techniques such as UV–visible and EPR spectroscopy to the analysis of protein Fe–NO complex formation in order to study how NO controls the activity of the bacterial transcriptional regulators NorR and NsrR. These studies have analysed NO-dependent biological activity both in vitro and in vivo using diverse biochemical, molecular and spectroscopic methods. Recently, we have applied ultrafast 2D-IR (two-dimensional IR) spectroscopy to the analysis of NO–protein interactions using Mb (myoglobin) and Cc (cytochrome c) as model haem proteins. The ultrafast fluctuations of Cc and Mb show marked differences, indicating altered flexibility of the haem pockets. We have extended this analysis to bacterial catalase enzymes that are known to play a role in the nitrosative stress response by detoxifying peroxynitrite. The first 2D-IR analysis of haem nitrosylation and perspectives for the future are discussed.

Electronic properties of the highly ruffled heme bound to the heme degrading enzyme IsdI

IsdI, a heme-degrading protein from Staphylococcus aureus, binds heme in a manner that distorts the normally planar heme prosthetic group to an extent greater than that observed so far for any other heme-binding protein. To understand better the relationship between this distinct structural characteristic and the functional properties of IsdI, spectroscopic, electrochemical, and crystallographic results are reported that provide evidence that this heme ruffling is essential to the catalytic activity of the protein and eliminates the need for the water cluster in the distal heme pocket that is essential for the activity of classical heme oxygenases. The lack of heme orientational disorder in (1)H-NMR spectra of the protein argues that the catalytic formation of β- and δ-biliverdin in nearly equal yield results from the ability of the protein to attack opposite sides of the heme ring rather than from binding of the heme substrate in two alternative orientations.

Nitrophorins: nitrite disproportionation reaction and other novel functionalities of insect heme-based nitric oxide transport proteins

Nitrophorins (NPs) comprise a unique class of heme proteins used by the blood-sucking insect Rhodnius prolixus to deliver the signaling gas molecule NO into the blood vessel of a host during feeding. Upon NO release, histamine can be scavenged by coordination to the heme iron. Although the protein is of similar size as the mammalian globin monomers and shares the same cofactor and proximal histidine coordination, nitrophorin structure, in contrast, is almost entirely composed of a β-barrel. Comparison of the NO and histamine association constants with the concentrations of both compounds invivo raises concerns about the very simple ligand release model in case of at least some of the NPs. Therefore, novel functionalities of the NPs were sought. As a result, catalysis of the nitrite disproportionation reaction was found, which leads to the formation of NO with nitrite as the sole substrate. This is the first example of a ferriheme protein that can perform this reaction. Furthermore, although NPs stabilize the ferriheme state, a peroxidase reactivity of the cofactor involving the higher oxidation state iron (Compound I/II) was studied with the potential to catalyze the oxidation of histamine and norepinephrine. In contrast to many other heme proteins including the globins, the ferroheme state was found to be extremely sensitive to O(2) , which is a consequence of the much lower reduction potential of the NPs, so that the 1-electron reduction of O(2) to O (•-)(2) becomes a thermodynamically favored process. Altogether, the detailed study of the NPs gives insight into the structure-function relationships required for the targeted delivery of diatomic gas molecules in biology. Moreover, the comparison of the structure-function relationships of the NPs (NO transporters) with those of the globins (O(2) transporters) will help to elucidate the architectural requirement for the respective tasks.

Relationship between protein structural fluctuations and rebinding dynamics in ferric haem nitrosyls

The interaction of nitric oxide (NO) with haem proteins is widespread in biology. In the current paper, we present the first ultrafast 2D-IR (two-dimensional infrared) spectroscopic analysis of haem nitrosylation, which has been combined with time-resolved IR pump–probe studies to investigate the relationship between equilibrium vibrational dynamics of the haem environment and ligand rebinding behaviour following photolysis of NO from the Fe(III)–NO site. Studies of two haem proteins, Mb (myoglobin) and Cc (cytochrome c), which play different physiological roles, reveal marked contrasts in the ultrafast fluctuations of the protein pockets containing the haem, showing that the Mb pocket is somewhat more flexible than that of Cc. This correlates strongly with slower observed photolysis rebinding kinetics of Mb–NO compared with Cc–NO, and indicates a direct link between ultrafast fluctuations and biological functionality. Furthermore, this indicates the validity of linear response theories in relation to protein ligand binding. Finally, 2D-IR shows that Cc–NO displays two distinct structural sub-sites at room temperature that do not exchange on the timescales accessible via the NO vibrational lifetime.