Next Article in Journal
Optimization of the Acid Cleavage of Proanthocyanidins and Other Polyphenols Extracted from Plant Matrices
Next Article in Special Issue
Concentrations of Transition Metal Ions in Rat Lungs after Tobacco Smoke Exposure and Treatment with His-Leu Dipeptide
Previous Article in Journal
One-Pot Solvent-Involved Synthesis of 5-O-Substituted 5H-Chromeno[2,3-b]pyridines
Previous Article in Special Issue
Sequence–Activity Relationship of ATCUN Peptides in the Context of Alzheimer’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase

Division of Theoretical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(1), 65; https://doi.org/10.3390/molecules28010065
Submission received: 18 November 2022 / Revised: 14 December 2022 / Accepted: 19 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue The Role of Metal Ions in Bio-Inorganic Chemistry)

Abstract

:
We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron–sulfur clusters with 1–4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400–500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is Mo III Fe 3 II Fe 4 III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0–E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4–E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.

1. Introduction

Gaseous nitrogen (N2) is the main component of our atmosphere, but nitrogen is still a limiting element for plant life and a major ingredient in synthetic fertilisers. The reason for this is that the triple bond in N≡N is extremely strong and inert, making N2 unavailable for most plants [1,2,3]. Industrially, N2 is converted to NH3 by the Born–Haber process, which was invited in the early 20th century and is considered as a major cause of the human population explosion [1]. It requires high temperature and pressure.
N2 + 8 e + 8 H+ + 16 ATP → 2 NH3 + H2 + 16 ADP + 16 Pi
However, a few bacteria and archaea can perform the same conversion at ambient temperature and pressure. This is accomplished by the enzyme nitrogenase (EC 1.18/19.6.1) [2,4,5,6,7,8,9,10]. The reaction is still quite demanding, requiring 16 molecules of ATP, eight electrons and eight protons [2,4,5]: X-ray crystallographic studies have shown that nitrogenase contains two unusual iron–sulfur clusters [6,7,8,9,10]. The P-cluster, Fe8S7Cys6, is essentially two merged [4Fe–4S] clusters with a central sulfide ligand coordinating to six iron ions (Figure 1a). It is employed for electron transfer from the Fe-protein, which donates the electrons, to the active site. The latter is the FeMo cluster, which is a complicated MoFe7S9C(homocitrate) cluster, with a central carbide ion (Figure 1b) and connected to the protein by one cysteine and one histidine residue at the opposite ends of the trigonal prismatic cluster. There exist also alternative nitrogenases, in which the Mo ion is replaced by V or Fe [11].
Nitrogenase has been extensively studied by both experimental [2,3,4,5,12,13,14,15,16,17] and computational methods [18,19,20,21,22]. The reaction is normally described by the Lowe–Thorneley cycle [23], in which eight intermediates are recognised, E0–E8, based on the number of added electrons and protons. It has been shown that the resting E0 state needs to be reduced by three or four electrons before the N2 substrate can bind [2,3]. It is believed that the binding is facilitated by the reductive elimination of H2 from the cluster, explaining why H2 is a compulsory byproduct of the reaction. However, many details of the reaction are still unknown and many conflicting mechanisms have been suggested [3,18].
The Lowe–Thorneley cycle emphasizes the importance of electron and proton transfer in the reaction cycle of nitrogenase. The driving force of electron transfer is the redox potential. Unfortunately, it is hard to measure redox potentials of the FeMo cluster in nitrogenase because the reaction cannot be arrested at certain En states [24]. The only certain redox potential is between the resting E0 state and a one-electron oxidised state, which is outside the Lowe–Thorneley cycle (we will denote it E−1), −0.042 V [5,24,25,26]. For the reduction of the resting state, a redox potential of −0.45 to −0.49 V have been reported, but they may represent a mixture of reduced states [24,25,27,28,29].
Redox potentials can also be calculated by computational methods. However, the accuracy is limited. For redox sites of the same type in variants of the same protein (i.e., from different organisms or mutations), calculations based on the Poisson–Boltzmann equation or similar methods may give mean errors of 0.03–0.11 V for relative redox potentials [30,31,32,33,34,35,36,37,38]. For absolute potentials and sites of different types, quantum mechanical (QM) calculations are needed and the accuracy is appreciably worse. Typical errors are 0.2–0.6 V [30,31,32,33,34] and a prediction of the potential of the FeMo cluster had an error of 1.3 V, leading to incorrect identification of the central carbide ion [39]. Even all-atom QM molecular dynamics and free-energy calculations did not give an accuracy better than 0.26 V [40].
Recently, we performed a comparison and calibration of various combined QM and molecular mechanics (QM/MM) methods to estimate redox potentials of 13 iron–sulfur clusters with 1–4 Fe ions [41]. We showed that the best results were obtained by QM-cluster calculations in a continuum solvent with a high dielectric constant, using a large QM model (~300 atoms), based on QM/MM structures. With such an approach, we obtained a mean absolute error of 0.17 V, after removal of a systematic error of 0.62 V. The maximum error among the 13 studied potentials was 0.44 V. However, even if the accuracy is rather mediocre, it is enough to make useful predictions, e.g., what redox couple is employed by the [4Fe–4S] ferredoxins.
In this study, we employ this calibrated approach to study redox potentials in nitrogenase. We study four issues: First, we examine whether the computational method works also for the more complicated P- and FeMo clusters in nitrogenase, with redox potentials that also involve proton transfer. Second, we obtain an independent check of the redox state and charge state of the FeMo cluster. Third, we examine the recent suggestion that the E0–E4 states should operate at a nearly constant potential, employing only a single redox couple [42,43]. Fourth, we study the redox potentials of different reaction mechanisms after the binding of the substrate (i.e., for the E4–E8 states).

2. Result and Discussion

2.1. Redox Potentials of the P-Cluster

We started the investigation with testing our methodology [41] on the P-cluster to see if the results are reliable also for the large iron–sulfur clusters in nitrogenase and for redox reactions involving protonation of the clusters. The P-cluster (Figure 1a) contains eight Fe ions and the resting PN state has been shown to be the fully reduced Fe 8 II state [44,45,46]. Three additional states have been experimentally observed, oxidised by 1–3 electrons [44,45,46]. They are denoted P1+–P3+. Only the first two states are believed to be involved in the catalytic mechanism, although there are some evidence that also the P2+ state may be used [46,47,48,49]. Crystallographic studies have shown that in the PN state, the cluster is essentially two [4Fe–4S] clusters merged by the S1 sulfide ion that coordinates to six Fe ions (cf. Figure 1a). In the P2+ state, Ser-188D becomes deprotonated and coordinates to Fe6 [48,50,51]. Likewise, the backbone N atom of Cys-88 (the sidechain of which is one of the ligands to the P-cluster) also becomes deprotonated and coordinates to Fe5 [48,50,51]. This leads to cleavage of the Fe5–S1 and Fe6–S1 bonds. The same structure is expected for the P3+ state. The structure of the P1+ state is more uncertain, because redox titrations indicated that only the P2+/P1+ redox-couple is pH-dependent, whereas no evidence was found for a coupled electron- and proton-transfer for the P1+/PN couple [52,53]. However, a recent crystal structure was interpreted to contain deprotonated Ser-188D and protonated Cys-88, although it is probably a mixture of the P1+ and the P2+ states [54,55,56].
Six states of the P-cluster were considered in this study as are described in Table 1. Each state was first QM/MM optimised with the small QM system and then subjected to a single-point energy calculation with the large QM system in a COSMO continuum solvent with a dielectric constant of 80. The QM/MM structures and the best BS states were taken from our previous publication [55]. The calculated redox potentials are listed in Table 2 and they are compared to experimentally measured redox potentials [5] in order to gauge the accuracy of the method when applied to the P-clusters of the nitrogenase. For the P1+ state, we tested three different protonation states.
It can be seen that the redox couples P2+ → P3+ and PNH2 → P1+H give errors of 0.05 and 0.33 V, respectively, compared to experimental potentials. These are within the range of errors observed in our previous study (maximum error 0.44 V) [41]. This is quite satisfying, especially considering that one of the considered redox potentials for the P-cluster involves coupled redox and protonation reaction, whereas the calibration study involved only pure redox reactions. This gives us confidence to apply the method also to the FeMo cluster.
For the P1+/PN redox potential, the calculated result is essentially independent on whether we include the proton transfer or not in the calculations (the PNH2 → P1+H2 and PNH2 → P1+H give calculated redox potentials that differ by only 0.02 V). This is in agreement with the experimental observation that the P1+/PN redox potential is pH-independent [52,53] and may solve the enigma why redox titrations did not observe any pH-dependence although crystal structures indicate that a deprotonation should be involved in this redox reaction [54,55,56].
For the P1+H → P2+ transition, the comparison with experimental results is somewhat harder, because the measured redox potential changes with pH, from −0.224 V at pH 6.0 to −0.348 V at pH 8.5 [52]. This has been interpreted to reflect the deprotonation of the backbone amide group of Cys-88: At very low pH, it is protonated in both the P1+ and P2+ states, whereas at high pH, it is deprotonated in both redox states. At intermediate pH (i.e., in the measured range), it is protonated in the P1+ state and deprotonated in the P2+ state [52]. Our calculations confirm these suggestions: We obtain a more positive potential for the P2+/P1+H couple ( E corr o = 0.1 V) than for the P2+/P1+ couple ( E corr o = −1.1 V). However, the restricted pH range of the measured potentials makes it hard to make a more detailed judgement of the calculated potentials. The experiments indicate that the redox potential of the pure P1+H → P2+ transition is larger than −0.224 V, indicating an error of less than 0.36 V for our calculated potential. Likewise, the redox potential of the pure P1+ → P2+ transition is more negative than −0.348 V, indicating an error of less than 0.72 V. Thus, our calculations confirm the experimental observation that the P1+ → P2+ transition involves a proton transfer [52,53], but we cannot obtain any detailed estimate of the error for this redox couple.

