Next Article in Journal
Is the ENaC Dysregulation in CF an Effect of Protein-Lipid Interaction in the Membranes?
Next Article in Special Issue
Occurrence, Evolution and Specificities of Iron-Sulfur Proteins and Maturation Factors in Chloroplasts from Algae
Previous Article in Journal
Different miRNA Profiles in Plasma Derived Small and Large Extracellular Vesicles from Patients with Neurodegenerative Diseases
Previous Article in Special Issue
Gene Expression Responses to Sequential Nutrient Deficiency Stresses in Soybean
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin

1
Dipartimento di Scienze Matematiche, Fisiche e Informatiche, Università di Parma, Parco Area delle Scienze 7/A, 43124 Parma, Italy
2
Center for Nano Science and Technology@PoliMi, Istituto Italiano di Tecnologia, via G. Pascoli 70/3, 20133 Milano, Italy
3
Departamento de Nutrición Vegetal, Estación Experimental de Aula Dei, Consejo Superior de Investigaciones Científicas (CSIC), Apartado 13034, 50080 Zaragoza, Spain
4
Dipartimento di Scienze degli Alimenti e del Farmaco, Università di Parma, Parco Area delle Scienze 27/A, 43124 Parma, Italy
5
IFN-CNR, Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(5), 2740; https://doi.org/10.3390/ijms22052740
Submission received: 31 December 2020 / Revised: 16 February 2021 / Accepted: 2 March 2021 / Published: 8 March 2021
(This article belongs to the Special Issue Iron and Sulfur in Plants 2.0)

Abstract

