Electronic Structure Investigation and Parametrization of Biologically Relevant Iron–Sulfur Clusters

Elucidation of factors shaping reactivity of 5'-deoxyadenosyl - a prominent organic radical in biology

This study investigates the factors modulating the reactivity of 5'-deoxyadenosyl (5'dAdo˙) radical, a potent hydrogen atom abstractor that forms in the active sites of radical SAM enzymes and that otherwise undergoes a rapid self-decay in aqueous solution. Here, we compare hydrogen atom abstraction (HAA) reactions between native substrates of radical SAM enzymes and 5'dAdo˙ in aqueous solution and in two enzymatic microenvironments. With that we reveal that HAA efficiency of 5'dAdo˙ is due to (i) the in situ formation of 5'dAdo˙ in a pre-ordered complex with a substrate, which attenuates the unfavorable effect of substrate:5'dAdo˙ complex formation, and (ii) the prevention of the conformational changes associated with self-decay by a tight active-site cavity. The enzymatic cavity, however, does not have a strong effect on the HAA activity of 5'dAdo˙. Thus, we performed an analysis of in-water HAA performed by 5'dAdo˙ based on a three-component thermodynamic model incorporating the diagonal effect of the free energy of reaction, and the off-diagonal effect of asynchronicity and frustration. To this aim, we took advantage of the straightforward relationship between the off-diagonal thermodynamic effects and the electronic-structure descriptor - the redistribution of charge between the reactants during the reaction. It allows to access HAA-competent redox and acidobasic properties of 5'dAdo˙ that are otherwise unavailable due to its instability upon one-electron reduction and protonation. The results show that all reactions feature a favourable thermodynamic driving force and tunneling, the latter of which lowers systematically barriers by ∼2 kcal mol-1. In addition, most of the reactions experience a favourable off-diagonal thermodynamic contribution. In HAA reactions, 5'dAdo˙ acts as a weak oxidant as well as a base, also 5'dAdo˙-promoted HAA reactions proceed with a quite low degree of asynchronicity of proton and electron transfer. Finally, the study elucidates the crucial and dual role of asynchronicity. It directly lowers the barrier as a part of the off-diagonal thermodynamic contribution, but also indirectly increases the non-thermodynamic part of the barrier by presumably controlling the adiabatic coupling between proton and electron transfer. The latter signals that the reaction proceeds as a hydrogen atom transfer rather than a proton-coupled electron transfer.