2.2. Accuracy of the Redox Potential of the FeMo Cluster and Oxidation Level of the Cluster

As mentioned above, only two experimental redox potentials have been reported for the FeMo cluster [24,25]. A potential of −0.042 V has been measured between the resting state and a one-electron oxidised state [26,45] (not involved in the normal reaction mechanism). Redox potentials for more reduced states are harder to measure, because reduction starts substrate or proton reduction. No firm measurement is available, but the potential between the resting state and a steady-state reduced state (which my represent more than one reduced state) has been estimated in four studies between −0.45 and −0.49 V [24,25,27,28,29]. In two cases, another potential of −0.30 to −0.32 V was also reported, but it may come from the P-cluster [24].
For the FeMo cluster, we wanted to study two related questions, viz. whether we can reproduce these potentials with our QM calculations, within the accuracy of the method, and whether we can identify the correct redox couple of the FeMo cluster. Recent Mössbauer, anomalous dispersion and QM investigations have suggested that the resting E0 state of the FeMo cluster is in the Mo III Fe 3 II Fe 4 III oxidation state [13,42,57,58]. However, this gives a large net negative charge of the cluster and its direct ligands (−5 e), which is only partly neutralised by two nearby arginine residues. Moreover, protonation energies of various groups of the FeMo cluster are unexpectedly large [59,60]. Therefore, an independent confirmation of the oxidation-state assignment is desirable.
To this end, twelve states of the FeMo cluster were considered, described in Table 3. We studied the resting E0 state, together with the one-electron oxidised and one-electron reduced states. For the latter, we considered both a structure with no protons added and a state with a proton added on the S2B µ2 bridging sulfide ion, as has been suggested by several QM investigations [18,59] and is also supported by experimental studies [61,62]. For the standard oxidation-state assignment ( Mo III Fe 3 II Fe 4 III for E0), these states are denoted E0, E−1, E1 and E1H. In addition, we considered two alternative charge states, either with two more or two electrons less (experimentally the resting state is a quartet [2,3], i.e., with an odd number of electrons, so electrons need to be added or removed in pairs). These states are called A0, A−1, A1 and A1H when two electrons were added (i.e., giving Mo III Fe 5 II Fe 2 III for the A0 state) and R0, R−1, R1 and R1H when two electrons were removed (i.e., giving Mo III Fe 1 II Fe 6 III for the R0 state; note that R1 = E−1 and A−1 = E1). The calculated redox potentials are listed in Table 4.
It can be seen that with the standard charge state ( Mo III Fe 3 II Fe 4 III for E0), our calculations reproduce the two experimental redox potentials with errors of 0.17 and 0.19 V, i.e., well within the error range observed in our previous study (maximum error 0.44 V) [41] and also for the P-cluster. However, the good results are observed only if it is assumed that the reduction of E0 is accompanied by the uptake of by a proton (i.e., E0 → E1H; for E0 → E1, the error is 0.6 V), showing that the calculations confirm that a proton transfer is involved in the redox reaction.
If we instead consider a FeMo cluster with two electrons less (i.e., with a Mo III Fe 1 II Fe 6 III assignment for the resting state, here called R0), the calculated potentials reproduce the experimental redox potentials worse: The R−1 → R0 transition gives and error of 2.0 V, much larger than the maximum error in our previous study [41], whereas R0 → R1H gives an error of 0.43 V.
Likewise, if we instead add two extra electrons ( Mo III Fe 5 II Fe 2 III for the resting state), we get a very large error (over 6 V) for the A−1 → A0 potential. For the A0 → A1H transition, the error is smaller, −0.48 V, but it is still somewhat larger than the maximum error observed in our calibration study. Thus, our calculations confirm that Mo III Fe 3 II Fe 4 III is the proper redox assignment for E0. Apparently, the E−1/E0 redox potential is more sensitive to the involved redox couple than the E0/E1H potential. The calculations also confirm that the E1 state is protonated.

2.3. Redox Potentials of the E0–E4 States of the FeMo Cluster

Next, we studied also the E0–E4 states of the FeMo cluster from Lowe–Thorneley reaction cycle. The aim was to examine the suggestion that all the E0–E4 states should have similar redox potentials, because they use only two formal redox states, Mo III Fe 3 II Fe 4 III and Mo III Fe 4 II Fe 3 III [42,43]. This is suggested to be accomplished by the added protons, which bind to Fe ions in the E2 and E4 states, thereby formally becoming hydride ions and changing the oxidation state of the Fe ions by two. Thus, the five states E0–E4 would formally be Fe 3 II Fe 4 III , Fe 4 II Fe 3 III H + , Fe 3 II Fe 4 III H + H , Fe 4 II Fe 3 III H 2 + H and Fe 3 II Fe 4 III H 2 + H 2 (leaving out Mo, which always is in the +III state).
We used mainly structures from previous studies [59,63,64,65,66,67] and included a few alternative structures for each state (except for E0 and E1H, for which there is a reasonable consensus), to see if we can discriminate between different possibilities using the redox potentials. The various structures are described in Table 5 (they are also shown in Figure S3) and the calculated redox potentials are listed in Table 6. We use the redox potential of the E0/E1H couple as a reference ( Δ E calc o = E calc o Δ E calc o E 0 / E 1 H ) to judge if all transitions have similar redox potentials.
In our previous studies, we found that with the TPSS functional, the best E2H2 structure has a proton on S2B and a hydride ion bridging Fe2 and Fe6 [59,64]. It can be seen that the estimated redox potential for the E1H → E2H2 transition is similar to that of the E0 → E1H transition, only 0.08 V less negative. This confirms the conjecture that the FeMo cluster in the early En states should have a constant redox potential.
If we instead use another structure for the E2H2 state with the extra two protons bound to S2B and as a hydride ion bridging Fe2 and Fe6, but with S2B dissociated from Fe2 (but not from Fe6; called E2H2′), the redox potential for the E1H → E2H2′ transition decreases by 0.25 V. This simply reflects that with TPSS-D4/def2-SV(P) and the 399-atom redox model, the E2H2′ structure is 0.25 eV (24 kJ/mol) less stable than the E2H2 structure. For the 184-atom QM/MM model, the difference is 15 kJ/mol. In our previous study, we showed that the relative stability of the structures with the protonated S2B group bridging Fe2 and Fe6 or dissociated from one of the two iron ions depends on what DFT functional is used [64]. For example, with the TPSSh functional, the E2H2′ structure is instead 11 kJ/mol more stable. Consequently, the redox potentials will also depend on the QM method used, with differences of ~0.3 V.
If we instead use a structure for the E2H2 state with the two protons on S2B and S5A (E2H2″), the calculated redox potential is 0.08 V more negative than with the most stable E2H2 state. Again, this reflects that this protonation state is 9 kJ/mol less stable. However, it also shows that the formal oxidation states in the FeMo cluster have only a minor influence on the redox potential. The E2H2″ state involves two protons on sulfide groups and therefore represents a doubly reduced formal Fe 5 II Fe 2 III H 2 + state in contrast to the Fe 3 II Fe 4 III H + H state for E2H2 and E2H2′. This shows that as long as the various protonation states have similar relative energies, they will also have similar redox potentials, showing that there is no major difference between different formal oxidation states of the Fe ions. This is in line with suggestions by Dance, who has pointed out that it is misleading to make a sharp contrast between protons on sulfides and hydride ions on Fe ions, because there is only a small difference in the charge populations on the H atom [18].
For the E3H3 state, we used a structure with the third proton bridging Fe3 and Fe7 (in addition a proton on S2B and a hydride ion bridging Fe2 and Fe6). The estimated redox potential for the E2H2 → E3H3 transition is 0.19 V more negative than that of the E0 → E1H transition. This is well within the maximum error in our calibration study 0.44 V [41], and therefore still in agreement with the expectation that the electrons can be donated to the FeMo cluster at a constant redox potential. We tested also another structure, taken from our previous systematic study [59], in which the third proton bound to Fe5 instead (E3H3′). It gave a 0.35 V more negative redox potential (i.e., further away from that of the E0 → E1H transition), reflecting that this structure is less stable. Interestingly, we found in contrast to our previous study that the broken-symmetry BS-14 state was more favourable for both these structures with the large QM model used in the redox calculations (but only for the best structure with the smaller QM/MM-optimised model).
For the E3H3 → E4H4 transition, our estimated redox potential is 0.41 V more positive than that of the E0 → E1H transition. This is within the maximum error in our calibration study [41], but considering that the redox potential of the E2H2/E3H3 couple was a bit to negative and this redox potential is a bit too positive, it might indicate that we have not yet found the best structure for the E3 state. For E4H4, we employed the best structure in our previous investigation of this state [63], viz. a structure with two protons on S2B and S5A and two hydride ions bridging Fe2/6 and Fe3/7. There are several possible conformations of such a structure. The best one has all H atoms pointing towards S3A, except the Fe3/7 hydride, which points towards S2B (all structures are shown in Figure S3). If we instead use a structure with the Fe2/6 hydride on the other side of S2B, i.e., the same face as the Fe2/6 hydride, the redox potential becomes 0.23 V more negative, reflecting that such a structure is 22 kJ/mol less stable. If we instead use the structure suggested by Hoffman and coworkers [2,21], i.e., with all four H atoms on the same face of the cluster (i.e., the two protons on S2B and S5A pointing in the opposite direction compared to E4H4 and E4H4′ structures), the redox potential becomes even more negative by 0.24 V. Likewise, if we use the best structure in our first investigation of the E4H4 state [59] (with H atoms on S2B, Fe2/6, Fe5 and Fe6), the redox potential becomes 0.15 V even more negative, reflecting that this structure is 59 kJ/mol less stable than the best state.
In conclusion, we find that for all four calculated redox potentials for the E0–E4 states are similar, within 0.41 V, i.e., within the accuracy of our method. This confirm the expectation that the redox potentials should be similar [42,43] so that they can accept electrons from the same source. However, the results are very sensitive to which structures are employed and which QM method and broken-symmetry state is used. On the other hand, we show that the formal oxidation states of the Fe ions (the number of protons on sulfide ions or hydride ions on Fe) is less important for the redox potentials, contrary to the suggestion that the redox potential depends on the formal oxidation state of the Fe ions [43,68].

