Charge Transfer at the Qo-Site of the Cytochrome bc1 Complex Leads to Superoxide Production

Long-Time Oxygen and Superoxide Localization in Arabidopsis thaliana Cryptochrome

Cryptochromes are proteins that are highly conserved across species and in many instances bind the flavin adenine dinucleotide (FAD) cofactor within their photolyase-homology region (PHR) domain. The FAD cofactor has multiple redox states that help catalyze reactions, and absorbs photons at about 450 nm, a feature linked to the light-related functions of cryptochrome proteins. Reactive oxygen species (ROS) are produced from redox reactions involving molecular oxygen and are involved in a myriad of biological processes. Superoxide O2•- is an exemplary ROS that may be formed through electron transfer from FAD to O2, generating an electron radical pair. Although the formation of a superoxide-FAD radical pair has been speculated, it is still unclear if the required process steps could be realized in cryptochrome. Here, we present results from molecular dynamics (MD) simulations of oxygen interacting with the PHR domain of Arabidopsis thaliana cryptochrome 1 (AtCRY1). Using MD simulation trajectories, oxygen binding locations are characterized through both the O2-FAD intermolecular distance and the local protein environment. Oxygen unbinding times are characterized through replica simulations of the bound oxygen. Simulations reveal that oxygen molecules can localize at certain sites within the cryptochrome protein for tens of nanoseconds, and superoxide molecules can localize for significantly longer. This relatively long-duration molecule binding suggests the possibility of an electron-transfer reaction leading to superoxide formation. Estimates of electron-transfer rates using the Marcus theory are performed for the identified potential binding sites. Molecular oxygen binding results are compared with recent results demonstrating long-time oxygen binding within the electron-transfer flavoprotein (ETF), another FAD binding protein.

Mitochondrial Cristae Morphology Reflecting Metabolism, Superoxide Formation, Redox Homeostasis, and Pathology

Significance: Mitochondrial (mt) reticulum network in the cell possesses amazing ultramorphology of parallel lamellar cristae, formed by the invaginated inner mitochondrial membrane. Its non-invaginated part, the inner boundary membrane (IBM) forms a cylindrical sandwich with the outer mitochondrial membrane (OMM). Crista membranes (CMs) meet IBM at crista junctions (CJs) of mt cristae organizing system (MICOS) complexes connected to OMM sorting and assembly machinery (SAM). Cristae dimensions, shape, and CJs have characteristic patterns for different metabolic regimes, physiological and pathological situations. Recent Advances: Cristae-shaping proteins were characterized, namely rows of ATP-synthase dimers forming the crista lamella edges, MICOS subunits, optic atrophy 1 (OPA1) isoforms and mitochondrial genome maintenance 1 (MGM1) filaments, prohibitins, and others. Detailed cristae ultramorphology changes were imaged by focused-ion beam/scanning electron microscopy. Dynamics of crista lamellae and mobile CJs were demonstrated by nanoscopy in living cells. With tBID-induced apoptosis a single entirely fused cristae reticulum was observed in a mitochondrial spheroid. Critical Issues: The mobility and composition of MICOS, OPA1, and ATP-synthase dimeric rows regulated by post-translational modifications might be exclusively responsible for cristae morphology changes, but ion fluxes across CM and resulting osmotic forces might be also involved. Inevitably, cristae ultramorphology should reflect also mitochondrial redox homeostasis, but details are unknown. Disordered cristae typically reflect higher superoxide formation. Future Directions: To link redox homeostasis to cristae ultramorphology and define markers, recent progress will help in uncovering mechanisms involved in proton-coupled electron transfer via the respiratory chain and in regulation of cristae architecture, leading to structural determination of superoxide formation sites and cristae ultramorphology changes in diseases. Antioxid. Redox Signal. 39, 635-683.

Charge Regulation in a Rieske Proton Pump Pinpoints Zero, One, and Two Proton-Coupled Electron Transfer

The degree to which redox-driven proton pumps regulate net charge during electron transfer (ΔZET) remains undetermined due to difficulties in measuring the net charge of solvated proteins. Values of ΔZET can reflect reorganization energies or redox potentials associated with ET and can be used to distinguish ET from proton(s)-coupled electron transfer (PCET). Here, we synthesized protein "charge ladders" of a Rieske [2Fe-2S] subunit from Thermus thermophilus (truncTtRp) and made 120 electrostatic measurements of ΔZET across pH. Across pH 5-10, truncTtRp is suspected of transitioning from ET to PCET, and then to two proton-coupled ET (2PCET). Upon reduction, we found that truncTtRp became more negative at pH 6.0 by one unit (ΔZET = -1.01 ± 0.14), consistent with single ET; was isoelectric at pH 8.8 (ΔZET = -0.01 ± 0.45), consistent with PCET; and became more positive at pH 10.6 (ΔZET = +1.37 ± 0.60), consistent with 2PCET. These ΔZET values are attributed to protonation of H154 and H134. Across pH, redox potentials of TtRp (measured previously) correlated with protonation energies of H154 and H134 and ΔZET for truncTtRp, supporting a discrete proton pumping mechanism for Rieske proteins at the Fe-coordinating histidines.

The mitochondrial respiratory chain

