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

The cytochrome b6f complex: plastoquinol oxidation and regulation of electron transport in chloroplasts

In oxygenic photosynthetic systems, the cytochrome b6f (Cytb6f) complex (plastoquinol:plastocyanin oxidoreductase) is a heart of the hub that provides connectivity between photosystems (PS) II and I. In this review, the structure and function of the Cytb6f complex are briefly outlined, being focused on the mechanisms of a bifurcated (two-electron) oxidation of plastoquinol (PQH2). In plant chloroplasts, under a wide range of experimental conditions (pH and temperature), a diffusion of PQH2 from PSII to the Cytb6f does not limit the intersystem electron transport. The overall rate of PQH2 turnover is determined mainly by the first step of the bifurcated oxidation of PQH2 at the catalytic site Qo, i.e., the reaction of electron transfer from PQH2 to the Fe2S2 cluster of the high-potential Rieske iron-sulfur protein (ISP). This point has been supported by the quantum chemical analysis of PQH2 oxidation within the framework of a model system including the Fe2S2 cluster of the ISP and surrounding amino acids, the low-potential heme b6L, Glu78 and 2,3,5-trimethylbenzoquinol (the tail-less analog of PQH2). Other structure-function relationships and mechanisms of electron transport regulation of oxygenic photosynthesis associated with the Cytb6f complex are briefly outlined: pH-dependent control of the intersystem electron transport and the regulatory balance between the operation of linear and cyclic electron transfer chains.

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.

Tunnel dynamics of quinone derivatives and its coupling to protein conformational rearrangements in respiratory complex I

Respiratory complex I in mitochondria and bacteria catalyzes the transfer of electrons from NADH to quinone (Q). The free energy available from the reaction is used to pump protons and to establish a membrane proton electrochemical gradient, which drives ATP synthesis. Even though several high-resolution structures of complex I have been resolved, how Q reduction is linked with proton pumping, remains unknown. Here, microsecond long molecular dynamics (MD) simulations were performed on Yarrowia lipolytica complex I structures where Q molecules have been resolved in the ~30 Å long Q tunnel. MD simulations of several different redox/protonation states of Q reveal the coupling between the Q dynamics and the restructuring of conserved loops and ion pairs. Oxidized quinone stabilizes towards the N2 FeS cluster, a binding mode not previously described in Yarrowia lipolytica complex I structures. On the other hand, reduced (and protonated) species tend to diffuse towards the Q binding sites closer to the tunnel entrance. Mechanistic and physiological relevance of these results are discussed.

Functional design of bacterial superoxide:quinone oxidoreductase

The superoxide anion - molecular oxygen reduced by a single electron - is produced in large amounts by enzymatic and adventitious reactions. It can perform a range of cellular functions, including bacterial warfare and iron uptake, signalling and host immune response in eukaryotes. However, it also serves as precursor for more deleterious species such as the hydroxyl anion or peroxynitrite and defense mechanisms to neutralize superoxide are important for cellular health. In addition to the soluble proteins superoxide dismutase and superoxide reductase, recently the membrane embedded diheme cytochrome b₅₆₁ (CybB) from E. coli has been proposed to act as a superoxide:quinone oxidoreductase. Here, we confirm superoxide and cellular ubiquinones or menaquinones as natural substrates and show that quinone binding to the enzyme accelerates the reaction with superoxide. The reactivity of the substrates is in accordance with the here determined midpoint potentials of the two b hemes (+48 and -23 mV / NHE). Our data suggest that the enzyme can work near the diffusion limit in the forward direction and can also catalyse the reverse reaction efficiently under physiological conditions. The data is discussed in the context of described cytochrome b₅₆₁ proteins and potential physiological roles of CybB.

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.

Multiscale modeling of stochastic dynamics processes with MBN Explorer

Stochastic dynamics describes processes in complex systems having the probabilistic nature. They can involve very different dynamical systems and occur on very different temporal and spatial scale. This paper discusses the concept of stochastic dynamics and its implementation in the popular program MBN Explorer. Stochastic dynamics in MBN Explorer relies on the Monte Carlo approach and permits simulations of physical, chemical, and biological processes. The paper presents the basic theoretical concepts underlying stochastic dynamics implementation and provides several examples highlighting its applicability to different systems, such as diffusing proteins seeking an anchor point on a cell membrane, deposition of nanoparticles on a surface leading to structures with fractal morphologies, and oscillations of compounds in an autocatalytic reaction. The chosen examples illustrate the diversity of applications that can be modeled by means of stochastic dynamics with MBN Explorer.