2.4. Redox Potentials of the E4–E8 States of the FeMo Cluster

Finally, we studied also redox potentials of the FeMo cluster in the later part of the reaction, after binding of N2 and its protonation to N2H2. In previous studies, we have suggested thermodynamically stable structures for the E4N2H2 to E8NH3 states, either with S2B bound or dissociated from the cluster [65,66,67]. We use these structures also in this study. They are described in Table 5 and are shown in Figure S3. In Table 6, the calculated redox potentials for both scenarios are presented.
It can be seen that both with and without S2B, the calculated redox potentials are all less negative than that for the E0/E1H couple, by 0.2–2.65 V. This reflects that once N2H2 has been formed, the following reactions are quite facile. Electrons tend to move towards sites with a more positive redox potential. Therefore, redox potentials that are more positive than those of the E0/E1H couple indicate that the electron transfer is more exothermic than in the E0 → E1 step. Thus, the electron transfers of the E4–E8 steps of the nitrogenase reaction are more downhill than those of the E0–E4 reactions. However, it might also indicate that the assumption that the bound N2 directly is protonated may be incorrect. In fact, the first protonation of N2 is the hardest step in the reduction of N2 to ammonia [69], and it is possible that it requires further reduction of the FeMo cluster before it is feasible (this part of the reaction was not studied in our previous studies). Our results indicate that this should be further studied.
The relative sizes of the four redox potentials for the E4–E8 states are also rather independent on whether S2B remains bound or is dissociated: The potentials of the first and third steps E4 → E5 and E6 → E7 are most positive, especially the latter, whereas the redox potentials of the other two steps are closer to that of the E0/E1H couple. For the four reductions with S2B still bound, our calculated redox potentials are 1.13, 0.41, 2.02 and 0.57 V, whereas with S2B dissociated, the four calculated redox potentials are 0.97, 0.20, 2.65 and 0.19 V (i.e., with a somewhat larger variation). The similarity of the trends for the two sets of potentials is conspicuous considering that the N–N bond is cleaved in the E5N2H3 → E6N2H4 transition with S2B bound, but in the E6N2H4 → E7N2H5 transition when S2B has dissociated.
In conclusion, our results show that later part of the reaction mechanism of Mo-nitrogenase give redox reactions that are more exothermic than that of the E0/E1H redox couple. Thus, we see no evidence that a stronger driving force is needed for the reaction, as has been suggested by Siegbahn [22]. Moreover, there is no large difference between the mechanisms with S2B bound or dissociated.

3. Methods

3.1. The Protein

The calculations were based on the 1.0-Å crystal structure of Mo nitrogenase from Azotobacter vinelandii (PDB code 3U7Q) [8]. The setup of the protein is identical to that of our previous studies [60,63,69,70]. The entire heterotetramer was considered in the calculations and the quantum mechanical (QM) calculations were concentrated on the FeMo clusters in the C subunit because there is a buried imidazole molecule from the solvent rather close to the active site (~11 Å) in the A subunit. The metal clusters not involved in the QM calculations were modelled by MM in the fully reduced and E0 resting states, respectively, using a QM charge model [70]. The protonation states of all residues were the same as before [70], and the homocitrate ligand was modelled in the singly protonated state with a proton shared between the hydroxyl group (O7 that coordinates to Mo) and the O1 carboxylate atom [57,70]. The protein was solvated in a sphere with a radius of 65 Å around the geometrical centre of the protein. Cl and Na+ ions were added to an ionic strength of 0.2 M [71]. The final system contained 133 915 atoms. For the protein, we used the Amber ff14SB force field [72] and water molecules were described by the TIP3P model [73]. The metal sites [70,74] were treated by a non-bonded model [75] and charges were obtained with the restrained electrostatic potential method [76].

3.2. QM Calculations