:
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 108 M−1s−1 range. Binding of the distal histidine to the heme competes with CO rebinding with extremely high rates (kh ~ 105 s−1). Histidine dissociation from the heme occurs with comparable rates, thus resulting in moderate equilibrium binding constants (KH ~ 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.

1. Introduction

Plant hemoglobins were first discovered in the root nodules of legumes. They were termed leghemoglobins and shown to be able to oxygenate reversibly. Subsequently, hemoglobins were also found in the nodules of actinorhizal plants and Parasponia (for a review of early references, see [1]). The only known function of leghemoglobins and other symbiotic hemoglobins of nodules is the transport and delivery of O2 to the nitrogen-fixing bacteria [1]. However, we know now that hemoglobins are widespread in nonsymbiotic tissues of all plants (for reviews, see [2,3,4,5,6,7]). These nonsymbiotic hemoglobins (now called phytoglobins (Glbs)) can be categorized into three classes according to phylogenetic analyses and biochemical properties [6,7,8,9,10]. Class 1 Glbs have an extremely high affinity for O2 and weak hexacoordination of the iron heme with a distal His (Hisd) residue, whereas class 2 Glbs have strong hexacoordination which results in O2 affinities that are ~100-fold lower than those of class 1 Glbs [10]. In both Glb classes, the presence of hydrophobic cavities and tunnels modulates ligand exchange between the active site on the heme distal site and the solvent, providing preferential pathways for the supply of reactants and egression of products [11,12,13,14]. Class 3 Glbs show high sequence similarity with bacterial truncated hemoglobins and relatively low O2 affinity [8,10,15]. The class 3 Glb of Arabidopsis thaliana (AtGlb3) is hexacoordinate when deoxygenated, but this state is transitory and reverts slowly (20 min) to the pentacoordinate form [15]. However, only a few Glbs of classes 2 and 3 have been characterized at the biophysical level.
The genes of the three Glb classes may respond differently to environmental conditions. Class 1 Glb genes are upregulated by hypoxia, osmotic and cold stress, nutrient deprivation, rhizobial and fungal infection and nitric oxide (NO) [6,16,17,18,19,20]. As shown for the class 1 Glbs of monocots and A. thaliana (AtGlb1), a function of class 1 Glbs is to modulate NO concentration, which allows the plant to adapt and survive under limiting O2 conditions [3]. This is thought to be due in part to the NO dioxygenase activity of the proteins, whereby the oxyferrous globin reacts with NO to yield nitrate and ferric globin [21,22]. Much less is known about the functions of class 2 and 3 Glbs. The class 2 Glb of A. thaliana (AtGlb2) is responsive to low temperatures and cytokinin treatment [23,24]. This protein also reacts with NO in vivo [22,25] and is involved in plant morphogenesis [25] and in the regulation of somatic embryogenesis [19]. Class 3 Glbs are unresponsive to several types of stress that affect the other two Glb classes, although they also respond to the supply of certain hormones in the roots [18]. A function for these Glbs in NO metabolism has been proposed in the symbiosis of plants with rhizobia and arbuscular mycorrhizal fungi [26]. In fact, a subclass of class 3 Glbs (Glb3-1) are highly expressed in nodules [18,27].
The recently identified MtGlb1-2 gene of the model legume Medicago truncatula is unique because it generates four alternative splice forms, designated MtGlb1-2.1 to MtGlb1-2.4, encoding proteins with one or two heme domains [28]. Thus, MtGlb1-2.1 and MtGlb1-2.2 have two heme domains, whereas MtGlb1-2.4 has one heme domain and it is unclear whether MtGlb1-2.3 has one or two heme domains. The MtGlb1-2 gene is transcriptionally activated by hypoxia and NO and is preferentially expressed in the meristematic region and vascular bundles of roots and root nodules [28]. In this work, we have used CO as a model ligand to investigate the binding properties of MtGlb1-2.1 and MtGlb1-2.4, as representative proteins with two hemes and one heme, respectively. These Glbs are extremely reactive toward the physiological ligands O2, NO and nitrite and show very high O2 affinities, NO dioxygenase activity (in the presence of O2) and nitrite reductase activity (in the absence of O2). A detailed kinetic characterization of the ligand rebinding processes required extension of temporal resolution to the short picosecond range, with the use of a hybrid femtosecond–nanosecond pump–probe setup. The results show unprecedented properties of these Glbs, with very large and fast geminate rebinding and very high bimolecular CO rebinding rates and Hisd binding and dissociation rates. Our studies with Hisd mutants confirmed the identity of the kinetic processes attributed to the Hisd binding and dissociation.

2. Results and Discussion

2.1. Unusually Fast CO Rebinding Kinetics to WT1 and WT4 after Nanosecond Laser Photolysis

The CO rebinding kinetics to wild-type MtGlb1-2.1 (WT1) and MtGlb1-2.4 (WT4) after nanosecond laser flash photolysis (LFP) are qualitatively very similar and are shown in Figure 1 for temperatures between 10 °C and 40 °C and for solutions equilibrated with 1 or 0.2 atm CO. As previously reported for other hexacoordinate hemoglobins, the overall progress curve is composed of two well-separated processes for both proteins. A fast, CO concentration-independent (unimolecular) rebinding is evident on the short nanosecond time scale, followed by a CO concentration-dependent (bimolecular) rebinding phase on the microsecond time scale. The effect of different CO concentrations at each temperature is easily appreciated by comparing the solid (1 atm CO) and the dotted (0.2 atm CO) curves.
The faster phase is commonly interpreted as unimolecular ligand rebinding from within the protein matrix and is called geminate rebinding. A minor amplitude decay is detectable in the long nanoseconds to short microseconds, whose CO concentration dependence is very weak. In contrast, the slower process occurring in the microsecond time scale shows a CO concentration dependence, indicating a bimolecular nature. Decreasing CO concentration results in slower and broader microsecond rebinding kinetics, best described by double exponential relaxations especially at higher temperatures (Figure S1, Table S1). Considering that the two proteins in their ferrous state show bis-histidyl hexacoordination [28], the heterogenous kinetics in the microsecond time scale is thus expected to reflect competitive distal histidine (Hisd) binding to, and dissociation from, the unliganded Fe2+ after photolysis, as previously reported for other plant and vertebrate hemoglobins [29,30,31]. The slowest portion of the kinetics is usually attributed to Hisd dissociation, and the amplitude of this process increases with temperature and upon decreasing CO concentration.
The qualitative effect of temperature on CO rebinding kinetics is easily appreciated by comparing rebinding curves at the same CO concentration in Figure 1. The first striking feature is that the amplitude of geminate rebinding increases dramatically upon decreasing the temperature, at the expense of the amplitude of the microsecond rebinding. Since the binding step to the heme is normally characterized by a very low enthalpic barrier, this finding is indicative of a large activation energy to exit from the protein matrix and reach the solvent. The microsecond process shows an increase in apparent rate at higher temperatures, directly revealing a thermal activation of the associated rate constant. Unlike other Glbs, the time required for Hisd dissociation and ligand rebinding completion is extraordinarily short (~1 ms). For comparison, other hexacoordinate Glbs from plants are characterized by Hisd dissociation rates spanning from 10 s−1 (class 2 Glbs) to 102 s−1 (class 1 Glbs) [10]. The Hisd dissociation rates from vertebrate hexacoordinate hemoglobins are much lower (~0.1–1 s−1) [5].
Importantly, the smaller than expected absorbance change at the end of the laser pulse and the steep drop of the signals in the short (nanoseconds) timescales (Figure 1) suggest that a relevant portion of the rebinding kinetics may fall within the 7 ns laser pulse width. Thus, in order to be able to perform a quantitative kinetic analysis, experiments with improved time resolution are necessary to complete the time course of the rebinding processes, accessing microscopic steps that occur in the subnanosecond time range.

2.2. Experiments with Subnanosecond Resolution Reveal a Large Geminate Rebinding

We thus applied a hybrid femtosecond–nanosecond pump–probe setup to collect the full time course of CO rebinding and gain access to the subnanosecond time frame. Besides allowing subnanosecond resolution, the setup allows collecting the kinetics up to 1 ms, seamlessly merging the ultrafast timescale with the time range normally covered in nanosecond LFP experiments. This feature of the apparatus represents a potentially major improvement over a more traditional approach, leaving a gap between the longest time accessible in femtosecond pump–probe (~1–2 ns) and the shortest time measured with nanosecond LFP (~20–30 ns) [32,33,34]. The proposed setup represents an alternative experimental approach that is simpler than previously reported instruments based on two independent and electronically synchronized femtosecond lasers [35,36].
Figure 2 compares the CO rebinding kinetics for WT1 at room temperature collected with the hybrid femtosecond–nanosecond pump–probe setup (gray dots) and with the nanosecond LFP setup (gray solid line). From simple visual inspection, it is evident that there is perfect agreement between the two curves in the time range where they overlap. Excellent agreement was also observed for MtGlb1-2.1 single and double mutants at the Hisd (see Section 2.4). Similar results were obtained for WT4 and its Hisd mutant.
The overall progress curve shows that a relevant fraction of CO rebinding occurs in the picoseconds, with a geminate phase accounting for 60% of the kinetics in WT1 (Figure 2). Similarly, the amplitude of geminate phase to WT4 is found to be ~70% of the overall rebinding. Thus, <50% of the photodissociated ligands in WT1 and WT4 escape geminate rebinding and migrate to the solvent, giving rise to the bimolecular phase.
Geminate rebinding to WT1 and WT4 is characterized by a biexponential kinetics (the two faster kinetic phases for each curve shown in Figure S2 and Table S2). The source of this heterogeneity may be related to either the presence of two heme domains with different binding rate constants (for WT1) or the presence of migration to nearby docking sites, from which the ligand is then rebound. Similar behavior was reported for geminate rebinding to truncated hemoglobins for which the identification of ligand egression and migration pathways from molecular modeling gave support to kinetic modeling of observed rebinding traces in some truncated globins from bacteria [32,33,37,38] and for nitrophorin 7 [39]. We favor the ligand migration interpretation for both WT1 and WT4 because the second-order rebinding is described by a single rate constant (see Supplementary Materials).

2.3. Description of the Rebinding Kinetics with a Microscopic Model

In the light of the above considerations, we modeled the observed kinetics using a reaction scheme (Scheme 1), where the unimolecular geminate phase comprises migration to a temporary docking site, accessible from the distal pocket, and the slower phase is the result of competitive binding to the heme of CO and the Hisd. The numerical method has been previously discussed [32,37,40].
The agreement between the experimental data and the fitting curves built with the numerical model is satisfactory, as can be appreciated in Figure 3a for WT1 and WT4. For simplicity, in this panel, we plotted the modeled rebinding kinetics on top of the measured kinetics. In Figure 3b–d we also report the time course of the reaction intermediates. The values of the microscopic rates obtained in the global analysis of data at 0.2 and 1 atm CO at room temperature are reported in Table 1.
The microscopic rate constant values retrieved from the fitting of pump–probe data using the minimal kinetic scheme proposed in Scheme 1 can be used to calculate the CO binding rate constant to the pentacoordinate species (kon = k−1kin/(kout + kin)), the observed CO binding rate (kobs = kon[CO]k−h/(k−h + kh + kon[CO])) and the Hisd equilibrium binding constant (KH = kh/k−h).
For both WT1 and WT4, the large amplitude kinetics in the picosecond–nanosecond scale (Figure 3b,c) is described through a significant direct CO rebinding to ferrous heme from a primary docking site in the distal cavity (k−1 ~108 s−1 for WT1 and ~109 s−1 for WT4), modulated by a transient population (kc and k−c ~108 s−1 for both proteins) of a secondary docking site, accessible from the distal pocket. The large amplitude of geminate rebinding in WT1 and WT4 suggests that the connection between the distal heme pocket and the solvent is inhibited in favor of the direct rebinding or internal migration. The innermost binding rate (k−1) shows values for WT1 and WT4 (Table 1) that are much higher than the ones reported for other hexacoordinate heme proteins, such as human neuroglobin (1.5 × 107 s−1 [41,42]) and the Glbs from A. thaliana (0.5 × 107 s−1 for AtGlb1 and 2.2 × 107 s−1 for AtGlb2 [12]) or from rice (Osat1-1; 7.6 × 107 s−1 [31]), and similar only to the ones typical of truncated HbO (3 × 108 s−1 for Thermobifida fusca [32], 8.37 × 108 s−1 for Pseudoalteromonas haloplanktis TAC125 [37], ~8 × 108 s−1 for truncated HbO from Mycobacterium tuberculosis [43]). Direct rebinding through the very high rate (k−1) is only partially reduced by the competitive migration to the inner docking site (kc = 2.3 × 108 s−1 for WT1 and 6.8 × 108 s−1 for WT4) and exit to the solvent (kout = 1.4 × 108 s−1 for WT1 and 4.9 × 108 s−1 for WT4). It is worth noting that the secondary docking site in WT4 provides a more stable kinetic trap to photodissociated ligands than in WT1 and results in a distinct kinetic step on the 10 ns time scale (Figure 3a, blue trace). The low fraction of photodissociated ligands reaching the solvent is mostly determined by the very high rates for direct rebinding and internal migration rather than by a low exit rate. Indeed, comparison with typical kout values for class 1 nonsymbiotic Glbs like AtGlb1 (9 × 107 s−1 [44]) and rice (OsGlb1-1, 7.33 × 107 s−1 [31]) shows that the exit rate for WT1 is 2-fold higher, whereas kout for WT4 is almost one order of magnitude larger. However, in those cases, the internal rates were much lower.
The structural and dynamical determinants for this behavior may be identified when structural information becomes available. Unfortunately, a topological description of cavities in the three-dimensional structure of WT1 and WT4 is not yet available, and it is not possible to identify the structural nature of the kinetic trap. However, the presence of cavities and migration pathways for ligands in hexacoordinate hemoglobins is well established in other plant Glbs, such as OsGlb1-1 [31] and AtGlb1 and AtGlb2 [12], as well as in vertebrate hemoglobins like human neuroglobin [41,45] and cytoglobin [46]. Stable or transient cavities and tunnels are usually lined by hydrophobic amino acids. Thus, they can accommodate diatomic gaseous ligands and function as kinetic traps [47]. A role of the kinetic traps in support of multisubstrate reactions, such as NO dioxygenase activity, has been suggested for myoglobin and hexacoordinate hemoglobins [48].

2.4. Rate Constants for Hisd Binding and Dissociation Are Extraordinarily High

A surprising result of the numerical analysis is the very high rate constants for association (kh = 7–8 × 105 s−1) and dissociation (k−h = 3–4 × 105 s−1) of the Hisd (Table 1), which are about three orders of magnitude larger than the reported values for any class 1 Glb (average kh = 130 s−1 and k−h = 75 s−1) or class 2 Glb (average values kh = 1500 s−1 and k−h = 25 s−1) [10], with unusually low energetic barriers (vide infra, Table 2). Moreover, Hisd dissociation rates from vertebrate hexacoordinate globins are much lower (~0.1–1 s−1) [5], and the highest kh value known so far is 1.4 × 104 s−1 and comes from a mollusk hexacoordinate hemoglobin [49]. Nevertheless, the resulting equilibrium binding constant of Hisd (KH) of ~2 (Table 1) to ferrous heme is in line with class 1 Glbs (average 1.7 ± 0.8) [5]. The microscopic rate constants enable us to calculate the binding rate constants to the pentacoordinate species as follows in Equation (1)
k =   k 1 k i n k i n +   k o u t
yielding kon = 1 × 108 M−1s−1 for WT1 and kon = 2.4 × 108 M−1s−1 for WT4. These values are at least one order of magnitude larger than for other class 1 nonsymbiotic plant Glbs (average value: (8.4 ± 14) × 106 M−1 s−1) [10]. A practical consequence of the very high dissociation rates is that the observed CO binding rate, Equation (2)
k = k o n C O k h k o n C O +   k h +   k h
falls in the 104 s−1 range (2.6 × 104 s−1 for WT1 and 6 × 104 s−1 for WT4) at 1 atm CO (Table 1).
In order to prove that the slowest kinetic step in the progress curve is associated with the formation of a bis-histidyl species, we obtained mutants in which the Hisd is replaced by Leu in one of the two heme domains of MtGlb1-2.1 (mutants 74 and 238) or in the two domains (mutant 74/238). Likewise, we replaced the Hisd with Leu in the single heme domain of MtGlb1-2.4 (mutant 109). Then, we examined the CO rebinding of the mutated proteins. The resulting progress curves of the 74/238 (Figure 3b) and 109 (Figure 3c) mutants show larger geminate rebinding. This indicates that the Hisd mutation interferes with ligand migration to internal cavities and to the solvent. On longer time scales, the bimolecular phase becomes monoexponential, confirming the identification of the Hisd dissociation with the slowest kinetic step in the WT proteins.
Comparisons of the best fit curves obtained for CO rebinding kinetics between WT1 and 74/238 (Figure 3b) or between WT4 and 109 (Figure 3c) show that the experimental data for the mutants can be described using the same kinetic scheme used for the WT proteins but without the kinetic step leading to the formation of the bis-histidyl species (GlbH, red traces). The amplitude of the geminate phase becomes larger in the mutants with respect to the WT proteins, particularly for the 109 mutant. This effect is the result of changes in k−1, kout and kin. The innermost binding rate k−1 increases ~2-fold in the 109 mutant and kout decreases by a similar factor, whereas kin increases ~3.5-fold and ~6.5-fold in the 74/238 and 109 mutants, respectively.
These results suggest a possible role of Hisd in modulating ligand exchange with the exterior and in modifying the plasticity of the protein structure. The kon rate remains practically unchanged in the 74/238 mutant but increases by a factor of 7 in the 109 mutant.

2.5. Dissecting the Contribution of the Hisd in Each Domain of MtGlb1-2.1

Figure 3d shows the effect of eliminating only one Hisd in MtGlb1-2.1. It is evident that the kinetics for single mutants are intermediate between those of WT1 and 74/238. The geminate amplitude becomes slightly larger in the 74 mutant and even more so in the 238 mutant. It is evident that the effect of the Hisd on the heme reactivity of each domain is similar but not identical. This observation is in line with the different nitrite reductase activity recently demonstrated by some of us for the two mutants [28].
The kinetic analyses, according to Scheme 1, are shown in Figure 4, and the fitting parameters are listed in Table 1. With respect to WT1, there are no substantial differences in the rate constants that describe ligand migration to the kinetic trap and to the solvent, whereas CO rebinding from the outside appears to be more influenced by mutations (~2.5-fold change in rate constant). According to these parameters, kon increases from 1 × 108 M−1 s−1 for WT1 to 1.5 × 108 M−1 s−1 for the 74 mutant and to 2 × 108 M−1 s−1 for the 238 mutant. The effect of these single mutations on kon confirms that accessibility of the distal cavity is influenced by Hisd, as already suggested for the 74/238 mutant. On the basis of these results, we propose that when the Hisd is replaced by a Leu residue, a more accessible exchange channel between the distal pocket and the solvent is opened, through which CO molecules can migrate in and out. As expected, due to the fact that MtGlb1-2.1 bears two hemes, the amplitude of the hexacoordinate species is partly reduced, but not removed, by single Hisd mutations (74 and 238). Moreover, the rate constants for ligation and deligation in single mutants are smaller than the corresponding WT1 values: change in kh ~2.5-fold for both mutants; change in k−h ~2.9-fold for the 74 mutant and 2.3-fold for the 238 mutant. This suggests that the single mutations affect the dynamics of the other domain and result in increased rate constants for Hisd binding and dissociation, with little consequences on the equilibrium binding constants KH = kh/k−h.

2.6. Thermal Activation of Kinetic Processes

Finally, we examined the temperature dependence of CO rebinding kinetics in the range of 10–40 °C (Figure 1). It is evident that, upon increasing temperature, the amplitude of the geminate phase decreases and the bimolecular phase speeds up. The former consideration indicates that protein fluctuations significantly assist CO migration inside the protein and the exchange between inside and outside, whereas the latter evidences the typical thermal activation behavior for diffusion-mediated processes. Kinetics acquired at the same temperature at 1 and 0.1 CO atm were simultaneously analyzed in order to improve the reliability of the best fit parameters. Each global analysis permits the estimation of the microscopic rate constants of the observed processes at a defined temperature. Based on previous studies on related Glbs and other hemoglobins [32,33,37], we assume in the analysis that the innermost rate constants k−1, kc and k−c are independent of temperature. From the temperature dependence of the microscopic rate constants, the activation energy of each kinetic step can be evaluated from the resulting linear Arrhenius plots (Figure S3). The values are reported in Table 2 for WT proteins and their mutants. Because only the temperature dependence of laser LFP traces can be used, only the thermal activation of the bimolecular phase and the hexacoordination process were studied. The exit rate to the solvent (kout) has an energy barrier similar for WT1 and WT4 (>10 kcal/mol), which is significantly higher than the values reported for AtGlb1 (2.7 kcal/mol [11]) and OsGlb1-1 (4.3 kcal/mol [31]). By contrast, the entry rate from the solvent (kin) shows an activation energy in line with the values reported for AtGlb1 (14.7 kcal/mol [11]) and OsGlb1-1 (12.4 kcal/mol [31]). The energy barrier associated with kout is only half as large for the mutants. The decrease of the thermal activation energy of kin is smaller in WT1 single mutants and is ~2-fold for the 109 and 74/238 mutants, where no Hisd residue is present. These results are in agreement with the ones retrieved from the analysis of pump–probe experiments on the mutants that suggest a more accessible pathway between the distal pocket and the solvent when the Hisd is substituted by a Leu residue. Interestingly, thermal activation associated with ligation (kh, 1.1 kcal/mol for WT1 and 1.4 kcal/mol for WT4) and deligation (k−h, 5.1 kcal/mol for WT1 and 2.6 kcal/mol for WT4) is heavily reduced with respect to AtGlb1 and OsGlb1-1 and other hexacoordinate globins in general. In fact, the energy barriers are 19 kcal/mol (AtGlb1) and 15.5 kcal/mol (OsGlb1-1) for kh, whereas values of 18 kcal/mol (AtGlb1) and 14 kcal/mol (OsGlb1-1) were retrieved for k−h. This low thermal activation is in line with the kinetic results that highlight very high reaction rates for Hisd binding and dissociation (kh and k−h).

3. Conclusions

The two splice forms, MtGlb1-2.1 and MtGlb1-2.4, of the recently identified MtGlb1-2 gene of the model legume M. truncatula encode proteins with two heme domains and one heme domain, respectively. Both are extremely reactive toward the physiological ligands O2, NO and nitrite and show very high O2 affinities, NO dioxygenase activity (in the presence of O2) and nitrite reductase activity (in the absence of O2). Ligand rebinding kinetics of the model diatomic ligand CO reveal unusual properties that may be fundamental to sustain the high reactivity toward multisubstrate reactions. Both proteins show unprecedentedly high rate constants for Hisd binding and dissociation, which results in moderate equilibrium binding constants, typical of class 1 Glbs. A very high rebinding rate to the heme Fe is responsible for an unusually large and fast geminate recombination, which mostly occurs in the subnanosecond time scale. The high binding rate to the heme results in very high bimolecular binding rates. As a result of the above properties, the observed binding rate of the diatomic ligands is extremely large and may be at the basis of the extreme reactivity of the proteins in the multisubstrate reactions they are involved in.

4. Materials and Methods

4.1. Protein Expression and Purification

The coding regions of MtGlb1-2.1 and MtGlb1-2.4 [28] were synthesized with codon optimization for expression in Escherichia coli. Gene synthesis and site-directed mutagenesis were performed by GenScript (Hong Kong). The MtGlb1-2.1 and MtGlb1-2.4 proteins and their mutant derivatives were cloned into the pET30a (+) expression vector (Novagene, Madison, WI, USA) with an N-terminal His-tag. The proteins were expressed in E. coli C41(DE3) cells (Lucigen, Middleton, WI, USA) and purified on nickel-affinity columns (HiTrap Chelating HP; GE Healthcare, Uppsala, Sweden) as previously described [28].

4.2. Nanosecond Laser Flash Photolysis of CO Rebinding Kinetics

Measurements were carried out using LFP instrumentation very similar to that described previously [12]. Briefly, the second harmonic of a nanosecond Nd:YAG laser (Surelite I-10, Continuum, San Jose, California, USA) was used to photodissociate CO molecules, and a continuous-wave probe beam from a 75 W Xe lamp was exploited to monitor transient absorbance changes at 436 nm. The time evolution of transmitted intensity was followed on a digital scope and then converted into absorbance change, ∆A (t). The temperature of the sample was precisely controlled with a Peltier element (Flash 100, Quantum Northwest, Inc. Liberty Lake, USA), allowing temperature stability better than 0.1 °C.

4.3. Ultrafast Spectroscopy of CO Rebinding Kinetics

Transient absorption experiments were performed in transmission geometry. An amplified Ti:sapphire laser (Integra-C, Quantronix, East Setauket, NY, USA) generated 100-fs pulses centered at 800 nm, at a repetition rate of 1 kHz. A broadband UV-Vis probe was generated by focusing the pulses into a thin CaF2 plate mounted on a moving stage to avoid laser-induced damage. For time delays <1 ns, the sample was excited with pump pulses centered at 530 nm (~100 fs FWHM), generated in a home-built noncollinear parametric amplifier [50] with pump–probe delay determined by an optical delay line. For delays >1 ns, pump pulses centered at 532 nm (700 ps FWHM) were obtained from a Q-switched Nd:YVO4 laser (Picolo, InnoLas, Krailling, Germany) which was electronically triggered and synchronized to the Ti:sapphire laser via an electronic delay [51]. After interaction with the sample, a spectrometer (Shamrock 303i Andor, Belfast, United Kingdom) dispersed the probe light onto a fast CCD array, enabling broadband shot-to-shot detection of differential absorption (∆A) spectra.

4.4. Data Analyses

The complete time evolution of CO rebinding kinetics was analyzed using a minimal kinetic model as previously described for other plant Glbs [11]. Numerically solving the differential equations associated with the scheme with the function ODE15s within MATLAB (MathWorks, Natick, MA, USA) and optimizing the parameters with a MATLAB version of the optimization package Minuit (CERN) allowed us to evaluate the microscopic rate constants of the observed processes.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/22/5/2740/s1. Figure S1: CO rebinding kinetics using nanosecond LFP in function of CO concentration at 40 °C. Figure S2: Analysis of the complete course of CO rebinding kinetics to WT1 and WT4. Figure S3: Arrhenius plots for the rate constants obtained from the analysis of nanosecond LFP data of WT1 and WT4. Table S1: Lifetimes from the fit of nanosecond LFP data at 40 °C at 1 atm and 0.1 atm CO for WT1 and WT4 and their mutants, using a sum of exponential decay functions. Table S2: Lifetimes from the fit of hybrid pump–probe data at room temperature at 1 CO atm for WT1 and WT4, using a sum of exponential decay functions.

Author Contributions

Conceptualization, S.A., C.V. and M.B.; methodology, S.A., S.B., C.V. and G.C.; validation, S.A., C.V. and M.B.; formal analysis, S.A. and A.J.B.; investigation, S.A., S.B., A.J.B., I.V. and C.P.-R.; resources, C.V., G.C. and M.B.; data curation, S.A. and C.V.; writing—original draft preparation, S.A., C.V., G.C. and M.B.; writing—review and editing, S.A., C.V. and M.B.; visualization, S.A., C.V. and M.B.; supervision, S.A., C.V. and M.B.; project administration, M.B.; funding acquisition, M.B. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Agencia Estatal de Investigación (AEI) grant AGL2017-85775-R, co-funded by the European Regional Development Fund, and by Government of Aragón (group A09_17R). The APC was funded by AEI and, in part, by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). G.C. acknowledges support from the PRIN 2017 Project 201795SBA3—HARVEST.