Structural and biophysical properties of a [4Fe4S] ferredoxin-like protein from Synechocystis sp. PCC 6803 with a unique two domain structure

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron‑sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2×[4Fe4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 ± 9 mV and - 516 ± 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

Structure of the human respiratory complex II

Human complex II is a key protein complex that links two essential energy-producing processes: the tricarboxylic acid cycle and oxidative phosphorylation. Deficiencies due to mutagenesis have been shown to cause mitochondrial disease and some types of cancers. However, the structure of this complex is yet to be resolved, hindering a comprehensive understanding of the functional aspects of this molecular machine. Here, we have determined the structure of human complex II in the presence of ubiquinone at 2.86 Å resolution by cryoelectron microscopy, showing it comprises two water-soluble subunits, SDHA and SDHB, and two membrane-spanning subunits, SDHC and SDHD. This structure allows us to propose a route for electron transfer. In addition, clinically relevant mutations are mapped onto the structure. This mapping provides a molecular understanding to explain why these variants have the potential to produce disease.

The nucleotide excision repair proteins through the lens of molecular dynamics simulations

Mutations that affect the proteins responsible for the nucleotide excision repair (NER) pathway can lead to diseases such as xeroderma pigmentosum, trichothiodystrophy, Cockayne syndrome, and Cerebro-oculo-facio-skeletal syndrome. Hence, understanding their molecular behavior is needed to elucidate these diseases' phenotypes and how the NER pathway is organized and coordinated. Molecular dynamics techniques enable the study of different protein conformations, adaptable to any research question, shedding light on the dynamics of biomolecules. However, as important as they are, molecular dynamics studies focused on DNA repair pathways are still becoming more widespread. Currently, there are no review articles compiling the advancements made in molecular dynamics approaches applied to NER and discussing: (i) how this technique is currently employed in the field of DNA repair, focusing on NER proteins; (ii) which technical setups are being employed, their strengths and limitations; (iii) which insights or information are they providing to understand the NER pathway or NER-associated proteins; (iv) which open questions would be suited for this technique to answer; and (v) where can we go from here. These questions become even more crucial considering the numerous 3D structures published regarding the NER pathway's proteins in recent years. In this work, we tackle each one of these questions, revising and critically discussing the results published in the context of the NER pathway.

Accurate prediction by AlphaFold2 for ligand binding in a reductive dehalogenase and implications for PFAS (per- and polyfluoroalkyl substance) biodegradation

Despite the success of AlphaFold2 (AF2), it is unclear how AF2 models accommodate for ligand binding. Here, we start with a protein sequence from Acidimicrobiaceae TMED77 (T7RdhA) with potential for catalyzing the degradation of per- and polyfluoroalkyl substances (PFASs). AF2 models and experiments identified T7RdhA as a corrinoid iron-sulfur protein (CoFeSP) which uses a norpseudo-cobalamin (BVQ) cofactor and two Fe₄S₄ iron-sulfur clusters for catalysis. Docking and molecular dynamics simulations suggest that T7RdhA uses perfluorooctanoic acetate (PFOA) as a substrate, supporting the reported defluorination activity of its homolog, A6RdhA. We showed that AF2 provides processual (dynamic) predictions for the binding pockets of ligands (cofactors and/or substrates). Because the pLDDT scores provided by AF2 reflect the protein native states in complex with ligands as the evolutionary constraints, the Evoformer network of AF2 predicts protein structures and residue flexibility in complex with the ligands, i.e., in their native states. Therefore, an apo-protein predicted by AF2 is actually a holo-protein awaiting ligands.

Redox sensitive human mitochondrial aconitase and its interaction with frataxin: In vitro and in silico studies confirm that it takes two to tango

Mitochondrial aconitase (ACO2) has been postulated as a redox sensor in the tricarboxylic acid cycle. Its high sensitivity towards reactive oxygen and nitrogen species is due to its particularly labile [4Fe-4S]²⁺ prosthetic group which yields an inactive [3Fe-4S]⁺ cluster upon oxidation. Moreover, ACO2 was found as a main oxidant target during aging and in pathologies where mitochondrial dysfunction is implied. Herein, we report the expression and characterization of recombinant human ACO2 and its interaction with frataxin (FXN), a protein that participates in the de novo biosynthesis of Fe-S clusters. A high yield of pure ACO2 (≥99%, 22 ± 2 U/mg) was obtained and kinetic parameters for citrate, isocitrate, and cis-aconitate were determined. Superoxide, carbonate radical, peroxynitrite, and hydrogen peroxide reacted with ACO2 with second-order rate constants of 10⁸, 10⁸, 10⁵, and 10² M^-1 s^-1, respectively. Temperature-induced unfolding assessed by tryptophan fluorescence of ACO2 resulted in apparent melting temperatures of 51.1 ± 0.5 and 43.6 ± 0.2 °C for [4Fe-4S]²⁺ and [3Fe-4S]⁺ states of ACO2, sustaining lower thermal stability upon cluster oxidation. Differences in protein dynamics produced by the Fe-S cluster redox state were addressed by molecular dynamics simulations. Reactivation of [3Fe-4S]⁺-ACO2 by FXN was verified by activation assays and direct iron-dependent interaction was confirmed by protein-protein interaction ELISA and fluorescence spectroscopic assays. Multimer modeling and protein-protein docking predicted an ACO2-FXN complex where the metal ion binding region of FXN approaches the [3Fe-4S]⁺ cluster, supporting that FXN is a partner for reactivation of ACO2 upon oxidative cluster inactivation.

Biomimetic Photo-Switches Softening Model Lipid Membranes

We report the synthesis and characterization of a novel photo-switch based on biomimetic cyclocurcumin analogous and interacting with the lipid bilayer, which can be used in the framework of oxygen-independent light-induced therapy. More specifically, by using molecular dynamics simulations and free energy techniques, we show that the inclusion of hydrophobic substituents is needed to allow insertion in the lipid membrane. After having confirmed experimentally that the substituents do not preclude the efficient photoisomerization, we show through UV-vis and dynamic light scattering measurements together with compression isotherms that the chromophore is internalized in both lipid vesicles and monomolecular film, respectively, inducing their fluidification. The irradiation of the chromophore-loaded lipid aggregates modifies their properties due to the different organization of the two diastereoisomers, E and Z. In particular, a competition between a fast structural reorganization and a slower expulsion of the chromophore after isomerization can be observed in the kinetic profiles recorded during E to Z photoisomerization. This report paves the way for future investigations in the optimization of biomimetic photoswitches potentially useful in modern light-induced therapeutic strategies.

Effects of Active-Center Reduction of Plant-Type Ferredoxin on Its Structure and Dynamics: Computational Analysis Using Molecular Dynamics Simulations

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

Unravelling the role of transient redox partner complexes in P450 electron transfer mechanics

The molecular evolution of cytochromes P450 and associated redox-driven oxidative catalysis remains a mystery in biology. It is widely believed that sterol 14α-demethylase (CYP51), an essential enzyme of sterol biosynthesis, is the ancestor of the whole P450 superfamily given its conservation across species in different biological kingdoms. Herein we have utilized X-ray crystallography, molecular dynamics simulations, phylogenetics and electron transfer measurements to interrogate the nature of P450-redox partner binding using the naturally occurring fusion protein, CYP51-ferredoxin found in the sterol-producing bacterium Methylococcus capsulatus. Our data advocates that the electron transfer mechanics in the M. capsulatus CYP51-ferredoxin fusion protein involves an ensemble of ferredoxin molecules in various orientations and the interactions are transient. Close proximity of ferredoxin, however, is required to complete the substrate-induced large-scale structural switch in the P450 domain that enables proton-coupled electron transfer and subsequent oxygen scission and catalysis. These results have fundamental implications regarding the early evolution of electron transfer proteins and for the redox reactions in the early steps of sterol biosynthesis. They also shed new light on redox protein mechanics and the subsequent diversification of the P450 electron transfer machinery in nature.

Activation of Glycyl Radical Enzymes─Multiscale Modeling Insights into Catalysis and Radical Control in a Pyruvate Formate-Lyase-Activating Enzyme

Pyruvate formate-lyase (PFL) is a glycyl radical enzyme (GRE) playing a pivotal role in the metabolism of strict and facultative anaerobes. Its activation is carried out by a PFL-activating enzyme, a member of the radical S-adenosylmethionine (rSAM) superfamily of metalloenzymes, which introduces a glycyl radical into the Gly radical domain of PFL. The activation mechanism is still not fully understood and is structurally based on a complex with a short model peptide of PFL. Here, we present extensive molecular dynamics simulations in combination with quantum mechanics/molecular mechanics (QM/MM)-based kinetic and thermodynamic reaction evaluations of a more complete activation model comprising the 49 amino acid long C-terminus region of PFL. We reveal the benefits and pitfalls of the current activation model, providing evidence that the bound peptide conformation does not resemble the bound protein-protein complex conformation with PFL, with implications for the activation process. Substitution of the central glycine with (S)- and (R)-alanine showed excellent binding of (R)-alanine over unstable binding of (S)-alanine. Radical stabilization calculations indicate that a higher radical stability of the glycyl radical might not be the sole origin of the evolutionary development of GREs. QM/MM-derived radical formation kinetics further demonstrate feasible activation barriers for both peptide and C-terminus activation, demonstrating why the crystalized model peptide system is an excellent inhibitory system for natural activation. This new evidence supports the theory that GREs converged on glycyl radical formation due to the better conformational accessibility of the glycine radical loop, rather than the highest radical stability of the formed peptide radicals.

Autophagy and evasion of the immune system by SARS-CoV-2. Structural features of the non-structural protein 6 from wild type and Omicron viral strains interacting with a model lipid bilayer

The viral cycle of SARS-CoV-2 is based on a complex interplay with the cellular machinery, which is mediated by specific proteins eluding or hijacking the cellular defense mechanisms. Among the complex pathways induced by the viral infection, autophagy is particularly crucial and is strongly influenced by the action of the non-structural protein 6 (Nsp6) interacting with the endoplasmic reticulum membrane. Importantly, differently from other non-structural proteins, Nsp6 is mutated in the recently emerged Omicron variant, suggesting a possible different role of autophagy. In this contribution we explore, for the first time, the structural properties of Nsp6 thanks to long-timescale molecular dynamics simulations and machine learning analysis, identifying the interaction patterns with the lipid membrane. We also show how the mutation brought by the Omicron variant may indeed modify some of the specific interactions, and more particularly help anchor the viral protein to the lipid bilayer interface.

The viral cycle of SARS-CoV-2 is based on a complex interplay with the cellular machinery, which is mediated by specific proteins eluding or hijacking the cellular defense mechanisms.

Computational Approach for Probing Redox Potential for Iron-Sulfur Clusters in Photosystem I

Photosystem I is a light-driven electron transfer device. Available X-ray crystal structure from Thermosynechococcus elongatus showed that electron transfer pathways consist of two nearly symmetric branches of cofactors converging at the first iron-sulfur cluster F_X, which is followed by two terminal iron-sulfur clusters F_A and F_B. Experiments have shown that F_X has lower oxidation potential than F_A and F_B, which facilitates the electron transfer reaction. Here, we use density functional theory and Multi-Conformer Continuum Electrostatics to explain the differences in the midpoint Em potentials of the F_X, F_A and F_B clusters. Our calculations show that F_X has the lowest oxidation potential compared to F_A and F_B due to strong pairwise electrostatic interactions with surrounding residues. These interactions are shown to be dominated by the bridging sulfurs and cysteine ligands, which may be attributed to the shorter average bond distances between the oxidized Fe ion and ligating sulfurs for F_X compared to F_A and F_B. Moreover, the electrostatic repulsion between the 4Fe-4S clusters and the positive potential of the backbone atoms is lowest for F_X compared to both F_A and F_B. These results agree with the experimental measurements from the redox titrations of low-temperature EPR signals and of room temperature recombination kinetics.

Effect of Oriented Electric Fields on Biologically Relevant Iron-Sulfur Clusters: Tuning Redox Reactivity for Catalysis

Enzyme-based iron-sulfur clusters, exemplified in families such as hydrogenases, nitrogenases, and radical S-adenosylmethionine enzymes, feature in many essential biological processes. The functionality of biological iron-sulfur clusters extends beyond simple electron transfer, relying primarily on the redox activity of the clusters, with a remarkable diversity for different enzymes. The active-site structure and the electrostatic environment in which the cluster resides direct this redox reactivity. Oriented electric fields in enzymatic active sites can be significantly strong, and understanding the extent of their effect on iron-sulfur cluster reactivity can inform first steps toward rationally engineering their reactivity. An extensive systematic density functional theory-based screening approach using OPBE/TZP has afforded a simple electric field-effect representation. The results demonstrate that the orientation of an external electric field of strength 28.8 MV cm^-1 at the center of the cluster can have a significant effect on its relative stability in the order of 35 kJ mol^-1. This shows clear implications for the reactivity of iron-sulfur clusters in enzymes. The results also demonstrate that the orientation of the electric field can alter the most stable broken-symmetry state, which further has implications on the directionality of initiated electron-transfer reactions. These insights open the path for manipulating the enzymatic redox reactivity of iron-sulfur cluster-containing enzymes by rationally engineering oriented electric fields within the enzymes.

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron-Sulfur Cubanes

Iron-sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. To provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III)₄S₄(SMe)₄]^-2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

