Comprehensive structural, infrared spectroscopic and kinetic investigations of the roles of the active-site arginine in bidirectional hydrogen activation by the [NiFe]-hydrogenase 'Hyd-2' from Escherichia coli

The active site of [NiFe]-hydrogenases contains a strictly-conserved pendant arginine, the guanidine head group of which is suspended immediately above the Ni and Fe atoms. Replacement of this arginine (R479) in hydrogenase-2 from E. coli results in an enzyme that is isolated with a very tightly-bound diatomic ligand attached end-on to the Ni and stabilised by hydrogen bonding to the Nζ atom of the pendant lysine and one of the three additional water molecules located in the active site of the variant. The diatomic ligand is bound under oxidising conditions and is removed only after a prolonged period of reduction with H2 and reduced methyl viologen. Once freed of the diatomic ligand, the R479K variant catalyses both H2 oxidation and evolution but with greatly decreased rates compared to the native enzyme. Key kinetic characteristics are revealed by protein film electrochemistry: most importantly, a very low activation energy for H2 oxidation that is not linked to an increased H/D isotope effect. Native electrocatalytic reversibility is retained. The results show that the sluggish kinetics observed for the lysine variant arise most obviously because the advantage of a more favourable low-energy pathway is massively offset by an extremely unfavourable activation entropy. Extensive efforts to establish the identity of the diatomic ligand, the tight binding of which is an unexpected further consequence of replacing the pendant arginine, prove inconclusive.

The E3 state of FeMoco: one hydride, two hydrides or dihydrogen?

Hydrides are present in the reduced states of the iron-molybdenum cofactor (FeMoco) of Mo nitrogenase and are believed to play a key mechanistic role in the dinitrogen reduction reaction catalyzed by the enzyme. Two hydrides are present in the E4 state according to 1H ENDOR and there is likely a single hydride in the E2 redox state. The 2-hydride E4 state has been experimentally observed to bind N2 and it has been speculated that E3 may bind N2 as well. However, the E3 state has not been directly observed and very little is known about its molecular and electronic structure or reactivity. In recent computational studies, we have explored the energy surfaces of the E2 and E4 by QM/MM modelling, and found that the most stable hydride isomers contain bridging or partially bridging hydrides with an open protonated belt sulfide-bridge. In this work we systematically explore the energy surface of the E3 redox state, comparing single hydride and two-hydride isomers with varying coordination and bridging vs. terminal sulfhydryl groups. We also include a model featuring a triply protonated carbide. The results are only mildly dependent on the QM-region size. The three most stable E3 isomers at the r2SCAN level of theory have in common: an open belt sulfide-bridge (terminal sulfhydryl group on Fe6) and either 2 bridging hydrides (between Fe2 and Fe6), 1 bridging-1-terminal hydride (around Fe2 and Fe6) or a dihydrogen ligand bound at the Fe2 site. Analyzing the functional dependency of the results, we find that functionals previously found to predict accurate structures of spin-coupled Fe/Mo dimers and FeMoco (TPSSh, B97-D3, r2SCAN, and B3LYP*) are in generally good agreement about the stability of these 3 E3 isomers. However, B3LYP*, similar to its parent B3LYP method, predicts a triply protonated carbide isomer as the most stable isomer, an unlikely scenario in view of the lack of experimental evidence for carbide protonation occurring in reduced FeMoco states. Distinguishing further between the 3 hydride isomers is difficult and this flexible coordination nature of hydrides suggests that multiple hydride isomers could be present during experimental conditions. N2 binding was explored and resulted in geometries with 2 bridging hydrides and N2 bound to either Fe2 or Fe6 with a local low-spin state on the Fe. N2 binding is predicted to be mildly endothermic, similar to the E2 state, and it seems unlikely that the E3 state is capable of binding N2.

Light-Driven Hydrogen Evolution Reaction Catalyzed by a Molybdenum-Copper Artificial Hydrogenase