Effect of Dipyridamole on Membrane Energization and Energy Transfer in Chromatophores of Rba. sphaeroides

Effect of dipyridamole (DIP) at concentrations up to 1 mM on fluorescent characteristics of light-harvesting complexes LH2 and LH1, as well as on conditions of photosynthetic electron transport chain in the bacterial chromatophores of Rba. sphaeroides was investigated. DIP was found to affect efficiency of energy transfer from the light-harvesting complex LH2 to the LH1-reaction center core complex and to produce the long-wavelength ("red") shift of the absorption band of light-harvesting bacteriochlorophyll molecules in the IR spectral region at 840-900 nm. This shift is associated with the membrane transition to the energized state. It was shown that DIP is able to reduce the photooxidized bacteriochlorophyll of the reaction center, which accelerated electron flow along the electron transport chain, thereby stimulating generation of the transmembrane potential on the chromatophore membrane. The results are important for clarifying possible mechanisms of DIP influence on the activity of membrane-bound functional proteins. In particular, they might be significant for interpreting numerous therapeutic effects of DIP.

Structure and Mechanism of Respiratory III-IV Supercomplexes in Bioenergetic Membranes

In the final steps of energy conservation in aerobic organisms, free energy from electron transfer through the respiratory chain is transduced into a proton electrochemical gradient across a membrane. In mitochondria and many bacteria, reduction of the dioxygen electron acceptor is catalyzed by cytochrome c oxidase (complex IV), which receives electrons from cytochrome bc1 (complex III), via membrane-bound or water-soluble cytochrome c. These complexes function independently, but in many organisms they associate to form supercomplexes. Here, we review the structural features and the functional significance of the nonobligate III2IV1/2 Saccharomyces cerevisiae mitochondrial supercomplex as well as the obligate III2IV2 supercomplex from actinobacteria. The analysis is centered around the Q-cycle of complex III, proton uptake by CytcO, as well as mechanistic and structural solutions to the electronic link between complexes III and IV.

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.

Catalytic Reactions and Energy Conservation in the Cytochrome bc1 and b6f Complexes of Energy-Transducing Membranes

This review focuses on key components of respiratory and photosynthetic energy-transduction systems: the cytochrome bc1 and b6f (Cytbc1/b6f) membranous multisubunit homodimeric complexes. These remarkable molecular machines catalyze electron transfer from membranous quinones to water-soluble electron carriers (such as cytochromes c or plastocyanin), coupling electron flow to proton translocation across the energy-transducing membrane and contributing to the generation of a transmembrane electrochemical potential gradient, which powers cellular metabolism in the majority of living organisms. Cytsbc1/b6f share many similarities but also have significant differences. While decades of research have provided extensive knowledge on these enzymes, several important aspects of their molecular mechanisms remain to be elucidated. We summarize a broad range of structural, mechanistic, and physiological aspects required for function of Cytbc1/b6f, combining textbook fundamentals with new intriguing concepts that have emerged from more recent studies. The discussion covers but is not limited to (i) mechanisms of energy-conserving bifurcation of electron pathway and energy-wasting superoxide generation at the quinol oxidation site, (ii) the mechanism by which semiquinone is stabilized at the quinone reduction site, (iii) interactions with substrates and specific inhibitors, (iv) intermonomer electron transfer and the role of a dimeric complex, and (v) higher levels of organization and regulation that involve Cytsbc1/b6f. In addressing these topics, we point out existing uncertainties and controversies, which, as suggested, will drive further research in this field.

Inhibition Mechanism of Antimalarial Drugs Targeting the Cytochrome bc1 Complex

Plasmodium falciparum (P. falciparum) is the main parasite known to cause malaria in humans. The antimalarial drug atovaquone is known to inhibit the Q_o-site of the cytochrome bc₁ complex of P. falciparum, which ultimately blocks ATP synthesis, leading to cell death. Through the years, mutations of the P. falciparum cytochrome bc₁ complex, causing resistance to atovaquone, have emerged. The present investigation applies molecular dynamics (MD) simulations to study how the specific mutations Y279S and L282V, known to cause atovaquone resistance in malarial parasites, affect the inhibition mechanism of two known inhibitors. Binding free energy estimates were obtained through free energy perturbation calculations but were unable to confidently resolve the effects of mutations due to the great complexity of the binding environment. Meanwhile, basic mechanistic considerations from the MD simulations provide a detailed characterization of inhibitor binding modes and indicate that the Y279S mutation weakens the natural binding of the inhibitors, while no conclusive effect of the L282V mutation could be observed.