We present a brief review of the mitochondrial respiratory chain with emphasis on complexes I, III and IV, which contribute to the generation of protonmotive force across the inner mitochondrial membrane, and drive the synthesis of ATP by the process called oxidative phosphorylation. The basic structural and functional details of these complexes are discussed. In addition, we briefly review the information on the so-called supercomplexes, aggregates of complexes I-IV, and summarize basic physiological aspects of cell respiration.

Long-Time Oxygen Localization in Electron Transfer Flavoprotein

Reactive oxygen species (ROS) exert a wide range of biological effects from beneficial regulatory function to deleterious oxidative stress. The electron transfer flavoprotein (ETF) is ubiquitous to life and is associated with aerobic metabolism and ROS production due to its location in the mitochondria. Quantifying oxygen localization within the ETF complex is critical for understanding the potential for electron transfer and radical pair formation between flavin adenine dinucleotide (FAD) cofactor and superoxide during ROS formation. Our study employed all-atom molecular dynamics simulations and identified several novel, long-lived oxygen binding sites within the ETF complex that appear near the FAD cofactor. Site locations, the local electrostatic environment, and characteristic oxygen binding times for each site were evaluated to establish factors that may lead to possible charge transfer reactions and superoxide formation within the ETF complex. The study revealed that some oxygen binding sites are naturally linked to protein domain features, suggesting opportunities to engineer and control ROS production and subsequent dynamics.

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.

Introducing VIKING: A Novel Online Platform for Multiscale Modeling

Various biochemical and biophysical processes, occurring on multiple time and length scales, can nowadays be studied using specialized software packages on supercomputer clusters. The complexity of such simulations often requires application of different methods in a single study and strong computational expertise. We have developed VIKING, a convenient web platform for carrying out multiscale computations on supercomputers. VIKING allows combining methods in standardized workflows, making complex simulations accessible to a broader biochemical and biophysical society.

Molecular Oxygen Binding in the Mitochondrial Electron Transfer Flavoprotein

Reactive oxygen species such as superoxide are potentially harmful byproducts of the aerobic metabolism in the inner mitochondrial membrane, and complexes I, II, III of the electron transport chain have been identified as primary sources. The mitochondrial fatty acid b-oxidation pathway may also play a yet uncharacterized role in reactive oxygen species generation, apparently at the level of the electron transfer flavoprotein:ubiquinone oxidoreductase (ETF:QO) and/or its redox partner electron-transfer flavoprotein (ETF). These enzymes comprise a key pathway through which electrons are sequentially shuttled from several dehydrogenases to the respiratory chain. The exact mechanisms of superoxide production have not been fully established, but a crucial starting point would be the binding of molecular oxygen within one of the protein complexes. The present investigation offers a comprehensive computational approach for the determination of binding modes and characteristic binding times of small molecules inside proteins, which is then used to reveal several O₂ binding sites near the flavin adenine dinucleotide cofactor of the ETF enzyme. The binding sites are further characterized to extract the necessary parameters for further studies of possible electron transfer between flavin and O₂ leading to radical pair formation and possible superoxide production.

Towards predicting intracellular radiofrequency radiation effects

Recent experiments have reported an effect of weak radiofrequency magnetic fields in the MHz-range on the concentrations of reactive oxygen species (ROS) in living cells. Since the energy that could possibly be deposited by the radiation is orders of magnitude smaller than the energy of molecular thermal motion, it was suggested that the effect was caused by the interaction of RF magnetic fields with transient radical pairs within the cells, affecting the ROS formation rates through the radical pair mechanism. It is, however, at present not entirely clear how to predict RF magnetic field effects at certain field frequency and intensity in nanoscale biomolecular systems. We suggest a possible recipe for interpreting the radiofrequency effects in cells by presenting a general workflow for calculation of the reactive perturbations inside a cell as a function of RF magnetic field strength and frequency. To justify the workflow, we discuss the effects of radiofrequency magnetic fields on generic spin systems to particularly illustrate how the reactive radicals could be affected by specific parameters of the experiment. We finally argue that the suggested workflow can be used to predict effects of radiofrequency magnetic fields on radical pairs in biological cells, which is specially important for wireless recharging technologies where one has to know of any harmful effects that exposure to such radiation might cause.

Atomistic Insights into Cryptochrome Interprotein Interactions

It is striking that the mechanism by which birds sense geomagnetic fields during the biannual migration seasons is not entirely understood. A protein believed to be responsible for avian magnetoreception is the flavoprotein cryptochrome (CRY), which fulfills many of the criteria for a magnetic field sensor. Some experiments, however, indicate that magnetoreception in birds may be disturbed by extremely weak radio frequency fields, an effect that likely cannot be described by an isolated CRY protein. An explanation can possibly be delivered if CRY binds to another protein inside a cell that would possess certain biochemical properties, and it is, therefore, important to identify possible intracellular CRY interaction partners. The goal of this study is to investigate a possible interaction between CRY4 and the iron-sulfur-containing assembly protein (ISCA1) from Erithacus rubecula (European robin), which has recently been proposed to be relevant for magnetic field sensing. The interaction between the proteins is established through classical molecular dynamics simulations for several possible protein-docking modes. The analysis of these simulations concludes that the ISCA1 complex and CRY4 are capable of binding; however, the peculiarities of this binding argue strongly against ISCA1 as relevant for magnetoreception.

On Atomistic Models for Molecular Oxygen

Molecular oxygen (O₂) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O₂ permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O₂ is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O₂ molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O₂ available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O₂, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O₂, in apolar solvents.