Orange protein (Orp) is a small bacterial metalloprotein of unknown function that harbors a unique molybdenum/copper (Mo/Cu) heterometallic cluster, [S₂MoS₂CuS₂MoS₂]^3-. In this paper, the performance of Orp as a catalyst for the photocatalytic reduction of protons into H₂ has been investigated under visible light irradiation. We report the complete biochemical and spectroscopic characterization of holo-Orp containing the [S₂MoS₂CuS₂MoS₂]^3- cluster, with docking and molecular dynamics simulations suggesting a positively charged Arg, Lys-containing pocket as the binding site. Holo-Orp exhibits excellent photocatalytic activity, in the presence of ascorbate as the sacrificial electron donor and [Ru(bpy)₃]Cl₂ as the photosensitizer, for hydrogen evolution with a maximum turnover number of 890 after 4 h irradiation. Density functional theory (DFT) calculations were used to propose a consistent reaction mechanism in which the terminal sulfur atoms are playing a key role in promoting H₂ formation. A series of dinuclear [S₂MS₂M'S₂MS₂]^(4n)- clusters, with M = Mo^VI, W^VI and M'⁽ⁿ⁺⁾ = Cu^I, Fe^I, Ni^I, Co^I, Zn^II, Cd^II were assembled in Orp, leading to different M/M'-Orp versions which are shown to display catalytic activity, with the Mo/Fe-Orp catalyst giving a remarkable turnover number (TON) of 1150 after 2.5 h reaction and an initial turnover frequency (TOF°) of 800 h^-1 establishing a record among previously reported artificial hydrogenases.

Understanding the Electronic Structure Basis for N2 Binding to FeMoco: A Systematic Quantum Mechanics/Molecular Mechanics Investigation

The FeMo cofactor (FeMoco) of Mo nitrogenase is responsible for reducing dinitrogen to ammonia, but it requires the addition of 3-4 e^-/H⁺ pairs before N₂ even binds. A binding site at the Fe2/Fe3/Fe6/Fe7 face of the cofactor has long been suggested based on mutation studies, with Fe2 or Fe6 nowadays being primarily discussed as possibilities. However, the nature of N₂ binding to the cofactor is enigmatic as the metal ions are coordinatively saturated in the resting state with no obvious binding site. Furthermore, the cofactor consists of high-spin Fe(II)/Fe(III) ions (antiferromagnetically coupled but also mixed-valence delocalized), which are not known to bind N₂. This suggests that an Fe binding site with a different molecular and electronic structure than the resting state must be responsible for the experimentally known N₂ binding in the E₄ state of FeMoco. We have systematically studied N₂ binding to Fe2 and Fe6 sites of FeMoco at the broken-symmetry QM/MM level as a function of the redox-, spin-, and protonation state of the cofactor. The local and global electronic structure changes to the cofactor taking place during redox events, protonation, Fe-S cleavage, hydride formation, and N₂ coordination are systematically analyzed. Localized orbital and quasi-restricted orbital analysis via diamagnetic substitution is used to get insights into the local single Fe ion electronic structure in various states of FeMoco. A few factors emerge as essential to N₂ binding in the calculations: spin-pairing of d_xz and d_yz orbitals of the N₂-binding Fe ion, a coordination change at the N₂-binding Fe ion aided by a hemilabile protonated sulfur, and finally hydride ligation. The results show that N₂ binding to E₀, E₁, and E₂ models is generally unfavorable, likely due to the high-energy cost of stabilizing the necessary spin-paired electronic structure of the N₂-binding Fe ion in a ligand environment that clearly favors high-spin states. The results for models of E₄, however, suggest a feasible model for why N₂ binding occurs experimentally in the E₄ state. E₄ models with two bridging hydrides between Fe2 and Fe6 show much more favorable N₂ binding than other models. When two hydrides coordinate to the same Fe ion, an increased ligand-field splitting due to octahedral coordination at either Fe2 or Fe6 is found. This altered ligand field makes it easier for the Fe ion to acquire a spin-paired electronic structure with doubly occupied d_xz and d_yz orbitals that backbond to a terminal N₂ ligand. Within this model for N₂ binding, both Fe2 and Fe6 emerge as possible binding site scenarios. Due to steric effects involving the His195 residue, affecting both the N₂ ligand and the terminal SH^- group, distinguishing between Fe2 and Fe6 is difficult; furthermore, the binding depends on the protonation state of His195.

Binding of exogenous cyanide reveals new active-site states in [FeFe] hydrogenases

[FeFe] hydrogenases are highly efficient metalloenyzmes for hydrogen conversion. Their active site cofactor (the H-cluster) is composed of a canonical [4Fe-4S] cluster ([4Fe-4S]_H) linked to a unique organometallic di-iron subcluster ([2Fe]_H). In [2Fe]_H the two Fe ions are coordinated by a bridging 2-azapropane-1,3-dithiolate (ADT) ligand, three CO and two CN^- ligands, leaving an open coordination site on one Fe where substrates (H₂ and H⁺) as well as inhibitors (e.g. O₂, CO, H₂S) may bind. Here, we investigate two new active site states that accumulate in [FeFe] hydrogenase variants where the cysteine (Cys) in the proton transfer pathway is mutated to alanine (Ala). Our experimental data, including atomic resolution crystal structures and supported by calculations, suggest that in these two states a third CN^- ligand is bound to the apical position of [2Fe]_H. These states can be generated both by "cannibalization" of CN^- from damaged [2Fe]_H subclusters as well as by addition of exogenous CN^-. This is the first detailed spectroscopic and computational characterisation of the interaction of exogenous CN^- with [FeFe] hydrogenases. Similar CN^--bound states can also be generated in wild-type hydrogenases, but do not form as readily as with the Cys to Ala variants. These results highlight how the interaction between the first amino acid in the proton transfer pathway and the active site tunes ligand binding to the open coordination site and affects the electronic structure of the H-cluster.

Ionization energies of metallocenes: a coupled cluster study of cobaltocene

Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T₁) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T₁) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.

Structural Correlations of Nitrogenase Active Sites using Nuclear Resonance Vibrational Spectroscopy and QM/MM Calculations

The biological conversion of N2 to NH3 is accomplished by the nitrogenase family, which is collectively comprised of three closely related but unique metalloenzymes. In the present study, we have...

Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco

The open-shell electronic structure of iron–sulfur clusters presents considerable challenges to quantum chemistry, with the complex iron–molybdenum cofactor (FeMoco) of nitrogenase representing perhaps the ultimate challenge for either wavefunction or density functional theory. While broken-symmetry density functional theory has seen some success in describing the electronic structure of such cofactors, there is a large exchange–correlation functional dependence in calculations that is not fully understood. In this work, we present a geometric benchmarking test set, FeMoD11, of synthetic spin-coupled Fe–Fe and Mo–Fe dimers, with relevance to the molecular and electronic structure of the Mo-nitrogenase FeMo cofactor. The reference data consists of high-resolution crystal structures of metal dimer compounds in different oxidation states. Multiple density functionals are tested on their ability to reproduce the local geometry, specifically the Fe–Fe/Mo–Fe distance, for both antiferromagnetically coupled and ferromagnetically coupled dimers via the broken-symmetry approach. The metal–metal distance is revealed not only to be highly sensitive to the amount of exact exchange in the functional but also to the specific exchange and correlation functionals. For the antiferromagnetically coupled dimers, the calculated metal–metal distance correlates well with the covalency of the bridging metal–ligand bonds, as revealed via the corresponding orbital analysis, Hirshfeld S/Fe charges, and Fe–S Mayer bond order. Superexchange via bridging ligands is expected to be the dominant interaction in these dimers, and our results suggest that functionals that predict accurate Fe–Fe and Mo–Fe distances describe the overall metal–ligand covalency more accurately and in turn the superexchange of these systems. The best performing density functionals of the 16 tested for the FeMoD11 test set are revealed to be either the nonhybrid functionals r²SCAN and B97-D3 or hybrid functionals with 10–15% exact exchange: TPSSh and B3LYP*. These same four functionals are furthermore found to reproduce the high-resolution X-ray structure of FeMoco well according to quantum mechanics/molecular mechanics (QM/MM) calculations. Almost all nonhybrid functionals systematically underestimate Fe–Fe and Mo–Fe distances (with r²SCAN and B97-D3 being the sole exceptions), while hybrid functionals with >15% exact exchange (including range-separated hybrid functionals) overestimate them. The results overall suggest r²SCAN, B97-D3, TPSSh, and B3LYP* as accurate density functionals for describing the electronic structure of iron–sulfur clusters in general, including the complex FeMoco cluster of nitrogenase.

Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