Modeling the effect of ion-induced shock waves and DNA breakage with the reactive CHARMM force field

Ion-induced DNA damage is an important effect underlying ion beam cancer therapy. This article introduces the methodology of modeling DNA damage induced by a shock wave caused by a projectile ion. Specifically it is demonstrated how single- and double strand breaks in a DNA molecule could be described by the reactive CHARMM (rCHARMM) force field implemented in the program MBN Explorer. The entire workflow of performing the shock wave simulations, including obtaining the crucial simulation parameters, is described in seven steps. Two exemplary analyses are provided for a case study simulation serving to: (a) quantify the shock wave propagation and (b) describe the dynamics of formation of DNA breaks. The article concludes by discussing the computational cost of the simulations and revealing the possible maximal computational time for different simulation set-ups.

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.

Bilirubin Links Heme Metabolism to Neuroprotection by Scavenging Superoxide

Bilirubin is one of the most frequently measured metabolites in medicine, yet its physiologic roles remain unclear. Bilirubin can act as an antioxidant in vitro, but whether its redox activity is physiologically relevant is unclear because many other antioxidants are far more abundant in vivo. Here, we report that depleting endogenous bilirubin renders mice hypersensitive to oxidative stress. We find that mice lacking bilirubin are particularly vulnerable to superoxide (O₂^⋅-) over other tested reactive oxidants and electrophiles. Whereas major antioxidants such as glutathione and cysteine exhibit little to no reactivity toward O₂^⋅-, bilirubin readily scavenges O₂^⋅-. We find that bilirubin's redox activity is particularly important in the brain, where it prevents excitotoxicity and neuronal death by scavenging O₂^⋅- during NMDA neurotransmission. Bilirubin's unique redox activity toward O₂^⋅- may underlie a prominent physiologic role despite being significantly less abundant than other endogenous and exogenous antioxidants.

Proliferating tumor cells mimick glucose metabolism of mature human erythrocytes

Mature human erythrocytes are dependent on anerobic glycolysis, i.e. catabolism (oxidation) of one glucose molecule to produce two ATP and two lactate molecules. Proliferating tumor cells mimick mature human erythrocytes to glycolytically generate two ATP molecules. They deliberately avoid or switch off their respiration, i.e. tricarboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS) machinery and consequently dispense with the production of additional 36 ATP molecules from one glucose molecule. This phenomenon is named aerobic glycolysis or Warburg effect. The present review deals with the fate of a glucose molecule after entering a mature human erythrocyte or a proliferating tumor cell and describes why it is useful for a proliferating tumor cell to imitate a mature erythrocyte. Blood consisting of plasma and cellular components (99% of the cells are erythrocytes) may be regarded as a mobile organ, constantly exercising a direct interaction with other organs. Therefore, the use of drugs, which influences the biological activity of erythrocytes, has an immediate effect on the entire organism.

