Multiscale description of avian migration: from chemical compass to behaviour modeling

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Importance of Polarizable Embedding for Absorption Spectrum Calculations of Arabidopsis thaliana Cryptochrome 1

Cryptochromes are essential flavoproteins for circadian rhythms and avian magnetoreception. Flavin adenine dinucleotide (FAD), a chromophore within cryptochromes, absorbs blue light, initiating electron transfer processes that lead to a biological signaling cascade. A key step in this cascade is the formation of the FAD semiquinone radical (FADH•), characterized through a specific red-light absorption. The absorption spectra of FADH• in cryptochromes are, however, significantly different from those recorded for the cofactor in solution, primarily due to protein-induced shifts in the absorption peaks. This study employs a multiscale approach, combining molecular dynamics (MD) simulations with quantum mechanical/molecular mechanical (QM/MM) methodologies, to investigate the influence of protein dynamics on embedded FADH• absorption. We emphasize the role of the protein's polarizable environment in the shaping of the absorption spectrum, crucial for accurate spectral predictions in cryptochromes. Our findings provide valuable insights into the absorption process, advancing our understanding of cryptochrome functioning.

Spin Dynamics of Flavoproteins

This short review reports the surprising phenomenon of nuclear hyperpolarization occurring in chemical reactions, which is called CIDNP (chemically induced dynamic nuclear polarization) or photo-CIDNP if the chemical reaction is light-driven. The phenomenon occurs in both liquid and solid-state, and electron transfer systems, often carrying flavins as electron acceptors, are involved. Here, we explain the physical and chemical properties of flavins, their occurrence in spin-correlated radical pairs (SCRP) and the possible involvement of flavin-carrying SCRPs in animal magneto-reception at earth's magnetic field.

Quantum-mechanical insights into the anisotropic response of the cryptochrome radical pair to a weak magnetic field

Cryptochrome photoreceptors contain a photochemically generated radical pair, which is thought to mediate sensing of the geomagnetic field direction in many living organisms. To gain insight into the response of the cryptochrome to a weak magnetic field, we have studied the quantum-mechanical hyperfine spin states of the radical pair. We identify quantum states responsible for the precise detection of the magnetic field direction, taking into account the strongly axial hyperfine interactions of each radical in the radical pair. The contribution of these states to the formation of the cryptochrome signaling state sharply increases when the magnetic field becomes orthogonal to the hyperfine axis of either radical. Due to such a response, the radical pair may be able to detect the particular field direction normal to the plane containing the hyperfine axes of the radicals.

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.

Computational Reconstruction and Analysis of Structural Models of Avian Cryptochrome 4

A recent study by Xu et al. (Nature, 2021, 594, 535-540) provided strong evidence that cryptochrome 4 (Cry4) is a key protein to endow migratory birds with the magnetic compass sense. The investigation compared the magnetic field response of Cry4 from migratory and nonmigratory bird species and suggested that a difference in magnetic sensitivity could exist. This finding prompted an in-depth investigation into Cry4 protein differences on the structural and dynamic levels. In the present study, the pigeon Cry4 (ClCry4) crystal structure was used to reconstruct the missing avian Cry4 protein structures via homology modeling for carefully selected bird species. The reconstructed Cry4 structure from European robin, Eurasian blackcap, zebra finch, chicken, and pigeon were subsequently simulated dynamically and analyzed. The studied avian Cry4 structures show flexibility in analogous regions pointing to similar activation mechanisms and/or signaling interaction partners. It can be concluded that the experimentally recorded difference in the magnetic field sensitivity of Cry4 from different birds is unlikely to be due to solely intrinsic dynamics of the proteins but requires additional factors that have not yet been identified.

Observations about utilitarian coherence in the avian compass