The E2 state of FeMoco: Hydride Formation versus Fe Reduction and a Mechanism for H2 Evolution

The iron‐molybdenum cofactor (FeMoco) is responsible for dinitrogen reduction in Mo nitrogenase. Unlike the resting state, E₀, reduced states of FeMoco are much less well characterized. The E₂ state has been proposed to contain a hydride but direct spectroscopic evidence is still lacking. The E₂ state can, however, relax back the E₀ state via a H₂ side‐reaction, implying a hydride intermediate prior to H₂ formation. This E₂→E₀ pathway is one of the primary mechanisms for H₂ formation under low‐electron flux conditions. In this study we present an exploration of the energy surface of the E₂ state. Utilizing both cluster‐continuum and QM/MM calculations, we explore various classes of E₂ models: including terminal hydrides, bridging hydrides with a closed or open sulfide‐bridge, as well as models without. Importantly, we find the hemilability of a protonated belt‐sulfide to strongly influence the stability of hydrides. Surprisingly, non‐hydride models are found to be almost equally favorable as hydride models. While the cluster‐continuum calculations suggest multiple possibilities, QM/MM suggests only two models as contenders for the E₂ state. These models feature either i) a bridging hydride between Fe₂ and Fe₆ and an open sulfide‐bridge with terminal SH on Fe₆ (E₂‐hyd) or ii) a double belt‐sulfide protonated, reduced cofactor without a hydride (E₂‐nonhyd). We suggest both models as contenders for the E₂ redox state and further calculate a mechanism for H₂ evolution. The changes in electronic structure of FeMoco during the proposed redox‐state cycle, E₀→E₁→E₂→E₀, are discussed.

Synthesis, Characterization, and Reaction Studies of Pd(II) Tripeptide Complexes

The aqueous synthesis of Pd(II) complexes with alkylated tripeptides led to the hydrolysis of the peptides at low pH values and mixtures of complexed peptides were formed. A non-aqueous synthetic route allowed the formation and isolation of single products and their characterization. Pd(II) complexes with α-Asp(OR)AlaGly(OR), β-Asp(OR)AlaGly(OR), and TrpAlaGly(OR) (R = H or alkyl) as tri and tetradentate chelates were characterized. The tridentate coordination mode was accompanied by a fourth monodentate ligand that was shown to participate in both ligand exchange reactions and a direct removal to form the tetradentate coordination mode. The tetradentate coordination revealed a rare a hemi labile carbonyl goup coordination mode to Pd(II). Reactivity with small molecules such as ethylene, acids, formate, and episulfide was investigated. Under acidic conditions and in the presence of ethylene; acetaldehyde was formed. The Pd(II) is a soft Lewis acid and thiophilic and the complexes abstract sulfur from episulfide at apparent modest catalytic rates. The complexes adopt a square planar geometry according to a spectroscopic analysis and DFT calculations that were employed to evaluate the most energetically favorable coordination geometry and compared with the observed infrared and NMR data.

Nudged Elastic Band Method for Molecular Reactions Using Energy-Weighted Springs Combined with Eigenvector Following

The climbing image nudged elastic band method (CI-NEB) is used to identify reaction coordinates and to find saddle points representing transition states of reactions. It can make efficient use of parallel computing as the calculations of the discretization points, the so-called images, can be carried out simultaneously. In typical implementations, the images are distributed evenly along the path by connecting adjacent images with equally stiff springs. However, for systems with a high degree of flexibility, this can lead to poor resolution near the saddle point. By making the spring constants increase with energy, the resolution near the saddle point is improved. To assess the performance of this energy-weighted CI-NEB method, calculations are carried out for a benchmark set of 121 molecular reactions. The performance of the method is analyzed with respect to the input parameters. Energy-weighted springs are found to greatly improve performance and result in successful location of the saddle points in less than a thousand energy and force evaluations on average (about a hundred per image) using the same set of parameter values for all of the reactions. Even better performance is obtained by stopping the calculation before full convergence and complete the saddle point search using an eigenvector following method starting from the location of the climbing image. This combination of methods, referred to as NEB-TS, turns out to be robust and highly efficient as it reduces the average number of energy and force evaluations down to a third, to 305. An efficient and flexible implementation of these methods has been made available in the ORCA software.

Localized and Delocalized States of a Diamine Cation: Resolution of a Controversy

Recent Rydberg spectroscopy measurements of a diamine molecule, N,N'-dimethylpiperazine (DMP), indicate the existence of a localized electronic state as well as a delocalized electronic state. This implies that the cation, DMP⁺, can similarly have its positive charge either localized on one of the N atoms or delocalized over both. This interpretation of the experiments has, however, been questioned based on coupled cluster calculations. In this article, results of high-level multireference configuration interaction calculations are presented where a localized state of DMP⁺ is indeed found to be present with an energy barrier separating it from the delocalized state. The energy difference between the two states is in excellent agreement with the experimental estimate. The results presented here, therefore, support the original interpretation of the experiments and illustrate a rare shortcoming of CCSD(T), the "gold standard" of quantum chemistry. These results have implications for the development of density functionals, as most functionals fail to produce the localized state.

