Upper bound on the biological effects of 50/60 Hz magnetic fields mediated by radical pairs

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: March 14, 2019 (This version)
Accepted Manuscript published: February 25, 2019 (Go to version)
Accepted: February 2, 2019
Received: December 6, 2018

1. Of interest
Substrate evaporation drives collective construction in termites

Giulio Facchini, Alann Rathery ... Andrea Perna

Research Article Apr 10, 2024
Further reading

Abstract
Introduction
Methods
Results
Discussion
Data availability
References
Article and author information
Metrics

Abstract

Prolonged exposure to weak (~1 μT) extremely-low-frequency (ELF, 50/60 Hz) magnetic fields has been associated with an increased risk of childhood leukaemia. One of the few biophysical mechanisms that might account for this link involves short-lived chemical reaction intermediates known as radical pairs. In this report, we use spin dynamics simulations to derive an upper bound of 10 parts per million on the effect of a 1 μT ELF magnetic field on the yield of a radical pair reaction. By comparing this figure with the corresponding effects of changes in the strength of the Earth’s magnetic field, we conclude that if exposure to such weak 50/60 Hz magnetic fields has any effect on human biology, and results from a radical pair mechanism, then the risk should be no greater than travelling a few kilometres towards or away from the geomagnetic north or south pole.

https://doi.org/10.7554/eLife.44179.001

Introduction

Residents in developed countries are almost continuously exposed to the 50 or 60 Hz electromagnetic fields generated by electrical appliances and power transmission lines. Although these fields are normally very weak in residential environments (electric component <100 V m⁻¹, magnetic component <1 μT (Swanson and Kheifets, 2006)), these extremely-low-frequency (ELF) fields have been mooted as a potential health hazard (Crumpton, 2005; Crumpton and Collins, 2004; Brocklehurst and McLauchlan, 1996). In relation to childhood leukaemia, ELF magnetic fields have been classified as ‘possibly carcinogenic to humans’ (Ahlbom et al., 2000; Greenland et al., 2000; IARC, 2002; Kheifets et al., 2010). Although there is scant evidence from animal experiments to support a causative link (IARC, 2002; WHO, 2007), numerous in vitro investigations have reported biological effects of ELF fields (reviewed in Juutilainen et al., 2006; Mattsson and Simkó, 2014; Juutilainen et al., 2018). Not many of these observations have been independently replicated (Lacy-Hulbert et al., 1998). To explain any genuine effects, there must exist biophysical mechanisms by which weak ELF magnetic fields interact with molecules in living organisms. More than 10 years ago, Swanson and Kheifets (Swanson and Kheifets, 2006) and Adair (Adair, 1991; Adair, 1999; Adair, 2000) assessed a range of potential mechanisms and concluded that none was likely to have biological consequences at magnetic flux densities of order 1 μT.

One of the mechanisms considered was the radical pair mechanism (Brocklehurst and McLauchlan, 1996). Unlike some of the others, this is unquestionably a genuine phenomenon supported by hundreds of laboratory studies of (mostly non-biological) chemical reactions in magnetic fields ranging from ~1 mT up to several Tesla (Steiner and Ulrich, 1989; Brocklehurst, 2002; Rodgers, 2009; Jones, 2016; Hore, 2012; Scaiano et al., 1994a). Radical pairs are short-lived reaction intermediates, typically formed in highly non-equilibrium electron-spin states. Their coherent spin dynamics, combined with spin-dependent reactivity, allow magnetic interactions a million times smaller than the thermal energy, k_BT, to have measurable, reproducible and predictable effects on chemical reaction yields. In the last few years, interest in the biological significance of <1 mT static magnetic fields has been kindled by new insights into the biophysical mechanism of the avian magnetic compass sense (Mouritsen, 2018; Hore and Mouritsen, 2016; Rodgers and Hore, 2009; Nordmann et al., 2017; Schulten et al., 1978). It now seems likely that migratory songbirds detect the direction of the geomagnetic field by means of photo-induced radical pairs generated within cryptochrome proteins located in their retinas (Hore and Mouritsen, 2016; Ritz et al., 2000; Liedvogel and Mouritsen, 2010). The growing support for this hypothesis includes an experimental demonstration that a radical pair reaction can be influenced by a static magnetic field of the same strength as the geomagnetic field (~50 μT) (Maeda et al., 2008). Whether this hypothesis proves to be correct or not, it is important to distinguish between functional responses, such as magnetoreception, which presumably have been optimised by evolutionary pressure, and fortuitous, non-functional magnetic effects which could in principle be damaging. Since humans seem to navigate without the benefit of in-built magnetic sensors, we focus here exclusively on the latter. However, one cannot exclude the possibility that other biological responses to the Earth’s magnetic field, for example at the cellular level, could have been useful at an early stage in human evolution and might still persist.

Prompted by their putative role in magnetoreception, cryptochromes have become the focus of recent discussions of fortuitous radical pair effects in biology (Bounds and Kuster, 2015; Lagroye et al., 2011; Vanderstraeten et al., 2015; Close, 2012; Close, 2014a; Close, 2014b; Landler and Keays, 2018; Krylov, 2017; Agliassa et al., 2018). Juutilainen et al., for example, have proposed a hypothesis to explain the link between environmental ELF magnetic fields and childhood leukaemia (‘magnetocarcinogenesis’) (Juutilainen et al., 2018). Because cryptochromes are key components of the circadian clock (Chaves et al., 2011; Kelleher et al., 2014; Hastings et al., 2014), and circadian systems are closely coupled to the regulation of DNA damage responses and defence against reactive oxygen species, it is possible that magnetic field effects on radical reactions in cryptochromes could lead to circadian dysregulation, genomic instability and ultimately cancer (Juutilainen et al., 2018).

Our purpose here is to extend Adair’s 1999 discussion (Adair, 1999) of radical pair magnetic field effects by modelling the spin dynamics of cryptochrome-based radical pairs using more advanced theoretical methods. The goal is to derive a realistic order-of-magnitude estimate of the largest likely effect of a 1 μT ELF field in the presence of the Earth’s magnetic field. This upper bound is compared with the predicted effects of modest changes in the strength of the geomagnetic field such as would be experienced by travelling around the globe and with the small changes in body temperature that occur naturally in healthy humans.

Methods

Outline

Request a detailed protocol

Our aim is to determine the change in the yield of a radical pair reaction caused by a 1 μT ELF field superimposed on the geomagnetic field. We make no attempt to link this change to any specific biological process; rather we seek to estimate the maximum primary magnetic field effect under chemically and physically plausible conditions. Calculations are based on the [FAD^•− TrpH^•+] radical pair that accounts for the observed effects of static magnetic fields on the photochemistry of purified cryptochromes (Hore and Mouritsen, 2016; Maeda et al., 2012; Sheppard et al., 2017; Kattnig et al., 2016a). This species is formed by the transfer of an electron from a tryptophan residue (TrpH) in the protein to the photo-excited, non-covalently bound, flavin adenine dinucleotide (FAD) chromophore, Figure 1(a). We consider a simplified spin system comprising the two unpaired electron spins, one on each radical, coupled to three nitrogen nuclei (¹⁴N, spin quantum number $I$ = 1) chosen for their large isotropic hyperfine coupling constants: N5 and N10 in FAD^•− ( $a$ = 523 μT and 189 μT) and N1 in TrpH^•+ ( $a$ = 322 μT), calculated using density functional theory (Lee et al., 2014). The radicals were assumed to have g-values equal to the free electron g-value, 2.0023. At the magnetic field strengths of interest here the difference in the Zeeman interactions of FAD^•− and TrpH^•+ is entirely negligible. We exclude the anisotropic components of the hyperfine interactions which are only relevant when the radicals are immobilised and aligned, as in the case of a magnetic compass sensor (Schulten et al., 1978). These hyperfine interactions condition the response of the radical pair to both static and ELF magnetic fields.

Figure 1

Download asset Open asset

The radical pair mechanism.

(a) Structures and atom numbers schemes for the FAD and TrpH molecules from which magnetically sensitive radical pairs are formed in purified cryptochromes. (b) Simple reaction scheme for a singlet-born radical pair able to react spin-selectively to form singlet and triplet reaction products. The red/blue arrows represent the coherent interconversion of the two forms of the radical pair. The reaction scheme is a simplified version of the cryptochrome photocycle (Hore and Mouritsen, 2016; Maeda et al., 2012). (c) Schematic dependence of the amounts of triplet radical pair and triplet product present as a function of time (in arbitrary units) after the formation of the radical pair in a singlet state. The triplet yield, $Φ_{T}$ , is the amount of triplet product formed once all radical pairs have reacted.

https://doi.org/10.7554/eLife.44179.002

The chemistry of this radical pair was modelled in a simplified fashion by means of competing, spin-selective, singlet and triplet reaction channels (Figure 1(b)). For the present purposes, the precise nature of these reactions is immaterial but to avoid unnecessary abstraction we briefly describe the reaction steps in the context of a cryptochrome-based magnetic compass sensor (Maeda et al., 2012). Formed in a spin-correlated singlet state (total electron spin quantum number, $S = 0$ ), the radical pair coherently interconverts with the corresponding triplet state ( $S = 1$ ) as a result of the internal hyperfine interactions and the Zeeman interactions of the electron spins with any external magnetic fields (Hore and Mouritsen, 2016). In cryptochrome, the singlet channel is a reverse electron transfer reaction that returns the charge-separated radical pair state to the diamagnetic ground state of the protein. The triplet channel is considered to produce the signalling state that initiates magnetic signal-transduction (Hore and Mouritsen, 2016). For simplicity, the two reaction channels were assigned identical first order rate constants, $k$ (the so-called ‘exponential model’ (Timmel et al., 1998)). Magnetic field effects arise from the competition between these two spin-conserving reactions together with the change in the extent and timing of the coherent singlet-triplet interconversion brought about by the external magnetic field(s). Although these calculations were performed for a highly simplified model of a radical pair state of cryptochrome, we will argue below (section titled 'Assumptions and approximations') that general conclusions can be drawn that are relevant to any chemically and physically realistic radical pair reaction.

Radical pairs can be influenced by magnetic interactions that are minuscule compared to $k_{B} T$ because they are formed in a spin-coherent state far removed from thermal equilibrium. To be sensitive to a weak static magnetic field, the coherence must persist for a time comparable to the period of the electron Larmor precession which, in a 50 μT magnetic field, is ~700 ns (Rodgers and Hore, 2009). Theoretical treatments often assume negligible spin relaxation during the radical pair lifetime; this is almost always unrealistic (Hore and Mouritsen, 2016; Kattnig et al., 2016b; Kattnig et al., 2016c; Worster et al., 2016). The reality is that random molecular motions modulate the local magnetic fields experienced by the electron spins causing irreversible loss of spin coherence. Spin relaxation was included here by means of a phenomenological exponential decay towards the equilibrium state (25% singlet, 75% triplet), with a rate constant $r$ = 10⁶ s⁻¹, corresponding to a spin-coherence lifetime, $r^{- 1}$ = 1 μs. This choice of r is discussed in 'Assumptions and approximations'.

Mutual exchange and dipolar interactions of the two unpaired electrons (Efimova and Hore, 2008) were ignored. The effects of external magnetic fields, both static and time-dependent, were quantified by calculating $Φ_{T}$ , the fraction of radical pairs that react via the triplet channel (‍ $0 \leq Φ_{T} \leq 1$ ‍, Figure 1(c)). We refer to this quantity as the triplet yield. The validity and implications of all the approximations and assumptions involved are discussed in 'Assumptions and approximations'.

