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

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.

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.

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.