Role of N-Oxide Moieties in Tuning Supramolecular Gel-State Properties

The role of specific interactions in the self-assembly process of low molecular weight gelators (LMWGs) was studied by altering the nonbonding interactions responsible for gel formation via structural modification of the gelator/nongelator. This was achieved by modifying pyridyl moieties of bis(pyridyl) urea-based hydrogelator (4–BPU) and the isomer (3–BPU) to pyridyl N–oxide compounds (L₁ and L₂, respectively). The modification of the functional groups resulted in the tuning of the gelation properties of the parent gelator, which induced/enhanced the gelation properties. The modified compounds displayed better mechanical and thermal stabilities and the introduction of the N–oxide moieties had a prominent effect on the morphologies of the gel network, which was evident from the scanning electron microscopy (SEM) images. The effect of various interactions due to the introduction of N–oxide moieties in the gel network formation was analyzed by comparing the solid-state interactions of the compounds using single crystal X-ray diffraction and computational studies, which were correlated with the enhanced gelation properties. This study shows the importance of specific nonbonding interactions and the spatial arrangement of the functional groups in the supramolecular gel network formation.

Caught in the Hinact : Crystal Structure and Spectroscopy Reveal a Sulfur Bound to the Active Site of an O2 -stable State of [FeFe] Hydrogenase

[FeFe] hydrogenases are the most active H2 converting catalysts in nature, but their extreme oxygen sensitivity limits their use in technological applications. The [FeFe] hydrogenases from sulfate reducing bacteria can be purified in an O2 -stable state called Hinact . To date, the structure and mechanism of formation of Hinact remain unknown. Our 1.65 Å crystal structure of this state reveals a sulfur ligand bound to the open coordination site. Furthermore, in-depth spectroscopic characterization by X-ray absorption spectroscopy (XAS), nuclear resonance vibrational spectroscopy (NRVS), resonance Raman (RR) spectroscopy and infrared (IR) spectroscopy, together with hybrid quantum mechanical and molecular mechanical (QM/MM) calculations, provide detailed chemical insight into the Hinact state and its mechanism of formation. This may facilitate the design of O2 -stable hydrogenases and molecular catalysts.

Scaffold diversity for enhanced activity of glycosylated inhibitors of fungal adhesion

Candida albicans is one of the most prevalent fungal pathogens involved in hospital acquired infections. It binds to glycans at the surface of epithelial cells and initiates infection. This process can be blocked by synthetic carbohydrates that mimic the structure of cell surface glycans. Herein we report the evaluation of a series of divalent glycosides featuring aromatic (benzene, squaramide) and bicyclic aliphatic (norbornene) scaffolds, with the latter being the first examples of their kind as small molecule anti-adhesion glycoconjugates. Galactosides 1 and 6, built on an aromatic core, were most efficient inhibitors of adhesion of C. albicans to buccal epithelial cells, displacing up to 36% and 48%, respectively, of yeast already attached to epithelial cells at 138 μM. Remarkably, cis-endo-norbornene 21 performed comparably to benzene-core derivatives. Conformational analysis reveals a preference for compounds 1 and 21 to adopt folded conformations. These results highlight the potential of norbornenes as a new class of aliphatic scaffolds for the synthesis of anti-adhesion compounds.

Quantum Mechanics/Molecular Mechanics Study of Resting-State Vanadium Nitrogenase: Molecular and Electronic Structure of the Iron-Vanadium Cofactor

The nitrogenase enzymes are responsible for all biological nitrogen reduction. How this is accomplished at the atomic level, however, has still not been established. The molybdenum-dependent nitrogenase has been extensively studied and is the most active catalyst for dinitrogen reduction of the nitrogenase enzymes. The vanadium-dependent form, on the other hand, displays different reactivity, being capable of CO and CO2 reduction to hydrocarbons. Only recently did a crystal structure of the VFe protein of vanadium nitrogenase become available, paving the way for detailed theoretical studies of the iron-vanadium cofactor (FeVco) within the protein matrix. The crystal structure revealed a bridging 4-atom ligand between two Fe atoms, proposed to be either a CO32- or NO3- ligand. Using a quantum mechanics/molecular mechanics model of the VFe protein, starting from the 1.35 Å crystal structure, we have systematically explored multiple computational models for FeVco, considering either a CO32- or NO3- ligand, three different redox states, and multiple broken-symmetry states. We find that only a [VFe7S8C(CO3)]2- model for FeVco reproduces the crystal structure of FeVco well, as seen in a comparison of the Fe-Fe and V-Fe distances in the computed models. Furthermore, a broken-symmetry solution with Fe2, Fe3, and Fe5 spin-down (BS7-235) is energetically preferred. The electronic structure of the [VFe7S8C(CO3)]2- BS7-235 model is compared to our [MoFe7S9C]- BS7-235 model of FeMoco via localized orbital analysis and is discussed in terms of local oxidation states and different degrees of delocalization. As previously found from Fe X-ray absorption spectroscopy studies, the Fe part of FeVco is reduced compared to FeMoco, and the calculations reveal Fe5 as locally ferrous. This suggests resting-state FeVco to be analogous to an unprotonated E1 state of FeMoco. Furthermore, V-Fe interactions in FeVco are not as strong compared to Mo-Fe interactions in FeMoco. These clear differences in the electronic structures of otherwise similar cofactors suggest an explanation for distinct differences in reactivity.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Insensitivity of Magnetic Coupling to Ligand Substitution in a Series of Tetraoxolene Radical-Bridged Fe2 Complexes