The Iron Maiden. Cytosolic Aconitase/IRP1 Conformational Transition in the Regulation of Ferritin Translation and Iron Hemostasis

Maintaining iron homeostasis is fundamental for almost all living beings, and its deregulation correlates with severe and debilitating pathologies. The process is made more complicated by the omnipresence of iron and by its role as a fundamental component of a number of crucial metallo proteins. The response to modifications in the amount of the free-iron pool is performed via the inhibition of ferritin translation by sequestering consensus messenger RNA (mRNA) sequences. In turn, this is regulated by the iron-sensitive conformational equilibrium between cytosolic aconitase and IRP1, mediated by the presence of an iron-sulfur cluster. In this contribution, we analyze by full-atom molecular dynamics simulation, the factors leading to both the interaction with mRNA and the conformational transition. Furthermore, the role of the iron-sulfur cluster in driving the conformational transition is assessed by obtaining the related free energy profile via enhanced sampling molecular dynamics simulations.

Two-Channel Model for Electron Transfer in a Dye-Catalyst-Dye Supramolecular Complex for Photocatalytic Water Splitting

To improve the performance of dye-sensitized photoelectrochemical cell (DS-PEC) devices for splitting water, the tailoring of the photocatalytic four-photon water oxidation half-reaction represents a principle challenge of fundamental significance. In this study, a Ru-based water oxidation catalyst (WOC) covalently bound to two 2,6-diethoxy-1,4,5,8-diimide-naphthalene (NDI) dye functionalities provides comparable driving forces and channels for electron transfer. Constrained ab initio molecular dynamics simulations are performed to investigate the photocatalytic cycle of this two-channel model for photocatalytic water splitting. The introduction of a second light-harvesting dye in the Ru-based dye-WOC-dye supramolecular complex enables two separate parallel electron-transfer channels, leading to a five-step catalytic cycle with three intermediates and two doubly oxidized states. The total spin S=1 is conserved during the catalytic process and the system with opposite spin on the oxidized NDI proceeds from the Ru=O intermediate to the final Ru-O2 intermediate with a triplet molecular 3 O2 ligand that is eventually released into the environment. The in-depth insight into the proposed photocatalytic cycle of the two-channel model provides a strategy for the development of novel high-efficiency supramolecular complexes for DS-PEC devices with buildup and conservation of spin multiplicity along the reaction coordinate as a design principle.