All QM calculations were performed with the Turbomole software (version 7.5) [77]. All calculations were performed with the TPSS [78] functional with the def2-SV(P) basis set [79], a combination that gave the best relative redox potentials in our previous study (lowest mean absolute deviation and maximum error, after removal of the mean signed error; B3LYP or larger basis sets gave more than twice as large mean absolute deviations and maximum errors, but a slightly smaller mean signed error) [41]. The calculations were sped up by expanding the Coulomb interactions in an auxiliary basis set, the resolution-of-identity (RI) approximation [80,81]. Empirical dispersion corrections were included with the DFT-D4 approach [82], as implemented in Turbomole. QM calculations were performed on both the FeMo cluster and the P-cluster, and two different sizes of the QM systems were employed, one smaller for QM/MM geometry optimisations and one larger for the redox-potential calculations.
In the QM/MM geometry optimisations, the P-cluster was modelled as Fe8S7Cys6, with the five of the Cys ligands modelled by CH3S, whereas Cys–88 was modelled by CH3CONHCH2CH2S, because the backbone amide group is deprotonated and coordinates to Fe5 in some the more oxidised states. Likewise, Ser-188D (i.e., belonging to subunit D, rather than C for all the other numbered residues) was included in the model as CH3OH because it is deprototated and coordinates to Fe6 in the oxidised states. The model contained 64 atoms for the fully protonated state and it is shown in Figure S1a in the Supporting Information. All QM/MM structures of the P-cluster were taken from our previous study [55].
The FeMo cluster was modelled by MoFe7S9C(homocitrate)(CH3S) (imidazole) in the QM/MM calculations, where the two last groups are models of Cys-275 and His-442. In addition, all groups that form hydrogen bonds to the FeMo cluster were also included in the QM model, viz. Arg-96, Gln-191 and His-195 (sidechains), Ser-278 and Arg-359 (both backbone and sidechain, including the CA and C and O atoms from Arg-277), Gly-356, Gly-357 and Leu-358 (backbone, including the CA and C and O atoms from Ile-355), as well as two water molecules. Finally, the sidechains of Glu-380, Val-70 and Phe-381 were also included. This QM system involved 191 atoms and is shown in Figure S1b.
For the redox-potential calculations, we used the largest QM system suggested in our previous investigation [41]. It included all functional groups in the proteins with any atom within 3.5 Å of a minimal QM system, consisting of all metal and sulfide ions together with all direct ligands (Cys, His and homocitrate). These QM systems were set up using our local program for BigQM calculations (changepdb) [83]. They contained ~400 atoms for the FeMo cluster and ~500 atoms for the P-cluster, and are shown in Figure S2. The calculations were based on the QM/MM-optimised structures.
In the redox calculations, the QM system was immersed into a continuum solvent, employing the conductor-like screening model (COSMO) [84,85] implemented in Turbomole. The default optimised COSMO atomic radii and a water solvent radius of 1.3 Å were employed to construct the solvent-accessible surface cavity [86], whereas a radius of 2.0 Å was used for Fe and Mo [87]. Structures for the QM + COSMO calculations were taken directly from the QM/MM calculations without further optimisation. The dielectric constant was 80, which gave the best redox potentials in our previous study [41].
Redox potentials (Eº) were calculated according to
Eº = E(ox) − E(red) − c
where E(ox) and E(red) are the energies of the oxidised and reduced states, and c is a correction factor (4.28 eV) to place the potentials on the scale of the standard hydrogen electrode [88]. The actual value of this factor has been much discussed and values between 4.05–4.44 eV have been suggested [88,89]. To avoid this problem, we use the method calibrated in our previous study to 13 different iron–sulfur clusters [41]. Therefore, we subtract from the c constant the mean signed error (MSE) obtained in this calibration study, MSE = −0.62 V or we subtract the calculated potential of the E0 → E1 transition from the redox potentials of the other states.
In some cases, it is known or assumed that the cluster takes up a proton during or after the electron transfer. In those cases, it was assumed that the proton comes from an imidazole molecule, studied with the same QM method in a COSMO continuum solvation solvent with the dielectric constant of water (80). Thus, we assume that the proton comes from a group with a pKa close to 7 (6.95). This changes the calculated redox potentials by −12.47 V.
The electronic structure of all QM systems was obtained with the broken-symmetry (BS) approach [90]: Each of the seven or eight Fe ions was modelled in the high-spin state, with either a surplus of α or β spin. Such a state can be selected in many different ways, giving rise to different BS states, which are specified by giving the number of the Fe ions with minority spin (the numbering of the Fe ions is shown in Figure 1) [55,60]. The various BS states were obtained either by swapping the coordinates of the Fe ions [91] or with the fragment approach by Szilagyi and Winslow [92].

3.3. QM/MM Calculations

QM/MM calculations were performed with the ComQum software [93,94]. In this approach, the protein and solvent are split into two subsystems: System 1 (the QM region) was relaxed by QM methods. System 2 contained the remaining part of the protein and the solvent, and it was kept fixed at the original coordinates (equilibrated crystal structure [70], to avoid the risk that different calculations end up in different local minima).
In the QM calculations, system 1 was represented by a wavefunction, whereas all the other atoms were represented by an array of partial point charges, one for each atom, taken from the MM setup. Thereby, the polarisation of the QM system by the surroundings is included in a self-consistent manner (electrostatic embedding). When there is a bond between systems 1 and 2 (a junction), the hydrogen link-atom approach was employed: The QM system was capped with hydrogen atoms, the positions of which are linearly related to the corresponding carbon atoms (carbon link atoms, CL) in the full system [93,95]. All atoms were included in the point-charge model, except the CL atoms [96]. ComQum employs a subtractive scheme with van der Waals link-atom corrections [97]. No cut-off is used for the QM and QM–MM interactions. The geometry optimisations were continued until the energy change between two iterations was less than 2.6 J/mol (10−6 a.u.) and the maximum norm of the Cartesian gradients was below 10−3 a.u.

4. Conclusions

In this study, we have investigated what information we can get from calculated redox potentials of the two metal clusters in Mo-nitrogenase. We employ our calibrated approach to calculate redox potentials for iron–sulfur clusters involving QM-cluster calculations in a continuum solvent with large QM models (400–500 atoms), based on structures from QM/MM optimisations [41]. We obtain several interesting results:
The calculations on the P-cluster show that our method gives approximately the same accuracy as for the simple iron–sulfur clusters with 1–4 Fe ions, with a maximum error of 0.33 V (0.44 V in our previous study [41]). This shows that the calculations are accurate also for the larger P-cluster and for redox reactions that include proton transfers.
The calculations confirm that the P1+ → P2+ transition involves a proton transfer (i.e., P1+H → P2+), as is also suggested by electrochemical and crystallographic studies [8,48,50,51,54,55,56].
The calculations show that the P1+H2/PNH2 and P1+H/PNH2 redox couples give very similar redox potentials, which may explain the experimental enigma that redox titrations do not show any pH dependence of the P1+/PN redox potential [52,53], although crystal structures indicate that also the PN → P1+ transition should involve a proton transfer [54,55,56].
For the FeMo-cluster, the calculations unambiguously identify Mo III Fe 3 II Fe 4 III as the proper formal oxidation state for the resting E0 state of the protein. This provides an independent confirmation of this oxidation state, also suggested by previous experimental and QM studies [13,42,57,58].
The calculations agree with experiments only if it is assumed that a proton is taken up together with the electron in the E0 → E1 reaction.
The calculations confirm that the E0 to E4H4 states all have similar redox potentials (within 0.41 V, i.e., lower than the estimated maximum error of the method), as expected for a site that should receive electrons from the same donor.
However, there is no major difference in the redox potentials between structures with protons on the µ2-bridging sulfide ions or hydride ions on the Fe ions (for E2H2 the difference is only 0.08 V). This shows that there is only minor differences between hydrogen atoms bound to S or Fe ions, as previously has been advocated by Dance [18] and that the formal oxidation states of the Fe ions are no good indicators of the redox potentials.
The redox potentials of the later steps of the reaction mechanism (E4N2H2 to E8) are more positive than that of the resting state (E0/E1H), showing that the reactions are more exothermic. The trends in the potentials do not change if S2B remains bound to the cluster or if it dissociates. This shows that there is no need of a significantly more negative potential for the Fe protein than measured.
The calculated redox potentials strongly depend on the structures used (the positions of the added protons), the broken-symmetry states and the QM method employed.
Thus, our calculations show that quite strong predictions can be provided by redox-potential calculations, even if the accuracy is rather poor compared to experimental measurements (a maximum error of 0.44 V), provided that the calculations are calibrated so that the expected errors are known. Such calculations can also identify possible problems in suggested reaction mechanisms. Our calculations also allow us to identify what redox reactions involve coupled electron and proton transfer, which is crucial to identify the detailed reaction mechanism of nitrogenase. Moreover, we have been able to explain the enigma why the P1+/PN redox potential is not pH dependent [52,53], although the crystal structures show that the a deprotonation reaction is involved [54,55,56].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28010065/s1, Figure S1: Structures of the P-cluster and the FeMo cluster, illustrating the QM systems used in the QM/MM geometry optimisations, as well as the names of the nearby residues; Figure S2: Models used for the redox calculations of the P-cluster and the FeMo cluster; Figure S3: Structures used for the E0–E8 states of the FeMo cluster.

Author Contributions

Conceptualization, U.R.; methodology, U.R.; formal analysis, U.R.; investigation, H.J. and O.K.G.S.; resources, U.R.; writing—original draft preparation, O.K.G.S.; writing—review and editing, H.J., O.K.G.S. and U.R.; visualization, H.J..; supervision, H.J. and U.R.; project administration, U.R.; funding acquisition, U.R. All authors have read and agreed to the published version of the manuscript.

Funding