The elucidation of magnetostructural correlations between bridging ligand substitution and strength of magnetic coupling is essential to the development of high-temperature molecule-based magnetic materials. Toward this end, we report the series of tetraoxolene-bridged Fe^II₂ complexes [(Me₃TPyA)₂Fe₂(^RL)]ⁿ⁺ (Me₃TPyA = tris(6-methyl-2-pyridylmethyl)amine; n = 2: ^OMeLH₂ = 3,6-dimethoxy-2,5-dihydroxo-1,4-benzoquinone, ^ClLH₂ = 3,6-dichloro-2,5-dihydroxo-1,4-benzoquinone, Na₂[^NO₂L] = sodium 3,6-dinitro-2,5-dihydroxo-1,4-benzoquinone; n = 4: ^SMe₂L = 3,6-bis(dimethylsulfonium)-2,5-dihydroxo-1,4-benzoquinone diylide) and their one-electron-reduced analogues. Variable-temperature dc magnetic susceptibility data reveal the presence of weak ferromagnetic superexchange between Fe^II centers in the oxidized species, with exchange constants of J = +1.2(2) (R = OMe, Cl) and +0.3(1) (R = NO₂, SMe₂) cm^-1. In contrast, X-ray diffraction, cyclic voltammetry, and Mössbauer spectroscopy establish a ligand-centered radical in the reduced complexes. Magnetic measurements for the radical-bridged species reveal the presence of strong antiferromagnetic metal-radical coupling, with J = -57(10), -60(7), -58(6), and -65(8) cm^-1 for R = OMe, Cl, NO₂, and SMe₂, respectively. The minimal effects of substituents in the 3- and 6-positions of ^RL^x-• on the magnetic coupling strength is understood through electronic structure calculations, which show negligible spin density on the substituents and associated C atoms of the ring. Finally, the radical-bridged complexes are single-molecule magnets, with relaxation barriers of U_eff = 50(1), 41(1), 38(1), and 33(1) cm^-1 for R = OMe, Cl, NO₂, and SMe₂, respectively. Taken together, these results provide the first examination of how bridging ligand substitution influences magnetic coupling in semiquinoid-bridged compounds, and they establish design criteria for the synthesis of semiquinoid-based molecules and materials.

Dissociation of the FEBID precursor cis-Pt(CO)2Cl2 driven by low-energy electrons

In this study, we present experimental and theoretical results on dissociative electron attachment and dissociative ionisation for the potential FEBID precursor cis-Pt(CO)₂Cl₂. UHV surface studies have shown that high purity platinum deposits can be obtained from cis-Pt(CO)₂Cl₂. The efficiency and energetics of ligand removal through these processes are discussed and experimental appearance energies are compared to calculated thermochemical thresholds. The present results demonstrate the potential effectiveness of electron-induced reactions in the deposition of this FEBID precursor, and these are discussed in conjunction with surface science studies on this precursor and the design of new FEBID precursors.

A model for dinitrogen binding in the E4 state of nitrogenase

Molybdenum nitrogenase is one of the most intriguing metalloenzymes in nature, featuring an exotic iron-molybdenum-sulfur cofactor, FeMoco, whose mode of action remains elusive. In particular, the molecular and electronic structure of the N2-binding E4 state is not known. In this study we present theoretical QM/MM calculations of new structural models of the E4 state of molybdenum-dependent nitrogenase and compare to previously suggested models for this enigmatic redox state. We propose two models as possible candidates for the E4 state. Both models feature two hydrides on the FeMo cofactor, bridging atoms Fe2 and Fe6 with a terminal sulfhydryl group on either Fe2 or Fe6 (derived from the S2B bridge) and the change in coordination results in local lower-spin electronic structure at Fe2 and Fe6. These structures appear consistent with the bridging hydride proposal put forward from ENDOR studies and are calculated to be lower in energy than other proposed models for E4 at the TPSSh-QM/MM level of theory. We critically analyze the DFT method dependency in calculations of FeMoco that has resulted in strikingly different proposals for this state. Importantly, dinitrogen binds exothermically to either Fe2 or Fe6 in our models, contrary to others, an effect rationalized via the unique ligand field (from the hydrides) at the Fe with an empty coordination site. A low-spin Fe site is proposed as being important to N2 binding. Furthermore, the geometries of these states suggest a feasible reductive elimination step that could follow, as experiments indicate. Via this step, two electrons are released, reducing the cofactor to yield a distorted 4-coordinate Fe2 or Fe6 that partially activates N2. We speculate that stabilization of an N2-bound Fe(i) at Fe6 (not found for Fe2 model) via reductive elimination is a crucial part of N2 activation in nitrogenases, possibly aided by the apical heterometal ion (Mo or V). By using protons from the sulfhydryl group (to regenerate the sulfide bridge between Fe2 and Fe6) and the nearby homocitrate hydroxy group, we calculate a plausible route to yield a diazene intermediate. This is found to be more favorable with the Fe6-bound model than the Fe2-bound model; however, this protonation is uphill in energy, suggesting protonation of N2 might occur later in the catalytic cycle or via another mechanism.

Resolving the structure of the E1 state of Mo nitrogenase through Mo and Fe K-edge EXAFS and QM/MM calculations

Biological nitrogen fixation is predominately accomplished through Mo nitrogenase, which utilizes a complex MoFe7S9C catalytic cluster to reduce N2 to NH3. This cluster requires the accumulation of three to four reducing equivalents prior to binding N2; however, despite decades of research, the intermediate states formed prior to N2 binding are still poorly understood. Herein, we use Mo and Fe K-edge X-ray absorption spectroscopy and QM/MM calculations to investigate the nature of the E1 state, which is formed following the addition of the first reducing equivalent to Mo nitrogenase. By analyzing the extended X-ray absorption fine structure (EXAFS) region, we provide structural insight into the changes that occur in the metal clusters of the protein when forming the E1 state, and use these metrics to assess a variety of possible models of the E1 state. The combination of our experimental and theoretical results supports that formation of E1 involves an Fe-centered reduction combined with the protonation of a belt-sulfide of the cluster. Hence, these results provide critical experiment and computational insight into the mechanism of this important enzyme.