Acknowledgments

S.A. and C.V. thank Alessandro Baistrocchi and Mattia Castellucchio for their contributions in performing LFP experiments during their Bachelor thesis work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

GlbPhytoglobin (or nonsymbiotic hemoglobin)
HisdDistal histidine
LFPLaser flash photolysis
NONitric oxide
O2Molecular oxygen
COCarbon monoxide
WTWild type

References

  1. Appleby, C.A. The origin and functions of hemoglobin in plants. Sci. Prog. 1992, 76, 365–398. [Google Scholar]
  2. Garrocho-Villegas, V.; Gopalasubramaniam, S.K.; Arredondo-Peter, R. Plant hemoglobins: What we know six decades after their discovery. Gene 2007, 398, 78–85. [Google Scholar] [CrossRef]
  3. Dordas, C. Nonsymbiotic hemoglobins and stress tolerance in plants. Plant Sci. 2009, 176, 433–440. [Google Scholar] [CrossRef] [PubMed]
  4. Jokipii-Lukkari, S.; Frey, A.D.; Kallio, P.T.; Häggman, H. Intrinsic non-symbiotic and truncated haemoglobins and heterologous Vitreoscilla haemoglobin expression in plants. J. Exp. Bot. 2009, 60, 409–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kakar, S.; Hoffman, F.G.; Storz, J.F.; Fabian, M.; Hargrove, M.S. Structure and reactivity of hexacoordinate hemoglobins. Biophys. Chem. 2010, 152, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hill, R.D. Non-symbiotic haemoglobins—What’s happening beyond nitric oxide scavenging? AoB Plants 2012, 2012, pls004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Becana, M.; Yruela, I.; Sarath, G.; Catalán, P.; Hargrove, M.S. Plant hemoglobins: A journey from unicellular green algae to vascular plants. New Phytol. 2020, 227, 1618–1635. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Wittenberg, J.B.; Bolognesi, M.; Wittenberg, B.A.; Guertin, M. Truncated hemoglobins: A new family of hemoglobins widely distributed in bacteria, unicellular eukaryotes, and plants. J. Biol. Chem. 2002, 277, 871–874. [Google Scholar] [CrossRef] [Green Version]
  9. Vinogradov, S.N.; Hoogewijs, D.; Bailly, X.; Arredondo-Peter, R.; Guertin, M.; Gough, J.; Dewilde, S.; Moens, L.; Vanfleteren, J.R. Three globin lineages belonging to two structural classes in genomes from the three kingdoms of life. Proc. Natl. Acad. Sci. USA 2005, 102, 11385–11389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Smagghe, B.J.; Hoy, J.A.; Percifield, R.; Kundu, S.; Hargrove, M.S.; Sarath, G.; Hilbert, J.L.; Watts, R.A.; Dennis, E.S.; Peacock, W.J.; et al. Correlations between oxygen affinity and sequence classifications of plant hemoglobins. Biopolymers 2009, 91, 1083–1096. [Google Scholar] [CrossRef]
  11. Bruno, S.; Faggiano, S.; Spyrakis, F.; Mozzarelli, A.; Cacciatori, E.; Dominici, P.; Grandi, E.; Abbruzzetti, S.; Viappiani, C. Different roles of protein dynamics and ligand migration in non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. Gene 2007, 398, 224–233. [Google Scholar] [CrossRef]
  12. Bruno, S.; Faggiano, S.; Spyrakis, F.; Mozzarelli, A.; Abbruzzetti, S.; Grandi, E.; Viappiani, C.; Feis, A.; Mackowiak, S.; Smulevich, G.; et al. The reactivity with CO of AHb1 and AHb2 from Arabidopsis thaliana is controlled by the distal His E7 and internal hydrophobic cavities. J. Am. Chem. Soc. 2007, 129, 2880–2889. [Google Scholar] [CrossRef]
  13. Spyrakis, F.; Luque, F.J.; Viappiani, C. Structural analysis in nonsymbiotic hemoglobins: What can we learn from inner cavities? Plant Sci. 2011, 181, 8–13. [Google Scholar] [CrossRef] [PubMed]
  14. Abbruzzetti, S.; Spyrakis, F.; Bidon-Chanal, A.; Luque, F.J.; Viappiani, C. Ligand migration through hemeprotein cavities: Insights from laser flash photolysis and molecular dynamics simulations. Phys. Chem. Chem. Phys. 2013, 15, 10686–10701. [Google Scholar] [CrossRef] [PubMed]
  15. Watts, R.A.; Hunt, P.W.; Hvitved, A.N.; Hargrove, M.S.; Peacock, W.J.; Dennis, E.S. A hemoglobin from plants homologous to truncated hemoglobins of microorganisms. Proc. Natl. Acad. Sci. USA 2001, 98, 10119–10124. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.H.; Kochian, L.V.; Doyle, J.J.; Garvin, D.F. Two tomato non-symbiotic haemoglobin genes are differentially expressed in response to diverse changes in mineral nutrient status. Plant Cell Environ. 2003, 26, 673–680. [Google Scholar] [CrossRef]
  17. Shimoda, Y.; Nagata, M.; Suzuki, A.; Abe, M.; Sato, S.; Kato, T.; Tabata, S.; Higashi, S.; Uchiumi, T. Symbiotic rhizobium and nitric oxide induce gene expression of non-symbiotic hemoglobin in Lotus japonicus. Plant Cell Physiol. 2005, 46, 99–107. [Google Scholar] [CrossRef] [PubMed]
  18. Bustos-Sanmamed, P.; Tovar-Méndez, A.; Crespi, M.; Sato, S.; Tabata, S.; Becana, M. Regulation of nonsymbiotic and truncated hemoglobin genes of Lotus japonicus in plant organs and in response to nitric oxide and hormones. New Phytol. 2011, 189, 765–776. [Google Scholar] [CrossRef] [Green Version]
  19. Mira, M.M.; Hill, R.D.; Stasolla, C. Phytoglobins improve hypoxic root growth by alleviating apical meristem cell death. Plant Physiol. 2016, 172, 2044–2056. [Google Scholar] [CrossRef] [Green Version]
  20. Berger, A.; Guinand, S.; Boscari, A.; Puppo, A.; Brouquisse, R. Medicago truncatula phytoglobin 1.1 controls symbiotic nodulation and nitrogen fixation via the regulation of nitric oxide concentration. New Phytol. 2020, 227, 84–98. [Google Scholar] [CrossRef] [Green Version]
  21. Smagghe, B.J.; Trent, J.T., 3rd; Hargrove, M.S. NO dioxygenase activity in hemoglobins is ubiquitous in vitro, but limited by reduction in vivo. PLoS ONE 2008, 3, e2039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Rubio, M.C.; Calvo-Begueria, L.; Díaz-Mendoza, M.; Elhiti, M.; Moore, M.; Matamoros, M.A.; James, E.K.; Díaz, I.; Pérez-Rontomé, C.; Villar, I.; et al. Phytoglobins in the nuclei, cytoplasm and chloroplasts modulate nitric oxide signaling and interact with abscisic acid. Plant J. 2019, 100, 38–54. [Google Scholar] [CrossRef] [PubMed]
  23. Trevaskis, B.; Watts, R.A.; Andersson, C.R.; Llewellyn, D.J.; Hargrove, M.S.; Olson, J.S.; Dennis, E.S.; Peacock, W.J. Two hemoglobin genes in Arabidopsis thaliana: The evolutionary origins of leghemoglobins. Proc. Natl. Acad. Sci. USA 1997, 94, 12230–12234. [Google Scholar] [CrossRef] [Green Version]
  24. Hunt, P.W.; Watts, R.A.; Trevaskis, B.; Llewelyn, D.J.; Burnell, J.; Dennis, E.S.; Peacock, W.J. Expression and evolution of functionally distinct haemoglobin genes in plants. Plant Mol. Biol. 2001, 47, 677–692. [Google Scholar] [CrossRef]
  25. Hebelstrup, K.H.; Shah, J.K.; Igamberdiev, A.U. The role of nitric oxide and hemoglobin in plant development and morphogenesis. Physiol. Plant. 2013, 148, 457–469. [Google Scholar] [CrossRef] [PubMed]
  26. Vieweg, M.F.; Hohnjec, N.; Küster, H. Two genes encoding different truncated hemoglobins are regulated during root nodule and arbuscular mycorrhiza symbioses of Medicago truncatula. Planta 2005, 220, 757–766. [Google Scholar] [CrossRef] [PubMed]
  27. Larrainzar, E.; Villar, I.; Rubio, M.C.; Pérez-Rontomé, C.; Huertas, R.; Sato, S.; Mun, J.-H.; Becana, M. Hemoglobins in the legume–Rhizobium symbiosis. New Phytol. 2020, 228, 472–484. [Google Scholar] [CrossRef]
  28. Villar, I.; Larrainzar, E.; Milazzo, L.; Pérez-Rontomé, C.; Rubio, M.C.; Smulevich, G.; Martínez, J.I.; Wilson, M.T.; Reeder, B.; Huertas, R.; et al. A plant gene encoding one-heme and two-heme hemoglobins with extreme reactivities toward diatomic gases and nitrite. Front. Plant Sci. 2020, 11, 1743. [Google Scholar] [CrossRef]
  29. Kriegl, J.M.; Bhattacharyya, A.J.; Nienhaus, K.; Deng, P.; Minkow, O.; Nienhaus, G.U. Ligand binding and protein dynamics in neuroglobin. Proc. Natl. Acad. Sci. USA 2002, 99, 7992–7997. [Google Scholar] [CrossRef] [Green Version]
  30. Abbruzzetti, S.; Bruno, S.; Faggiano, S.; Grandi, E.; Mozzarelli, A.; Viappiani, C. Time-resolved methods in Biophysics. 2. Monitoring haem proteins at work with nanosecond laser flash photolysis. Photochem. Photobiol. Sci. 2006, 5, 1109–1120. [Google Scholar] [CrossRef] [PubMed]
  31. Bisht, N.K.; Abbruzzetti, S.; Uppal, S.; Bruno, S.; Spyrakis, F.; Mozzarelli, A.; Viappiani, C.; Kundu, S. Ligand migration and hexacoordination in type 1 non-symbiotic rice hemoglobin. Biochim. Biophys. Acta-Proteins Proteom. 2011, 1814, 1042–1053. [Google Scholar] [CrossRef] [PubMed]
  32. Marcelli, A.; Abbruzzetti, S.; Bustamante, J.P.; Feis, A.; Bonamore, A.; Boffi, A.; Gellini, C.; Salvi, P.R.; Estrin, D.A.; Bruno, S.; et al. Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca. PLoS ONE 2012, 7, e39884. [Google Scholar]
  33. Bustamante, J.P.; Abbruzzetti, S.; Marcelli, A.; Gauto, D.; Boechi, L.; Bonamore, A.; Boffi, A.; Bruno, S.; Feis, A.; Foggi, P.; et al. Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins. J. Phys. Chem. B 2014, 118, 1234–1245. [Google Scholar] [CrossRef]
  34. Knipp, M.; Ogata, H.; Soavi, G.; Cerullo, G.; Allegri, A.; Abbruzzetti, S.; Bruno, S.; Viappiani, C.; Bidon-Chanal, A.; Luque, F.J. Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7. F100 Research 2015, 4, 45. [Google Scholar] [CrossRef] [Green Version]
  35. Benabbas, A.; Ye, X.; Kubo, M.; Zhang, Z.; Maes, E.M.; Montfort, W.R.; Champion, P.M. Ultrafast dynamics of diatomic ligand binding to nitrophorin 4. J. Am. Chem. Soc. 2010, 132, 2811–2820. [Google Scholar] [CrossRef] [Green Version]
  36. Yu, A.; Ye, X.; Ionascu, D.; Cao, W.; Champion, P.M. Two-color pump-probe laser spectroscopy instrument with picosecond time-resolved electronic delay and extended scan range. Rev. Sci. Instrum. 2005, 76, 114301. [Google Scholar] [CrossRef]
  37. Boubeta, F.M.; Boechi, L.; Estrin, D.; Patrizi, B.; Di Donato, M.; Iagatti, A.; Giordano, D.; Verde, C.; Bruno, S.; Abbruzzetti, S.; et al. Cold-adaptation signatures in the ligand rebinding kinetics to the truncated hemoglobin of the antarctic bacterium Pseudoalteromonas haloplanktis TAC125. J. Phys. Chem. B 2018, 122, 11649–11661. [Google Scholar] [CrossRef]
  38. Giordano, D.; Boubeta, F.; Prisco, G.D.; Estrin, D.; Smulevich, G.; Viappiani, C.; Verde, C. Conformational flexibility drives cold adaptation in Pseudoalteromonas haloplanktis TAC125 globins. Antiox. Redox Signal. 2020, 32, 396–411. [Google Scholar] [CrossRef] [PubMed]
  39. Abbruzzetti, S.; Allegri, A.; Bidon-Chanal, A.; Ogata, H.; Soavi, G.; Cerullo, G.; Bruno, S.; Montali, C.; Luque, F.J.; Viappiani, C. Electrostatic tuning of the ligand binding mechanism by Glu27 in nitrophorin 7. Sci. Rep. 2018, 8, 10855. [Google Scholar] [CrossRef]
  40. Sottini, S.; Abbruzzetti, S.; Viappiani, C.; Ronda, L.; Mozzarelli, A. Determination of microscopic rate constants for CO binding and migration in myoglobin encapsulated in silica gels. J. Phys. Chem. B 2005, 109, 19523–19528. [Google Scholar] [CrossRef]
  41. Abbruzzetti, S.; Faggiano, S.; Bruno, S.; Spyrakis, F.; Mozzarelli, A.; Dewilde, S.; Moens, L.; Viappiani, C. Ligand migration through the internal hydrophobic cavities in human neuroglobin. Proc. Natl. Acad. Sci. USA 2009, 106, 18984–18989. [Google Scholar] [CrossRef] [Green Version]
  42. Uzan, J.; Dewilde, S.; Burmester, T.; Hankeln, T.; Moens, L.; Hamdane, D.; Marden, M.C.; Kiger, L. Neuroglobin and other hexacoordinated hemoglobins show a weak temperature dependence of oxygen binding. Biophys. J. 2004, 87, 1196–1204. [Google Scholar] [CrossRef] [Green Version]
  43. Jasaitis, A.; Ouellet, H.; Lambry, J.C.; Martin, J.L.; Friedman, J.M.; Guertin, M.; Vos, M.H. Ultrafast heme-ligand recombination in truncated hemoglobin HbO from Mycobacterium tuberculosis: A ligand cage. Chem. Phys. 2012, 396, 10–16. [Google Scholar] [CrossRef]
  44. Faggiano, S.; Abbruzzetti, S.; Spyrakis, F.; Grandi, E.; Viappiani, C.; Bruno, S.; Mozzarelli, A.; Cozzini, P.; Astegno, A.; Dominici, P.; et al. Structural plasticity and functional implications of internal cavities in distal mutants of type 1 non-symbiotic hemoglobin AHb1 from Arabidopsis thaliana. J. Phys. Chem. B 2009, 113, 16028–16038. [Google Scholar] [CrossRef] [PubMed]
  45. Anselmi, M.; DiNola, A.; Amadei, A. Kinetics of carbon monoxide migration and binding in solvated neuroglobin as revealed by molecular dynamics simulations and quantum mechanical calculations. J. Phys. Chem. B 2011, 115, 2436–2446. [Google Scholar] [CrossRef] [PubMed]
  46. Gabba, M.; Abbruzzetti, S.; Bruno, F.S.S.; Mozzarelli, A.; Luque, F.J.; Viappiani, C.; Cozzini, P.; Nardini, M.; Germani, F.; Bolognesi, M.; et al. CO rebinding kinetics and molecular dynamics simulations highlight dynamic regulation of internal cavities in human cytoglobin. PLoS ONE 2013, 8, e49770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Spyrakis, F.; Lucas, F.; Bidon-Chanal, A.; Viappiani, C.; Guallar, V.; Luque, F.J. Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2. Biochim. Biophys. Acta-Proteins Proteom. 2013, 1834, 1957–1967. [Google Scholar] [CrossRef]
  48. Brunori, M.; Giuffrè, A.; Nienhaus, K.; Nienhaus, G.U.; Scandurra, F.M.; Vallone, B. Neuroglobin, nitric oxide, and oxygen: Functional pathways and conformational changes. Proc. Natl. Acad. Sci. USA 2005, 102, 8483–8488. [Google Scholar] [CrossRef] [Green Version]
  49. Dewilde, S.; Ebner, B.; Vinck, E.; Gilany, K.; Hankeln, T.; Burmester, T.; Kreiling, J.; Reinisch, C.; Vanfleteren, J.R.; Kiger, L.; et al. The nerve hemoglobin of the bivalve mollusc Spisula solidissima: Molecular cloning, ligand binding studies, and phylogenetic analysis. J. Biol. Chem. 2006, 281, 5364–5372. [Google Scholar] [CrossRef] [Green Version]
  50. Manzoni, C.; Polli, D.; Cerullo, G. Two-color pump-probe system broadly tunable over the visible and the near infrared with sub-30 fs temporal resolution. Rev. Sci. Instrum. 2006, 77, 023103. [Google Scholar] [CrossRef]
  51. Barker, A.J.; Chen, K.; Hodgkiss, J.M. Distance distributions of photogenerated charge pairs in organic photovoltaic cells. J. Am. Chem. Soc. 2014, 136, 12018–12026. [Google Scholar] [CrossRef]