Abbreviations: TCA: tricarboxylic acid cycle; OXPHOS: oxidative phosphorylation; GSH: reduced state of glutathione; NFκB: Nuclear factor of kappa B; PKB (Akt): protein kinase B; NOS: nitric oxide synthase; IgG: immune globulin G; H₂S: hydrogen sulfide; slanDCs: Human 6-sulfo LacNAc-expressing dendritic cells; IL-8: interleukin-8; LPS: lipopolysaccharide; ROS: reactive oxygen species; PPP: pentose phosphate pathway; NADPH: nicotinamide adenine dinucleotide phosphate hydrogen; R5P: ribose-5-phophate; NAD: nicotinamide adenine dinucleotide; FAD: flavin adenine dinucleotide; O₂^●−: superoxide anion; G6P: glucose 6-phosphate; HbO₂: Oxyhemoglobin; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GAP: glyceraldehyde-3-phosphate; 1,3-BPG: 1,3-bis-phosphoglycerate; 2,3-BPG: 2,3-bisphosphoglycerte; PGAM1: phosphoglycerate mutase 1; 3-PG: 3-phosphoglycerate; 2-PG: 2-phosphoglycerate; MIPP1: Multiple inositol polyphosphate phosphatase; mTORC1: mammalian target of rapamycin complex 1; Ru5P: ribulose 5-phosphate; ox-PPP: oxidative branch of pentose phosphate pathway; PGK: phosphoglycerate kinase; IFN-γ: interferon-γ; LDH: lactate dehydrogenase; STAT3: signal transducer and activator of transcription 3; Rheb: Ras homolog enriched in Brain; H₂O₂: hydrogen peroxide; ROOH: lipid peroxide; SOD: superoxide dismutase; MRC: mitochondrial respiratory chain; MbFe²⁺-O₂: methmyoglobin; RNR: ribonucleotide reductase; PRPP: phosphoribosylpyrophosphate; PP_i: pyrophosphate; GSSG: oxidized state of glutathione; non-ox-PPP: non-oxidative branch of pentose phosphate pathway; RPI: ribose-5-phosphate isomerase; RPE: ribulose 5-phosphate 3-epimerase; X5P: xylulose 5-phosphate; TK: transketolase; TA: transaldolase; F6P: fructose-6-phosphate; AR2: aldose reductase 2; SD: sorbitol dehydrogenase; HK: hexokinase; MG: mehtylglyoxal; DHAP: dihydroxyacetone phosphate; TILs: tumor-infiltrating lymphocytes; MCTs: monocarboxylate transporters; pHi: intracellular pH; Hif-1α: hypoxia-induced factor 1; NHE1: sodium/H⁺ (Na⁺/H⁺) antiporter 1; V-ATPase: vacuolar-type proton ATPase; CAIX: carbonic anhydrase; CO₂: carbon dioxide; HCO₃⁻: bicarbonate; NBC: sodium/bicarbonate (Na⁺/HCO₃⁻) symporter; pHe: extracellular pH; GLUT-1: glucose transporter 1; PGK-1: phosphoglycerate kinase 1

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.

Scavenging of superoxide by a membrane-bound superoxide oxidase

Superoxide is a reactive oxygen species produced during aerobic metabolism in mitochondria and prokaryotes. It causes damage to lipids, proteins and DNA and is implicated in cancer, cardiovascular disease, neurodegenerative disorders and aging. As protection, cells express soluble superoxide dismutases, disproportionating superoxide to oxygen and hydrogen peroxide. Here, we describe a membrane-bound enzyme that directly oxidizes superoxide and funnels the sequestered electrons to ubiquinone in a diffusion-limited reaction. Experiments in proteoliposomes and inverted membranes show that the protein is capable of efficiently quenching superoxide generated at the membrane in vitro. The 2.0 Å crystal structure shows an integral membrane di-heme cytochrome b poised for electron transfer from the P-side and proton uptake from the N-side. This suggests that the reaction is electrogenic and contributes to the membrane potential while also conserving energy by reducing the quinone pool. Based on this enzymatic activity, we propose that the enzyme family be denoted superoxide oxidase (SOO).

Molecular mechanism of DNA association with single-stranded DNA binding protein

During DNA replication, the single-stranded DNA binding protein (SSB) wraps single-stranded DNA (ssDNA) with high affinity to protect it from degradation and prevent secondary structure formation. Although SSB binds ssDNA tightly, it can be repositioned along ssDNA to follow the advancement of the replication fork. Using all-atom molecular dynamics simulations, we characterized the molecular mechanism of ssDNA association with SSB. Placed in solution, ssDNA–SSB assemblies were observed to change their structure spontaneously; such structural changes were suppressed in the crystallographic environment. Repeat simulations of the SSB–ssDNA complex under mechanical tension revealed a multitude of possible pathways for ssDNA to come off SSB punctuated by prolonged arrests at reproducible sites at the SSB surface. Ensemble simulations of spontaneous association of short ssDNA fragments with SSB detailed a three-dimensional map of local affinity to DNA; the equilibrium amount of ssDNA bound to SSB was found to depend on the electrolyte concentration but not on the presence of the acidic tips of the SSB tails. Spontaneous formation of ssDNA bulges and their diffusive motion along SSB surface was directly observed in multiple 10-µs-long simulations. Such reptation-like motion was confined by DNA binding to high-affinity spots, suggesting a two-step mechanism for SSB diffusion.

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.

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.