The role of the dihedral angle and excited cation states in ionization and dissociation of mono-halogenated biphenyls; a combined experimental and theoretical coupled cluster study

We present a combined theoretical and experimental study on the ionization and primary fragmentation channels of the mono-halogenated biphenyls; 2-chlorobiphenyl, 2-bromobiphenyl and 2-iodobiphenyl. The ionization energies (IEs) of the 2-halobiphenyls and the appearance energies (AEs) of the principal fragments are determined through electron impact ionization, while quantum mechanical calculations at the coupled cluster level of theory are used to elucidate the observed processes and the associated dynamics. The primary fragmentation channels are the direct loss of the halogen upon ionization, the loss of the respective hydrogen halides (HX) as well as loss of the hydrogen halide and an additional hydrogen. We find that the dihedral angle strongly influences the relative potential energy of the neutral and the cation on their respective ground state surfaces, an effect caused by the strong influence of the nuclear motion on the conjugation between the phenyl rings. For the principal dissociative ionization channels from the mono-halogenated biphenyls we reason that these can not be described as statistical decay from the ground state cation, but must rather be understood as direct, state-selective processes from specific excited cationic states characterized through local ionization of either the halogenated or the non-substituted phenyl ring.

Multistep Explicit Solvation Protocol for Calculation of Redox Potentials

The calculation of molecular redox potentials in aqueous solution presents a challenge to quantum chemistry due to the need to calculate charged, open-shell species experiencing large solvent effects. Traditionally, redox potentials are calculated via the use of density functional theory and continuum solvation methods, but such protocols have been found to often suffer from large errors, particularly in the case of aqueous solution. While explicit solvation models hold promise of higher accuracy to describe solvent effects in general, their complicated use and lack of well-defined, reliable protocols has hindered their adoption. In this study, we present an explicit-solvation-based approach for the calculation of molecular redox potentials. We combine the use of affordable semiempirical QM/MM molecular dynamics (making use of the recently proposed GFN-xTB method by Grimme et al.) for both redox states and use the linear response approximation to relate vertical ionization energies to the adiabatic redox potential. Simulation length, averaging over snapshots, and accounting for bulk and polarization effects are systematically evaluated using phenol as a working example. We find that it is crucial to reliably account for bulk solvation effects in these calculations, as well as polarization effects which we divide up into short-range and long-range contributions. The short-range polarization contribution is accounted for via QM-region expansion, while the long-range contribution is accounted for via Drude-polarizable QM/MM. Our multistep protocol has been coded to be used in a fully automatic way in a local version of Chemshell. It has been evaluated on a test set of oxidation potentials of organic molecules and found to give gas-solution redox shifts with a mean absolute error of 0.13 eV with respect to experiment, compared to mean absolute errors of 0.26 and 0.21 eV with CPCM and SMD continuum models, respectively.

Electron interactions with the heteronuclear carbonyl precursor H2FeRu3(CO)13 and comparison with HFeCo3(CO)12: from fundamental gas phase and surface science studies to focused electron beam induced deposition

In the current contribution we present a comprehensive study on the heteronuclear carbonyl complex H2FeRu3(CO)13 covering its low energy electron induced fragmentation in the gas phase through dissociative electron attachment (DEA) and dissociative ionization (DI), its decomposition when adsorbed on a surface under controlled ultrahigh vacuum (UHV) conditions and exposed to irradiation with 500 eV electrons, and its performance in focused electron beam induced deposition (FEBID) at room temperature under HV conditions. The performance of this precursor in FEBID is poor, resulting in maximum metal content of 26 atom % under optimized conditions. Furthermore, the Ru/Fe ratio in the FEBID deposit (≈3.5) is higher than the 3:1 ratio predicted. This is somewhat surprising as in recent FEBID studies on a structurally similar bimetallic precursor, HFeCo3(CO)12, metal contents of about 80 atom % is achievable on a routine basis and the deposits are found to maintain the initial Co/Fe ratio. Low temperature (≈213 K) surface science studies on thin films of H2FeRu3(CO)13 demonstrate that electron stimulated decomposition leads to significant CO desorption (average of 8-9 CO groups per molecule) to form partially decarbonylated intermediates. However, once formed these intermediates are largely unaffected by either further electron irradiation or annealing to room temperature, with a predicted metal content similar to what is observed in FEBID. Furthermore, gas phase experiments indicate formation of Fe(CO)4 from H2FeRu3(CO)13 upon low energy electron interaction. This fragment could desorb at room temperature under high vacuum conditions, which may explain the slight increase in the Ru/Fe ratio of deposits in FEBID. With the combination of gas phase experiments, surface science studies and actual FEBID experiments, we can offer new insights into the low energy electron induced decomposition of this precursor and how this is reflected in the relatively poor performance of H2FeRu3(CO)13 as compared to the structurally similar HFeCo3(CO)12.

QM/MM calculations reveal a bridging hydroxo group in a vanadium nitrogenase crystal structure

A recent crystal structure of VFe protein is revealed via QM/MM calculations to contain an hydroxo group, likely derived from water.

Low energy electron-induced decomposition of (η5-Cp)Fe(CO)2Mn(CO)5, a potential bimetallic precursor for focused electron beam induced deposition of alloy structures

Low energy electron-induced decomposition of a potential bimetallic nanofabrication precursor is studied in gas-phase, at surfaces and by quantum chemical calculations.

Formation and decay of negative ion states up to 11 eV above the ionization energy of the nanofabrication precursor HFeCo3(CO)12

Dissociative electron attachment, 11 eV above the ionization energy of the focused electron beam induced deposition (FEBID) precursor HFeCo₃(CO)₁₂. A unique observation with potential significance for FEBID precursor design.

