Noncovalent interactions that tune the reactivities of the flavins in bifurcating electron transferring flavoprotein.

Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize noncovalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin (ETflavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin (Bfflavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of ETflavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of ETflavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement, and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the Bfflavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the Bfflavin O4. This would stabilize the anionic Bfflavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.

Ultrafast protein response in the Pfr state of Cph1 phytochrome.

Photoisomerization is a fundamental process in several classes of photoreceptors. Phytochromes sense red and far-red light in their Pr and Pfr states, respectively. Upon light absorption, these states react via individual photoreactions to the other state. Cph1 phytochrome shows a photoisomerization of its phycocyanobilin (PCB) chromophore in the Pfr state with a time constant of 0.7 ps. The dynamics of the PCB chromophore has been described, but whether or not the apoprotein exhibits an ultrafast response too, is not known. Here, we compare the photoreaction of 13C/15N labeled apoprotein with unlabeled apoprotein to unravel ultrafast apoprotein dynamics in Cph1. In the spectral range from 1750 to 1620 cm-1 we assigned several signals due to ultrafast apoprotein dynamics. A bleaching signal at 1724 cm-1 is tentatively assigned to deprotonation of a carboxylic acid, probably Asp207, and signals around 1670 cm-1 are assigned to amide I vibrations of the capping helix close to the chromophore. These signals remain after photoisomerization. The apoprotein dynamics appear upon photoexcitation or concomitant with chromophore isomerization. Thus, apoprotein dynamics occur prior to and after photoisomerization on an ultrafast time-scale. We discuss the origin of the ultrafast apoprotein response with the 'Coulomb hammer' mechanism, i.e. an impulsive change of electric field and Coulombic force around the chromophore upon excitation.

Infrared Spectroscopy Elucidates the Inhibitor Binding Sites in a Metal-Dependent Formate Dehydrogenase.

Biological carbon dioxide (CO2 ) reduction is an important step by which organisms form valuable energy-richer molecules required for further metabolic processes. The Mo-dependent formate dehydrogenase (FDH) from Rhodobacter capsulatus catalyzes reversible formate oxidation to CO2 at a bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor. To elucidate potential substrate binding sites relevant for the mechanism, we studied herein the interaction with the inhibitory molecules azide and cyanate, which are isoelectronic to CO2 and charged as formate. We employed infrared (IR) spectroscopy in combination with density functional theory (DFT) and inhibition kinetics. One distinct inhibitory molecule was found to bind to either a non-competitive or a competitive binding site in the secondary coordination sphere of the active site. Site-directed mutagenesis of key amino acid residues in the vicinity of the bis-MGD cofactor revealed changes in both non-competitive and competitive binding, whereby the inhibitor is in case of the latter interaction presumably bound between the cofactor and the adjacent Arg587.

Photoinduced reaction mechanisms in prototypical and bathy phytochromes.

Phytochromes, found in plants, fungi, and bacteria, exploit light as a source of information to control physiological processes via photoswitching between two states of different physiological activity, i.e. a red-absorbing Pr and a far-red-absorbing Pfr state. Depending on the relative stability in the dark, bacterial phytochromes are divided into prototypical and bathy phytochromes, where the stable state is Pr and Pfr, respectively. In this work we studied representatives of these groups (prototypical Agp1 and bathy Agp2 from Agrobacterium fabrum) together with the bathy-like phytochrome XccBphP from Xanthomonas campestris by resonance Raman and IR difference spectroscopy. In all three phytochromes, the photoinduced conversions display the same mechanistic pattern as reflected by the chromophore structures in the various intermediate states. We also observed in each case the secondary structure transition of the tongue, which is presumably crucial for the function of phytochrome. The three phytochromes differ in details of the chromophore conformation in the various intermediates and the energetic barrier of their respective decay reactions. The specific protein environment in the chromophore pocket, which is most likely the origin for these small differences, also controls the proton transfer processes concomitant to the photoconversions. These proton translocations, which are tightly coupled to the structural transition of the tongue, presumably proceed via the same mechanism along the Pr → Pfr conversion whereas the reverse Pfr → Pr photoconversion includes different proton transfer pathways. Finally, classification of phytochromes in prototypical and bathy (or bathy-like) phytochromes is discussed in terms of molecular structure and mechanistic properties.

Local Electric Field Changes during the Photoconversion of the Bathy Phytochrome Agp2.

Phytochromes switch between a physiologically inactive and active state via a light-induced reaction cascade, which is initiated by isomerization of the tetrapyrrole chromophore and leads to the functionally relevant secondary structure transition of a protein segment (tongue). Although details of the underlying cause-effect relationships are not known, electrostatic fields are likely to play a crucial role in coupling chromophores and protein structural changes. Here, we studied local electric field changes during the photoconversion of the dark state Pfr to the photoactivated state Pr of the bathy phytochrome Agp2. Substituting Tyr165 and Phe192 in the chromophore pocket by para-cyanophenylalanine (pCNF), we monitored the respective nitrile stretching modes in the various states of photoconversion (vibrational Stark effect). Resonance Raman and IR spectroscopic analyses revealed that both pCNF-substituted variants undergo the same photoinduced structural changes as wild-type Agp2. Based on a structural model for the Pfr state of F192pCNF, a molecular mechanical-quantum mechanical approach was employed to calculate the electric field at the nitrile group and the respective stretching frequency, in excellent agreement with the experiment. These calculations serve as a reference for determining the electric field changes in the photoinduced states of F192pCNF. Unlike F192pCNF, the nitrile group in Y165pCNF is strongly hydrogen bonded such that the theoretical approach is not applicable. However, in both variants, the largest changes of the nitrile stretching modes occur in the last step of the photoconversion, supporting the view that the proton-coupled restructuring of the tongue is accompanied by a change of the electric field.

A Resonance Raman Marker Band Characterizes the Slow and Fast Form of Cytochrome c Oxidase.

Cytochrome c oxidase (CcO) in its as-isolated form is known to exist in a slow and fast form, which differ drastically in their ability to bind oxygen and other ligands. While preparation methods have been established that yield either the fast or the slow form of the protein, the underlying structural differences have not been identified yet. Here, we have performed surface enhanced resonance Raman (SERR) spectroscopy of CcO immobilized on electrodes in both forms. SERR spectra obtained in resonance with the heme a3 metal-to-ligand charge transfer (MLCT) transition at 650 nm displayed a sharp vibrational band at 748 or 750 cm-1 when the protein was in its slow or fast form, respectively. DFT calculations identified the band as a mode of the His-419 ligand that is sensitive to the oxygen ligand and the protonation state of Tyr-288 within the binuclear complex. Potential-dependent SERR spectroscopy showed a redox-induced change of this band around 525 mV versus Ag/AgCl exclusively for the fast form, which coincides with the redox potential of the Tyr-O/Tyr-O- transition. Our data points to a peroxide ligand in the resting state of CcO for both forms. The observed frequencies and redox sensitivities of the Raman marker band suggest that a radical Tyr-288 is present in the fast form and a protonated Tyr-288 in the slow form.

Hidden in Plain Sight: Natural Products of Commensal Microbiota as an Environmental Selection Pressure for the Rise of New Variants of SARS-CoV-2.

Since the emergence of SARS-CoV-2, little attention has been paid to the interplay between the interaction of virus and commensal microbiota. Here, we used molecular docking and dynamics simulations to study the interaction of some of the known metabolites and natural products (NPs) produced by commensal microbiota with the receptor binding domain (RBD) of the spike glycoprotein of SARS-CoV-2. The results predict that NPs of commensal microbiota such as bile acids and non-ribosomal peptides (NRPs), of which some are siderophores, bind to the wild-type RBD and interfere with its binding to the ACE2 receptor. N501Y mutation, which is present in many of the emerging variants of the virus, abolishes the predicted binding pocket of bile acids and NRPs. Based on these findings, available experimental data showing that bile acids reduce the binding affinity of wild-type RBD to the ACE2 receptor, and the data suggesting that the respiratory tract microbiota affect viral infection we put forward the following proposal: mutations such as N501Y enable the RBD to bind to the ACE2 receptor more effectively in the presence of NPs produced by the respiratory tract bacteria thereby, increasing the infectivity rate of the virus. We hope our data stimulate future works to better understand the interactions of NPs produced by commensal microbiota with respiratory viruses like SARS-CoV-2.

Ultrafast proton release reaction and primary photochemistry of phycocyanobilin in solution observed with fs-time-resolved mid-IR and UV/Vis spectroscopy.

Deactivation processes of photoexcited (λex = 580 nm) phycocyanobilin (PCB) in methanol were investigated by means of UV/Vis and mid-IR femtosecond (fs) transient absorption (TA) as well as static fluorescence spectroscopy, supported by density-functional-theory calculations of three relevant ground state conformers, PCBA, PCBB and PCBC, their relative electronic state energies and normal mode vibrational analysis. UV/Vis fs-TA reveals time constants of 2.0, 18 and 67 ps, describing decay of PCBB*, of PCBA* and thermal re-equilibration of PCBA, PCBB and PCBC, respectively, in line with the model by Dietzek et al. (Chem Phys Lett 515:163, 2011) and predecessors. Significant substantiation and extension of this model is achieved first via mid-IR fs-TA, i.e. identification of molecular structures and their dynamics, with time constants of 2.6, 21 and 40 ps, respectively. Second, transient IR continuum absorption (CA) is observed in the region above 1755 cm-1 (CA1) and between 1550 and 1450 cm-1 (CA2), indicative for the IR absorption of highly polarizable protons in hydrogen bonding networks (X-H…Y). This allows to characterize chromophore protonation/deprotonation processes, associated with the electronic and structural dynamics, on a molecular level. The PCB photocycle is suggested to be closed via a long living (> 1 ns), PCBC-like (i.e. deprotonated), fluorescent species.

Tracking the route of molecular oxygen in O2-tolerant membrane-bound [NiFe] hydrogenase.

[NiFe] hydrogenases catalyze the reversible splitting of H2 into protons and electrons at a deeply buried active site. The catalytic center can be accessed by gas molecules through a hydrophobic tunnel network. While most [NiFe] hydrogenases are inactivated by O2, a small subgroup, including the membrane-bound [NiFe] hydrogenase (MBH) of Ralstonia eutropha, is able to overcome aerobic inactivation by catalytic reduction of O2 to water. This O2 tolerance relies on a special [4Fe3S] cluster that is capable of releasing two electrons upon O2 attack. Here, the O2 accessibility of the MBH gas tunnel network has been probed experimentally using a "soak-and-freeze" derivatization method, accompanied by protein X-ray crystallography and computational studies. This combined approach revealed several sites of O2 molecules within a hydrophobic tunnel network leading, via two tunnel entrances, to the catalytic center of MBH. The corresponding site occupancies were related to the O2 concentrations used for MBH crystal derivatization. The examination of the O2-derivatized data furthermore uncovered two unexpected structural alterations at the [4Fe3S] cluster, which might be related to the O2 tolerance of the enzyme.

Red, Orange, Green: Light- and Temperature-Dependent Color Tuning in a Cyanobacteriochrome.

Cyanobacteriochromes (CBCRs) are photoreceptor proteins that photoconvert between two parent states and thereby regulate various biological processes. An intriguing property is their variable ultraviolet-visible (UV-vis) absorption that covers the entire spectral range from the far-red to the near-UV region and thus makes CBCRs promising candidates for optogenetic applications. Here, we have studied Slr1393, a CBCR that photoswitches between red- and green-absorbing states (Pr and Pg, respectively). Using UV-vis absorption, fluorescence, and resonance Raman (RR) spectroscopy, a further orange-absorbing state O600 that is in thermal equilibrium with Pr was identified. The different absorption properties of the three states were attributed to the different lengths of the conjugated π-electron system of the phycocyanobilin chromophore. In agreement with available crystal structures and supported by quantum mechanics/molecular mechanics (QM/MM) calculations, the most extended conjugation holds for Pr whereas it is substantially reduced in Pg. Here, the two outer pyrrole rings D and A are twisted out of the plane defined by inner pyrrole rings B and C. For the O600 state, the comparison of the experimental RR spectra with QM/MM-calculated spectra indicates a partially distorted ZZZssa geometry in which ring A is twisted while ring D and the adjacent methine bridge display essentially the same geometry as Pr. The quantitative analysis of temperature-dependent spectra yields an enthalpy barrier of ∼30 kJ/mol for the transition from Pr to O600. This reaction is associated with the movement of a conserved tryptophan residue from the chromophore binding pocket to a solvent-exposed position.

Intramolecular Proton Transfer Controls Protein Structural Changes in Phytochrome.

Phytochromes are biological photoswitches that interconvert between two parent states (Pr and Pfr). The transformation is initiated by photoisomerization of the tetrapyrrole chromophore, followed by a sequence of chromophore and protein structural changes. In the last step, a phytochrome-specific peptide segment (tongue) undergoes a secondary structure change, which in prokaryotic phytochromes is associated with the (de)activation of the output module. The focus of this work is the Pfr-to-Pr photoconversion of the bathy bacteriophytochrome Agp2 in which Pfr is the thermodynamically stable state. Using spectroscopic techniques, we studied the structural and functional consequences of substituting Arg211, Tyr165, His278, and Phe192 close to the biliverdin (BV) chromophore. In Pfr, substitutions of these residues do not affect the BV structure. The characteristic Pfr properties of bathy phytochromes, including the protonated propionic side chain of ring C (propC) of BV, are preserved. However, replacing Arg211 or Tyr165 blocks the photoconversion in the Meta-F state, prior to the secondary structure transition of the tongue and without deprotonation of propC. The Meta-F state of these variants displays low photochemical activity, but electronic excitation causes ultrafast alterations of the hydrogen bond network surrounding the chromophore. In all variants studied here, thermal back conversion from the photoproducts to Pfr is decelerated but substitution of His278 or Phe192 is not critical for the Pfr-to-Pr photoconversion. These variants do not impair deprotonation of propC or the α-helix/ß-sheet transformation of the tongue during the Meta-F-to-Pr decay. Thus, we conclude that propC deprotonation is essential for restructuring of the tongue.

The C-Terminal VPRTES Tail of LL-37 Influences the Mode of Attachment to a Lipid Bilayer and Antimicrobial Activity.

Cathelicidins are a family of host defense antimicrobial peptides in mammalian species. Among them, LL-37 is the only peptide of this family found in humans. Although LL-37 has been intensively investigated in the past, the mode of exerting its bactericidal activity through the specific interactions with bacterial membranes remains elusive. In this work, we combined microbiological and computational approaches with a tool box of experimental biophysical techniques, including conventional and surface-enhanced infrared absorption spectroscopy as well as fluorescence spectroscopy to characterize the structural and dynamic properties of LL-37 and shorter variants adsorbed on POPC/POPG (9:1) lipid bilayers as mimics of bacterial membranes. First, microbiological assays demonstrate that, while LL-32 and, in a lesser degree, LL-37 show hemolysis and antimicrobial activity, LL-20 remains practically inactive. Second, by comparing experimental and computational data of LL-37 with LL-20, we explained the bactericidal activity of the active peptide core as a consequence of an increased flexibility of the peptide structure, leading to reactive dangling charged side chains. Third, permeabilization assays showed a concentration-dependent membrane disruption activity of LL-37 and LL-32: at high peptide concentrations, LL-32 shows higher activity than LL-37, while, at low peptide concentrations, both peptides show similar activities. Responsible for this behavior is the C-terminal VPRTES tail (Ct-VPRTES tail), which, according to atomistic simulations, is able to promote the insertion of the peptide in the membrane and plays an essential role in controlling ordered peptide oligomerization on the surface of the membrane.

Role of the Propionic Side Chains for the Photoconversion of Bacterial Phytochromes.

Bacteriophytochromes harboring a biliverdin IXα (BV) chromophore undergo photoinduced reaction cascades to switch between physiologically inactive and active states. Employing vibrational spectroscopic and computational methods, we analyzed the role of propionic substituents of BV in the transformations between parent states Pr and Pfr in prototypical (Agp1) and bathy (Agp2) phytochromes from Agrobacterium fabrum. Both proteins form adducts with BV monoesters (BVM), esterified at propionic side chain B (PsB) or C (PsC), but in each case, only one monoester adduct is reactive. In the reactive Agp2-BVM-B complex (esterified at ring B), the Pfr dark state displays the structural properties characteristic of bathy phytochromes, including a protonated PsC. As in native Agp2, PsC is deprotonated in the final step of the Pfr phototransformation. However, the concomitant α-helix/ß-sheet secondary structure change of the tongue is blocked at the stage of unfolding of the coiled loop region. This finding and the shift of the tautomeric equilibrium of BVM toward the enol form are attributed to the drastic changes in the electrostatic potential. The calculations further suggest that deprotonation of PsC and the protonation state of His278 control the reactivity of the enol tautomer, thereby accounting for the extraordinarily slow thermal reversion. Although strong perturbations of the electrostatic potential are also found for Agp1-BVM, the consequences for the Pr-to-Pfr phototransformation are less severe. Specifically, the structural transition of the tongue is not impaired and thermal reversion is even accelerated. The different response of Agp1 and Agp2 to monoesterification of BV points to different photoconversion mechanisms.

Monitoring the Orientational Changes of Alamethicin during Incorporation into Bilayer Lipid Membranes.

Antimicrobial peptides (AMPs) are the first line of defense after contact of an infectious invader, for example, bacterium or virus, with a host and an integral part of the innate immune system of humans. Their broad spectrum of biological functions ranges from cell membrane disruption over facilitation of chemotaxis to interaction with membrane-bound or intracellular receptors, thus providing novel strategies to overcome bacterial resistances. Especially, the clarification of the mechanisms and dynamics of AMP incorporation into bacterial membranes is of high interest, and different mechanistic models are still under discussion. In this work, we studied the incorporation of the peptaibol alamethicin (ALM) into tethered bilayer lipid membranes on electrodes in combination with surface-enhanced infrared absorption (SEIRA) spectroscopy. This approach allows monitoring the spontaneous and potential-induced ion channel formation of ALM in situ. The complex incorporation kinetics revealed a multistep mechanism that points to peptide-peptide interactions prior to penetrating the membrane and adopting the transmembrane configuration. On the basis of the anisotropy of the backbone amide I and II infrared absorptions determined by density functional theory calculations, we employed a mathematical model to evaluate ALM reorientations monitored by SEIRA spectroscopy. Accordingly, ALM was found to adopt inclination angles of ca. 69°-78° and 21° in its interfacially adsorbed and transmembrane incorporated states, respectively. These orientations can be stabilized efficiently by the dipolar interaction with lipid head groups or by the application of a potential gradient. The presented potential-controlled mechanistic study suggests an N-terminal integration of ALM into membranes as monomers or parallel oligomers to form ion channels composed of parallel-oriented helices, whereas antiparallel oligomers are barred from intrusion.

Solvation of diclofenac in water from atomistic molecular dynamics simulations - interplay between solute-solute and solute-solvent interactions.

The solubility-permeability relationship of active pharmaceutical ingredients determines the efficacy of their usage. Diclofenac (DCL), which is a widely used nonsteroidal anti-inflammatory drug, is characterized by extremely good membrane permeability, but low water solubility limiting drug effectiveness. The present research focuses on the fundamental explanation of this limitation using the combination of ab initio and classical molecular dynamics simulations of different ionic forms of DCL in water, namely, ionized, un-ionized and the mixture of them both. The analysis of diclofenac solvation in an aqueous environment is used to understand the origin of drug precipitation, especially in gastric pH. The used computational approach reveals the formation of micelle-like self-associated aggregates of diclofenac in water as the result of intermolecular π-π interactions and C-Hπ hydrogen bonds. The DCL aggregation in water is shown to depend mostly on drug concentration, protonation and temperature of the aqueous environment. The detected self-association properties of the drug in water are likely to be of great importance during the development of new drug formulations and fabrication of drug adsorbents for wastewater.

Redox induced protonation of heme propionates in cytochrome c oxidase: Insights from surface enhanced resonance Raman spectroscopy and QM/MM calculations.

Understanding the coupling between heme reduction and proton translocation in cytochrome c oxidase (CcO) is still an open problem. The propionic acids of heme a3 have been proposed to act as a proton loading site (PLS) in the proton pumping pathway, yet this proposal could not be verified by experimental data so far. We have set up an experiment where the redox states of the two hemes in CcO can be controlled via external electrical potential. Surface enhanced resonance Raman (SERR) spectroscopy was applied to simultaneously monitor the redox state of the hemes and the protonation state of the heme propionates. Simulated spectra based on QM/MM calculations were used to assign the resonant enhanced CH2 twisting modes of the propionates to the protonation state of the individual heme a and heme a3 propionates respectively. The comparison between calculated and measured H2OD2O difference spectra allowed a sound band assignment. In the fully reduced enzyme at least three of the four heme propionates were found to be protonated whereas in the presence of a reduced heme a and an oxidized heme a3 only protonation of one heme a3 propionates was observed. Our data supports the postulated scenario where the heme a3 propionates are involved in the proton pathway.

Carbon Monoxide Dehydrogenase Reduces Cyanate to Cyanide.

The biocatalytic function of carbon monoxide dehydrogenase (CODH) has a high environmental relevance owing to its ability to reduce CO2 . Despite numerous studies on CODH over the past decades, its catalytic mechanism is not yet fully understood. In the present combined spectroscopic and theoretical study, we report first evidences for a cyanate (NCO- ) to cyanide (CN- ) reduction at the C-cluster. The adduct remains bound to the catalytic center to form the so-called CN- -inhibited state. Notably, this conversion does not occur in crystals of the Carboxydothermus hydrogenoformans CODH enzyme (CODHIICh ), as indicated by the lack of the corresponding CN- stretching mode. The transformation of NCO- , which also acts as an inhibitor of the two-electron-reduced Cred2 state of CODH, could thus mimic CO2 turnover and open new perspectives for elucidation of the detailed catalytic mechanism of CODH.

The Molybdenum Active Site of Formate Dehydrogenase Is Capable of Catalyzing C-H Bond Cleavage and Oxygen Atom Transfer Reactions.

Formate dehydrogenases (FDHs) are capable of performing the reversible oxidation of formate and are enzymes of great interest for fuel cell applications and for the production of reduced carbon compounds as energy sources from CO2. Metal-containing FDHs in general contain a highly conserved active site, comprising a molybdenum (or tungsten) center coordinated by two molybdopterin guanine dinucleotide molecules, a sulfido and a (seleno-)cysteine ligand, in addition to a histidine and arginine residue in the second coordination sphere. So far, the role of these amino acids in catalysis has not been studied in detail, because of the lack of suitable expression systems and the lability or oxygen sensitivity of the enzymes. Here, the roles of these active site residues is revealed using the Mo-containing FDH from Rhodobacter capsulatus. Our results show that the cysteine ligand at the Mo ion is displaced by the formate substrate during the reaction, the arginine has a direct role in substrate binding and stabilization, and the histidine elevates the pKa of the active site cysteine. We further found that in addition to reversible formate oxidation, the enzyme is further capable of reducing nitrate to nitrite. We propose a mechanistic scheme that combines both functionalities and provides important insights into the distinct mechanisms of C-H bond cleavage and oxygen atom transfer catalyzed by formate dehydrogenase.

Domain motions and electron transfer dynamics in 2Fe-superoxide reductase.

Superoxide reductases are non-heme iron enzymes that represent valuable model systems for the reductive detoxification of reactive oxygen species. In the present study, we applied different theoretical methods to study the structural dynamics of a prototypical 2Fe-superoxide reductase and its influence on electron transfer towards the active site. Using normal mode and essential dynamics analyses, we could show that enzymes of this type are capable of well-defined, electrostatically triggered domain movements, which may allow conformational proofreading for cellular redox partners involved in intermolecular electron transfer. Moreover, these global modes of motion were found to enable access to molecular configurations with decreased tunnelling distances between the active site and the enzyme's second iron centre. Using all-atom classical molecular dynamics simulations and the tunnelling pathway model, however, we found that electron transfer between the two metal sites is not accelerated under these conditions. This unexpected finding suggests that the unperturbed enzymatic structure is optimized for intramolecular electron transfer, which provides an indirect indication of the biological relevance of such a mechanism. Consistently, efficient electron transfer was found to depend on a distinct route, which is accessible via the equilibrium geometry and characterized by a quasi conserved tyrosine that could enable multistep-tunnelling (hopping). Besides these explicit findings, the present study demonstrates the importance of considering both global and local protein dynamics, and a generalized approach for the functional analysis of these aspects is provided.