Calculation of $Φ_{T}$

Request a detailed protocol

In the absence of spin relaxation and chemical reactions, the probability that a radical pair is in a triplet state a time $t$ after being created in a singlet state (compare equation (4) of Timmel et al., 1998) is

p_{T}^{'} (t) = 1 - \frac{1}{M} \sum_{m}^{} \sum_{n}^{} | ⟨ m | {\hat{P}}^{S} | n ⟩ |^{2} \cos ([ω_{m} - ω_{n}] t),

in which $ℏ ω_{i} = ⟨ i | \hat{H} | i ⟩$ , where $| i ⟩$ and $ℏ ω_{i}$ are, respectively, the eigenstates and eigenvalues of the radical pair spin Hamiltonian, $\hat{H}$ , which accounts for the hyperfine and Zeeman interactions of the radicals. ${\hat{P}}^{S}$ is the singlet projection operator and $M$ is the total number of nuclear spin configurations (Timmel et al., 1998). The oscillatory time-dependence of the triplet probability reflects the formation of the radical pair in a coherent, non-stationary, superposition state.

The spin Hamiltonian of the model [FAD^•− TrpH^•+] radical pair was

\hat{H} = ω ({\hat{S}}_{A z} + {\hat{S}}_{B z}) + a_{1} {\hat{S}}_{A} . {\hat{I}}_{1} + a_{2} {\hat{S}}_{A} . {\hat{I}}_{2} + a_{3} {\hat{S}}_{B} . {\hat{I}}_{3}

with the following spin angular momentum operators: ${\hat{S}}_{A}$ and ${\hat{S}}_{A z}$ , electron in FAD^•−; ${\hat{S}}_{B}$ and ${\hat{S}}_{B z}$ , electron in TrpH^•+ (or Z^•, see later); ${\hat{I}}_{1}$ and ${\hat{I}}_{2}$ nitrogens N5 and N10 in FAD^•−; ${\hat{I}}_{3}$ ‍, nitrogen N1 in TrpH^•+. $ω = γ_{e} B$ is the electron Larmor frequency in a static magnetic field of strength $B$ and $a_{n}$ ( $n = 1 - 3$ ) are the hyperfine coupling constants.

Spin relaxation was introduced phenomenologically (compare equation (19) of Bagryansky et al., 2007) as:

p_{T} (t) = \frac{3}{4} + (p_{T}^{'} (t) - \frac{3}{4}) \exp (- r t),

where $r$ is the relaxation rate constant.

Following (Timmel et al. (1998), the chemical fate of the radical pair was modelled by means of separate spin-selective reactions of the singlet and triplet pairs (Figure 1(b)). To simplify the calculation, these two processes were assigned identical first order rate constants, $k$ . The ultimate yield of the product formed spin-selectively from the triplet radical pair state at a time $t ≫ k^{- 1}$ is therefore:

\begin{matrix} Φ_{T} = k \int_{0}^{\infty} p_{T} (t) e^{- k t} d t \\ = \frac{3}{4} + \frac{k}{4 (k + r)} - \frac{1}{M} \sum_{m}^{} \sum_{n}^{} | ⟨ m | {\hat{P}}^{S} | n ⟩ |^{2} \frac{k (k + r)}{{(k + r)}^{2} + {(ω_{m} - ω_{n})}^{2}} . \end{matrix}

Equation (4) was used to calculate $Φ_{T}$ and hence $m f e_{ELF}$ and $m f e_{GMF}$ (see below).

Calculation of $m f e_{ELF}$ and $m f e_{GMF}$

Request a detailed protocol

The same model radical pair was modelled with a weak ELF field superimposed on the geomagnetic field (GMF, $B_{0}$ = 50 μT). Epidemiological studies of childhood leukaemia (Ahlbom et al., 2000; Greenland et al., 2000; Kheifets et al., 2010) found a two-fold increase in risk for chronic exposure to ELF magnetic fields of root-mean-square magnetic flux density ≥0.4 μT with an average strength of ~0.7 μT. The latter corresponds to a peak intensity of 0.7 × √2 ≈ 1.0 μT. We therefore take the ELF field to have peak magnetic flux density $B_{1}$ = 1.0 μT and frequency 50 Hz. This field is assumed to be linearly polarised and aligned parallel to the GMF such that the total field experienced by the radicals varies between 49 and 51 μT. Any other relative orientation of the two fields would result in a smaller range of total field strengths and smaller ELF effects (see 'Assumptions and approximations').

Generally requiring lifetimes between 100 ns and 10 μs for a significant low field effect (see 'Static magnetic fields'), magnetically sensitive radical pairs have a fleeting existence compared to the 20 ms period of a 50 Hz electromagnetic wave. Consequently the ELF field can be treated as effectively static during the lifetime of a radical pair (Scaiano et al., 1994b). We suppose that radical pairs are formed continuously, for example by photo-excitation of cryptochrome. Each radical pair therefore experiences a different, effectively static, magnetic field whose intensity, $B$ , depends on the phase of the ELF field, α, which barely changes during the lifetime of the pair,

B = B_{0} + B_{1} \cos α,

with α randomly distributed in the range ( $0, π$ ). The net effect of the ELF field on an ensemble of independently created radical pairs is an average over α,

\bar{Φ_{T} (B_{0}, B_{1})} = \frac{1}{π} \int_{0}^{π} Φ_{T} (B) d α .

We define the magnetic field effect, $m f e_{ELF}$ , as the relative difference between $\bar{Φ_{T} (B_{0}, B_{1})}$ and $Φ_{T} (B_{0})$ , the triplet yield in the absence of the ELF field:

m f e_{ELF} = \frac{\bar{Φ_{T} (B_{0}, B_{1})} - Φ_{T} (B_{0})}{Φ_{T} (B_{0})} .

$m f e_{ELF}$ may be evaluated by using $B_{1} << B_{0}$ to expand $Φ_{T} (B)$ as a Taylor series around $B = B_{0}$ , to second order in $B - B_{0}$ :

Φ_{T} (B) \approx Φ_{T} (B_{0}) + (B - B_{0}) Φ_{T}^{(1)} (B_{0}) + \frac{1}{2} {(B - B_{0})}^{2} Φ_{T}^{(2)} (B_{0}),

where $Φ_{T}^{(n)} (B_{0})$ is the $n$ -th derivative of $Φ_{T} (B)$ evaluated at $B = B_{0}$ . Combining Equations (5), (6) and (8) gives the average triplet yield as a sum of zero and second order terms in $B_{1}$ :

\bar{Φ_{T} (B_{0}, B_{1})} \approx Φ_{T} (B_{0}) + \frac{1}{4} B_{1}^{2} Φ_{T}^{(2)} (B_{0}) .

Combining Equations (7) and (9) gives:

m f e_{ELF} \approx \frac{1}{4} B_{1}^{2} \frac{Φ_{T}^{(2)} (B_{0})}{Φ_{T} (B_{0})} .

When $B_{1} << B_{0}$ , the magnetic field effect, $m f e_{ELF}$ , is thus proportional to the second derivative (i.e. the curvature) of $Φ_{T} (B)$ at $B = B_{0}$ and to the square of the intensity of the ELF field (Adair, 1999; Adair, 1994).

By analogy with Equation (7) , we also define a geomagnetic field effect, $m f e_{GMF}$ , as the fractional change in the triplet yield in the absence of an ELF field, when the static magnetic field is changed from $B_{0}$ to $B_{0} + Δ B_{0}$ where $| Δ B_{0} | ≪ B_{0}$ :

m f e_{GMF} = \frac{Φ_{T} (B_{0} + Δ B_{0}) - Φ_{T} (B_{0})}{Φ_{T} (B_{0})} .

Expansion of $Φ_{T} (B)$ as a Taylor series around $B = B_{0}$ , to first order in $Δ B_{0}$ , gives

m f e_{GMF} \approx Δ B_{0} \frac{Φ_{T}^{(1)} (B_{0})}{Φ_{T} (B_{0})},

which should be compared with Equation (10). In contrast to $m f e_{ELF}$ , $m f e_{GMF}$ depends on the first derivative (i.e. the gradient) of $Φ_{T}$ and is linear in $Δ B_{0}$ .

Results

Static magnetic fields

We first look at the dependence of the triplet yield, $Φ_{T}$ , on the strength of an external static magnetic field, $B_{0}$ , to provide a basis for the discussion of ELF magnetic field effects. Figure 2(a) shows $Φ_{T}$ as a function of $B_{0}$ in the range 0–5 mT for seven values of the radical pair lifetime, $τ = k^{- 1}$ , between 30 ns and 100 μs. The salient features are as follows. (a) For lifetimes greatly in excess of 1 μs, $Φ_{T}$ is almost independent of $B_{0}$ and approximately equal to 0.75. (b) For intermediate lifetimes ( $τ$ ≈ 1 μs), the magnetic field effect is bi-phasic: a small increase in $Φ_{T}$ for $B_{0}$ < 1 mT is followed by a larger decrease which levels out at fields in excess of 5 mT. The initial ‘bump’, known as the ‘low field effect’ (LFE), has been observed in experiments on a variety of radical pair reactions (Maeda et al., 2012; Kattnig et al., 2016a; Timmel et al., 1998; Brocklehurst, 1976; Eveson et al., 2000). (c) When the lifetime of the radical pair is much shorter than 1 μs, the LFE vanishes and the magnetic field effect becomes mono-phasic. (d) Compared to $B_{0} = 0$ , a typical geomagnetic field ( $B_{0}$ = 50 μT) produces the largest change in $Φ_{T}$ when the LFE is at its most prominent, that is for lifetimes, $τ$ ≈ 1 μs. This can be seen more clearly in Figure 2(b).

Figure 2 with 1 supplement see all

Download asset Open asset

Effects of a static magnetic field on a model of the [FAD^•− TrpH^•+] radical pair in cryptochrome.

(a) Dependence of the reaction yield $Φ_{T}$ on the static magnetic field strength, $B_{0}$ , for different values of the radical pair lifetime, $τ$ . The spin relaxation time is 1 μs. (b) An expanded view of the low field region of (a) with $B_{0}$ = 50 μT indicated by the dashed line. The difference between $Φ_{T}$ at zero field and at 50 μT is largest when $τ$ ≈ 1 μs. The corresponding calculations for a spin relaxation time of 0.1 μs are shown in Figure 2—figure supplement 1.

https://doi.org/10.7554/eLife.44179.003

All the features of Figure 2 that relate to Earth-strength magnetic fields have been observed experimentally and can readily be understood (Rodgers, 2009; Timmel et al., 1998). Briefly, if the radical pair exists for much less than 1 μs, there is no time for significant Larmor precession (frequency, 1.4 MHz) in a 50 μT field and hence no additional singlet−triplet interconversion. If the radical pair lives much longer than 1 μs, spin relaxation destroys the spin coherence and establishes a 1:3 singlet:triplet ratio before the radicals can react. In both cases the effect of a 50 μT magnetic field effect is minimal (Kattnig et al., 2016b; Kattnig et al., 2016c; Worster et al., 2016). The bi-phasic character of $Φ_{T}$ arises from two distinct mechanisms: oscillation of low frequency zero-quantum coherences at small $B_{0}$ and energetic isolation of two of the three triplet sub-levels at high $B_{0}$ (Timmel et al., 1998; Till et al., 1998). The slight irregularities in $Φ_{T}$ , visible in some of the traces in Figure 2(a), arise from avoided energy-level crossings (Timmel et al., 2001; Hiscock et al., 2016). Such features are usually imperceptible for more realistic spin systems containing a larger number of nuclear spins.

Two other features of Figure 2 will be relevant for our discussion of ELF effects. At $B_{0}$ ≈ 50 μT, the gradient, ‍ $d Φ_{T} / d B_{0}$ , is positive and the curvature, $d^{2} Φ_{T} / d B_{0}^{2}$ , mostly negative (i.e. concave downward).

ELF magnetic fields

Before showing simulations of $m f e_{ELF}$ , we anticipate some of its properties. The dependence of $m f e_{ELF}$ on the curvature of $Φ_{T}$ in Equation (10) is rationalised in Figure 3(a) which shows that $m f e_{ELF}$ = 0 if $Φ_{T}$ depends linearly on $B$ in the range $(B_{0} - B_{1}) \leq B \leq (B_{0} + B_{1})$ . In that case, the average effect of static and ELF fields together (blue arrow) is the same as for the static field alone (green arrow). The change in $Φ_{T}$ when 0 ≤ α ≤ $\frac{1}{2} π$ exactly cancels that for $\frac{1}{2} π$ ≤ α ≤ π. This is not the case when $Φ_{T}$ has a non-linear dependence on $B$ , Figure 3(b).

Figure 3

Download asset Open asset

Schematic representations (thick black lines) of the dependence of the reaction yield, $Φ_{T}$ , on the strength of a static magnetic field.

The orange arrows indicate the yields for the maximum and minimum values of $B$ in Equation (5). The green arrows show the yields when $B = B_{0}$ . The blue arrows mark the free radical yields averaged over the phase of the ELF field, Equation (6) . (a) When $Φ_{T}$ is linear in $B$ , the effect of the GMF and the ELF field *together* equals that of the GMF *alone*. (b) When $Φ_{T}$ is non-linear, the effects of GMF plus ELF and GMF-alone differ. $Φ_{T}$ in (b) has been drawn with a negative curvature (concave downward), with the result that $\bar{Φ_{T} (B_{0}, B_{1})} < Φ_{T} (B_{0})$ and $m f e_{ELF} < 0$ (Equation (7) and (10)). A positive curvature (concave upward) would give the opposite signs.

https://doi.org/10.7554/eLife.44179.005

From Figures 2 and 3, we can anticipate that for $B_{0}$ = 50 μT and $B_{1}$ = 1.0 μT, $m f e_{ELF}$ will be largest when the LFE is strongest, that is when τ ≈ 1 μs. Figure 4(a) shows that this is indeed the case. For lifetimes in the range 1 ns−1 ms, the largest change in the reaction yield caused by the ELF magnetic field is −1.2 parts per million (ppm); $m f e_{ELF}$ is negative because $Φ_{T}^{(2)} (B_{0})$ < 0 for most values of τ. The significance of the sign of $m f e_{ELF}$ is discussed in 'Signs of $m f e_{ELF}$ and $m f e_{GMF}$ '.

Figure 4 with 3 supplements see all

Download asset Open asset

Magnetic field effects (in ppm) on model radical pairs as a function of the lifetime of the radical pair, τ, in the range 1 ns−1 ms.

The spin relaxation time is 1 μs. (a) and (c) A model of the [FAD^•− TrpH^•+] radical pair in cryptochrome, containing two ¹⁴N nuclei in the FAD^•− radical and one in the TrpH^•+ radical. (b) and (d) A model of the [FAD^•− Z^•] radical pair in cryptochrome in which the Z^• radical has no hyperfine interactions. (a) and (b) $m f e_{ELF}$ (Equation (10)) for $B_{0}$ = 50 μT, $B_{1}$ = 1.0 μT. (c) and (d) $m f e_{GMF}$ (Equation (12)) for $B_{0}$ = 50 μT, $Δ B_{0}$ = −1.0 μT. Note the different vertical scales of the four panels. The validity of the Taylor series expansion (Equation (9)) was confirmed by numerical integration over α, the phase of the ELF field (Figure 4—figure supplement 1). Corresponding calculations for a spin relaxation time of 0.1 μs are shown in Figure 4—figure supplement 2. Corresponding versions of (a) and (c) when 25 μT ≤ $B_{0}$ ≤ 65 μT are shown in Figure 4—figure supplement 3.

https://doi.org/10.7554/eLife.44179.006

Although magnetic field effects on purified cryptochromes arise from [FAD^•− TrpH^•+], there is some evidence that, in vivo, FAD^•− is partnered by a radical with fewer, smaller hyperfine interactions than TrpH^•+ (Lee et al., 2014; Ritz et al., 2009). The simulation in Figure 4(a) was therefore repeated for a radical pair, [FAD^•− Z^•], identical to [FAD^•− TrpH^•+] except that the second radical, Z^•, has no hyperfine interactions. Previous studies, both experimental and theoretical, have shown that such highly asymmetric radical pairs show larger low field effects than when the hyperfine interactions are more evenly distributed between the two radicals (Lee et al., 2014; Stass et al., 1995; Rodgers et al., 2007). As Figure 4(b) shows, this is also the case for ELF fields. The maximum ELF effect for the model [FAD^•− Z^•] pair is −14 ppm, roughly an order of magnitude larger than the largest $m f e_{ELF}$ for [FAD^•− TrpH^•+].

Comparison of ELF and GMF effects

Figure 4(a) and (b) suggest that the effect of a 1 μT ELF magnetic field is likely to be no more than a few parts per million. How should such numbers be interpreted in terms of any potential biological significance? The simplest approach would be to assert that changes in the yield of a chemical reaction smaller than, say, 100 ppm (i.e. 0.01%) can be dismissed as insignificant. Although arguably reasonable, the choice of any such threshold is necessarily arbitrary, especially as we have deliberately refrained from considering specific biological processes. An alternative would be to calculate a ‘signal-to-noise’ ratio in which the ‘signal’ would be a change in concentration (e.g. of a reactive oxygen species) induced by the ELF field and the ‘noise’ would be the average background fluctuation in that concentration. Then, with some confidence, one could reject potential ELF effects that had a predicted signal-to-noise ratio less than 1.0. Clearly, this too is unsatisfactory due to lack of knowledge of both the ‘noise’ and the link between the radical pair reaction and a biological end-point. We have therefore chosen not to rely solely on the absolute values of $m f e_{ELF}$ but instead to compare predicted ELF effects with the changes in $Φ_{T}$ that would arise from small variations in the strength of the geomagnetic field such as might be experienced by moving to a different point on the Earth’s surface. Performing the two calculations using the same model and the same parameters allows ELF effects to be put into a more readily appreciated perspective. Furthermore, by taking the ratio of $m f e_{ELF}$ (Equation (10)) and $m f e_{GMF}$ (Equation (12)) under identical conditions, any errors arising from the approximations and assumptions of the model will tend to cancel.

To facilitate comparison with the ELF effects shown in Figure 4(a) and (b) (for which $B_{1}$ = 1.0 μT), we keep $B_{0}$ fixed at 50 μT and choose $Δ B_{0}$ = −1.0 μT (negative, so that $m f e_{GMF}$ has, for the most part, the same sign as $m f e_{ELF}$ ). The results are shown in Figure 4(c) and (d) for the model [FAD^•− TrpH^•+] and [FAD^•− Z^•] radical pairs respectively. In both cases, the largest $m f e_{GMF}$ exceeds the largest $m f e_{ELF}$ by two orders of magnitude. In the following we denote these maximum magnetic field effects $m f e_{GMF}^{max}$ and $m f e_{ELF}^{max}$ .

Equation (12) and (10) show that $m f e_{GMF}$ and $m f e_{ELF}$ scale, respectively, linearly in $Δ B_{0}$ and quadratically in $B_{1}$ . We can therefore use Figure 4 to predict $m f e_{GMF}$ , $m f e_{ELF}$ and their ratio for different values of $Δ B_{0}$ and $B_{1}$ , provided $| Δ B_{0} | ≪ B_{0}$ and $B_{1} ≪ B_{0}$ . For example, for Figures 4(a) and 4(c), $m f e_{GMF}^{max}$ is approximately 280 times larger than $m f e_{ELF}^{max}$ . $m f e_{GMF}$ and $m f e_{ELF}$ should therefore have similar maximum amplitudes when $Δ B_{0}$ is reduced by a factor of 280. Thus, for the model of [FAD^•− TrpH^•+], a −3.6 nT change in the geomagnetic field should have roughly the same effect on $Φ_{T}$ as would a 1 μT ELF field. For the simpler [FAD^•− Z^•] radical pair (Figures 4(b) and 4(d), where the ratio of $m f e_{GMF}^{max}$ to $m f e_{ELF}^{max}$ is 150), the corresponding $Δ B_{0}$ would be ≈ −6.7 nT.

Comparison with temperature effects

A different comparison, which may also help to put the predicted ELF effects into context, relates to the daily variation in body temperature which, for a healthy adult, is typically ±0.5°C. The two parameters in our calculation that are most likely to be temperature-sensitive are the rate constants for the reactions ( $k$ ) and the relaxation ( $r$ ) of the radicals. By analogy with Equation (10) and (12), we define $T e_{Δ k, Δ r}$ ( $T e$ = temperature effect) as the change in $Φ_{T}$ resulting from small increases ( $Δ k$ and $Δ r$ ) in $k$ and $r$ , respectively,

T e_{Δ k, Δ r} = \frac{Φ_{T}^{(k + Δ k, r + Δ r)} - Φ_{T}^{(k, r)}}{Φ_{T}^{(k, r)}} .

To obtain crude estimates of the temperature dependence of $Φ_{T}$ , we assume that both $k$ and $r$ increase/decrease by 0.1% for a temperature rise/fall of 1°C. Assuming Arrhenius behaviour (rate $\propto \exp [- E_{a} / R T]$ ) and $T$ = 37°C, this corresponds to a small activation energy, $E_{a}$ = 0.80 kJ mol⁻¹ (≈ 0.31 $R T$ ). Larger activation energies would lead to more dramatic changes in these rate constants.

Taking $k$ = $r$ = 10⁶ s⁻¹, $B_{0}$ = 50 μT and $Δ B_{0}$ = $B_{1}$ = 0, we find $T e_{Δ k, 0}$ = −36 ppm and $T e_{0, Δ r}$ =+21 ppm for [FAD^•− TrpH^•+] and $T e_{Δ k, 0}$ = −54 ppm and $T e_{0, Δ r}$ = +27 ppm for [FAD^•− Z^•]. The two temperature effects are additive: increasing both rate constants by 0.1% simultaneously results in $T e_{Δ k, Δ r}$ = −15 ppm for [FAD^•− TrpH^•+] and −27 ppm for [FAD^•− Z^•]. These numbers should be compared with the corresponding values of $m f e_{ELF}$ in Figure 4(a) and (b).

For convenience, all the numerical results reported in this section are collected in Table 1, with $Ω_{ELF}^{GMF}$ defined as the ratio $m f e_{GMF}^{max} / m f e_{ELF}^{max}$ .

Table 1

Changes in reaction product yields (in ppm) arising from small changes in magnetic field strengths and temperature

https://doi.org/10.7554/eLife.44179.010

		[FAD^•− TrpH^•+]	[FAD^•−Z^•]
$m f e_{ELF}^{max}$ ^*	$B_{1}$ =1.0 μT	−1.2	−14
$m f e_{GMF}^{max}$ ^*	$Δ B_{0}$ = −1.0 μT	−330	−2100
$Ω_{ELF}^{GMF}$ ^†		280	150
$T e_{Δ k, 0}$ ^‡	$B_{1} = Δ B_{0} = 0$	−36	−54
$T e_{0, Δ r}$ ^§	$B_{1} = Δ B_{0} = 0$	+21	+27
$T e_{Δ k, Δ r}$ ^¶	$B_{1} = Δ B_{0} = 0$	−15	−27

* $r$ = 10⁶ s⁻¹, $B_{0}$ = 50 μT

^† $Ω_{ELF}^{GMF} = m f e_{GMF}^{max} / m f e_{ELF}^{max}$
^‡ $k$ = $r$ = 10⁶ s⁻¹, $Δ k$ = 10³ s⁻¹, $Δ r$ = 0, $B_{0}$ = 50 μT

^§ $k$ = $r$ = 10⁶ s⁻¹, $Δ k$ = 0, $Δ r$ = 10³ s⁻¹, $B_{0}$ = 50 μT
^¶ $k$ = $r$ = 10⁶ s⁻¹, $Δ k$ = $Δ r$ = 10³ s⁻¹, $B_{0}$ = 50 μT

Discussion

Signs of $m f e_{ELF}$ and $m f e_{GMF}$

Equations (10) and (12) predict that $m f e_{ELF}$ and $m f e_{GMF}$ depend, respectively, on the curvature and gradient of $Φ_{T}$ at $B_{0}$ = 50 μT and Figure 2 shows that these quantities are, respectively, negative and positive for radical pairs that have lifetimes between ~100 ns and 10 μs. Our simulations (Figure 4) bear out this expectation. The presence of a weak ELF field or a small decrease in the strength of the geomagnetic field both reduce the yield of the triplet product and correspondingly increase $Φ_{S}$ , the singlet yield (‍ $Φ_{S} + Φ_{T} = 1$ , Figure 1(b) and (c)). The signs of $m f e_{ELF}$ and $m f e_{GMF}$ are both reversed if the radical pair starts out as a triplet rather than a singlet. Given our wish to be agnostic about the chemical and biological details, it makes no sense to interpret the absolute signs of $m f e_{ELF}$ and $m f e_{GMF}$ .

Nor are the relative signs of $m f e_{ELF}$ and $m f e_{GMF}$ unconditionally useful for interpreting experimental data. Although $m f e_{ELF}$ and $m f e_{GMF}$ have mostly negative values for both radical pairs in Figure 4 (where $r$ = 10⁶ s⁻¹), positive values of both quantities are predicted when $r$ = 10⁷ s⁻¹ (Figure 4—figure supplement 2) as a result of the different dependence of $Φ_{T}$ on $B_{0}$ (Figure 2—figure supplement 1).

Assumptions and approximations

Before attempting to draw conclusions from the data in Table 1, we first assess the validity of the assumptions and approximations used to obtain them. The most important of these concern (a) spin relaxation, (b) hyperfine interactions, (c) electron-electron interactions, (d) relative motion of the radicals, and (e) the ELF field.

(a) Spin relaxation of organic radicals at physiological temperatures is likely to be faster than the 1 μs we have assumed here and unlikely to be slower (Kattnig et al., 2016b). Although one could imagine longer relaxation times (e.g. ~10 μs) in an avian protein (e.g. cryptochrome) optimised for exquisitely sensitive detection of the geomagnetic field, it is highly improbable that these properties could arise by chance. It would require that the random thermal motions of the radicals had very low amplitude or very high frequency or both (Hore and Mouritsen, 2016). Neither seems probable in a biological environment at physiological temperatures. If spin relaxation is faster than 1 μs, all magnetic field effects will be smaller than those in Table 1. For example, $m f e_{ELF}$ and $m f e_{GMF}$ are both reduced by roughly an order of magnitude when $r^{- 1}$ = 0.1 μs, although their ratio is not greatly changed (see Figure 4—figure supplement 2). By contrast, changing the relaxation time from 1 μs to 10 μs increases $m f e_{ELF}^{max}$ and $m f e_{GMF}^{max}$ by factors of 2.7 and 1.7 respectively, and consequently reduces $Ω_{ELF}^{GMF}$ by 40%. There are only minor changes in $m f e_{ELF}^{max}$ and $m f e_{GMF}^{max}$ for relaxation times longer than 10 μs. This is because 10 μs is already an order of magnitude longer than the reciprocal of the electron Larmor frequency (~700 ns) in a static magnetic field of 50 μT. Although increasing the relaxation time beyond 10 μs would allow the radical pair to respond more sensitively to static fields weaker than 50 μT, the magnetic field effect at 50 μT field is already at its maximum.

(b) Organic radicals with the properties needed to show magnetic field effects almost always have more hyperfine interactions than the one or two that have been included here (Gerson and Huber, 2003). The semi-occupied molecular orbital that contains the unpaired electron in an organic radical is often delocalised so that the electron interacts with several nearby hydrogen and/or nitrogen nuclei. Our experience of spin dynamics simulations has been that the more hyperfine interactions included in both radicals, the smaller the effects of weak magnetic fields. Conversely, the magnetic sensitivity is generally larger when one of the radicals has just a few small hyperfine interactions or none at all (Lee et al., 2014). As we have seen with [FAD^•− Z^•], removing the single hyperfine interaction from the model TrpH^•+ radical, increases $m f e_{ELF}$ and $m f e_{GMF}$ by about an order of magnitude (and reduces $Ω_{ELF}^{GMF}$ by a factor of 2, from 280 to 150, Table 1). Biologically relevant radicals with no or very few hydrogen and nitrogen atoms near the unpaired electron are scarce. Superoxide, $O_{2}^{• -}$ , is one and nitric oxide, ${NO}^{•}$ , is another but both normally have such fast spin relaxation (probably ~1 ns (Karogodina et al., 2009; Karogodina et al., 2011)) that they can only participate in magnetically sensitive reactions under rather special circumstances (Hogben et al., 2009; Kattnig, 2017). The only other radical discussed in the context of magnetoreception is that formed by oxidation of ascorbic acid, Asc^•−. [FAD^•−Asc^•−] is predicted to show magnetic field effects larger than [FAD^•− TrpH^•+] but smaller than [FAD^•− Z^•] (Kattnig et al., 2016a; Lee et al., 2014; Nielsen et al., 2017). In short, it is difficult to imagine a biologically plausible radical pair whose hyperfine couplings make it significantly more sensitive to weak magnetic fields than the simplified model of [FAD^•− TrpH^•+] we have considered here.

(c) In all of the simulations discussed here, we have ignored the possibility that the unpaired electrons in the two radicals have a mutual spin-spin interaction. Although this has been common practice in theoretical treatments of radical pair magnetoreception (Hore and Mouritsen, 2016), it only starts to be a good approximation when the exchange and dipolar interactions, which tend to block singlet-triplet interconversion, are smaller than the geomagnetic field (~50 μT). This in turn would require the radicals to be separated by more than 3.8 nm (Efimova and Hore, 2008). At such large distances, it is extremely doubtful whether spin-selective recombination would be fast enough to compete with spin relaxation. Inclusion of realistic spin-spin interactions can easily attenuate the predicted magnetic field effects by an order of magnitude or more (Efimova and Hore, 2008).

(d) We have treated the radicals as fixed in space partly because that is the case in cryptochrome but more importantly because the effects of weak fields are more pronounced when the radicals are unable to move freely. Translationally diffusing radicals generally show stronger LFEs when their motions are restricted, for example by confinement within a micelle (Eveson et al., 2000; Turro and Weed, 1983; Sakaguchi and Hayashi, 1984). One of the reasons for this is that without such constraints, translational motion is an efficient source of spin relaxation via modulation of exchange interactions (Shushin, 1991). Another is that a fraction of the radical pairs inevitably diffuse apart without ever having the opportunity to undergo the spin-selective reaction(s) essential for a magnetic field effect.

(e) We have taken the ELF field to be linearly polarised and aligned parallel to the GMF such that the total field experienced by the radicals varies between 49 and 51 μT. Any other choice would result in a smaller range of total field strengths and therefore smaller ELF effects. If a 1 μT ELF field is randomly oriented with respect to a 50 μT static field, then the average total field fluctuates between 49.37 and 50.64 μT, that is a variation of ~0.64 μT instead of 1 μT. Given the quadratic dependence on the amplitude of the ELF field (Equation (10)), our choice of parallel fields overestimates the ELF effect by a factor of ~(0.64)⁻² ≈ 2.4 compared to random orientations.

To summarise, in all five cases (a-e) we have deliberately designed the model and chosen its parameters in such a way as to maximise and/or overestimate the sensitivity to weak magnetic fields.

However, as we now discuss, we have also made a few assumptions that could in practice, or at least in principle, cause $m f e_{ELF}$ and/or $m f e_{GMF}$ to exceed the values in Table 1. First we deal with three factors related to the reaction scheme in Figure 1(b). (f) Rather than starting as a singlet, the radical pair might be formed in a triplet state or arise from the encounter of radicals with uncorrelated electron spins (‘F-pairs’) (Steiner and Ulrich, 1989; Salikhov et al., 1984). (g) The singlet and triplet states of the radical pair could react with different rate constants (‘asymmetric recombination’) (Lewis et al., 2014). (h) Only one of the competing reaction channels needs to be spin-selective: the other can be independent of spin and proceed with identical singlet and triplet rate constants, as is the case for [FAD^•− TrpH^•+] in cryptochrome (‘semi-spin-selective recombination’) (Maeda et al., 2012). Simulations of the model [FAD^•− TrpH^•+] radical pair show that relaxing these three assumptions can increase $m f e_{ELF}$ by up to a factor of 8 (Blunsden, 2018). For example, $m f e_{ELF}^{max}$ increases by 40% if $k_{S}$ and $k_{T}$ are no longer constrained to be equal but is 20% smaller for semi-spin-selective recombination. Somewhat larger changes in $m f e_{ELF}^{max}$ are found when the radical pair starts in a triplet state: an 8-fold increase for asymmetric recombination and 7-fold when only the singlet recombination step is spin-selective. However, in all cases $m f e_{GMF}^{max}$ shows similar changes such that $Ω_{ELF}^{GMF}$ is no more than 20% larger than that for symmetric recombination of a singlet-born radical pair. We stress that these increases in magnetic sensitivity are maxima, corresponding to particular choices of rate constants, initial spin state and reaction scheme. Although evolution could have exploited such conditions, for example to make a more efficient magnetic compass, it seems improbable that they would occur by chance.

(i) In calculating the effects of a 1 μT ELF magnetic field we have used a geomagnetic field strength (50 μT) appropriate for northern Europe. Repeating the [FAD^•− TrpH^•+] simulations (Figure 4) with $B_{0}$ = 25 μT and 65 μT (the extreme values of the Earth’s field) give $m f e_{ELF}$ values respectively 2.2 times bigger and 1.4 times smaller than for $B_{0}$ = 50 μT (see Figure 4—figure supplement 3). Similar effects were found for $m f e_{GMF}$ which was 1.4 times bigger when $B_{0}$ = 25 μT and 1.2 times smaller when $B_{0}$ = 65 μT. We note that Swanson and Kheifets considered the possibility that ELF magnetic fields could have different consequences at different locations on the Earth’s surface due to variations in the geomagnetic field. To test this, they analysed 15 epidemiological studies and found ‘some, but rather limited and not statistically significant, evidence’ for a correlation between ELF exposure and incidence of childhood leukaemia (Swanson and Kheifets, 2012).

Over the years, various mechanisms that could amplify small radical pair magnetic field effects have been suggested, none of which we have so far considered. Four are mentioned here. First (j), is the possibility that a superparamagnetic nanoparticle could boost the strength of the ELF field experienced by a nearby radical pair (Binhi, 2008). Briefly, the idea is that the 50/60 Hz field could align the instantaneous, fluctuating magnetic moment of the nanoparticle such that the magnetic field close to its surface would be much stronger than the external ELF field but would still oscillate at 50/60 Hz. We can assess the likely importance of this effect using a thermodynamic argument based on ferritin, a naturally occurring superparamagnetic protein, which has an instantaneous magnetic moment, $m \approx 300 μ_{B}$ ( $μ_{B}$ is the Bohr magneton) (Kilcoyne and Cywinski, 1995; Vohralik and Lam, 2009; Worster and Hore, 2018). The energy of the interaction of ferritin with a $B_{1}$ = 1 μT magnetic field, $\sim m B_{1}$ , is a million times smaller than the thermal energy, $k_{B} T$ , at physiological temperature. The alignment of the nanoparticle’s magnetic moment, and therefore the amplification of the ELF field, will consequently be negligible.

In the context of cryptochrome-based magnetoreception, two amplification mechanisms have been proposed, and in one case demonstrated experimentally. They rely on (k) slow radical termination reactions (Kattnig et al., 2016a) and (l) reactions of the radicals with paramagnetic scavengers (Kattnig, 2017; Kattnig and Hore, 2017). The latter mechanism has the interesting and potentially important property that it seems to be immune to fast spin relaxation in one of the radicals (Kattnig, 2017), opening the possibility that superoxide and other reactive oxygen species could be involved in responses to weak magnetic fields. Until the mechanism is confirmed experimentally, it is difficult to know how much amplification to expect.

The final amplification mechanism (m) is more speculative (giving, perhaps, greater scope for disproportionate enhancements of ELF magnetic field effects). It requires a reaction scheme involving chemical feedback in which a radical pair is one of the key reaction intermediates (Eichwald and Walleczek, 1996; Eichwald and Walleczek, 1997; Eichwald and Walleczek, 1998; Møller and Olsen, 1999; Møller et al., 2000; Møller and Olsen, 2000; Purtov, 2004). The inherent non-linearity of such reactions could at least in principle allow small magnetically induced changes in radical pair lifetimes to have a disproportionately large effect on, for example, the amplitude of chemical oscillations. Further, one could imagine an oscillating reaction with a cycle time that coincidentally matched the period of the ELF field which might then be able to ‘pump’ the oscillations to higher amplitudes. However, this would require some degree of phase-coherence between the ELF field and the biochemical oscillator. Unless there is a mechanism by which an environmental ELF field could entrain a biological oscillator, coherent pumping and therefore selective amplification of 50/60 Hz magnetic field effects would seem rather unlikely. It is difficult to think of any mechanism which could explain a specific and disproportionate response of a radical pair reaction to an ELF field.

To sum up, we have discussed 13 assumptions and approximations in the simulations on which this article is based. In our judgement, the first five (a-e) are the most important and almost certainly lead to overestimates of the effects of both ELF and static fields. The next three (f-h) could, if operating simultaneously, enhance $m f e_{ELF}$ by up to an order of magnitude, but only for specific combinations of rate constants. Using a value other than 50 μT for the Earth’s magnetic field (i) would increase or decrease $m f e_{ELF}$ and $m f e_{GMF}$ , but only by a factor of ~2. We judge the next three (j-l) unlikely to have any substantial consequences for biological radical pair reactions that are fortuitously sensitivity to magnetic fields. All of (a-l) should increase or decrease $m f e_{ELF}^{max}$ and $m f e_{GMF}^{max}$ to approximately the same degree making their ratio, $Ω_{ELF}^{GMF}$ , much less sensitive to the exact conditions of the simulations. Finally, it is difficult to assess the importance of amplification via chemical feedback (m) given the highly speculative nature of this mechanism. Certainly, a set of coupled reactions would require some very unusual (and unknown) properties to be able to boost magnetic field effects preferentially at 50/60 Hz. It is not inconceivable that there exists in human biology systems at the cellular level that have evolved sensitivity to the Earth’s magnetic field. If so, then some of the above assumptions may be less secure than presented leading to underestimates of magnetic field effects. However, it is still difficult to imagine situations that would selectively enhance responses at 50/60 Hz.

Additional aspects

Before concluding, we want to mention briefly five additional aspects of radical pair chemistry/physics.

First, when a radical pair is formed in a singlet state, but not when it is formed as an unpolarised triplet, its electron spins are entangled (Hogben et al., 2012). The possibility that entanglement arises naturally in the ‘warm, wet and noisy’ milieu of a living cell has attracted a certain amount of attention from theoretical physicists (Cai et al., 2010; Gauger et al., 2011; Kominis, 2012; Pauls et al., 2013; Tiersch et al., 2014; Zhang et al., 2014). Nevertheless, there is no apparent way in which this entanglement can enhance magnetic responses (Hore and Mouritsen, 2016).

Second, humans have only been exposed to widespread anthropogenic ELF fields since the early days of electrification in the late 19th century. If there is a mechanism by which radical pair reactions can be unusually sensitive to 50/60 Hz magnetic fields, it cannot have evolved in only ~100 years and would have to be a chance consequence of, for example, an unknown cellular response to the Earth’s static magnetic field that had evolved at a much earlier stage in human development.

Third, radical pairs are well known to show resonant responses to radiofrequency magnetic fields in the frequency range 1–100 MHz (Henbest et al., 2004; Frankevich and Kubarev, 1982; Woodward et al., 2001). Similar effects cannot occur at ELF frequencies because spin relaxation will destroy all spin coherence on a timescale much faster than the 20 ms period of a 50 Hz oscillation. To put it another way: with a ~1 μs coherence lifetime, any resonance would be ~1 MHz wide, thus precluding any possibility that the sensitivity to a 50 Hz field could be greater than that for a static field.

Fourth, it appears that radical pairs can only be formed in cryptochromes when the FAD cofactor is correctly bound. Without the flavin, there is very limited scope for the intra-protein electron transfer reactions that could produce magnetically responsive radical pairs. Current evidence suggests that Type 2 vertebrate cryptochromes may be ‘vestigial’ flavoproteins, unable to bind FAD (Kutta et al., 2017). If true, then human cryptochromes (exclusively Type 2) are unlikely to form radical pairs and therefore improbable as agents of biological radical pair effects. By contrast, Type 4 cryptochromes, which are found in birds, fish and reptiles (Kobayashi et al., 2000; Kubo et al., 2006; Kubo et al., 2010), do bind FAD (Ozturk et al., 2009; Mitsui et al., 2015; Qin et al., 2016; Wang et al., 2018) and appear to be fit for purpose as magnetoreceptors in migratory songbirds (Günther et al., 2018).

Fifth, not all radical pair reactions are magnetically sensitive. For example, several flavin-dependent enzymes known or thought to proceed via radical pair intermediates show no magnetic field effects because catalysis is not rate-limited by spin-selective reaction steps (Messiha et al., 2015).

Conclusions

We believe the values of $m f e_{ELF}^{max}$ in Table 1 (−1.2 ppm and −14 ppm for the two model radical pairs considered) provide a good estimate of the maximum fortuitous effect of a 1 μT 50/60 Hz magnetic field on any radical pair reaction. Given the arguments in 'Assumptions and approximations' and the scarcity of biological radicals devoid of significant hyperfine interactions, we propose that 10 ppm should be seen as an upper limit on the change in the yield of a radical pair reaction produced by a 1 μT ELF field in the presence of the Earth’s magnetic field (25–65 μT).

To put this number in context, we have estimated the changes in the yields of radical pair reactions that might result from ±0.5°C variations in temperature. These changes (Table 1) are of the same order of magnitude (~10 ppm) as, or larger than, $m f e_{ELF}^{max}$ .

Table 1 also contains estimates of the maximum effect ( $m f e_{GMF}^{max}$ ) of a 1 μT reduction in the strength of the geomagnetic field in the absence of an ELF field (−330 ppm and −2100 ppm). For similar reasons, we believe these values should also be appropriate for radical pair reactions that have not been optimised for a role in magnetic sensing. The ratios of the two maximum magnetic field effects, $Ω_{ELF}^{GMF} = m f e_{GMF}^{max} / m f e_{ELF}^{max}$ , in Table 1 (280 and 150 for the two model radical pairs) suggest that $Ω_{ELF}^{GMF}$ is likely to be in the range 100–500 for any radical pair reaction. In other words, a 1 μT decrease in the strength of the geomagnetic field should have an effect 100–500 times larger than would a 1 μT ELF field in the presence of the geomagnetic field. Or, equivalently, a 2–10 nT change in the strength of the geomagnetic field should have a similar effect to that of a 1 μT ELF field.

The last result can be put into context by considering the variation in the strength of the geomagnetic field over the surface of the Earth: ~65 μT at the magnetic poles and ~25 μT at the magnetic equator. Given the circumference of the Earth (~40,000 km), the average variation in the geomagnetic field along a line of longitude is roughly 4 nT km⁻¹. Therefore, the maximum effect of a 1 μT ELF field on a radical pair reaction (10 ppm) should be comparable to the maximum effect of travelling 0.5–2.5 km along a north-south axis.

A related comparison, on a much smaller length scale, is provided by measurements of the local distortions in the Earth’s magnetic field in homes in the UK caused by the proximity of ferrous objects (Swanson, 1994). Variations in static field strength of order 1 μT can be experienced by movement from room to room, corresponding to much larger effects on radical pair reactions than could be expected from a 1 μT ELF field.

It may also be relevant to note that the Earth’s magnetic field has a component of amplitude 25–50 nT that oscillates with a 24 hr period (caused by the tidal effect of the sun’s gravity on the Earth’s atmosphere (Liboff, 2014)). Potentially, therefore, the natural diurnal variation in the geomagnetic field could have a larger effect on radical pair chemistry than a 1 μT 50/60 Hz field.

To conclude, the predicted effects of 1 μT ELF magnetic fields in the presence of the Earth’s magnetic field are small (<10 ppm) and similar to or smaller than effects on the same reactions resulting either from travelling a few kilometres or from natural fluctuations in body temperature neither of which would normally be considered as potentially harmful to human health.

Implicit in everything, we have written here is the assumption that the current theory of the radical pair mechanism is complete. We cannot exclude the possibility that, despite more than 40 years of research, there is some completely unknown aspect of radical pair spin chemistry that confers a disproportionate sensitivity to ELF magnetic fields. We cannot begin to imagine what that aspect might be except to note a possibly related observation. Migratory birds are prevented from using their magnetic compass by extraordinarily weak broadband radiofrequency noise (Ritz et al., 2009; Ritz et al., 2004; Thalau et al., 2005; Winklhofer et al., 2013; Engels et al., 2014; Kavokin et al., 2014; Wiltschko et al., 2015; Schwarze et al., 2016) the predicted effects of which are far too small to be consistent with our current understanding of radical pair spin dynamics (Kavokin, 2009; Hiscock et al., 2017). It is possible, therefore, that a deeper understanding of the mechanism of avian magnetoreception will bring new insights into the risks associated with exposure to weak environmental 50/60 Hz magnetic fields.

Data availability

All data generated during this study are included in the manuscript and supporting files. Code files have been provided for Figures 2 and 4.

References

1. Adair RK
(1991) Constraints on biological effects of weak extremely-low-frequency electromagnetic fields
Physical Review A 43:1039–1048.
https://doi.org/10.1103/PhysRevA.43.1039
- PubMed
- Google Scholar
1. Adair RK
(1994) Biological responses to weak 60-Hz electric and magnetic fields must vary as the square of the field strength
PNAS 91:9422–9425.
https://doi.org/10.1073/pnas.91.20.9422
- PubMed
- Google Scholar
1. Adair RK
(1999) Effects of very weak magnetic fields on radical pair reformation
Bioelectromagnetics 20:255–263.
https://doi.org/10.1002/(SICI)1521-186X(1999)20:4<255::AID-BEM6>3.0.CO;2-W
- Google Scholar
1. Adair RK
(2000) Static and low-frequency magnetic field effects: health risks and therapies
Reports on Progress in Physics 63:415–454.
https://doi.org/10.1088/0034-4885/63/3/204
- Google Scholar
(2018) Geomagnetic field impacts on cryptochrome and phytochrome signaling
Journal of Photochemistry and Photobiology B: Biology 185:32–40.
https://doi.org/10.1016/j.jphotobiol.2018.05.027
- Google Scholar
1. Ahlbom A
2. Day N
3. Feychting M
4. Roman E
5. Skinner J
6. Dockerty J
7. Linet M
8. McBride M
9. Michaelis J
10. Olsen JH
11. Tynes T
12. Verkasalo PK
(2000) A pooled analysis of magnetic fields and childhood leukaemia
British Journal of Cancer 83:692–698.
https://doi.org/10.1054/bjoc.2000.1376
- PubMed
- Google Scholar
(2007) Quantum beats in radical pairs. russ
Chemical Reviews 76:493–506.
https://doi.org/10.1070/RC2007v076n06ABEH003715
- Google Scholar
1. Binhi V
(2008) Do naturally occurring magnetic nanoparticles in the human body mediate increased risk of childhood leukaemia with EMF exposure?
International Journal of Radiation Biology 84:569–579.
https://doi.org/10.1080/09553000802195323
- PubMed
- Google Scholar
Thesis
1. Blunsden J
(2018)
Effects of Weak Magnetic Fields on Radical Pair Reactions Part II Thesis

University of Oxford.
- Google Scholar
1. Bounds PL
2. Kuster N
(2015) Is cryptochrome a primary sensor of extremely low frequency magnetic fields in childhood leukemia?
Biophysical Journal 108:562a.
https://doi.org/10.1016/j.bpj.2014.11.3081
- Google Scholar
1. Brocklehurst B
(1976) Spin correlation in the geminate recombination of radical ions in hydrocarbons. Part 1.—Theory of the magnetic field effect
J. Chem. Soc., Faraday Trans. 2 72:1869–1884.
https://doi.org/10.1039/F29767201869
- Google Scholar
1. Brocklehurst B
(2002) Magnetic fields and radical reactions: recent developments and their role in nature
Chemical Society Reviews 31:301–311.
https://doi.org/10.1039/b107250c
- PubMed
- Google Scholar
1. Brocklehurst B
2. McLauchlan KA
(1996) Free radical mechanism for the effects of environmental electromagnetic fields on biological systems
International Journal of Radiation Biology 69:3–24.
https://doi.org/10.1080/095530096146147
- PubMed
- Google Scholar
(2010) Quantum control and entanglement in a chemical compass
Physical Review Letters 104:220502.
https://doi.org/10.1103/PhysRevLett.104.220502
- PubMed
- Google Scholar
1. Chaves I
2. Pokorny R
3. Byrdin M
4. Hoang N
5. Ritz T
6. Brettel K
7. Essen LO
8. van der Horst GT
9. Batschauer A
10. Ahmad M
(2011) The cryptochromes: blue light photoreceptors in plants and animals
Annual Review of Plant Biology 62:335–364.
https://doi.org/10.1146/annurev-arplant-042110-103759
- PubMed
- Google Scholar
1. Close J
(2012) Are stress responses to geomagnetic storms mediated by the cryptochrome compass system?
Proceedings of the Royal Society B: Biological Sciences 279:2081–2090.
https://doi.org/10.1098/rspb.2012.0324
- Google Scholar
1. Close JP
(2014a) The compass within the clock - Part 1: the hypothesis of magnetic fields as secondary zeitgebers to the circadian system—logical and scientific objections
Hypothesis 12:e1.
https://doi.org/10.5779/hypothesis.v12i1.359
- Google Scholar
1. Close JP
(2014b) The compass within the clock - Part 2: does cryptochrome radical-pair based signalling contribute to the temperature-robustness of circadian systems?
Hypothesis 12:e2.
https://doi.org/10.5779/hypothesis.v12i1.360
- Google Scholar
1. Crumpton MJ
(2005) The Bernal lecture 2004 are low-frequency electromagnetic fields a health hazard?
Philosophical Transactions of the Royal Society B: Biological Sciences 360:1223–1230.
https://doi.org/10.1098/rstb.2005.1663
- Google Scholar
1. Crumpton MJ
2. Collins AR
(2004) Are environmental electromagnetic fields genotoxic?
DNA Repair 3:1385–1387.
https://doi.org/10.1016/j.dnarep.2004.05.006
- PubMed
- Google Scholar
1. Efimova O
2. Hore PJ
(2008) Role of exchange and dipolar interactions in the radical pair model of the avian magnetic compass
Biophysical Journal 94:1565–1574.
https://doi.org/10.1529/biophysj.107.119362
- PubMed
- Google Scholar
1. Eichwald C
2. Walleczek J
(1996) Model for magnetic field effects on radical pair recombination in enzyme kinetics
Biophysical Journal 71:623–631.
https://doi.org/10.1016/S0006-3495(96)79263-9
- PubMed
- Google Scholar
1. Eichwald C
2. Walleczek J
(1997) Low-frequency-dependent effects of oscillating magnetic fields on radical pair recombination in enzyme kinetics
The Journal of Chemical Physics 107:4943–4950.
https://doi.org/10.1063/1.474858
- Google Scholar
1. Eichwald C
2. Walleczek J
(1998) Magnetic field perturbations as a tool for controlling enzyme-regulated and oscillatory biochemical reactions
Biophysical Chemistry 74:209–224.
https://doi.org/10.1016/S0301-4622(98)00180-X
- PubMed
- Google Scholar
1. Engels S
2. Schneider NL
3. Lefeldt N
4. Hein CM
5. Zapka M
6. Michalik A
7. Elbers D
8. Kittel A
9. Hore PJ
10. Mouritsen H
(2014) Anthropogenic electromagnetic noise disrupts magnetic compass orientation in a migratory bird
Nature 509:353–356.
https://doi.org/10.1038/nature13290
- PubMed
- Google Scholar
(2000) The effects of weak magnetic fields on radical recombination reactions in micelles
International Journal of Radiation Biology 76:1509–1522.
https://doi.org/10.1080/09553000050176270
- PubMed
- Google Scholar
Book
1. Frankevich EL
2. Kubarev SI
(1982)
Spectroscopy of reaction yield detected magnetic resonance

In: Clarke R. H, editors. Triplet State ODMR Spectroscopy. New York: Wiley. pp. 137–183.
- Google Scholar
1. Gauger EM
2. Rieper E
3. Morton JJ
4. Benjamin SC
5. Vedral V
(2011) Sustained quantum coherence and entanglement in the avian compass
Physical Review Letters 106:040503.
https://doi.org/10.1103/PhysRevLett.106.040503
- PubMed
- Google Scholar
Book
1. Gerson F
2. Huber W
(2003) Electron Spin Resonance Spectroscopy of Organic Radicals
Weinheim: Wiley-VCH.
https://doi.org/10.1002/3527601627
- Google Scholar
1. Greenland S
2. Sheppard AR
3. Kaune WT
4. Poole C
5. Kelsh MA
(2000) A pooled analysis of magnetic fields, wire codes, and childhood leukemia
Epidemiology 11:624–634.
https://doi.org/10.1097/00001648-200011000-00003
- Google Scholar
1. Günther A
2. Einwich A
3. Sjulstok E
4. Feederle R
5. Bolte P
6. Koch KW
7. Solov'yov IA
8. Mouritsen H
(2018) Double-cone localization and seasonal expression pattern suggest a role in magnetoreception for European robin cryptochrome 4
Current Biology 28:211–223.
https://doi.org/10.1016/j.cub.2017.12.003
- PubMed
- Google Scholar
(2014) Circadian pacemaking in cells and circuits of the suprachiasmatic nucleus
Journal of Neuroendocrinology 26:2–10.
https://doi.org/10.1111/jne.12125
- PubMed
- Google Scholar
1. Henbest KB
2. Kukura P
3. Rodgers CT
4. Hore PJ
5. Timmel CR
(2004) Radio frequency magnetic field effects on a radical recombination reaction: a diagnostic test for the radical pair mechanism
Journal of the American Chemical Society 126:8102–8103.
https://doi.org/10.1021/ja048220q
- PubMed
- Google Scholar
1. Hiscock HG
2. Worster S
3. Kattnig DR
4. Steers C
5. Jin Y
6. Manolopoulos DE
7. Mouritsen H
8. Hore PJ
(2016) The quantum needle of the avian magnetic compass
PNAS 113:4634–4639.
https://doi.org/10.1073/pnas.1600341113
- PubMed
- Google Scholar
(2017) Disruption of magnetic compass orientation in migratory birds by radiofrequency electromagnetic fields
Biophysical Journal 113:1475–1484.
https://doi.org/10.1016/j.bpj.2017.07.031
- PubMed
- Google Scholar
(2009) Possible involvement of superoxide and dioxygen with cryptochrome in avian magnetoreception: origin of Zeeman resonances observed by in vivo EPR spectroscopy
Chemical Physics Letters 480:118–122.
https://doi.org/10.1016/j.cplett.2009.08.051
- Google Scholar
(2012) Entanglement and sources of magnetic anisotropy in radical pair-based avian magnetoreceptors
Physical Review Letters 109:220501.
https://doi.org/10.1103/PhysRevLett.109.220501
- PubMed
- Google Scholar
1. Hore PJ
(2012) Are biochemical reactions affected by weak magnetic fields?
PNAS 109:1357–1358.
https://doi.org/10.1073/pnas.1120531109
- PubMed
- Google Scholar
1. Hore PJ
2. Mouritsen H
(2016) The radical-pair mechanism of magnetoreception
Annual Review of Biophysics 45:299–344.
https://doi.org/10.1146/annurev-biophys-032116-094545
- PubMed
- Google Scholar
Book
1. IARC
(2002)
Static and Extremely Low-Frequency (ELF) Electric and Magnetic Fields, 80

Lyon, France: International Agency for Research on Cancer.
- Google Scholar
1. Jones AR
(2016) Magnetic field effects in proteins
Molecular Physics 114:1691–1702.
https://doi.org/10.1080/00268976.2016.1149631
- Google Scholar
(2006) Do extremely low frequency magnetic fields enhance the effects of environmental carcinogens? A meta-analysis of experimental studies
International Journal of Radiation Biology 82:1–12.
https://doi.org/10.1080/09553000600577839
- PubMed
- Google Scholar
(2018) Magnetocarcinogenesis: is there a mechanism for carcinogenic effects of weak magnetic fields?
Proceedings of the Royal Society B: Biological Sciences 285:20180590.
https://doi.org/10.1098/rspb.2018.0590
- Google Scholar
(2009) Magnetic field effect in the reaction of recombination of nitric oxide and superoxide anion
Applied Magnetic Resonance 36:195–208.
https://doi.org/10.1007/s00723-009-0018-2
- Google Scholar
(2011) Kinetic magnetic-field effect involving the small biologically relevant inorganic radicals NO and O₂^•−
ChemPhysChem 12:1714–1728.
https://doi.org/10.1002/cphc.201100178
- PubMed
- Google Scholar
1. Kattnig DR
2. Evans EW
3. Déjean V
4. Dodson CA
5. Wallace MI
6. Mackenzie SR
7. Timmel CR
8. Hore PJ
(2016a) Chemical amplification of magnetic field effects relevant to avian magnetoreception
Nature Chemistry 8:384–391.
https://doi.org/10.1038/nchem.2447
- PubMed
- Google Scholar
(2016b) Electron spin relaxation in cryptochrome-based magnetoreception
Physical Chemistry Chemical Physics 18:12443–12456.
https://doi.org/10.1039/C5CP06731F
- PubMed
- Google Scholar
(2016c) Electron spin relaxation can enhance the performance of a cryptochrome-based magnetic compass sensor
New Journal of Physics 18:063007.
https://doi.org/10.1088/1367-2630/18/6/063007
- Google Scholar
1. Kattnig DR
(2017) Radical-pair-based magnetoreception amplified by radical scavenging: resilience to spin relaxation
The Journal of Physical Chemistry B 121:10215–10227.
https://doi.org/10.1021/acs.jpcb.7b07672
- PubMed
- Google Scholar
1. Kattnig DR
2. Hore PJ
(2017) The sensitivity of a radical pair compass magnetoreceptor can be significantly amplified by radical scavengers
Scientific Reports 7:11640.
https://doi.org/10.1038/s41598-017-09914-7
- PubMed
- Google Scholar
1. Kavokin KV
(2009) The puzzle of magnetic resonance effect on the magnetic compass of migratory birds
Bioelectromagnetics 30:402–410.
https://doi.org/10.1002/bem.20485
- PubMed
- Google Scholar
(2014) Magnetic orientation of garden warblers (Sylvia borin) under 1.4 MHz radiofrequency magnetic field
Journal of the Royal Society Interface 11:20140451.
https://doi.org/10.1098/rsif.2014.0451
- Google Scholar
(2014) Circadian molecular clocks and cancer
Cancer Letters 342:9–18.
https://doi.org/10.1016/j.canlet.2013.09.040
- PubMed
- Google Scholar
1. Kheifets L
2. Ahlbom A
3. Crespi CM
4. Draper G
5. Hagihara J
6. Lowenthal RM
7. Mezei G
8. Oksuzyan S
9. Schüz J
10. Swanson J
11. Tittarelli A
12. Vinceti M
13. Wunsch Filho V
(2010) Pooled analysis of recent studies on magnetic fields and childhood leukaemia
British Journal of Cancer 103:1128–1135.
https://doi.org/10.1038/sj.bjc.6605838
- PubMed
- Google Scholar
1. Kilcoyne SH
2. Cywinski R
(1995) Ferritin: a model superparamagnet
Journal of Magnetism and Magnetic Materials 140-144:1466–1467.
https://doi.org/10.1016/0304-8853(94)00626-1
- Google Scholar
1. Kobayashi Y
2. Ishikawa T
3. Hirayama J
4. Daiyasu H
5. Kanai S
6. Toh H
7. Fukuda I
8. Tsujimura T
9. Terada N
10. Kamei Y
11. Yuba S
12. Iwai S
13. Todo T
(2000) Molecular analysis of zebrafish photolyase/cryptochrome family: two types of cryptochromes present in zebrafish
Genes to Cells 5:725–738.
https://doi.org/10.1046/j.1365-2443.2000.00364.x
- PubMed
- Google Scholar
1. Kominis IK
(2012) Magnetic sensitivity and entanglement dynamics of the chemical compass
Chemical Physics Letters 542:143–146.
https://doi.org/10.1016/j.cplett.2012.06.014
- Google Scholar
1. Krylov VV
(2017) Biological effects related to geomagnetic activity and possible mechanisms
Bioelectromagnetics 38:497–510.
https://doi.org/10.1002/bem.22062
- PubMed
- Google Scholar
1. Kubo Y
2. Akiyama M
3. Fukada Y
4. Okano T
(2006) Molecular cloning, mRNA expression, and immunocytochemical localization of a putative blue-light photoreceptor CRY4 in the chicken pineal gland
Journal of Neurochemistry 97:1155–1165.
https://doi.org/10.1111/j.1471-4159.2006.03826.x
- PubMed
- Google Scholar
1. Kubo Y
2. Takeuchi T
3. Okano K
4. Okano T
(2010) Cryptochrome genes are highly expressed in the ovary of the african clawed frog, Xenopus tropicalis
PLOS ONE 5:e9273.
https://doi.org/10.1371/journal.pone.0009273
- PubMed
- Google Scholar
(2017) Vertebrate cryptochromes are vestigial flavoproteins
Scientific Reports 7:44906.
https://doi.org/10.1038/srep44906
- PubMed
- Google Scholar
(1998) Biological responses to electromagnetic fields
The FASEB Journal 12:395–420.
https://doi.org/10.1096/fasebj.12.6.395
- PubMed
- Google Scholar
(2011) ELF magnetic fields: animal studies, mechanisms of action
Progress in Biophysics and Molecular Biology 107:369–373.
https://doi.org/10.1016/j.pbiomolbio.2011.09.003
- PubMed
- Google Scholar
1. Landler L
2. Keays DA
(2018) Cryptochrome: the magnetosensor with a sinister side?
PLOS Biology 16:e3000018.
https://doi.org/10.1371/journal.pbio.3000018
- PubMed
- Google Scholar
1. Lee AA
2. Lau JCS
3. Hogben HJ
4. Biskup T
5. Kattnig DR
6. Hore PJ
(2014) Alternative radical pairs for cryptochrome-based magnetoreception
Journal of the Royal Society Interface 11:20131063.
https://doi.org/10.1098/rsif.2013.1063
- Google Scholar
(2014) Asymmetric recombination and electron spin relaxation in the semiclassical theory of radical pair reactions
The Journal of Chemical Physics 141:044111.
https://doi.org/10.1063/1.4890659
- PubMed
- Google Scholar
1. Liboff AR
(2014) Why are living things sensitive to weak magnetic fields?
Electromagnetic Biology and Medicine 33:241–245.
https://doi.org/10.3109/15368378.2013.809579
- PubMed
- Google Scholar
1. Liedvogel M
2. Mouritsen H
(2010) Cryptochromes—a potential magnetoreceptor: what do we know and what do we want to know?
Journal of the Royal Society Interface 7:S147–S162.
https://doi.org/10.1098/rsif.2009.0411.focus
- Google Scholar
1. Maeda K
2. Henbest KB
3. Cintolesi F
4. Kuprov I
5. Rodgers CT
6. Liddell PA
7. Gust D
8. Timmel CR
9. Hore PJ
(2008) Chemical compass model of avian magnetoreception
Nature 453:387–390.
https://doi.org/10.1038/nature06834
- PubMed
- Google Scholar
1. Maeda K
2. Robinson AJ
3. Henbest KB
4. Hogben HJ
5. Biskup T
6. Ahmad M
7. Schleicher E
8. Weber S
9. Timmel CR
10. Hore PJ
(2012) Magnetically sensitive light-induced reactions in cryptochrome are consistent with its proposed role as a magnetoreceptor
PNAS 109:4774–4779.
https://doi.org/10.1073/pnas.1118959109
- PubMed
- Google Scholar
1. Mattsson MO
2. Simkó M
(2014) Grouping of experimental conditions as an approach to evaluate effects of extremely low-frequency magnetic fields on oxidative response in in vitro studies
Frontiers in Public Health 2:132.
https://doi.org/10.3389/fpubh.2014.00132
- PubMed
- Google Scholar
(2015) Magnetic field effects as a result of the radical pair mechanism are unlikely in redox enzymes
Journal of the Royal Society Interface 12:20141155.
https://doi.org/10.1098/rsif.2014.1155
- Google Scholar
1. Mitsui H
2. Maeda T
3. Yamaguchi C
4. Tsuji Y
5. Watari R
6. Kubo Y
7. Okano K
8. Okano T
(2015) Overexpression in yeast, photocycle, and in vitro structural change of an avian putative magnetoreceptor cryptochrome 4
Biochemistry 54:1908–1917.
https://doi.org/10.1021/bi501441u
- PubMed
- Google Scholar
1. Møller AC
2. Olsen LF
(1999) Effect of magnetic fields on an oscillating enzyme reaction
Journal of the American Chemical Society 121:6351–6354.
https://doi.org/10.1021/ja990834l
- Google Scholar
(2000) Further studies of the effect of magnetic fields on the oscillating peroxidase–oxidase reaction
Physical Chemistry Chemical Physics 2:3443–3446.
https://doi.org/10.1039/b003641m
- Google Scholar
1. Møller AC
2. Olsen LF
(2000) Perturbations of simple oscillations and complex dynamics in the peroxidase−oxidase reaction using magnetic fields
The Journal of Physical Chemistry B 104:140–146.
https://doi.org/10.1021/jp993284m
- Google Scholar
1. Mouritsen H
(2018) Long-distance navigation and magnetoreception in migratory animals
Nature 558:50–59.
https://doi.org/10.1038/s41586-018-0176-1
- PubMed
- Google Scholar
(2017) Ascorbic acid may not be involved in cryptochrome-based magnetoreception
Journal of the Royal Society Interface 14:20170657.
https://doi.org/10.1098/rsif.2017.0657
- Google Scholar
(2017) Magnetoreception-a sense without a receptor
PLOS Biology 15:e2003234.
https://doi.org/10.1371/journal.pbio.2003234
- PubMed
- Google Scholar
1. Ozturk N
2. Selby CP
3. Song SH
4. Ye R
5. Tan C
6. Kao YT
7. Zhong D
8. Sancar A
(2009) Comparative photochemistry of animal Type 1 and Type 4 cryptochromes
Biochemistry 48:8585–8593.
https://doi.org/10.1021/bi901043s
- PubMed
- Google Scholar
1. Pauls JA
2. Zhang Y
3. Berman GP
4. Kais S
(2013) Quantum coherence and entanglement in the avian compass
Physical Review E 87:062704.
https://doi.org/10.1103/PhysRevE.87.062704
- Google Scholar
1. Purtov PA
(2004) External magnetic fields as a possible cause of stability disturbance of stationary states far from equilibrium in reactions involving radical pairs
Applied Magnetic Resonance 26:83–97.
https://doi.org/10.1007/BF03166564
- Google Scholar
1. Qin S
2. Yin H
3. Yang C
4. Dou Y
5. Liu Z
6. Zhang P
7. Yu H
8. Huang Y
9. Feng J
10. Hao J
11. Hao J
12. Deng L
13. Yan X
14. Dong X
15. Zhao Z
16. Jiang T
17. Wang HW
18. Luo SJ
19. Xie C
(2016) A magnetic protein biocompass
Nature Materials 15:217–226.
https://doi.org/10.1038/nmat4484
- PubMed
- Google Scholar
1. Ritz T
2. Adem S
3. Schulten K
(2000) A model for photoreceptor-based magnetoreception in birds
Biophysical Journal 78:707–718.
https://doi.org/10.1016/S0006-3495(00)76629-X
- PubMed
- Google Scholar
(2004) Resonance effects indicate a radical-pair mechanism for avian magnetic compass
Nature 429:177–180.
https://doi.org/10.1038/nature02534
- PubMed
- Google Scholar
1. Ritz T
2. Wiltschko R
3. Hore PJ
4. Rodgers CT
5. Stapput K
6. Thalau P
7. Timmel CR
8. Wiltschko W
(2009) Magnetic compass of birds is based on a molecule with optimal directional sensitivity
Biophysical Journal 96:3451–3457.
https://doi.org/10.1016/j.bpj.2008.11.072
- PubMed
- Google Scholar
1. Rodgers CT
2. Norman SA
3. Henbest KB
4. Timmel CR
5. Hore PJ
(2007) Determination of radical re-encounter probability distributions from magnetic field effects on reaction yields
Journal of the American Chemical Society 129:6746–6755.
https://doi.org/10.1021/ja068209l
- PubMed
- Google Scholar
1. Rodgers CT
(2009) Magnetic field effects in chemical systems
Pure and Applied Chemistry 81:19–43.
https://doi.org/10.1351/PAC-CON-08-10-18
- Google Scholar
1. Rodgers CT
2. Hore PJ
(2009) Chemical magnetoreception in birds: the radical pair mechanism
PNAS 106:353–360.
https://doi.org/10.1073/pnas.0711968106
- PubMed
- Google Scholar
1. Sakaguchi Y
2. Hayashi H
(1984) Laser-photolysis study of the photochemical reactions of naphthoquinones in a sodium dodecyl sulfate micelle under high magnetic fields
The Journal of Physical Chemistry 88:1437–1440.
https://doi.org/10.1021/j150651a040
- Google Scholar
Book
(1984)
Spin Polarization and Magnetic Field Effects in Radical Reactions

New York: Elsevier.
- Google Scholar
(1994a)
Model for the rationalization of magnetic field effects in vivo. Application of the radical pair mechanism to biological systems

Photochemistry and Photobiology 59:585–589.
- PubMed
- Google Scholar
(1994b) Application of the radical pair mechanism to free radicals in organized systems: can the effects of 60 Hz be predicted from studies under static fields?
Bioelectromagnetics 15:549–554.
https://doi.org/10.1002/bem.2250150608
- Google Scholar
(1978) A biomagnetic sensory mechanism based on magnetic field modulated coherent electron spin motion
Zeitschrift Für Physikalische Chemie 111:1–5.
https://doi.org/10.1524/zpch.1978.111.1.001
- Google Scholar
1. Schwarze S
2. Schneider NL
3. Reichl T
4. Dreyer D
5. Lefeldt N
6. Engels S
7. Baker N
8. Hore PJ
9. Mouritsen H
(2016) Weak broadband electromagnetic fields are more disruptive to magnetic compass orientation in a night-migratory songbird (Erithacus rubecula) than strong narrow-band fields
Frontiers in Behavioral Neuroscience 10:55.
https://doi.org/10.3389/fnbeh.2016.00055
- PubMed
- Google Scholar
1. Sheppard DMW
2. Li J
3. Henbest KB
4. Neil SRT
5. Maeda K
6. Storey J
7. Schleicher E
8. Biskup T
9. Rodriguez R
10. Weber S
11. Hore PJ
12. Timmel CR
13. Mackenzie SR
(2017) Millitesla magnetic field effects on the photocycle of Drosophila melanogaster cryptochrome
Scientific Reports 7:42228.
https://doi.org/10.1038/srep42228
- Google Scholar
1. Shushin AI
(1991) The effect of the spin exchange interaction on SNP and RYDMR spectra of geminate radical pairs
Chemical Physics Letters 181:274–278.
https://doi.org/10.1016/0009-2614(91)90366-H
- Google Scholar
(1995) Manifestation of quantum coherence upon recombination of radical ion pairs in weak magnetic fields. systems with non-equivalent nuclei
Chemical Physics Letters 233:444–450.
https://doi.org/10.1016/0009-2614(94)01489-I
- Google Scholar
1. Steiner UE
2. Ulrich T
(1989) Magnetic field effects in chemical kinetics and related phenomena
Chemical Reviews 89:51–147.
https://doi.org/10.1021/cr00091a003
- Google Scholar
1. Swanson J
(1994) Measurements of static magnetic fields in homes in the UK and their implication for epidemiological studies of exposure to alternating magnetic fields
Journal of Radiological Protection 14:67–75.
https://doi.org/10.1088/0952-4746/14/1/007
- Google Scholar
1. Swanson J
2. Kheifets L
(2006) Biophysical mechanisms: a component in the weight of evidence for health effects of power-frequency electric and magnetic fields
Radiation Research 165:470–478.
https://doi.org/10.1667/RR3522.1
- PubMed
- Google Scholar
1. Swanson J
2. Kheifets L
(2012) Could the geomagnetic field be an effect modifier for studies of power-frequency magnetic fields and childhood leukaemia?
Journal of Radiological Protection 32:413–418.
https://doi.org/10.1088/0952-4746/32/4/413
- PubMed
- Google Scholar
1. Thalau P
2. Ritz T
3. Stapput K
4. Wiltschko R
5. Wiltschko W
(2005) Magnetic compass orientation of migratory birds in the presence of a 1.315 MHz oscillating field
Naturwissenschaften 92:86–90.
https://doi.org/10.1007/s00114-004-0595-8
- PubMed
- Google Scholar
(2014) Approaches to measuring entanglement in chemical magnetometers
The Journal of Physical Chemistry A 118:13–20.
https://doi.org/10.1021/jp408569d
- PubMed
- Google Scholar
(1998) The influence of very small magnetic fields on radical recombination reactions in the limit of slow recombination
Chemical Physics Letters 298:7–14.
https://doi.org/10.1016/S0009-2614(98)01158-0
- Google Scholar
(1998) Effects of weak magnetic fields on free radical recombination reactions
Molecular Physics 95:71–89.
https://doi.org/10.1080/00268979809483134
- Google Scholar
(2001) Model calculations of magnetic field effects on the recombination reactions of radicals with anisotropic hyperfine interactions
Chemical Physics Letters 334:387–395.
https://doi.org/10.1016/S0009-2614(00)01436-6
- Google Scholar
1. Turro NJ
2. Weed GC
(1983) Micellar systems as supercages for reactions of geminate radical pairs. magnetic effects
Journal of the American Chemical Society 105:1861–1868.
https://doi.org/10.1021/ja00345a031
- Google Scholar
(2015) Could magnetic fields affect the circadian clock function of cryptochromes? Testing the basic premise of the cryptochrome hypothesis (ELF magnetic fields)
Health Physics 109:84–89.
https://doi.org/10.1097/HP.0000000000000292
- PubMed
- Google Scholar
1. Vohralik PF
2. Lam SKH
(2009) NanoSQUID detection of magnetization from ferritin nanoparticles
Superconductor Science and Technology 22:064007.
https://doi.org/10.1088/0953-2048/22/6/064007
- Google Scholar
1. Wang X
2. Jing C
3. Selby CP
4. Chiou YY
5. Yang Y
6. Wu W
7. Sancar A
8. Wang J
(2018) Comparative properties and functions of Type 2 and Type 4 pigeon cryptochromes
Cellular and Molecular Life Sciences 75:4629–4641.
https://doi.org/10.1007/s00018-018-2920-y
- PubMed
- Google Scholar
Book
1. WHO
(2007)
Extremely Low Frequency Fields

Geneva, Switzerland: World Health Organization.
- Google Scholar
1. Wiltschko R
2. Thalau P
3. Gehring D
4. Niessner C
5. Ritz T
6. Wiltschko W
(2015) Magnetoreception in birds: the effect of radio-frequency fields
Journal of the Royal Society Interface 12:20141103.
https://doi.org/10.1098/rsif.2014.1103
- Google Scholar
(2013) Avian magnetic compass can be tuned to anomalously low magnetic intensities
Proceedings of the Royal Society B: Biological Sciences 280:20130853.
https://doi.org/10.1098/rspb.2013.0853
- Google Scholar
(2001) Radio frequency magnetic field effects on electron-hole recombination
Physical Review Letters 87:077602.
https://doi.org/10.1103/PhysRevLett.87.077602
- PubMed
- Google Scholar
(2016) Spin relaxation of radicals in cryptochrome and its role in avian magnetoreception
The Journal of Chemical Physics 145:035104.
https://doi.org/10.1063/1.4958624
- PubMed
- Google Scholar
1. Worster SB
2. Hore PJ
(2018) Proposal to use superparamagnetic nanoparticles to test the role of cryptochrome in magnetoreception
Journal of the Royal Society Interface 15:20180587.
https://doi.org/10.1098/rsif.2018.0587
- Google Scholar
1. Zhang Y
2. Berman GP
3. Kais S
(2014) Sensitivity and entanglement in the avian chemical compass
Physical Review E 90:042707.
https://doi.org/10.1103/PhysRevE.90.042707
- Google Scholar

Article and author information

Author details

PJ Hore

Department of Chemistry, Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford, United Kingdom

Contribution
conducted all aspects of the work, from conceptualisation through to writing, and revised the final article

For correspondence
peter.hore@chem.ox.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8863-570X

Funding

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

I am grateful to Jukka Juutilainen, Jim Metcalfe, and John Swanson for helpful comments on the manuscript and to the Electromagnetic Fields Biological Research Trust and its Scientific Advisory Committee for financial support and stimulating discussions over a period of many years.

Version history

Received: December 6, 2018
Accepted: February 2, 2019
Accepted Manuscript published: February 25, 2019 (version 1)
Version of Record published: March 14, 2019 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

2,646

views
292

downloads
35

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

PJ Hore

(2019)

Upper bound on the biological effects of 50/60 Hz magnetic fields mediated by radical pairs

eLife 8:e44179.

https://doi.org/10.7554/eLife.44179

Share this article

Cite this article

The radical pair mechanism.

Effects of a static magnetic field on a model of the [FAD•− TrpH•+] radical pair in cryptochrome.

Schematic representations (thick black lines) of the dependence of the reaction yield, ΦT, on the strength of a static magnetic field.

Magnetic field effects (in ppm) on model radical pairs as a function of the lifetime of the radical pair, τ, in the range 1 ns−1 ms.

Changes in reaction product yields (in ppm) arising from small changes in magnetic field strengths and temperature

Author details

PJ Hore

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Effects of a static magnetic field on a model of the [FAD^•− TrpH^•+] radical pair in cryptochrome.

Schematic representations (thick black lines) of the dependence of the reaction yield, $Φ_{T}$ , on the strength of a static magnetic field.