In single electron collisions with the heteronuclear metal carbonyl compound HFeCo₃(CO)₁₂ we observe the formation of long-lived negative ion states up to about 20 eV, 11 eV above its ionization energy. These transient negative ions (TNIs) relax through dissociation (dissociative electron attachment, DEA), losing up to all 12 CO ligands, demonstrating their resilience towards reemission of the captured electron – even at such very high energies. This is unique in DEA and we hypothesize that this phenomenon is rooted in the orbital structure enabling a scaffold of multi-particle, electronically excited resonances. We support this with calculated MO-diagrams revealing dense bands of energy levels near the HOMO–LUMO gap. HFeCo₃(CO)₁₂ is a promising focused electron beam induced deposition (FEBID) precursor and we argue that its unusual DEA behavior relates to its exceptional performance in FEBID. This may be general to a class of molecules with high potential for nano-fabrication by FEBID.

Proton Shuttling and Reaction Paths in Dissociative Electron Attachment to o- and p-Tetrafluorohydroquinone, an Experimental and Theoretical Study

Here we present a combined experimental and theoretical study on the fragmentation of o- and p-tetrafluorohydroquinone upon low energy electron attachment. Despite an identical ring-skeleton and identical functional groups in these constitutional isomers, they show distinctly different fragmentation patterns, a phenomenon that cannot be explained by distinct resonances or different thermochemistry. Using high-level quantum chemical calculations with the computationally affordable domain localized pair natural orbital approach, DLPNO-CCSD(T), we are able to provide a complete and accurate description of the respective reaction dynamics, revealing proton shuttling and transition states for competing channels as the explanation for the different behavior of these isomers. The results represent a "schoolbook example" of how the combination of experiment and modern high-level theory may today provide a thorough understanding of complex reaction dynamics by computationally affordable means.

Comparative electronic structures of nitrogenase FeMoco and FeVco

An investigation of the active site cofactors of the molybdenum and vanadium nitrogenases (FeMoco and FeVco) was performed using high-resolution X-ray spectroscopy. Synthetic heterometallic iron-sulfur cluster models and density functional theory calculations complement the study of the MoFe and VFe holoproteins using both non-resonant and resonant X-ray emission spectroscopy. Spectroscopic data show the presence of direct iron-heterometal bonds, which are found to be weaker in FeVco. Furthermore, the interstitial carbide is found to perturb the electronic structures of the cofactors through highly covalent Fe-C bonding. The implications of these conclusions are discussed in light of the differential reactivity of the molybdenum and vanadium nitrogenases towards various substrates. Possible functional roles for both the heterometal and the interstitial carbide are detailed.

X-ray Absorption and Emission Spectroscopic Studies of [L2Fe2S2](n) Model Complexes: Implications for the Experimental Evaluation of Redox States in Iron-Sulfur Clusters

Herein, a systematic study of [L2Fe2S2](n) model complexes (where L = bis(benzimidazolato) and n = 2-, 3-, 4-) has been carried out using iron and sulfur K-edge X-ray absorption (XAS) and iron Kβ and valence-to-core X-ray emission spectroscopies (XES). These data are used as a test set to evaluate the relative strengths and weaknesses of X-ray core level spectroscopies in assessing redox changes in iron-sulfur clusters. The results are correlated to density functional theory (DFT) calculations of the spectra in order to further support the quantitative information that can be extracted from the experimental data. It is demonstrated that due to canceling effects of covalency and spin state, the information that can be extracted from Fe Kβ XES mainlines is limited. However, a careful analysis of the Fe K-edge XAS data shows that localized valence vs delocalized valence species may be differentiated on the basis of the pre-edge and K-edge energies. These findings are then applied to existing literature Fe K-edge XAS data on the iron protein, P-cluster, and FeMoco sites of nitrogenase. The ability to assess the extent of delocalization in the iron protein vs the P-cluster is highlighted. In addition, possible charge states for FeMoco on the basis of Fe K-edge XAS data are discussed. This study provides an important reference for future X-ray spectroscopic studies of iron-sulfur clusters.

The Fe-V Cofactor of Vanadium Nitrogenase Contains an Interstitial Carbon Atom

The first direct evidence is provided for the presence of an interstitial carbide in the Fe-V cofactor of Azotobacter vinelandii vanadium nitrogenase. As for our identification of the central carbide in the Fe-Mo cofactor, we employed Fe Kβ valence-to-core X-ray emission spectroscopy and density functional theory calculations, and herein report the highly similar spectra of both variants of the cofactor-containing protein. The identification of an analogous carbide, and thus an atomically homologous active site in vanadium nitrogenase, highlights the importance and influence of both the interstitial carbide and the identity of the heteroatom on the electronic structure and catalytic activity of the enzyme.

Molecular structure of 1,2-bis(trifluoromethyl)-1,1,2,2-tetramethyldisilane in the gas, liquid, and solid phases: unusual conformational changes between phases

The molecular structure of 1,2-bis(trifluoromethyl)-1,1,2,2-tetramethyldisilane has been determined in three different phases (solid, liquid, and gas) using various spectroscopic and diffraction techniques. Both the solid-state and gas-phase investigations revealed only one conformer to be present in the sample analyzed, whereas the liquid phase revealed the presence of three conformers. The data have been reproduced using computational methods and a rationale is presented for the observation of three conformers in the liquid state.

The discovery of Mo(III) in FeMoco: reuniting enzyme and model chemistry

Biological nitrogen fixation is enabled by molybdenum-dependent nitrogenase enzymes, which effect the reduction of dinitrogen to ammonia using an Fe₇MoS₉C active site, referred to as the iron molybdenum cofactor or FeMoco. In this mini-review, we summarize the current understanding of the molecular and electronic structure of FeMoco. The advances in our understanding of the active site structure are placed in context with the parallel evolution of synthetic model studies. The recent discovery of Mo(III) in the FeMoco active site is highlighted with an emphasis placed on the important role that model studies have played in this finding. In addition, the reactivities of synthetic models are discussed in terms of their relevance to the enzymatic system.