From Synthetic to Biological Fe4 S4 Complexes: Redox Properties Correlated to Function of Radical S-Adenosylmethionine Enzymes

By employing the computational protocol for calculation of reduction potentials of the Fe₄ S₄ -containing species validated using a representative series of well-defined synthetic complexes, we focused on redox properties of two prototypical radical SAM enzymes to reveal how they transform SAM into the reactive 5'-deoxyadenosyl radical, and how they tune this radical for its proper biological function. We found the reduction potential of SAM is indeed elevated by 0.3-0.4 V upon coordination to Fe₄ S₄ , which was previously speculated in the literature. This makes a generation of 5'-deoxyadenosyl radical from SAM less endergonic (by ca. 7-9 kcal mol^-1 ) and hence more feasible in both enzymes as compared to the identical process in water. Furthermore, our calculations indicate that the enzyme-bound 5'-deoxyadenosyl radical has a significantly lower reduction potential than in referential aqueous solution, which may help the enzymes to suppress potential side redox reactions and simultaneously elevate its proton-philic character, which may, in turn, promote the radical hydrogen-atom abstraction ability.

Development of parameters compatible with the CHARMM36 force field for [Fe4S4]2+ clusters and molecular dynamics simulations of adenosine-5'-phosphosulfate reductase in GROMACS 2019

DFT calculations were used to obtain parameters compatible with the CHARMM36 force field for iron-sulfur clusters (Fe-S) of the type [Fe₄S₄]²⁺ that are coordinated to dissimilatory adenosine-5'-phosphosulfate reductase (APSrAB). Classical molecular dynamics (MD) simulations were performed on two APSrAB systems to validate the parameters and verify the stability of the studied systems. The time analysis of the parameters inserted into the force field was in reasonable agreement with the experimental X-ray diffraction data. The analysis of the time evolution of the studied systems indicated that these systems and, in particular, the clusters in their respective cavities had a good stability and were in agreement with what was observed in previous works. The parameters obtained provide the basis for the study of APSrAB as well as other systems that contain [Fe₄S₄]²⁺ through the CHARMM36 force field.

PRE-driven protein NMR structures: an alternative approach in highly paramagnetic systems

Metalloproteins play key roles across biology, and knowledge of their structure is essential to understand their physiological role. For those metalloproteins containing paramagnetic states, the enhanced relaxation caused by the unpaired electrons often makes signal detection unfeasible near the metal center, precluding adequate structural characterization right where it is more biochemically relevant. Here, we report a protein structure determination by NMR where two different sets of restraints, one containing Nuclear Overhauser Enhancements (NOEs) and another containing Paramagnetic Relaxation Enhancements (PREs), are used separately and eventually together. The protein PioC from Rhodopseudomonas palustris TIE-1 is a High Potential Iron-Sulfur Protein (HiPIP) where the [4Fe-4S] cluster is paramagnetic in both oxidation states at room temperature providing the source of PREs used as alternative distance restraints. Comparison of the family of structures obtained using NOEs only, PREs only, and the combination of both reveals that the pairwise root-mean-square deviation (RMSD) between them is similar and comparable with the precision within each family. This demonstrates that, under favorable conditions in terms of protein size and paramagnetic effects, PREs can efficiently complement and eventually replace NOEs for the structural characterization of small paramagnetic metalloproteins and de novo-designed metalloproteins by NMR. DATABASES: The 20 conformers with the lowest target function constituting the final family obtained using the full set of NMR restraints were deposited to the Protein Data Bank (PDB ID: 6XYV). The 20 conformers with the lowest target function obtained using NOEs only (PDB ID: 7A58) and PREs only (PDB ID: 7A4L) were also deposited to the Protein Data Bank. The chemical shift assignments were deposited to the BMRB (code 34487).

Monte Carlo simulation and thermodynamic integration applied to protein charge transfer

We introduce a combination of Monte Carlo simulation and thermodynamic integration methods to address a model problem in free energy computations, electron transfer in proteins. The feasibility of this approach is tested using the ferredoxin protein from Clostridium acidurici. The results are compared to numerical solutions of the Poisson-Boltzmann equation and data from recent molecular dynamics simulations on charge transfer in a protein complex, the NrfHA nitrite reductase of Desulfovibrio vulgaris. Despite the conceptual and computational simplicity of the Monte Carlo approach, the data agree well with those obtained by other methods. A link to experiments is established via the cytochrome subunit of the bacterial photosynthetic reaction center of Rhodopseudomonas viridis.

Mutation effects on charge transport through the p58c iron-sulfur protein

Growing experimental evidence indicates that iron–sulfur proteins play key roles in DNA repair and replication. In particular, charge transport between [Fe₄S₄] clusters, mediated by proteins and DNA, may convey signals to coordinate enzyme action. Human primase is a well studied [Fe₄S₄] protein, and its p58c domain (which contains an [Fe₄S₄] cluster) plays a role in the initiation of DNA replication. The Y345C mutation in p58c is linked to gastric tumors and may influence the protein-mediated charge transport. The complexity of protein–DNA systems, and the intricate electronic structure of [Fe₄S₄] clusters, have impeded progress into understanding functional charge transport in these systems. In this study, we built force fields to describe the high potential [Fe₄S₄] cluster in both oxidation states. The parameterization is compatible with AMBER force fields and enabled well-balanced molecular dynamics simulations of the p58c–RNA/DNA complex relevant to the initiation of DNA replication. Using the molecular mechanics Poisson–Boltzmann and surface area solvation method on the molecular dynamics trajectories, we find that the p58c mutation induces a modest change in the p58c–duplex binding free energy in agreement with recent experiments. Through kinetic modeling and analysis, we identify key features of the main charge transport pathways in p58c. In particular, we find that the Y345C mutation partially changes the composition and frequency of the most efficient (and potentially relevant to the biological function) charge transport pathways between the [Fe₄S₄] cluster and the duplex. Moreover, our approach sets the stage for a deeper understanding of functional charge transfer in [Fe₄S₄] protein–DNA complexes.

Functional electron transfer between the [Fe₄S₄] cluster and the nucleic acid is impacted by a Y345C mutation in the p58c subunit of human primase.

Radical Stabilization Energies for Enzyme Engineering: Tackling the Substrate Scope of the Radical Enzyme QueE

Experimental assessment of catalytic reaction mechanisms and profiles of radical enzymes can be severely challenging due to the reactive nature of the intermediates and sensitivity of cofactors such as iron-sulfur clusters. Here, we present an enzyme-directed computational methodology for the assessment of thermodynamic reaction profiles and screening for radical stabilization energies (RSEs) for the assessment of catalytic turnovers in radical enzymes. We have applied this new screening method to the radical S-adenosylmethione enzyme 7-carboxy-7-deazaguanine synthase (QueE), following a detailed molecular dynamics (MD) analysis that clarifies the role of both specific enzyme residues and bound Mg²⁺, Ca²⁺, or Na⁺. The MD simulations provided the basis for a statistical approach to sample different conformational outcomes. RSE calculation at the M06-2X/6-31+G* level of theory provided the most computationally cost-effective assessment of enzyme-based energies, facilitated by an initial triage using semiempirical methods. The impact of intermolecular interactions on RSE was clearly established, and application to the assessment of potential alternative substrates (focusing on radical clock type rearrangements) proposes a selection of carbon-substituted analogues that would react to afford cyclopropylcarbinyl radical intermediates as candidates for catalytic turnover by QueE.

Bona fide choline monoxygenases evolved in Amaranthaceae plants from oxygenases of unknown function: Evidence from phylogenetics, homology modeling and docking studies

Few land plants can synthesize and accumulate the osmoprotectant glycine betaine (GB) even though this metabolic trait has major adaptive importance given the prevalence of drought, hypersaline soils or cold. GB is synthesized from choline in two reactions catalyzed by choline monooxygenases (CMOs) and enzymes of the family 10 of aldehyde dehydrogenases (ALDH10s) that gained betaine aldehyde dehydrogenase activity (BADH). Homolog genes encoding CMO and ALDH10 enzymes are present in all known land plant genomes, but since GB-non-accumulators plants lack the BADH-type ALDH10 isozyme, they would be expected to also lack the CMO activity to avoid accumulation of the toxic betaine aldehyde. To explore CMOs substrate specificity, we performed amino acid sequence alignments, phylogenetic analysis, homology modeling and docking simulations. We found that plant CMOs form a monophyletic subfamily within the Rieske/mononuclear non-heme oxygenases family with two clades: CMO1 and CMO2, the latter diverging from CMO1 after gene duplication. CMO1 enzymes are present in all plants; CMO2s only in the Amaranthaceae high-GB-accumulators plants. CMO2s, and particularly their mononuclear non-heme iron domain where the active site is located, evolved at a faster rate than CMO1s, which suggests positive selection. The homology model and docking simulations of the spinach CMO2 enzyme showed at the active site three aromatic residues forming a box with which the trimethylammonium group of choline could interact through cation-π interactions, and a glutamate, which also may interact with the trimethylammonium group through a charge-charge interaction. The aromatic box and the carboxylate have been shown to be critical for choline binding in other proteins. Interestingly, these residues are conserved in CMO2 proteins but not in CMO1 proteins, where two of these aromatic residues are leucine and the glutamate is asparagine. These findings reinforce our proposal that the CMO1s physiological substrate is not choline but a still unknown metabolite.

Charge transfer through a fragment of the respiratory complex I and its regulation: an atomistic simulation approach

We simulate electron transfer within a fragment of the NADH:ubiquinone oxidoreductase (respiratory complex I) of the hyperthermophilic bacterium Aquifex aeolicus. We apply molecular dynamics simulations, thermodynamic integration, and a thermodynamic network least squares analysis to compute two key parameters of Marcus' theory of charge transfer, the thermodynamic driving force and the reorganization energy. Intramolecular contributions to the Gibbs free energy differences of electron and hydrogen transfer processes, ΔG, are accessed by calibrating against experimental redox titration data. This approach permits the computation of the interactions between the species NAD+, FMNH2, N1a-, and N3-, and the construction of a free energy surface for the flow of electrons within the fragment. We find NAD+ to be a strong candidate for the regulation of charge transfer.

Why the Flavin Adenine Dinucleotide (FAD) Cofactor Needs To Be Covalently Linked to Complex II of the Electron-Transport Chain for the Conversion of FADH2 into FAD

A covalently bound flavin cofactor is predominant in the succinate‐ubiquinone oxidoreductase (SQR; Complex II), an essential component of aerobic electron transport, and in the menaquinol‐fumarate oxidoreductase (QFR), the anaerobic counterpart, although it is only present in approximately 10 % of the known flavoenzymes. This work investigates the role of this 8α‐N3‐histidyl linkage between the flavin adenine dinucleotide (FAD) cofactor and the respiratory Complex II. After parameterization with DFT calculations, classical molecular‐dynamics simulations and quantum‐mechanics calculations for Complex II:FAD and Complex II:FADH₂, with and without the covalent bond, were performed. It was observed that the covalent bond is essential for the active‐center arrangement of the FADH₂/FAD cofactor. Removal of this bond causes a displacement of the isoalloxazine group, which influences interactions with the protein, flavin solvation, and possible proton‐transfer pathways. Specifically, for the noncovalently bound FADH₂ cofactor, the N1 atom moves away from the His‐A365 and His‐A254 residues and the N5 atom moves away from the glutamine‐62A residue. Both of the histidine and glutamine residues interact with a chain of water molecules that cross the enzyme, which is most likely involved in proton transfer. Breaking this chain of water molecules could thereby compromise proton transfer across the two active sites of Complex II.

Rubidium chalcogenido diferrates(III) containing dimers [Fe2 Q 6]6− of edge-sharing tetrahedra (Q=O, S, Se)

The two isotypic rubidium chalcogenido diferrates Rb12[Fe2 Q 6](Q 2)3 (Q=S/Se), which both form needles with green-metallic lustre, were synthesized from Rb2S, elemental iron, rubidium and sulfur (Q=S) or from the pure elements (Q=Se) at maximum temperatures of 500–800°C. Their triclinic crystal structures were determined by means of X-ray single crystal data (space group P1̅, a=863.960(10)/903.2(3), b=942.790(10)/982.1(3), c=1182.70(2)/1227.4(4) pm, α=77.4740(10)/77.262(6), β=71.5250(10)/71.462(6), γ=63.7560(10)/63.462(5)°, Z=1, R1=0.0308/0.0658 for Q=S/Se). The structures contain isolated dinuclear anions [FeIII 2 Q 6]6− composed of two edge-sharing [FeQ 4] tetrahedra (dFe −Q =223.4–232.3/236.2–244.8 pm), which are also found in the two polymorphs of the pure alkali diferrates Rb6[Fe2 Q 6]. The diferrate ions are arranged in layers running in the a/b plane around z=0. Inbetween (around z ≈ 1 2 $z \approx {1 \over 2}$ ), two crystallographically different disulfide/diselenide ions Q 2 2 − $Q_2^{2 - }$ (dQ −Q =211.1–213.4/237.9–241.1 pm), which are arranged in slightly puckered 36 nets, are intercalated. The intra-anionic distances and angles, the Rb coordination numbers and the molar volumes of these two ‘double-salts’ are in accordance with their corresponding reference compounds, Rb6[Fe2 Q 6] and Rb2 Q 2. In addition, the two polymorphs of Rb6[Fe2Se6], which are both isotypic with the sulfido analogous (Cs6[Ga2Se6]-type, monoclinic, space group P21 /c, a=827.84(5), b=1329.51(7), c=1074.10(6) pm, β=127.130(5)°, R1=0.0443 and Ba6[Al2Sb6]-type, orthorhombic, space group Cmce, a=1963.70(3), b=718.98(3), c=1348.40(7) pm, R1=0.0264) were prepared and characterized to complete the series of alkali diferrates(III) with oxido, sulfido and selenido ligands. The electronic band structures of the three Rb salts Rb6[Fe2 Q 6], which have been calculated within the GGA+U approach applying an AFM spin ordering in the dimers and appropriate Hubbard parameters, allow a comparison of the chemical bonding characteristics (e.g. covalency) and the magnetic properties (magnetic moments) within the series of chalcogenido ligands. An analysis of the spin densities enables a comparative consideration of the mechanisms crucial for the magnetic ordering in chalcogenido ferrates. Ultimately, the electronic structure of the new compound Rb12[Fe2S6](S2)3 nicely compares with those of the S2-free reference compound Rb6[Fe2S6].