Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

Reduction of metmyoglobin by inorganic disulfide species

The mechanism of the metal centered reduction of metmyoglobin (MbFeIII) by inorganic disulfide species has been studied by combined spectroscopic and kinetic analyses, under argon atmosphere. The process is kinetically characterized by biexponential time traces, for variable ratios of excess disulfide to protein, in the pH interval 6.6-8.0. Using UV-vis and resonance Raman spectroscopies, we observed that MbFeIII is converted into a low spin hexacoordinated ferric complex, tentatively assigned as MbFeIII(HSS-)/MbFeIII(SS2-), in an initial fast step. The complex is slowly converted into a pentacoordinated ferrous form, assigned as MbFeII according to the resonance Raman records. The reduction is a pH-dependent process, but independent of the initial disulfide concentration, suggesting the unimolecular decomposition of the intermediate complex following a reductive homolysis. We estimated the rate of the fast formation of the complex at pH 7.4 (kon = 3.7 × 103 M-1 s-1), and a pKa2 = 7.5 for the equilibrium MbFeIII(HSS-)/MbFeIII(SS2-). Also, we estimated the rate for the slow reduction at the same pH (kred = 10-2 s-1). A reaction mechanism compliant with the experimental results is proposed. This mechanistic study provides a differential kinetic signature for the reactions of disulfide compared to sulfide species on metmyoglobin, which may be considered in other hemeprotein systems.