Molybdenum L-Edge XAS Spectra of MoFe Nitrogenase

A molybdenum L-edge X-ray absorption spectroscopy (XAS) study is presented for native and oxidized MoFe protein of nitrogenase as well as Mo-Fe model compounds. Recently collected data on MoFe protein (in oxidized and reduced forms) is compared to previously published Mo XAS data on the isolated FeMo cofactor in NMF solution and put in context of the recent Mo K-edge XAS study, which showed a Mo^III assignment for the molybdenum atom in FeMoco. The L₃-edge data are interpreted within a simple ligand-field model, from which a time-dependent density functional theory (TDDFT) approach is proposed as a way to provide further insights into the analysis of the molybdenum L₃-edges. The calculated results reproduce well the relative spectral trends that are observed experimentally. Ultimately, these results give further support for the Mo^III assignment in protein-bound FeMoco, as well as isolated FeMoco.

High-resolution molybdenum K-edge X-ray absorption spectroscopy analyzed with time-dependent density functional theory

X-ray absorption spectroscopy (XAS) is a widely used experimental technique capable of selectively probing the local structure around an absorbing atomic species in molecules and materials. When applied to heavy elements, however, the quantitative interpretation can be challenging due to the intrinsic spectral broadening arising from the decrease in the core-hole lifetime. In this work we have used high-energy resolution fluorescence detected XAS (HERFD-XAS) to investigate a series of molybdenum complexes. The sharper spectral features obtained by HERFD-XAS measurements enable a clear assignment of the features present in the pre-edge region. Time-dependent density functional theory (TDDFT) has been previously shown to predict K-pre-edge XAS spectra of first row transition metal compounds with a reasonable degree of accuracy. Here we extend this approach to molybdenum K-edge HERFD-XAS and present the necessary calibration. Modern pure and hybrid functionals are utilized and relativistic effects are accounted for using either the Zeroth Order Regular Approximation (ZORA) or the second order Douglas-Kroll-Hess (DKH2) scalar relativistic approximations. We have found that both the predicted energies and intensities are in excellent agreement with experiment, independent of the functional used. The model chosen to account for relativistic effects also has little impact on the calculated spectra. This study provides an important calibration set for future applications of molybdenum HERFD-XAS to complex catalytic systems.

Identification of a spin-coupled Mo(iii) in the nitrogenase iron–molybdenum cofactor

The molybdenum atom in FeMoco is imperative to the high activity of the enzyme and has been proposed to be Mo(iv). We demonstrate that only Mo(iii) fits Mo HERFD XAS data, the first example of Mo(iii) in biology. Theoretical calculations further reveal an unusual spin-coupled Mo(iii).

Conformational Properties of 1-Halogenated-1-Silacyclohexanes, C5H10SiHX (X = Cl, Br, I): Gas Electron Diffraction, Low-Temperature NMR, Temperature-Dependent Raman Spectroscopy, and Quantum-Chemical Calculations

The molecular structures of axial and equatorial conformers of cyclo-C₅H₁₀SiHX (X = Cl, Br, I) as well as the thermodynamic equilibrium between these species was investigated by means of gas electron diffraction, dynamic nuclear magnetic resonance, temperature-dependent Raman spectroscopy, and quantum-chemical calculations applying CCSD(T), MP2, and DFT methods. According to the experimental and calculated results, all three compounds exist as a mixture of two chair conformers of the six-membered ring. The two chair forms of C_s symmetry differ in the axial or equatorial position of the X atom. In all cases, the axial conformer is preferred over the equatorial one. When the experimental uncertainties are taken into account, all of the experimental and theoretical results for the conformational energy (E_axial - E_equatorial) fit into a remarkably narrow range of -0.50 ± 0.15 kcal mol^-1. It was found by NBO analysis that the axial conformers are unfavorable in terms of steric energy and conjugation effects and that they are stabilized mainly by electrostatic interactions. The conformational energies for C₆H₁₁X and cyclo-C₅H₁₀SiHX (X = F, Cl, Br, I, At) were compared using CCSD(T) calculations. In both series, fluorine is predicted to have a lower conformational preference (cyclohexane equatorial, silacyclohexane axial) than Cl, Br, and I. It is predicted that astatine would behave very similarly to Cl, Br, and I within each series.

Modelling zwitterions in solution: 3-fluoro-γ-aminobutyric acid (3F-GABA)

The conformations and relative stabilities of folded and extended 3-fluoro-γ-aminobutyric acid (3F-GABA) conformers were studied using explicit solvation models. Geometry optimisations in the gas phase with one or two explicit water molecules favour folded and neutral structures containing intramolecular NH···O-C hydrogen bonds. With three or five explicit water molecules zwitterionic minima are obtained, with folded structures being preferred over extended conformers. The stability of folded versus extended zwitterionic conformers increases on going from a PCM continuum solvation model to the microsolvated complexes, though extended structures become less disfavoured with the inclusion of more water molecules. Full explicit solvation was studied with a hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme and molecular dynamics simulations, including more than 6000 TIP3P water molecules. According to free energies obtained from thermodynamic integration at the PM3/MM level and corrected for B3LYP/MM total energies, the fully extended conformer is more stable than folded ones by about -4.5 kJ mol(-1). B3LYP-computed (3)J(F,H) NMR spin-spin coupling constants, averaged over PM3/MM-MD trajectories, agree best with experiment for this fully extended form, in accordance with the original NMR analysis. The seeming discrepancy between static PCM calculations and experiment noted previously is now resolved. That the inexpensive semiempirical PM3 method performs so well for this archetypical zwitterion is encouraging for further QM/MM studies of biomolecular systems.