It is hypothesised that the avian compass relies on spin dynamics in a recombining radical pair. Quantum coherence has been suggested as a resource to this process that nature may utilise to achieve increased compass sensitivity. To date, the true functional role of coherence in these natural systems has remained speculative, lacking insights from sufficiently complex models. Here, we investigate realistically large radical pair models with up to 21 nuclear spins, inspired by the putative magnetosensory protein cryptochrome. By varying relative radical orientations, we reveal correlations of several coherence measures with compass fidelity. Whilst electronic coherence is found to be an ineffective predictor of compass sensitivity, a robust correlation of compass sensitivity and a global coherence measure is established. The results demonstrate the importance of realistic models, and appropriate choice of coherence measure, in elucidating the quantum nature of the avian compass.

Anisotropic magnetic field effects in the re-oxidation of cryptochrome in the presence of scavenger radicals

The avian compass and many other of nature's magnetoreceptive traits are widely ascribed to the protein cryptochrome. There, magnetosensitivity is thought to emerge as the spin dynamics of radicals in the applied magnetic field enters in competition with their recombination. The first and dominant model makes use of a radical pair. However, recent studies have suggested that magnetosensitivity could be markedly enhanced for a radical triad, the primary radical pair of which undergoes a spin-selective recombination reaction with a third radical. Here, we test the practicality of this supposition for the reoxidation reaction of the reduced FAD cofactor in cryptochrome, which has been implicated with light-independent magnetoreception but appears irreconcilable with the classical radical pair mechanism (RPM). Based on the available realistic cryptochrome structures, we predict the magnetosensitivity of radical triad systems comprising the flavin semiquinone, the superoxide, and a tyrosine or ascorbyl scavenger radical. We consider many hyperfine-coupled nuclear spins, the relative orientation and placement of the radicals, their coupling by the electron-electron dipolar interaction, and spin relaxation in the superoxide radical in the limit of instantaneous decoherence, which have not been comprehensively considered before. We demonstrate that these systems can provide superior magnetosensitivity under realistic conditions, with implications for dark-state cryptochrome magnetoreception and other biological magneto- and isotope-sensitive radical recombination reactions.

Network explanations and explanatory directionality

Network explanations raise foundational questions about the nature of scientific explanation. The challenge discussed in this article comes from the fact that network explanations are often thought to be non-causal, i.e. they do not describe the dynamical or mechanistic interactions responsible for some behaviour, instead they appeal to topological properties of network models describing the system. These non-causal features are often thought to be valuable precisely because they do not invoke mechanistic or dynamical interactions and provide insights that are not available through causal explanations. Here, I address a central difficulty facing attempts to move away from causal models of explanation; namely, how to recover the directionality of explanation. Within causal models, the directionality of explanation is identified with the direction of causation. This solution is no longer available once we move to non-causal accounts of explanation. I will suggest a solution to this problem that emphasizes the role of conditions of application. In doing so, I will challenge the idea that sui generis mathematical dependencies are the key to understand non-causal explanations. The upshot is a conceptual account of explanation that accommodates the possibility of non-causal network explanations. It also provides guidance for how to evaluate such explanations. This article is part of the theme issue 'Unifying the essential concepts of biological networks: biological insights and philosophical foundations'.

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.

Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome

Cryptochrome proteins are activated by the absorption of blue light, leading to the formation of radical pairs through electron transfer in the active site. Recent experimental studies have shown that once some of the amino acid residues in the active site of Xenopus laevis cryptochrome DASH are mutated, radical-pair formation is still observed. In this study, we computationally investigate electron-transfer pathways in the X. laevis cryptochrome DASH by extensively equilibrating a previously established homology model using molecular dynamics simulations and then mutating key amino acids involved in the electron transfer. The electron-transfer pathways are then probed by using tight-binding density-functional theory. We report the alternative electron-transfer pathways resolved at the molecular level and, through comparison of amino acid sequences for cryptochromes from different species, we demonstrate that one of these alternative electron-transfer pathways could be general for all cryptochrome DASH proteins.

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.

Modeling MesoBioNano systems with MBN Studio made easy

This paper introduces MesoBioNano (MBN) Studio - a graphical user interface for a popular multiscale simulation package MBN Explorer. MBN Studio has been developed to facilitate setting up and starting MBN Explorer calculations, monitoring their progress and examining the calculation results. It is tailored for any calculations that are supported by MBN Explorer, such as for example the single-point energy calculations, structure optimization, molecular dynamics, and kinetic Monte Carlo simulations. Apart from that MBN Studio has built-in tools allowing the calculation and analysis of specific characteristics that are determined by the output of the simulations, such as the diffusion coefficients of molecular species, melting temperatures and associated heat capacities, radial distribution function; a dedicated modeling plug-in allows constructing molecular systems in a quick and efficient manner. Employing this plug-in, one can easily construct molecular systems of different geometries (e.g., spherical or ellipsoidal nanoparticles, cubic crystalline samples) with various atomic composition. The paper presents the first public release of MBN Studio and provides an overview of its significant capabilities, as well as the reference point for further extensions.

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.

Computational reconstruction reveals a candidate magnetic biocompass to be likely irrelevant for magnetoreception

Birds use the magnetic field of the Earth to navigate during their annual migratory travel. The possible mechanism to explain the biophysics of this compass sense involves electron transfers within the photoreceptive protein cryptochrome. The magnetoreceptive functioning of cryptochromes is supposedly facilitated through an iron rich polymer complex which couples to multiple cryptochromes. The present investigation aims to independently reconstruct this complex and describe its interaction with Drosophila melanogaster cryptochromes. The polymer complex consists of ISCA1 protein monomers with internally bound iron sulphur clusters and simultaneously binds ten cryptochromes. Through molecular dynamics we have analysed the stability of the ISCA1-cryptochrome complex and characterized the interaction at the binding sites between individual cryptochrome and ISCA1. It is found that the cryptochrome binding to the ISCA1 polymer is not uniform and that the binding affinity depends on its placement along the ISCA1 polymer. This finding supports the claim that the individual ISCA1 monomer acts as possible intracellular interaction partner of cryptochrome, but the proposed existence of an elongated ISCA1 polymer with multiple attached cryptochromes appears to be questionable.

Magnetic field effect in natural cryptochrome explored with model compound

Many animals sense the Earth's magnetic-field and use it for navigation. It is proposed that a light-dependent quantum effect in cryptochrome proteins, residing in the retina, allows for such an iron-free spin-chemical compass. The photochemical processes, spin-dynamics and its magnetic field dependence in natural cryptochrome are not fully understood by the in vivo and in vitro studies. For a deeper insight into these biophysical mechanisms in cryptochrome, we had introduced a flavin-tryptophan dyad (F10T). Here we present the magnetic field dependence of 1H photo-CIDNP NMR on F10T and a theoretical model for low-field photo-CIDNP of F10T. This model provides mixing mechanism of energy-levels and spin-dynamics at low magnetic fields. Photo-CIDNP has been observed even at Earth's magnetic field (~0.05 mT). These experiments prove F10T to be an excellent model compound establishing the key mechanism of avian-magnetoreception and provide insight into the optimal behaviour of cryptochrome at Earth's magnetic field.

The sensitivity of a radical pair compass magnetoreceptor can be significantly amplified by radical scavengers

Birds have a remarkable ability to obtain navigational information from the Earth's magnetic field. The primary detection mechanism of this compass sense is uncertain but appears to involve the quantum spin dynamics of radical pairs formed transiently in cryptochrome proteins. We propose here a new version of the current model in which spin-selective recombination of the radical pair is not essential. One of the two radicals is imagined to react with a paramagnetic scavenger via spin-selective electron transfer. By means of simulations of the spin dynamics of cryptochrome-inspired radical pairs, we show that the new scheme offers two clear and important benefits. The sensitivity to a 50 μT magnetic field is greatly enhanced and, unlike the current model, the radicals can be more than 2 nm apart in the magnetoreceptor protein. The latter means that animal cryptochromes that have a tetrad (rather than a triad) of tryptophan electron donors can still be expected to be viable as magnetic compass sensors. Lifting the restriction on the rate of the spin-selective recombination reaction also means that the detrimental effects of inter-radical exchange and dipolar interactions can be minimised by placing the radicals much further apart than in the current model.