This investigation has been supported by grants from the Swedish Research Council (project 2018-05003) and from China Scholarship Council. The computations were performed on computer resources provided by the Swedish National Infrastructure for Computing (SNIC) at Lunarc at Lund University, NSC at Linköping University and HPC2N at Umeå University, partially funded by the Swedish Research Council (grant 2018-05973).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Smith, B.E. Nitrogenase Reveals Its Inner Secrets. Science 2002, 297, 1654–1655. [Google Scholar] [CrossRef] [PubMed]
  2. Hoffman, B.M.; Lukoyanov, D.; Yang, Z.-Y.; Dean, D.R.; Seefeldt, L.C. Mechanism of Nitrogen Fixation by Nitrogenase: The Next Stage. Chem. Rev. 2014, 114, 4041–4062. [Google Scholar] [CrossRef]
  3. Seefeldt, L.C.; Yang, Z.-Y.; Lukoyanov, D.A.; Harris, D.F.; Dean, D.R.; Raugei, S.; Hoffman, B.M. Reduction of Substrates by Nitrogenases. Chem. Rev. 2020, 120, 5082–5106. [Google Scholar] [CrossRef] [PubMed]
  4. Burgess, B.K.; Lowe, D.J. Mechanism of Molybdenum Nitrogenase. Chem. Rev. 1996, 96, 2983–3012. [Google Scholar] [CrossRef] [PubMed]
  5. Schmid, B.; Chiu, H.-J.; Ramakrishnan, V.; Howard, J.B.; Rees, D.C. Nitrogenase. In Handbook of Metalloproteins; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2006; pp. 1025–1036. ISBN 9780470028636. [Google Scholar]
  6. Kim, J.; Rees, D.C. Structural Models for the Metal Centers in the Nitrogenase Molybdenum-Iron Protein. Science 1992, 257, 1677–1682. [Google Scholar] [CrossRef] [PubMed]
  7. Einsle, O.; Tezcan, F.A.; Andrade, S.L.A.; Schmid, B.; Yoshida, M.; Howard, J.B.; Rees, D.C. Nitrogenase MoFe-Protein at 1.16 Å Resolution: A Central Ligand in the FeMo-Cofactor. Science 2002, 297, 1696. [Google Scholar] [CrossRef] [Green Version]
  8. Spatzal, T.; Aksoyoglu, M.; Zhang, L.; Andrade, S.L.A.; Schleicher, E.; Weber, S.; Rees, D.C.; Einsle, O. Evidence for Interstitial Carbon in Nitrogenase FeMo Cofactor. Science 2011, 334, 940. [Google Scholar] [CrossRef] [Green Version]
  9. Spatzal, T.; Perez, K.A.; Einsle, O.; Howard, J.B.; Rees, D.C. Ligand Binding to the FeMo-Cofactor: Structures of CO-Bound and Reactivated Nitrogenase. Science 2014, 345, 1620–1623. [Google Scholar] [CrossRef] [Green Version]
  10. Einsle, O. Nitrogenase FeMo Cofactor: An Atomic Structure in Three Simple Steps. J. Biol. Inorg. Chem. 2014, 19, 737–745. [Google Scholar] [CrossRef]
  11. Eady, R.R. Structure−Function Relationships of Alternative Nitrogenases. Chem. Rev. 1996, 96, 3013–3030. [Google Scholar] [CrossRef]
  12. Lancaster, K.M.; Roemelt, M.; Ettenhuber, P.; Hu, Y.; Ribbe, M.W.; Neese, F.; Bergmann, U.; DeBeer, S. X-Ray Emission Spectroscopy Evidences a Central Carbon in the Nitrogenase Iron-Molybdenum Cofactor. Science 2011, 334, 974–977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Bjornsson, R.; Lima, F.A.; Spatzal, T.; Weyhermüller, T.; Glatzel, P.; Bill, E.; Einsle, O.; Neese, F.; DeBeer, S. Identification of a Spin-Coupled Mo(III) in the Nitrogenase Iron–Molybdenum Cofactor. Chem. Sci. 2014, 5, 3096–3103. [Google Scholar] [CrossRef] [Green Version]
  14. Igarashi, R.Y.; Laryukhin, M.; Dos Santos, P.C.; Lee, H.-I.; Dean, D.R.; Seefeldt, L.C.; Hoffman, B.M. Trapping H Bound to the Nitrogenase FeMo-Cofactor Active Site during H2 Evolution: Characterization by ENDOR Spectroscopy. J. Am. Chem. Soc. 2005, 127, 6231–6241. [Google Scholar] [CrossRef] [PubMed]
  15. Hoeke, V.; Tociu, L.; Case, D.A.; Seefeldt, L.C.; Raugei, S.; Hoffman, B.M. High-Resolution ENDOR Spectroscopy Combined with Quantum Chemical Calculations Reveals the Structure of Nitrogenase Janus Intermediate E4(4H). J. Am. Chem. Soc. 2019, 141, 11984–11996. [Google Scholar] [CrossRef]
  16. Roth, L.E.; Tezcan, F.A. X-ray Crystallography. Methods Mol. Biol. 2011, 766, 147–164. [Google Scholar] [CrossRef]
  17. Van Stappen, C.; Decamps, L.; Cutsail, G.E.; Bjornsson, R.; Henthorn, J.T.; Birrell, J.A.; DeBeer, S. The Spectroscopy of Nitrogenases. Chem. Rev. 2020, 120, 5005–5081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Dance, I. Computational Investigations of the Chemical Mechanism of the Enzyme Nitrogenase. ChemBioChem 2020, 21, 1671–1709. [Google Scholar] [CrossRef]
  19. Kästner, J.; Blöchl, P.E. Ammonia Production at the FeMo Cofactor of Nitrogenase: Results from Density Functional Theory. J. Am. Chem. Soc. 2007, 129, 2998–3006. [Google Scholar] [CrossRef]
  20. Varley, J.B.; Wang, Y.; Chan, K.; Studt, F.; Nørskov, J.K. Mechanistic Insights into Nitrogen Fixation by Nitrogenase Enzymes. Phys. Chem. Chem. Phys. 2015, 17, 29541–29547. [Google Scholar] [CrossRef]
  21. Raugei, S.; Seefeldt, L.C.; Hoffman, B.M. Critical Computational Analysis Illuminates the Reductive-Elimination Mechanism That Activates Nitrogenase for N2 Reduction. Proc. Natl. Acad. Sci. USA 2018, 115, 10521–10530. [Google Scholar] [CrossRef]
  22. Siegbahn, P.E.M. The Mechanism for Nitrogenase Including All Steps. Phys. Chem. Chem. Phys. 2019, 21, 15747–15759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Thorneley, R.N.F.; Lowe, D.J. Kinetics and Mechanism of the Nitrogenase Enzyme System. In Molybdenum Enzymes; Spiro, T.G., Ed.; Wiley: New York, NY, USA, 1985; pp. 221–284. [Google Scholar]
  24. Rutledge, H.L.; Tezcan, F.A. Electron Transfer in Nitrogenase. Chem. Rev. 2020, 120, 5158–5193. [Google Scholar] [CrossRef] [PubMed]
  25. Igarashi, R.Y.; Seefeldt, L.C. Nitrogen Fixation: The Mechanism of the Mo-Dependent Nitrogenase. Crit. Rev. Biochem. Mol. Biol. 2003, 38, 351–384. [Google Scholar] [CrossRef] [PubMed]
  26. O’Donnell, M.J.; Smith, B.E. Electron-Paramagnetic-Resonance Studies on the Redox Properties of the Molybdenum-Iron Protein of Nitrogenase between +50 and −450 MV. Biochem. J. 1978, 173, 831–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Watt, G.D.; Burns, A.; Lough, S.; Tennent, D.L. Redox and Spectroscopic Properties of Oxidized MoFe Protein from Azotobacter Vinelandii. Biochemistry 1980, 19, 4926–4932. [Google Scholar] [CrossRef]
  28. Watt, G.D.; Bulen, W.A. Redox Potentials of FeMo Cofactor in Nitrogenase. In Proceedings of the First International Symposium on Nitrogen Fixation; Newton, W.E., Nymans, C.J., Eds.; Washington State University Press: Washington, DC, USA, 1976; Volume 1, p. 123. [Google Scholar]
  29. Lough, S.; Burns, A.; Watt, G.D. Redox Reactions of the Nitrogenase Complex from Azotobacter Vinelandii. Biochemistry 1983, 22, 4062–4066. [Google Scholar] [CrossRef]
  30. Torres, R.A.; Lovell, T.; Noodleman, L.; Case, D.A. Density Functional and Reduction Potential Calculations of Fe4S4 Clusters. J. Am. Chem. Soc. 2003, 125, 1923–1936. [Google Scholar] [CrossRef]
  31. Ullmann, M.G.; Noodleman, L.; Case, D.A. Density Functional Calculation of PKa Values and Redox Potentials in the Bovine Rieske Iron-Sulfur Protein. JBIC J. Biol. Inorg. Chem. 2002, 7, 632–639. [Google Scholar] [CrossRef]
  32. Mouesca, J.-M.; Chen, J.L.; Noodleman, L.; Bashford, D.; Case, D.A. Density Functional/Poisson-Boltzmann Calculations of Redox Potentials for Iron-Sulfur Clusters. J. Am. Chem. Soc. 1994, 116, 11898–11914. [Google Scholar] [CrossRef]
  33. Noodleman, L.; Lovell, T.; Liu, T.; Himo, F.; Torres, R.A. Insights into Properties and Energetics of Iron–Sulfur Proteins from Simple Clusters to Nitrogenase. Curr. Opin. Chem. Biol. 2002, 6, 259–273. [Google Scholar] [CrossRef]
  34. Noodleman, L.; Han, W.-G. Structure, Redox, PKa, Spin. A Golden Tetrad for Understanding Metalloenzyme Energetics and Reaction Pathways. J. Biol. Inorg. Chem. 2006, 11, 674–694. [Google Scholar] [CrossRef]
  35. Kim, C.-H.; Newton, W.E.; Dean, D.R. Role of the MoFe Protein.Alpha.-Subunit Histidine-195 Residue in FeMo-Cofactor Binding and Nitrogenase Catalysis. Biochemistry 1995, 34, 2798–2808. [Google Scholar] [CrossRef] [PubMed]
  36. Perrin, B.S., Jr.; Niu, S.; Ichiye, T. Calculating Standard Reduction Potentials of [4Fe–4S] Proteins. J. Comput. Chem. 2013, 34, 576–582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Perrin, B.S., Jr.; Miller, B.T.; Schalk, V.; Woodcock, H.L.; Brooks, B.R.; Ichiye, T. Web-Based Computational Chemistry Education with CHARMMing III: Reduction Potentials of Electron Transfer Proteins. PLOS Comput. Biol. 2014, 10, e1003739. [Google Scholar] [CrossRef] [Green Version]
  38. Stephens, P.J.; Jollie, D.R.; Warshel, A. Protein Control of Redox Potentials of Iron−Sulfur Proteins. Chem. Rev. 1996, 96, 2491–2514. [Google Scholar] [CrossRef]
  39. Lovell, T.; Liu, T.; Case, D.A.; Noodleman, L. Structural, Spectroscopic, and Redox Consequences of a Central Ligand in the FeMoco of Nitrogenase: A Density Functional Theoretical Study. J. Am. Chem. Soc. 2003, 125, 8377–8383. [Google Scholar] [CrossRef]
  40. Cheng, J.; Liu, X.; VandeVondele, J.; Sulpizi, M.; Sprik, M. Redox Potentials and Acidity Constants from Density Functional Theory Based Molecular Dynamics. Acc. Chem. Res. 2014, 47, 3522–3529. [Google Scholar] [CrossRef] [Green Version]
  41. Jafari, S.; Tavares Santos, Y.A.; Bergmann, J.; Irani, M.; Ryde, U. Benchmark Study of Redox Potentials Calculation for Iron–Sulfur Clusters in Proteins. Inorg. Chem. 2022, 61, 5991–6007. [Google Scholar] [CrossRef]
  42. Spatzal, T.; Schlesier, J.; Burger, E.-M.; Sippel, D.; Zhang, L.; Andrade, S.L.A.; Rees, D.C.; Einsle, O. Nitrogenase FeMoco Investigated by Spatially Resolved Anomalous Dispersion Refinement. Nat. Commun. 2016, 7, 10902. [Google Scholar] [CrossRef] [Green Version]
  43. Doan, P.E.; Telser, J.; Barney, B.M.; Igarashi, R.Y.; Dean, D.R.; Seefeldt, L.C.; Hoffman, B.M. 57Fe ENDOR Spectroscopy and ‘Electron Inventory’ Analysis of the Nitrogenase E4 Intermediate Suggest the Metal-Ion Core of FeMo-Cofactor Cycles Through Only One Redox Couple. J. Am. Chem. Soc. 2011, 133, 17329–17340. [Google Scholar] [CrossRef]
  44. Yang, Z.-Y.; Danyal, K.; Seefeldt, L.C. Nitrogen Fixation. Methods Mol. Biol. 2011, 766, 267–292. [Google Scholar]
  45. Pierik, A.J.; Wassink, H.; Haaker, H.; Hagen, W.R. Redox Properties and EPR Spectroscopy of the P Clusters of Azotobacter Vinelandii MoFe Protein. Eur. J. Biochem. 1993, 212, 51–61. [Google Scholar] [CrossRef] [PubMed]
  46. Rupnik, K.; Hu, Y.; Lee, C.C.; Wiig, J.A.; Ribbe, M.W.; Hales, B.J. P+ State of Nitrogenase P-Cluster Exhibits Electronic Structure of a [Fe4S4]+ Cluster. J. Am. Chem. Soc. 2012, 134, 13749–13754. [Google Scholar] [CrossRef] [Green Version]
  47. Katz, F.E.H.; Owens, C.P.; Tezcan, F.A. Electron Transfer Reactions in Biological Nitrogen Fixation. Isr. J. Chem. 2016, 56, 682–692. [Google Scholar] [CrossRef]
  48. Owens, C.P.; Katz, F.E.H.; Carter, C.H.; Oswald, V.F.; Tezcan, F.A. Tyrosine-Coordinated P-Cluster in G. Diazotrophicus Nitrogenase: Evidence for the Importance of O-Based Ligands in Conformationally Gated Electron Transfer. J. Am. Chem. Soc. 2016, 138, 10124–10127. [Google Scholar] [CrossRef] [PubMed]
  49. Seefeldt, L.C.; Hoffman, B.M.; Dean, D.R. Electron Transfer in Nitrogenase Catalysis. Curr. Opin. Chem. Biol. 2012, 16, 19–25. [Google Scholar] [CrossRef] [Green Version]
  50. Peters, J.W.; Stowell, M.H.B.; Soltis, S.M.; Finnegan, M.G.; Johnson, M.K.; Rees, D.C. Redox-Dependent Structural Changes in the Nitrogenase P-Cluster. Biochemistry 1997, 36, 1181–1187. [Google Scholar] [CrossRef]
  51. Rupnik, K.; Lee, C.C.; Wiig, J.A.; Hu, Y.; Ribbe, M.W.; Hales, B.J. Nonenzymatic Synthesis of the P-Cluster in the Nitrogenase MoFe Protein: Evidence of the Involvement of All-Ferrous [Fe4S4]0 Intermediates. Biochemistry 2014, 53, 1108–1116. [Google Scholar] [CrossRef]
  52. Lanzilotta, W.N.; Christiansen, J.; Dean, D.R.; Seefeldt, L.C. Evidence for Coupled Electron and Proton Transfer in the [8Fe-7S] Cluster of Nitrogenase. Biochemistry 1998, 37, 11376–11384. [Google Scholar] [CrossRef]
  53. Chan, J.M.; Christiansen, J.; Dean, D.R.; Seefeldt, L.C. Spectroscopic Evidence for Changes in the Redox State of the Nitrogenase P-Cluster during Turnover. Biochemistry 1999, 38, 5779–5785. [Google Scholar] [CrossRef]
  54. Keable, S.M.; Zadvornyy, O.A.; Johnson, L.E.; Ginovska, B.; Rasmussen, A.J.; Danyal, K.; Eilers, B.J.; Prussia, G.A.; LeVan, A.X.; Raugei, S.; et al. Structural Characterization of the P1+ Intermediate State of the P-Cluster of Nitrogenase. J. Biol. Chem. 2018, 293, 9629–9635. [Google Scholar] [CrossRef] [PubMed]
  55. Cao, L.; Börner, M.C.; Bergmann, J.; Caldararu, O.; Ryde, U. Geometry and Electronic Structure of the P-Cluster in Nitrogenase Studied by Combined Quantum Mechanical and Molecular Mechanical Calculations and Quantum Refinement. Inorg. Chem. 2019, 58, 9672–9690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Cao, L.; Ryde, U. Quantum Refinement with Multiple Conformations: Application to the P-Cluster in Nitrogenase. Acta Crystallogr. Sect. D 2020, 76, 1145–1156. [Google Scholar] [CrossRef] [PubMed]
  57. Bjornsson, R.; Neese, F.; DeBeer, S. Revisiting the Mössbauer Isomer Shifts of the FeMoco Cluster of Nitrogenase and the Cofactor Charge. Inorg. Chem. 2017, 56, 1470–1477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Siegbahn, P.E.M. Model Calculations Suggest That the Central Carbon in the FeMo-Cofactor of Nitrogenase Becomes Protonated in the Process of Nitrogen Fixation. J. Am. Chem. Soc. 2016, 138, 10485–10495. [Google Scholar] [CrossRef]
  59. Cao, L.; Caldararu, O.; Ryde, U. Protonation and Reduction of the FeMo Cluster in Nitrogenase Studied by Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations. J. Chem. Theory Comput. 2018, 14, 6653–6678. [Google Scholar] [CrossRef] [Green Version]
  60. Cao, L.; Ryde, U. Influence of the Protein and DFT Method on the Broken-Symmetry and Spin States in Nitrogenase. Int. J. Quantum Chem. 2018, 118, e25627. [Google Scholar] [CrossRef]
  61. Van Stappen, C.; Thorhallsson, A.T.; Decamps, L.; Bjornsson, R.; DeBeer, S. Resolving the Structure of the E1 State of Mo Nitrogenase through Mo and Fe K-Edge EXAFS and QM/MM Calculations. Chem. Sci. 2019, 10, 9807–9821. [Google Scholar] [CrossRef] [Green Version]
  62. Van Stappen, C.; Davydov, R.; Yang, Z.-Y.; Fan, R.; Guo, Y.; Bill, E.; Seefeldt, L.C.; Hoffman, B.M.; DeBeer, S. Spectroscopic Description of the E1 State of Mo Nitrogenase Based on Mo and Fe X-Ray Absorption and Mössbauer Studies. Inorg. Chem. 2019, 58, 12365–12376. [Google Scholar] [CrossRef] [Green Version]
  63. Cao, L.; Ryde, U. What Is the Structure of the E4 Intermediate in Nitrogenase? J. Chem. Theory Comput. 2020, 16, 1936–1952. [Google Scholar] [CrossRef] [Green Version]
  64. Jiang, H.; Svensson, O.K.G.; Ryde, U. QM/MM Study of Partial Dissociation of S2B for the E2 Intermediate of Nitrogenase. Inorg. Chem. 2022, 61, 18067–18076. [Google Scholar] [CrossRef] [PubMed]
  65. Cao, L.; Ryde, U. Putative Reaction Mechanism of Nitrogenase after Dissociation of a Sulfide Ligand. J. Catal. 2020, 391, 247–259. [Google Scholar] [CrossRef]
  66. Jiang, H.; Ryde, U. Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of Any Sulfide Ligand. Chem.-Eur. J. 2022, 28, e202103933. [Google Scholar] [CrossRef] [PubMed]
  67. Jiang, H.; Svensson, O.K.G.; Cao, L.; Ryde, U. Proton Transfer Pathways in Nitrogenase with and without Dissociated S2B. Angew. Chemie Int. Ed. 2022, 61, e202208544. [Google Scholar] [CrossRef] [PubMed]
  68. Rohde, M.; Sippel, D.; Trncik, C.; Andrade, S.L.A.; Einsle, O. The Critical E4 State of Nitrogenase Catalysis. Biochemistry 2018, 57, 5497–5504. [Google Scholar] [CrossRef]
  69. Cao, L.; Ryde, U. Extremely Large Differences in DFT Energies for Nitrogenase Models. Phys. Chem. Chem. Phys. 2019, 21, 2480–2488. [Google Scholar] [CrossRef] [Green Version]
  70. Cao, L.; Caldararu, O.; Ryde, U. Protonation States of Homocitrate and Nearby Residues in Nitrogenase Studied by Computational Methods and Quantum Refinement. J. Phys. Chem. B 2017, 121, 8242–8262. [Google Scholar] [CrossRef] [Green Version]
  71. Barney, B.M.; McClead, J.; Lukoyanov, D.; Laryukhin, M.; Yang, T.C.; Dean, D.R.; Hoffman, B.M.; Seefeldt, L.C. Diazene (HN=NH) Is a Substrate for Nitrogenase: Insights into the Pathway of N2 Reduction. Biochemistry 2007, 46, 6784–6794. [Google Scholar] [CrossRef] [Green Version]
  72. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. Ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from Ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  73. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  74. Cao, L.; Ryde, U. N2H2 Binding to the Nitrogenase FeMo Cluster, Studied by QM/MM Methods. J. Biol. Inorg. Chem. 2020, 25, 521–540. [Google Scholar] [CrossRef] [PubMed]
  75. Hu, L.; Ryde, U. Comparison of Methods to Obtain Force-Field Parameters for Metal Sites. J. Chem. Theory Comput. 2011, 7, 2452–2463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Bayly, C.I.; Cieplak, P.; Cornell, W.D.; Kollman, P.A. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges: The RESP Model. J. Phys. Chem. 1993, 97, 10269–10280. [Google Scholar] [CrossRef]
  77. Furche, F.; Ahlrichs, R.; Hättig, C.; Klopper, W.; Sierka, M.; Weigend, F. Turbomole. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2014, 4, 91–100. [Google Scholar] [CrossRef]
  78. Tao, J.; Perdew, J.P.; Staroverov, V.N.; Scuseria, G.E. Climbing the Density Functional Ladder: Non-Empirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 146401. [Google Scholar] [CrossRef] [Green Version]
  79. Schäfer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian Basis Sets for Atoms Li to Kr. J. Chem. Phys. 1992, 97, 2571–2577. [Google Scholar] [CrossRef] [Green Version]
  80. Eichkorn, K.; Treutler, O.; Öhm, H.; Häser, M.; Ahlrichs, R. Auxiliary Basis Sets to Approximate Coulomb Potentials. Chem. Phys. Lett. 1995, 240, 283–289. [Google Scholar] [CrossRef]
  81. Eichkorn, K.; Weigend, F.; Treutler, O.; Ahlrichs, R. Auxiliary Basis Sets for Main Row Atoms and Transition Metals and Their Use to Approximate Coulomb Potentials. Theor. Chem. Acc. 1997, 97, 119–124. [Google Scholar] [CrossRef]
  82. Caldeweyher, E.; Ehlert, S.; Hansen, A.; Neugebauer, H.; Spicher, S.; Bannwarth, C.; Grimme, S. A Generally Applicable Atomic-Charge Dependent London Dispersion Correction. J. Chem. Phys. 2019, 150, 154122. [Google Scholar] [CrossRef]
  83. Hu, L.; Söderhjelm, P.; Ryde, U. Accurate Reaction Energies in Proteins Obtained by Combining QM/MM and Large QM Calculations. J. Chem. Theory Comput. 2013, 9, 640–649. [Google Scholar] [CrossRef] [Green Version]
  84. Klamt, A.; Schüürmann, G. Cosmo—A New Approach To Dielectric Screening in Solvents With Explicit Expressions for the Screening Energy and Its Gradient. J. Chem. Soc.-Perkin Trans. 2 1993, 1993, 799–805. [Google Scholar] [CrossRef]
  85. Schäfer, A.; Klamt, A.; Sattel, D.; Lohrenz, J.C.W.; Eckert, F. COSMO Implementation in TURBOMOLE: Extension of an Efficient Quantum Chemical Code towards Liquid Systems. Phys. Chem. Chem. Phys. 2000, 2, 2187–2193. [Google Scholar] [CrossRef] [Green Version]
  86. Klamt, A.; Jonas, V.; Bürger, T.; Lohrenz, J.C.W. Refinement and Parametrization of COSMO-RS. J. Phys. Chem. A 1998, 102, 5074–5085. [Google Scholar] [CrossRef]
  87. Sigfridsson, E.; Ryde, U. Comparison of Methods for Deriving Atomic Charges from the Electrostatic Potential and Moments. J. Comput. Chem. 1998, 19, 377–395. [Google Scholar] [CrossRef]
  88. Kelly, C.P.; Cramer, C.J.; Truhlar, D.G. Aqueous Solvation Free Energies of Ions and Ion-Water Clusters Based on an Accurate Value for the Absolute Aqueous Solvation Free Energy of the Proton. J. Phys. Chem. B 2006, 110, 16066–16081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Cheng, C.; Hayashi, S. Ab Initio Evaluation of the Redox Potential of Cytochrome C. J. Chem. Theory Comput. 2021, 17, 1194–1207. [Google Scholar] [CrossRef]
  90. Lovell, T.; Li, J.; Liu, T.; Case, D.A.; Noodleman, L. FeMo Cofactor of Nitrogenase: A Density Functional Study of States MN, MOX, MR, and MI. J. Am. Chem. Soc. 2001, 123, 12392–12410. [Google Scholar] [CrossRef]
  91. Greco, C.; Fantucci, P.; Ryde, U.; de Gioia, L. Fast Generation of Broken-Symmetry States in a Large System Including Multiple Iron–Sulfur Assemblies: Investigation of QM/MM Energies, Clusters Charges, and Spin Populations. Int. J. Quantum Chem. 2011, 111, 3949–3960. [Google Scholar] [CrossRef] [Green Version]
  92. Szilagyi, R.K.; Winslow, M.A. On the Accuracy of Density Functional Theory for Iron-Sulfur Clusters. J. Comput. Chem. 2006, 27, 1385–1397. [Google Scholar] [CrossRef]
  93. Ryde, U. The Coordination of the Catalytic Zinc in Alcohol Dehydrogenase Studied by Combined Quantum-Chemical and Molecular Mechanics Calculations. J. Comput. Aided. Mol. Des. 1996, 10, 153–164. [Google Scholar] [CrossRef] [Green Version]
  94. Ryde, U.; Olsson, M.H.M. Structure, Strain, and Reorganization Energy of Blue Copper Models in the Protein. Int. J. Quantum Chem. 2001, 81, 335–347. [Google Scholar] [CrossRef]
  95. Reuter, N.; Dejaegere, A.; Maigret, B.; Karplus, M. Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches. J. Phys. Chem. A 2000, 104, 1720–1735. [Google Scholar] [CrossRef]
  96. Hu, L.; Söderhjelm, P.; Ryde, U. On the Convergence of QM/MM Energies. J. Chem. Theory Comput. 2011, 7, 761–777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Cao, L.; Ryde, U. On the Difference between Additive and Subtractive QM/MM Calculations. Front. Chem. 2018, 6, 89. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The cores of (a) the P-cluster (PN state) and (b) the FeMo cluster (E0 state) with atom names indicated.