Figure 1. CO rebinding kinetics after nanosecond laser flash photolysis (LFP) of (a) WT1 and (b) WT4 solutions equilibrated with 1 atm CO (solid lines) and 0.2 atm CO (dotted lines) at 10 °C (black), 20 °C (red), 30 °C (green) and 40 °C (blue). Data are reported as the progress curve representing the fraction of deoxy molecules, N(t), as a function of time after photolysis. The absorbance change was monitored at 436 nm, following laser photolysis at 532 nm.
Figure 1. CO rebinding kinetics after nanosecond laser flash photolysis (LFP) of (a) WT1 and (b) WT4 solutions equilibrated with 1 atm CO (solid lines) and 0.2 atm CO (dotted lines) at 10 °C (black), 20 °C (red), 30 °C (green) and 40 °C (blue). Data are reported as the progress curve representing the fraction of deoxy molecules, N(t), as a function of time after photolysis. The absorbance change was monitored at 436 nm, following laser photolysis at 532 nm.
Ijms 22 02740 g001
Figure 2. CO rebinding kinetics for WT1 (gray dots) solution equilibrated with 1 atm CO and at room temperature collected with the hybrid femtosecond–nanosecond pump–probe setup. The corresponding kinetics measured after nanosecond LFP is reported as a gray solid line. For comparison, the rebinding kinetics to single and double Hisd mutants are reported (orange, 238; magenta, 74; olive, 74/238; see Section 2.4 for definitions).
Figure 2. CO rebinding kinetics for WT1 (gray dots) solution equilibrated with 1 atm CO and at room temperature collected with the hybrid femtosecond–nanosecond pump–probe setup. The corresponding kinetics measured after nanosecond LFP is reported as a gray solid line. For comparison, the rebinding kinetics to single and double Hisd mutants are reported (orange, 238; magenta, 74; olive, 74/238; see Section 2.4 for definitions).
Ijms 22 02740 g002
Scheme 1. Minimal kinetic scheme proposed for the reaction of CO with Glbs. Deoxy (pentacoordinate) and bis-histidyl (hexacoordinate) species are indicated by GlbP and GlbH, respectively. Immediately after photodissociation of GlbCO, the ligand is in a primary docking site (Glb:CO) within the distal pocket. From this site, the ligand can be rebound (rate k−1), migrate to a nearby temporary docking site (tr)1 with rate kc (reverse rate k−c), or exit to the solvent with rate kout. Ligands are rebound by GlbP from the solvent with rate kin, whereas binding and dissociation of Hisd occur with rates kh and k−h.
Scheme 1. Minimal kinetic scheme proposed for the reaction of CO with Glbs. Deoxy (pentacoordinate) and bis-histidyl (hexacoordinate) species are indicated by GlbP and GlbH, respectively. Immediately after photodissociation of GlbCO, the ligand is in a primary docking site (Glb:CO) within the distal pocket. From this site, the ligand can be rebound (rate k−1), migrate to a nearby temporary docking site (tr)1 with rate kc (reverse rate k−c), or exit to the solvent with rate kout. Ligands are rebound by GlbP from the solvent with rate kin, whereas binding and dissociation of Hisd occur with rates kh and k−h.
Ijms 22 02740 sch001
Figure 3. Analysis of the complete course of CO rebinding kinetics to WT1, WT4 and their mutants. The progress curve for the CO rebinding reaction after the hybrid pump–probe experiment is plotted as concentration of deoxy heme versus time. (a) Comparison between WT1 and WT4. The experimental data are superimposed on the fitted curves obtained using Scheme 1. (b,c) Replacement of Hisd by Leu affects the time evolution of the CO rebinding kinetics, as demonstrated by the comparison between WT1 and its mutant 74/238 in (b) and between WT4 and its mutant 109 in c. The fitting curves are superimposed on the experimental data (circles). In (b,c), the time courses of the other relevant species are shown using the same color code as for the kinetic model reported in Scheme 1. (d) Comparison between the kinetics of WT1 and its 74, 238 and 74/238 mutants.
Figure 3. Analysis of the complete course of CO rebinding kinetics to WT1, WT4 and their mutants. The progress curve for the CO rebinding reaction after the hybrid pump–probe experiment is plotted as concentration of deoxy heme versus time. (a) Comparison between WT1 and WT4. The experimental data are superimposed on the fitted curves obtained using Scheme 1. (b,c) Replacement of Hisd by Leu affects the time evolution of the CO rebinding kinetics, as demonstrated by the comparison between WT1 and its mutant 74/238 in (b) and between WT4 and its mutant 109 in c. The fitting curves are superimposed on the experimental data (circles). In (b,c), the time courses of the other relevant species are shown using the same color code as for the kinetic model reported in Scheme 1. (d) Comparison between the kinetics of WT1 and its 74, 238 and 74/238 mutants.
Ijms 22 02740 g003
Figure 4. Analysis of CO rebinding kinetics to the (a) 238 and (b) 74 mutants at room temperature, 1 atm CO. The fitting curves are superimposed on the experimental data (open circles). The time course of the relevant species shown in the kinetic scheme (Scheme 1) are reported: Glb:CO (violet), (tr)1 (light blue), GlbH (red) and GlbP (green).
Figure 4. Analysis of CO rebinding kinetics to the (a) 238 and (b) 74 mutants at room temperature, 1 atm CO. The fitting curves are superimposed on the experimental data (open circles). The time course of the relevant species shown in the kinetic scheme (Scheme 1) are reported: Glb:CO (violet), (tr)1 (light blue), GlbH (red) and GlbP (green).
Ijms 22 02740 g004
Table 1. Microscopic rate constants from the fit of pump–probe data at room temperature.
Table 1. Microscopic rate constants from the fit of pump–probe data at room temperature.
Rate
Constants
UnitsMtGlb1-2.1MtGlb1-2.4
WT17423874/238WT4109
k−1×108 s−12.0 ± 0.41.9 ± 0.42.6 ± 0.51.6 ± 0.314 ± 323 ± 4
kh×105 s−18 ± 13.8 ± 0.74.2 ± 0.8-7 ± 1-
k−h×105 s−13.8 ± 0.71.3 ± 0.31.6 ± 0.3-3.4 ± 0.7-
kout×108 s−11.4 ± 0.31.1 ± 0.20.9 ± 0.20.7 ± 0.15 ± 12.1 ± 0.4
kin×108 M−1s−11.4 ± 0.33.7 ± 0.73.4 ± 0.74.8 ± 0.91.0 ± 0.26.5 ± 0.9
kc×108 s−12.3 ± 0.42.9 ± 0.62.5 ± 0.52.3 ± 0.46.8 ± 0.93.3 ± 0.7
k−c×108 s−11.2 ± 0.21.8 ± 0.42.7 ± 0.52.5 ± 0.50.8 ± 0.21.7 ± 0.3
kon×108 M−1s−11.0 ± 0.21.5 ± 0.32.0 ± 0.41.4 ± 0.32.4 ± 0.417 ± 3
KH 2.2 ± 0.42.9 ± 0.62.5 ± 0.5-2.0 ± 0.4-
kobs×104 s−12.6 ± 0.52.7 ± 0.53.7 ± 0.7-6 ± 1-
Table 2. Activation energies (kcal/mol) determined from Arrhenius plots in the temperature range 10–40 °C. Microscopic rate constants are retrieved from the fit of LFP data.
Table 2. Activation energies (kcal/mol) determined from Arrhenius plots in the temperature range 10–40 °C. Microscopic rate constants are retrieved from the fit of LFP data.
Rate ConstantsMtGlb1-2.1MtGlb1-2.4
WT17423874/238WT4109
kout10.7 ± 0.48.2 ± 0.27.0 ± 0.56.1 ± 0.211 ± 1.04.1 ± 10.7
kin13.6 ± 0.911.7 ± 0.111 ± 1.07.8 ± 0.314 ± 2.0-
kh1.1 ± 0.22.7 ± 1.44.8 ± 0.5-1.4 ± 0.4-
k−h5.1 ± 0.24.3 ± 0.63.4 ± 0.8-2.6 ± 0.2-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abbruzzetti, S.; Barker, A.J.; Villar, I.; Pérez-Rontomé, C.; Bruno, S.; Cerullo, G.; Viappiani, C.; Becana, M. Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin. Int. J. Mol. Sci. 2021, 22, 2740. https://doi.org/10.3390/ijms22052740

AMA Style

Abbruzzetti S, Barker AJ, Villar I, Pérez-Rontomé C, Bruno S, Cerullo G, Viappiani C, Becana M. Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin. International Journal of Molecular Sciences. 2021; 22(5):2740. https://doi.org/10.3390/ijms22052740

Chicago/Turabian Style

Abbruzzetti, Stefania, Alex J. Barker, Irene Villar, Carmen Pérez-Rontomé, Stefano Bruno, Giulio Cerullo, Cristiano Viappiani, and Manuel Becana. 2021. "Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin" International Journal of Molecular Sciences 22, no. 5: 2740. https://doi.org/10.3390/ijms22052740

APA Style

Abbruzzetti, S., Barker, A. J., Villar, I., Pérez-Rontomé, C., Bruno, S., Cerullo, G., Viappiani, C., & Becana, M. (2021). Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin. International Journal of Molecular Sciences, 22(5), 2740. https://doi.org/10.3390/ijms22052740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop