Parameterization of the prosthetic redox centers of the bacterial cytochrome bc 1 complex for atomistic molecular dynamics simulations

Structural basis of respiratory complex adaptation to cold temperatures

In response to cold, mammals activate brown fat for respiratory-dependent thermogenesis reliant on the electron transport chain. Yet, the structural basis of respiratory complex adaptation upon cold exposure remains elusive. Herein, we combined thermoregulatory physiology and cryoelectron microscopy (cryo-EM) to study endogenous respiratory supercomplexes from mice exposed to different temperatures. A cold-induced conformation of CI:III2 (termed type 2) supercomplex was identified with a ∼25° rotation of CIII2 around its inter-dimer axis, shortening inter-complex Q exchange space, and exhibiting catalytic states that favor electron transfer. Large-scale supercomplex simulations in mitochondrial membranes reveal how lipid-protein arrangements stabilize type 2 complexes to enhance catalytic activity. Together, our cryo-EM studies, multiscale simulations, and biochemical analyses unveil the thermoregulatory mechanisms and dynamics of increased respiratory capacity in brown fat at the structural and energetic level.

Structural basis of respiratory complexes adaptation to cold temperatures

In response to cold, mammals activate brown fat for respiratory-dependent thermogenesis reliant on the electron transport chain (1, 2). Yet, the structural basis of respiratory complex adaptation to cold remains elusive. Herein we combined thermoregulatory physiology and cryo-EM to study endogenous respiratory supercomplexes exposed to different temperatures. A cold-induced conformation of CI:III 2 (termed type 2) was identified with a ∼25° rotation of CIII 2 around its inter-dimer axis, shortening inter-complex Q exchange space, and exhibiting different catalytic states which favor electron transfer. Large-scale supercomplex simulations in lipid membrane reveal how unique lipid-protein arrangements stabilize type 2 complexes to enhance catalytic activity. Together, our cryo-EM studies, multiscale simulations and biochemical analyses unveil the mechanisms and dynamics of respiratory adaptation at the structural and energetic level.

Modulating the pro-apoptotic activity of cytochrome c at a biomimetic electrified interface

Programmed cell death via apoptosis is a natural defence against excessive cell division, crucial for fetal development to maintenance of homeostasis and elimination of precancerous and senescent cells. Here, we demonstrate an electrified liquid biointerface that replicates the molecular machinery of the inner mitochondrial membrane at the onset of apoptosis. By mimicking in vivo cytochrome c (Cyt c) interactions with cell membranes, our platform allows us to modulate the conformational plasticity of the protein by simply varying the electrochemical environment at an aqueous-organic interface. We observe interfacial electron transfer between an organic electron donor decamethylferrocene and O₂, electrocatalyzed by Cyt c. This interfacial reaction requires partial Cyt c unfolding, mimicking Cyt c in vivo peroxidase activity. As proof of concept, we use our electrified liquid biointerface to identify drug molecules, such as bifonazole, that can potentially down-regulate Cyt c and protect against uncontrolled neuronal cell death in neurodegenerative disorders.

Hydrogen bonding rearrangement by a mitochondrial disease mutation in cytochrome bc 1 perturbs heme b H redox potential and spin state

To perform their specific electron-transfer relay functions, hemes commonly adopt low spin states with fine-tuned redox potentials. Understanding molecular elements controlling these properties is crucial for the description of natural proteins and engineering redox-active systems. We describe unusual effects of mitochondrial disease-related mutation in cytochrome bc₁, based on which we identify a dual role of hydrogen bonding to the propionate group of heme b_H. We observe that stabilization of the hydrogen bond in mutant enhances the redox potential but destabilizes the low spin state of oxidized heme. This demonstrates a critical role of the hydrogen bonding, and heme-protein interactions in general, to secure a suitable redox potential and spin state, a notion that might be universal for other heme proteins.

Hemes are common elements of biological redox cofactor chains involved in rapid electron transfer. While the redox properties of hemes and the stability of the spin state are recognized as key determinants of their function, understanding the molecular basis of control of these properties is challenging. Here, benefiting from the effects of one mitochondrial disease–related point mutation in cytochrome b, we identify a dual role of hydrogen bonding (H-bond) to the propionate group of heme b_H of cytochrome bc₁, a common component of energy-conserving systems. We found that replacing conserved glycine with serine in the vicinity of heme b_H caused stabilization of this bond, which not only increased the redox potential of the heme but also induced structural and energetic changes in interactions between Fe ion and axial histidine ligands. The latter led to a reversible spin conversion of the oxidized Fe from 1/2 to 5/2, an effect that potentially reduces the electron transfer rate between the heme and its redox partners. We thus propose that H-bond to the propionate group and heme-protein packing contribute to the fine-tuning of the redox potential of heme and maintaining its proper spin state. A subtle balance is needed between these two contributions: While increasing the H-bond stability raises the heme potential, the extent of increase must be limited to maintain the low spin and diamagnetic form of heme. This principle might apply to other native heme proteins and can be exploited in engineering of artificial heme-containing protein maquettes.

Modeling the Energy Landscape of Side Reactions in the Cytochrome bc1 Complex

Much of the metabolic molecular machinery responsible for energy transduction processes in living organisms revolves around a series of electron and proton transfer processes. The highly redox active enzymes can, however, also pose a risk of unwanted side reactions leading to reactive oxygen species, which are harmful to cells and are a factor in aging and age-related diseases. Using extensive quantum and classical computational modeling, we here show evidence of a particular superoxide production mechanism through stray reactions between molecular oxygen and a semiquinone reaction intermediate bound in the mitochondrial complex III of the electron transport chain, also known as the cytochrome b c 1 complex. Free energy calculations indicate a favorable electron transfer from semiquinone occurring at low rates under normal circumstances. Furthermore, simulations of the product state reveal that superoxide formed at the Q _o -site exclusively leaves the b c 1 complex at the positive side of the membrane and escapes into the intermembrane space of mitochondria, providing a critical clue in further studies of the harmful effects of mitochondrial superoxide production.

Ergodicity-Breaking in Thermal Biological Electron Transfer? Cytochrome C Revisited

It was recently suggested that certain redox proteins operate in an ergodicity-breaking regime to facilitate biological electron transfer (ET). A signature for this is a large variance reorganization free energy (several electronvolts) but a significantly smaller Stokes reorganization free energy due to incomplete protein relaxation on the time scale of the ET event. Here we investigate whether this picture holds for oxidation of cytochrome c in aqueous solution, at various levels of theory including classical molecular dynamics with two additive and one electronically polarizable force field, and QM/MM calculations with the QM region treated by full electrostatic DFT embedding and by the perturbed matrix method. Sampling the protein and energy gap dynamics over more than 250 ns, we find no evidence for ergodicity-breaking effects. In particular, the inclusion of electronic polarizability of the heme group at QM/MM levels did not induce nonergodic effects, contrary to previous reports by Matyushov et al. The well-known problem of overestimation of reorganization free energies with additive force fields is cured when the protein and solvent are treated as electronically polarizable. Ergodicity-breaking effects may occur in other redox proteins, and our results suggest that long simulations, ideally on the ET time scale, with electronically polarizable force fields are required to obtain strong numerical evidence for them.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Mutations at the Qo Site of the Cytochrome bc1 Complex Strongly Affect Oxygen Binding

The homodimeric bc₁ protein complex is embedded in membranes of mitochondria and photosynthetic bacteria, where it transports protons across the membrane to maintain an electrostatic potential used to drive ATP synthesis as part of the respiratory or photosynthetic pathways. The reaction cycle of the bc₁ complex is driven by series of redox processes involving substrate molecules from the membrane, but occasional side reactions between an intermediate semiquinone substrate and molecular oxygen are suspected to be a source of toxic superoxide, which is believed to be a factor in aging. The present investigation employs molecular dynamics simulations to study the effect of mutations in the Q_o binding sites of the bc₁ complex on the ability of oxygen molecules to migrate to and bind at various locations within the complex. It is found that the mutations strongly affect the ability of oxygen to bind at the Q_o sites, and moreover, different behavior of the two monomers of the bc₁ complex is observed. The conformational differences at the Q_o sites of the two monomers are studied in detail and discussed. The anionic form of semiquinone was identified as leading to the greatest opportunity for side reactions with oxygen.

Mechanism of the Primary Charge Transfer Reaction in the Cytochrome bc1 Complex

The bc₁ complex is a critical enzyme for the ATP production in photosynthesis and cellular respiration. Its biochemical function relies on the so-called Q-cycle, which is well established and operates via quinol substrates that bind inside the protein complex. Despite decades of research, the quinol-protein interaction, which initiates the Q-cycle, has not yet been completely described. Furthermore, the initial charge transfer reactions of the Q-cycle lack a physical description. The present investigation utilizes classical molecular dynamics simulations in tandem with quantum density functional theory calculations, to provide a complete and consistent quantitative description of the primary events that occur within the bc₁ complex upon quinol binding. In particular, the electron and proton transfer reactions that trigger the Q-cycle in the bc₁ complex from Rhodobacter capsulatus are studied. The coupled nature of these charge transfer reactions was revealed by obtaining the transition energy path connecting configurations of the Q_o-site prior and after the transfers. The analysis of orbitals and partial charge distribution of the different states of the Q_o-site has further supported the conclusion. Finally, key structural elements of the bc₁ complex that trigger the charge transfer reactions were established, manifesting the importance of the environment in the process, which is furthermore evidenced by free energy calculations.

Atomistic determinants of co-enzyme Q reduction at the Qi-site of the cytochrome bc1 complex

The cytochrome (cyt) bc₁ complex is an integral component of the respiratory electron transfer chain sustaining the energy needs of organisms ranging from humans to bacteria. Due to its ubiquitous role in the energy metabolism, both the oxidation and reduction of the enzyme’s substrate co-enzyme Q has been studied vigorously. Here, this vast amount of data is reassessed after probing the substrate reduction steps at the Q_i-site of the cyt bc₁ complex of Rhodobacter capsulatus using atomistic molecular dynamics simulations. The simulations suggest that the Lys251 side chain could rotate into the Q_i-site to facilitate binding of half-protonated semiquinone – a reaction intermediate that is potentially formed during substrate reduction. At this bent pose, the Lys251 forms a salt bridge with the Asp252, thus making direct proton transfer possible. In the neutral state, the lysine side chain stays close to the conserved binding location of cardiolipin (CL). This back-and-forth motion between the CL and Asp252 indicates that Lys251 functions as a proton shuttle controlled by pH-dependent negative feedback. The CL/K/D switching, which represents a refinement to the previously described CL/K pathway, fine-tunes the proton transfer process. Lastly, the simulation data was used to formulate a mechanism for reducing the substrate at the Q_i-site.

Spontaneous Binding of Molecular Oxygen at the Qo-Site of the bc1 Complex Could Stimulate Superoxide Formation

A key part of the respiratory and photosynthetic pathways is the bc1 protein complex embedded in the inner membrane of mitochondria and the plasma membrane of photosynthetic bacteria. The protein complex pumps protons across the membrane to maintain an electrostatic potential, which is in turn used to drive ATP synthesis. This molecular machinery, however, is suspected to be a source of superoxide, which is toxic to the cell, even in minuscular quantities, and believed to be a factor in aging. Through molecular dynamics simulations, we investigate here the migration of molecular oxygen in the bc1 complex in order to identify possible reaction sites that could lead to superoxide formation. It is found, in particular, that oxygen penetrates spontaneously the Qo binding site of the bc1 complex in the presence of an intermediate semiquinone radical, thus making the Qo-site a strong candidate for being a center of superoxide production.

Duplications of an iron–sulphur tripeptide leads to the formation of a protoferredoxin

Based on UV-Vis, NMR, and EPR spectroscopies and DFT and molecular dynamics calculations, a model prebiotic [2Fe–2S] tripeptide was shown to accept and donate electrons.

Identification of ubiquinol binding motifs at the Qo-site of the cytochrome bc1 complex

Enzymes of the bc₁ complex family power the biosphere through their central role in respiration and photosynthesis. These enzymes couple the oxidation of quinol molecules by cytochrome c to the transfer of protons across the membrane, to generate a proton-motive force that drives ATP synthesis. Key for the function of the bc₁ complex is the initial redox process that involves a bifurcated electron transfer in which the two electrons from a quinol substrate are passed to different electron acceptors in the bc₁ complex. The electron transfer is coupled to proton transfer. The overall mechanism of quinol oxidation by the bc₁ complex is well enough characterized to allow exploration at the atomistic level, but details are still highly controversial. The controversy stems from the uncertain binding motifs of quinol at the so-called Q_o active site of the bc₁ complex. Here we employ a combination of classical all atom molecular dynamics and quantum chemical calculations to reveal the binding modes of quinol at the Q_o-site of the bc₁ complex from Rhodobacter capsulatus. The calculations suggest a novel configuration of amino acid residues responsible for quinol binding and support a mechanism for proton-coupled electron transfer from quinol to iron–sulfur cluster through a bridging hydrogen bond from histidine that stabilizes the reaction complex.

Molecular Dynamics Simulation of Nitrobenzene Dioxygenase Using AMBER Force Field

Molecular dynamics simulation of the oxygenase component of nitrobenzene dioxygenase (NBDO) system, a member of the naphthalene family of Rieske nonheme iron dioxygenases, has been carried out using the AMBER force field combined with a new set of parameters for the description of the mononuclear nonheme iron center and iron-sulfur Rieske cluster. Simulation results provide information on the structure and dynamics of nitrobenzene dioxygenase in an aqueous environment and shed light on specific interactions that occur in its catalytic center. The results suggest that the architecture of the active site is stabilized by key hydrogen bonds, and Asn258 positions the substrate for oxidation. Analysis of protein-water interactions reveal the presence of a network of solvent molecules at the entrance to the active site, which could be of potential catalytic importance.