Figure 1. The cores of (a) the P-cluster (PN state) and (b) the FeMo cluster (E0 state) with atom names indicated.
Molecules 28 00065 g001
Table 1. The various states studied for the P-cluster. The table lists an abbreviation (Abb.) for each state, the protonation status for Cys-88 and Ser-188D (P = protonated, D = deprotonated), the spin state (S), the BS state (specifying Fe ions with minority spin) and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
Table 1. The various states studied for the P-cluster. The table lists an abbreviation (Abb.) for each state, the protonation status for Cys-88 and Ser-188D (P = protonated, D = deprotonated), the spin state (S), the BS state (specifying Fe ions with minority spin) and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
StateAbb.Cys-88Ser-188SBSch
reducedPNH2PP01247−4
one-electron oxidisedP1+H2PP½1247−3
P1+HPD½1247−4
P1+DD½1247−5
two-electron oxidisedP2+DD4358−4
three-electron oxidisedP3+DD7/2358−3
Table 2. Redox potentials of five redox couples for the P-cluster in V. E calc o is the raw redox potentials from the COSMO calculations, E corr o are calculated redox potentials corrected by the mean signed error (MSE = −0.62 V) from our previous study [41] (i.e., E corr o = E calc o + 0.62 V), E exp o is the experimental redox potentials [5,24,25,45] and E err o is the error for the various calculations. Results with E err o < 0.44 V are marked in bold face.
Table 2. Redox potentials of five redox couples for the P-cluster in V. E calc o is the raw redox potentials from the COSMO calculations, E corr o are calculated redox potentials corrected by the mean signed error (MSE = −0.62 V) from our previous study [41] (i.e., E corr o = E calc o + 0.62 V), E exp o is the experimental redox potentials [5,24,25,45] and E err o is the error for the various calculations. Results with E err o < 0.44 V are marked in bold face.
Redox Couple       E calc o E corr o E exp o   E err o
PNH2 → P1+H2−0.620.00−0.3090.31
PNH2 → P1+H−0.600.02−0.3090.33
P1+H → P2+−0.480.14<−0.224 a<0.36
P1+ → P2+−1.69−1.07<−0.348 a<0.72
P2+ → P3+−0.580.040.09−0.05
a The measured redox potential for the P2+/P1+ couple is pH dependent, decreasing from −0.224 V at pH 6.0 to −0.348 V at pH 8.5 (cf. the text) [52].
Table 3. States studied for the resting, one-electron reduced or one-electron oxidised states of the FeMo cluster. The table lists the protonation status of the S2B atom (P = protonated, D = deprotonated), the spin state (S), the BS state (specifying Fe ions with minority spin), and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
Table 3. States studied for the resting, one-electron reduced or one-electron oxidised states of the FeMo cluster. The table lists the protonation status of the S2B atom (P = protonated, D = deprotonated), the spin state (S), the BS state (specifying Fe ions with minority spin), and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
StateS2B SBSch
R−1D12350
R0D3/2235−1
R1D2235−2
R1HP2235−1
E−1D1235−2
E0D3/2235−3
E1D2235−4
E1HP2235−3
A−1D2235−4
A0D3/2235−5
A1D2235−6
A1HP2235−5
Table 4. Redox potentials for various redox couples of the FeMo cluster in V. E calc o is the raw redox potentials, E corr o = E calc o + 0.62 V (the MSE in our previous study [41]), E exp o is the experimental redox potentials (for the E0/E1 couple, we used −0.47 V, which is in the middle of the range of reported values) [5,24,25,27,28,29] and E err o is the error for the various calculations. Results with   E err o < 0.44 V are marked in bold face.
Table 4. Redox potentials for various redox couples of the FeMo cluster in V. E calc o is the raw redox potentials, E corr o = E calc o + 0.62 V (the MSE in our previous study [41]), E exp o is the experimental redox potentials (for the E0/E1 couple, we used −0.47 V, which is in the middle of the range of reported values) [5,24,25,27,28,29] and E err o is the error for the various calculations. Results with   E err o < 0.44 V are marked in bold face.
Transition E calc o E corr o E exp o   E err o
R−1 → R01.311.93−0.0421.98
R0 → R10.421.04(−0.47)1.51
R0 → R1H−0.66−0.04−0.470.43
E−1 → E0−0.490.13−0.0420.17
E0 → E1−1.69−1.07(−0.47)−0.60
E0 → E1H−1.28−0.66−0.47−0.19
A−1 → A0−6.88−6.29−0.042−6.21
A0 → A1−3.46−2.84(−0.47)−2.37
A0 → A1H−1.57−0.95−0.47−0.48
Table 5. Structures studied for the E0–E8 states of the FeMo cluster. The table lists the protonated atoms or the added substrate, the spin state (S), the BS state (specifying the Fe ions with minority spin), and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
Table 5. Structures studied for the E0–E8 states of the FeMo cluster. The table lists the protonated atoms or the added substrate, the spin state (S), the BS state (specifying the Fe ions with minority spin), and the net charge of the QM systems (ch; same for both for QM/MM and for the redox calculations).
StateProtonated Atoms/SubstrateSBSch
E03/2235−3
E1HS2B(3)2235−3
E2H2S2B(3), Fe2/6(5)3/2247−3
E2H2S2B (dissoc from Fe2), Fe2/6(5)3/2247−3
E2H2S2B(3), S5A(3)3/2247−3
E3H3S2B(3), Fe2/6(5), Fe3/7(2)114−3
E3H3S2B(3), Fe2/6(5), Fe5114−3
E4H4S2B(3), Fe2/6(3), Fe3/7(2), S5A(3)½14−3
E4H4S2B(3), Fe2/6(5), Fe3/7(2), S5A(3)½14−3
E4H4S2B(5), Fe2/6(5), Fe3/7(2), S5A(2)½14−3
E4H4″′S2B(3), Fe2/6(3), Fe6, Fe5½14−3
With S2B still bound
E4N2H2HNNH2 (proton from HCA)½147−3
E5N2H3H2NNH2 (proton from HCA)1147−3
E6N2H4NH2 + NH3 (proton from HCA)½147−3
E6NHNH2 (proton from HCA)½147−3
E7NH2NH3 (proton from HCA)1147−3
E8NH3NH3½147−3
With S2B dissociated
E4N2H2NNH2½147−1
E5N2H3HNNH21147−1
E6N2H4H2NNH2½147−1
E7N2H5NH2 + NH30147−1
E7NH2NH21147−1
E8NH3NH3½147−1
Table 6. Calculated redox potentials for the E0–E8 states of the FeMo cluster. The last column ( Δ E calc o ) reports the difference in the calculated redox potential compared to that of the E0 → E1H transition. Redox potentials for the most favourable structures of the first four transitions are shown in bold face.
Table 6. Calculated redox potentials for the E0–E8 states of the FeMo cluster. The last column ( Δ E calc o ) reports the difference in the calculated redox potential compared to that of the E0 → E1H transition. Redox potentials for the most favourable structures of the first four transitions are shown in bold face.
Transition E calc o Δ E calc o
E0 → E1H−1.280.00
E1H → E2H2−1.200.08
E1H → E2H2−1.45−0.17
E1H → E2H2−1.290.00
E2H2 → E3H3−1.47−0.19
E2H2 → E3H3−1.81−0.53
E3H3 → E4H4−0.870.41
E3H3 → E4H4−1.100.18
E3H3 → E4H4−1.34−0.06
E3H3 → E4H4″′−1.49−0.21
With S2B
E4N2H2 → E5N2H3−0.151.13
E5N2H3 → E6N2H4−0.870.41
E6NH → E7NH20.742.02
E7NH2 → E8NH3−0.710.57
Without S2B
E4N2H2′ → E5N2H3−0.310.97
E5N2H3′ → E6N2H4−1.070.20
E6N2H4′ → E7N2H51.372.65
E7NH2′ → E8NH3−1.090.19
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiang, H.; Svensson, O.K.G.; Ryde, U. Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase. Molecules 2023, 28, 65. https://doi.org/10.3390/molecules28010065

AMA Style

Jiang H, Svensson OKG, Ryde U. Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase. Molecules. 2023; 28(1):65. https://doi.org/10.3390/molecules28010065

Chicago/Turabian Style

Jiang, Hao, Oskar K. G. Svensson, and Ulf Ryde. 2023. "Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase" Molecules 28, no. 1: 65. https://doi.org/10.3390/molecules28010065

Article Metrics

Back to TopTop