Radical pair model for magnetic field effects on NMDA receptor activity

The N-methyl-D-aspartate receptor is a prominent player in brain development and functioning. Perturbations to its functioning through external stimuli like magnetic fields can potentially affect the brain in numerous ways. Various studies have shown that magnetic fields of varying strengths affect these receptors. We propose that the radical pair mechanism, a quantum mechanical process, could explain some of these field effects. Radicals of the form \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[\text{RO}^\cdot \text{Mg(H}_2\text{O})_n^{+\cdot }]$$\end{document}[RO·Mg(H2O)n+·], where R is a protein residue that can be Serine or Tyrosine, are considered for this study. The variation in the singlet fractional yield of the radical pairs, as a function of magnetic field strength, is calculated to understand how the magnetic field affects the products of the radical pair reactions. Based on the results, the radical pair mechanism is a likely candidate for explaining the magnetic field effects observed on the receptor activity. The model predicts changes in the behaviour of the system as magnetic field strength is varied and also predicts certain isotope effects. The results further suggest that similar effects on radical pairs could be a plausible explanation for various magnetic field effects within the brain.


Introduction
Magnetic field effects (MFEs) in biological systems have garnered considerable attention from the academic community.This has been studied in the context of numerous systems, such as magnetoreception 1 , the working of the circadian clock 2 , genetics 3 , and various others [4][5][6][7][8][9][10][11][12][13][14][15] .In the last few decades, there has been significant progress in modelling different MFEs in these systems based on quantum mechanical phenomena 16 .One model that has been quite promising in explaining some of these effects is the radical pair mechanism (RPM) 17 .This mechanism explains the behavior of radical pair electrons, which are pairs of correlated electrons formed on separate molecules from the rupture of chemical bonds or through electron hopping.The RPM model has been most commonly applied in explaining avian magnetoreception in migratory birds 18 .It has also been used to model magnetoreception in other migratory organisms 19 , as well.
There are numerous phenomena involving MFEs in biological systems that have yet to be explained by any mechanisms.One of the most fascinating areas to study such phenomena is the brain.It has been found that exposure to magnetic fields affects various functions in the brain.The frequency of neuron firing, pain sensitivity and inhibition, production of melatonin, and functioning of the pineal gland are some examples of these effects 6 .
Understanding the effects of magnetic fields on the brain is especially important since magnetic stimulations are already a well-known and well-used method for non-invasive brain stimulation used as both diagnostic and therapeutic tools [20][21][22] .Since magnetic stimulations can modulate various brain functions, they also have the potential to be used as a treatment for numerous neuropsychiatric diseases 23,24 .Understanding the mechanism by which magnetic fields affect important parts of the brain could potentially lead to further advances in such non-invasive treatment tools.
Despite the existence of a large number of studies on these phenomena, there still needs to be more understanding of the mechanisms by which most of these effects occur.However, the RPM model has recently been applied to studying various phenomena within the brain with promising results.According to this mechanism, the electrons on the RP involved evolves under the influence of various interactions, like the Zeeman and hyperfine interactions, into a superposition of singlet and triplet states.The rate and yield of the final products formed from these states are affected by external magnetic fields, which is the cause of the various experimental observations.Lithium effects on hyperactivity 25 , magnetic field and lithium effects on the working of the circadian clock 26 , hypomagnetic field effects on neurogenesis 27 , hypomagnetic field effects on microtubule reorganization 28 , and Xenon-induced general anesthesia 29 are examples of the various phenomena that have been modelled using this mechanism.These results encourage us to apply the mechanism to understand other phenomena within the brain.
One of the examples of such phenomena is the MFEs on the functioning of the N-Methyl-D-Aspartate receptors (NMDARs) within the brain.Several experimental studies successfully show that both static and oscillating magnetic fields of different strenghts can affect the activity of NMDARs 30,31 .These receptors are critical in all stages of development in higher organisms and are heavily involved in various brain functions, such as neuronal development and synaptic plasticity 32 .The wide range of roles this receptor plays in the brain makes it a rich subject for research, with numerous studies dedicated to understanding its structure, function, and physiological roles over the years 33,34 .The variations in the different subunits that form this receptor increase the complexity and functionality of the various forms of the receptor 35,36 .As these receptors are deeply involved in the development and functioning of the brain, they are also crucial players in numerous neurological and psychological disorders.NMDARs have been shown to be involved in ischemic strokes, schizophrenia, excitotoxic brain injury, memory and learning impairments, and various other disorders [37][38][39][40] , thus making it extremely important that we understand the effects of external MFEs on these receptors and how they affect different processes within the brain.
Hirai et al. demonstrated in two separate studies 41,42 that immature rat hippocampal cells cultured in the presence of a 100 mT sustained or repetitive static magnetic field exhibit an increase in intracellular free Ca 2+ ion concentration upon increasing NMDA concentration.In both experiments, upregulation of certain mRNAs that encode some of the NMDAR subunits was observed.
We propose that an RPM model could provide an explanation for the observed effects.The MFE observed in these experiments may be occurring at any level, ranging from the transcription of the NMDAR subunits mRNAs to a direct effect on the functioning of the receptor.We attempt to identify the possible radical pairs involved in the various reaction pathways in these levels that may be affected by the static magnetic field (SMF).The potential RPs are identified based on the experimental evidence presented previously and an understanding of the mechanisms involved in the formation and function of the NMDAR.We study the effects of the magnetic field by calculating the product yield of the possible radical pair reactions involved in the system.On the basis of our findings, we conclude that an RPM model could provide a plausible explanation for the magnetic field effects observed on the NMDAR activity.This also provides further indication that radical pairs could be a key player in various magnetic field effects that occur in the brain.We further predict variations in the behaviour of the system as the magnetic field strength varies, which could be experimentally studied.We also show that there may be certain isotope dependent effects that could also be of experimental interest.

Results from Previous Experiments
Hirai et al. 41 conducted an experiment in which immature rat hippocampal cells were cultured for three days in vitro (DIV) in the presence of 100 mT SMF.Following the exposure, an increase in the intake of Ca 2+ ions was observed as the NMDA concentration within the cells was increased.of sustained exposure to SMF on the extent of Ca 2+ intake into rat hippocampal cells.Immature rat hippocampal cells were cultured for three days in vitro (DIV) in the presence of 100 mT SMF for this experiment.The values in the table are expressed as percentages over the maximal reading found at the end of the experiments from an addition of 10 µM of Ca 2+ ionophore(A23187).

NMDA
In a separate study conducted by Hirai et al. 42 , immature rat hippocampal cells were repeatedly exposed to 100 mT SMF for 15 minutes per day for eight days.Upon harvesting the cells 24 hours after the final exposure, it was found that there was once again a significant increase in intracellular Ca 2+ ion concentration.An upregulation in the mRNA expression of just the NR2B subunit of NMDAR by a factor (approximate) of 1.4 was also reported under the effect of such an exposure.
Tables 1 and 2 42 of repeated exposure to SMF on the extend of Ca 2+ intake into rat hippocampal cells.Immature rat hippocampal cells were repeatedly exposed to 100 mT SMF for 15 minutes per day for eight days for this experiment.The values in the table are expressed as percentages over the maximal reading found at the end of the experiments from an addition of 10 µM of Ca 2+ ionophore(A23187).
From the same experiments, Hirai et al. also reported an increase in the binding activity of Activator Protein-1(AP-1), which is a transcription factor that has been shown to be an active regulator of NR2B transcription 43 .

Possible Role of Phosphorylation in MFE on NMDARs
The first possible explanation for the MFE is an increase in NMDAR activity as a result of a rise in the number of mRNAs that form these receptors.Prybylowski et al. 44 demonstrated that an increase in the number of NR2 subunits within cerebellar granule cells causes an increase in the concentration of NMDARs within the cells.The same is most likely true for hippocampal cells.
As stated previously, an increase in the NR2B mRNA concentration is observed within the cells in the presence of an SMF 41,42 .It is highly likely that an increase in mRNA concentration would also result in an increase in NR2B subunit concentration within the cells.According to the findings of Prybylowski et al., the increase in NR2B subunits should result in an increase in the number of NMDARs in the cell.An increase in the number of NMDARs provides more channels for the Ca 2+ ions to enter the cells, thus causing an increase in the concentration of these ions within the cells, as the concentration of NMDA in the culture increases.
The upregulation in NR2B mRNA expression observed in the experiments could have occurred as the result of an MFE occurring at any stage of the mRNA formation.The increase in AP-1 binding activity reported by Hirai et al. could be a possible explanation for the same.AP-1 transcription factors typically consist of heterodimers of Jun and Fos proteins or homodimers of Jun proteins.The DNA binding activity of AP-1 is regulated by phosphorylation and dephosphorylation processes of the c-Jun protein, which belongs to the Jun protein family.Phosphorylation of Ser 63 and Ser 73 in the N-terminal domain of c-Jun stimulates the transcriptional activity of the protein 45,46 .Phosphorylation has been reported to be affected by MFEs [47][48][49] , making this process a likely candidate for the SMF effects.
As mentioned previously, it is also possible that the MFE occurs directly on the NMDAR or its subunits.Phosphorylation of serine and tyrosine sites in NMDAR subunits has been shown to regulate NMDAR activity 50,51 .For instance, phosphorylation of Ser 896 and Ser 897 together increases NMDAR surface expression.Receptor currents increase in NMDARs containing NR2A subunits when Ser 1291 and Ser 1312 are phosphorylated.Tyrosine phosphorylation of NR2A also potentiates NMDAR currents.These phosphorylation processes can also be influenced by MFEs.
We propose that the MFEs reported on NMDAR activity may be modeled by an RPM model based on the effects of SMF on one of the aforementioned processes.The MFE could occur via either a Serine phosphorylation, which will be referred to as the "Ser Pathway," or a tyrosine phosphorylation, which will be referred to as the "Tyr Pathway."

RPM Model of Oxyradicals and Hydrated Magnesium Cations
Buchachenko et al. 52 theorized that the transfer of phosphate groups to proteins induced by protein kinases may involve an ion-radical mechanism.To confirm this, they studied the catalysis of phosphorylation of prothrombin by prothrombin kinases with 24 Mg 2+ and 25 Mg 2+ ions.They successfully demonstrated an isotope dependence in the kinase's efficiency, which could be due to the nuclear spin of the Mg 2+ isotopes.It was proposed that this indicates that the ion-radical mechanism could be a plausible mechanism for enzymatic phosphorylation processes.The mechanism involves the following reactions: According to this mechanism, phosphorylation could involve an oxyradical, RO • (where R represents a protein residue), and a hydrated magnesium cation radical, Mg(H 2 O) n +• .
Depending on the pathway involved, the protein residue for our system could be either Serine or Tyrosine.Based on this, we propose that the MFEs can be modeled by an RPM model involving the RP [RO • Mg(H 2 O) n +• ] formed during this reaction.
The reactions involved in the transcriptional level phosphorylation processes involve only the Ser Pathway, whereas the Ser or Tyr Pathways may be involved in the MFEs on NMDAR subunits.We model the MFEs on the RPs using a simplified Hamiltonian consisting only of Zeeman and hyperfine interactions (HFIs),which is a good approximation for sufficiently distant radicals 18,53 .We assume that radical electrons have the same g-value as free electrons, which is a good approximation at the field values we are interested in.We also consider only the isotropic Fermi contact contribution for the HFIs, as the molecular arrangement is likely to be random.Therefore, the Hamiltonian would have the following form: Here, ŜA and ŜB are the electron spin operators of the RPs labelled A and B, respectively.The Larmour frequency of the electrons due to Zeeman interaction is given by ω.ÎA represents the nuclear spin operator associated with the oxyradical nucleus that contributes most to the HFI.
We perform our calculations considering the contributions from the natural abundance of the different magnesium isotopes.It is well known that 10% of naturally occurring magnesium is composed of the isotope 25 Mg, which is a spinful isotope.In the equation above, ÎB represents the nuclear spin operator of any naturally occurring 25 Mg nucleus present in the system.The terms a A and a B are the hyperfine coupling constants for each of the nuclei.
We use density functional theory (DFT) calculations to determine the hyperfine coupling constant (HFCC) values for all spinful nuclei involved.The details for the same have been given in the Methods section.For our calculations, we consider the HFI contribution from only the nucleus with the highest HFCC, which is a commonly used approximation.In the case of the Serine oxyradical, we find that the highest HF contribution is from a hydrogen nucleus with a A = 7.45 mT.For the Tyrosine oxyradical, the highest HFCC is of one of the hydrogen nuclei with a A 1 = 1.86 mT.The HFCC of the 25 Mg isotope is a B = −11.22mT.

Singlet Yield Calculations
The fractional singlet yield (FSY) of the reaction with and without the magnetic field indicates how the SMF affects the rate of formation of the different reaction products.This could lead to further changes in the reaction rates or product concentrations of various reactions involved in the pathways leading to the experimentally observed effects 26 .
The FSY of the reaction is calculated by tracking the spin state of the RP during the reaction.This calculation can be carried out by solving the Liouville-von Neumann equation, which describes the evolution of the density matrix over time 17 .For a general singlet-born RP under the effect of a weak magnetic field, the fractional yield for periods larger than the lifetime of the RPs is found from the eigenvalues and eigenvectors of the Hamiltonian as: where is the nuclear spin multiplicity, |m⟩ and |n⟩ are the eigenstates of the hamiltonian, Ĥ, with the eigenenergies given by ω m and ω n , respectively, and PS is the singlet projection operator acting on the electron spins.Here, we assume that the singlet and triplet reaction rates are equal and are denoted by k.Finally, r is the relaxation rate or spin-coherence lifetime of the radical pairs 54 .
As mentioned above, the FSY provides information about the yield of the different reaction products with and without the effect of the SMF.We connect this change in the yield of the products to the change in the intracellular concentration of Ca 2+ ions, in order to understand the MFE on the cells.
Due to the lack of experimental data for the exact values of the reaction rate, k, and relaxation rate, r, of the RP we consider, we explore the k/r space for an approximate range of potential values.We consider the ratio of the FSY at the control magnetic field strength(B 0 ) to the FSY at 100 mT(B exp ), which is the experimental field strength used: We plot this quantity against a range of possible k and r values.The magnitude of the ratio reflects the extent of the magnetic field effects on radical pair electrons, which we compare to the magnitude of the ratio of Ca 2+ concentrations within the cells, as shown in Table 1 & Table 2.

Ser Pathway
Fig. 1a depicts a plot of S versus k and r for the radical pair containing the serine oxyradical.In this calculation, we consider the contribution of all Magnesium isotopes according to their natural abundance.As the magnetic field strength for the control setup of the experiment conducted by Hirai et al. 41,42 ranged from B 0 = 0 − 0.3 mT, we take B 0 = 0.15 mT as our control value.Here, the lighter yellow colored areas of the plot within the black outline represents the values of k and r for which there is the greatest variation in the FSY between B 0 and B exp .These points have the potential to capture the behaviour of the system, as they match the magnitude of the ratio of Ca 2+ concentration in the control vs. experimental setups.The values of k and r are consistent with typical values considered in the literature on the RPM.

Tyr Pathway
Similar to what we have done for the case of the Ser Pathway, the quantity S was calculated and plotted against phyically feasible ranges of k and r for the radical pair involving the tyrosine oxyradical.Fig. 2a shows the same.The FSY was found to exhibit a few spikes at certain values of the MF strength, B, for both the Ser and Tyr Pathway.These spikes are associated with certain eigenstates overlapping in energy due to the effect of the spin of the 25 Mg isotope.A dip in the FSY values near zero was also found at higher values of k, as seen from the examples given in Figs. 1 & 2. The extent of the drop increases as k increases for both Ser and Tyr Pathways.The drop in value is seen in the behaviour of the system under the effect of both isotopes.More results regarding the contribution from the individual isotopes are presented in the Appendix of the paper.

Discussion
The primary objective of this study was to determine whether the radical pair mechanism can provide a plausible explanation for the magnetic field effects on NMDAR activity reported by Hirai et al. 41,42 .We find that a simple RP model is capable of capturing the MFEs on the system, as the magnitude of the MFE on the singlet yield corresponds well to the magnitude of the experimental MFE observed on the Ca 2+ ion concentration within the cells.It has been shown above that both Serine and Tyrosine phosphorylation could be involved in the upregulation of NMDAR activity in the presence of a static magnetic field.Serine phosphorylation is involved in activating the AP1 transcription factor, which is a major regulator of the NR2B subunit.As described previously, an increase in NR2B can increase the number of NMDARs in the system, resulting in an increase in the calcium ion influx into the cells.We have also seen how phosphorylation of certain NMDAR Subunits at Serine and Tyrosine sites may also cause an increase in NMDAR activity, which again raises the possibility that magnetic field effects occur through the phosphorylation processes.We thus propose that the radical pairs involved in the ion-radical mechanism of phosphorylation could be a likely candidate for the RPM model.We thus consider radical pairs of the form where the protein residue, R could be a Tyrosine or Serine residue.
Based on the experimental evidence found by Hirai et al., there is an increase in the concentration of Ca 2+ ions within the hippocampal cells in the presence of magnetic fields.For cells exposed to a sustained or repetitive SMF, an increase in concentration of the ions by an average factor of about 1.5 and 1.4, respectively, has been observed.So far, our results indicate that this variation in the behavior of the system in the presence of the magnetic field could indeed be explained using an RPM model.
The radical pairs involving both the Serine and Tyrosine pathway models exhibit significant variation in singlet yield percentage in the presence of the magnetic field.We find that both models capture a variation in FSY by a factor of over 1.3 in the presence of the 100 mT magnetic field for certain plausible values of k and r.Based on these results, we propose that the increase in the yield of the singlet product of the reaction could be leading to the increase in the Ca 2+ ion concentration through various biochemical reactions within the system.Though the increase of the singlet yield of the model is not exactly equal to the experimentally observed increase in the Ca 2+ concentration, it is plausible that an amplification of the effects could occur through the various processes involved in the Ca 2+ ion intake into the cells.Multiple studies have shown that such amplifications are a promising possibility in numerous systems involving the RPM model [55][56][57][58] .Therefore, we propose that both the Ser and Tyr Pathways could potentially explain the MFE observed in experiments.However, there is limited evidence of the formation of serine radicals, particularly in biological systems [59][60][61] .It is worth noting that the radical pairs under consideration could also be involved in various other reactions within the biological system, which could also play a role in the observed magnetic field effects.
Based on the results we have presented, we propose that the radical pair mechanism is a promising model to explain the effects of the magnetic field effects on NMDAR activity.New experiments studying the effect of lower magnetic field strength on NMDA activity could give further insight into the behaviour of the system and the accuracy of the results from our model.Further experiments could also be conducted to study the likelihood of Serine and Tyrosine radical formation in the system and the role they could play in the pathways leading to the MFE on NMDA activity.The most common method for radical detection in biological systems is EPR spectroscopy.Direct EPR spectroscopy could be used to detect the presence of the Tyrosine radical as they are long-lived due to stabilization by electron delocalization.There is some evidence of EPR and time-resolved EPR being successfully used to study Serine radicals 59,61 .These could potentially be used to study if these radicals are present in a given system.Experiments that study the effects of the magnetic field in a system with a higher concentration of the 25 Mg isotope could also be interesting, as we find certain isotope dependent effects in the calculation of FSY, which is discussed in detail in the Appendix section.
It is possible that similar models based on the RPM could be used to understand magnetic field effects on various other biological processes.Furthering our understanding of how external stimuli like magnetic fields affect the different processes within the brain could lead to a much deeper understanding of its working.It could also provide significant insights into finding ways to diagnose and treat various neuropsychiatric and psychological disorders, as well, thus making this simultaneously an intriguing and practically relevant field of study.

DFT Analysis
The DFT calculations for our radical pairs were performed using the ORCA package 62 .The optimization of the molecular structures was performed using the dispersion-corrected PBE0 functional and the def2-TZVP basis set.
The hyperfine coupling constants were calculated from the orbitals obtained from the optimization calculations.The B3LYP functional and def2-TZVP basis set were used for the calculation of both the HFCCs, a A , and a B .The solvent effects in every calculation were considered using the conductor-like polarizable continuum model (CPCM) with a dielectric constant of 2. shows the behaviour when the calculations consider the spinless isotopes, 24 Mg and 26 Mg.Fig. Ib shows the case where the HFI contribution from the spinful isotope, 25 Mg, is considered.We find that higher values of S are found when only the spinful isotope is considered for the calculation.An example of the variation of fractional singlet yield percentage with magnetic field strength from 0.15 mT to 100 mT for a choice of k and r from the feasible region is given in Figs.IIa and IIb.We choose k = 2 × 10 6 s −1 and r = 3 × 10 5 s −1 for the same.We find that the behaviour of the system is significantly different between the spinless and spinful cases.The singlet yield at low field are substantially lower for the spinful case.There are sudden spikes in the value at certain MF strengths, as well, which can be associated with ethe igenenergies of certain states crossing over.The differences that one observes between the spinful and spinless isotopes for this case are similar to those for the Ser pathway.

B. Isotope Contributions at Low Field
We plot examples in Fig. V from both the Ser and Tyr pathway by considering the contribution from the spinless and spinful isotopes separately to understand the contribution of the same to the drop in FSY at MF strength close to 0 mT.We find that the drop in FSY values occur with or without the contribution from the 25 Mg isotope.

C. HFCCs of Other Nuclei
For the calculations in this paper, we consider only the HFI contribution of the nuclei of the oxyradical with the highest HFCC.This is a common approximation that is used in such calculations for the sake of simplicity.
The next highest HFCC of a nucleus in the Serine oxyradical was 5.73 mT.Similarly, in the case of the Tyrosine oxyradical, the next highest HFCCs were -0.64 mT and -0.53 mT. 12/13

Figure 1 .
Figure 1.(a) Ratio of the FSY percentages, S plotted in the k − r plane for the Ser Oxyradical case.The lighter yellow areas are closest in value to the ratios seen in Tables 1 & 2. (b) Variation of the FSY percentage with magnetic field strength ranging from 0.15 mT to 100 mT.The figure is plotted keeping r = 2 × 10 5 s −1 .The reaction rates are k = 2 × 10 6 s −1 and k = 7.5 × 10 6 s −1 , which correspond to S = 1.32 and S = 1.17.All the plots are plotted considering the HFI contribution from the hydrogen atom of the Serine oxyradical with the highest HFCC, a A = 7.45 mT.We also consider the contribution of the 25 Mg isotope at its natural abundance of 10%.The HFI contribution from the same is a B = −11.22mT.

Fig. 2
also shows examples of the variation of FSY percentage with magnetic field strength from 0.15 mT to 100 mT for two different values of k and r.

Figure 2 .
Figure 2. (a) Ratio of the FSY percentages, S plotted in the k − r plane.The lighter yellow areas are closer in value to the ratios seen in Tables 1 & 2. (b) Variation of the FSY percentage with magnetic field strength ranging from 0.15 mT to 100 mT.The figure is plotted keeping r = 2 × 10 5 s −1 .The reaction rates are k = 2 × 10 6 s −1 and k = 7.5 × 10 6 s −1 , which correspond to S = 1.32 and S = 1.17.All the figures are plotted considering the HFI contribution from the hydrogen atom of the Tyrosine oxyradical with the highest HFCC, a A = 1.86 mT.We also consider the contribution of the 25 Mg isotope at its natural abundance of 10%.The HFI contribution from the same is a B = −11.22mT.

Fig.I shows
Fig.I shows the quantity S for the radical pair with the serine oxyradical plotted against the reaction rate, k, and relaxation rate, r.Both calculations consider the contribution from the nucleus of Serine oxyradical that has the highest HFCC, .Fig.Iashows the behaviour when the calculations consider the spinless isotopes,24 Mg and 26 Mg.Fig.Ibshows the case where the HFI contribution from the spinful isotope,25 Mg, is considered.We find that higher values of S are found when only the spinful isotope is considered for the calculation.

Figure I .
Figure I. Ratio of the fractional singlet yield percentages, S plotted against k(s −1 ) and r(s −1 ) for Ser radical case: (a) System containing only 24 Mg and 26 Mg isotopes is considered; (b) System containing 25 Mg isotope with the HFCC, a B = −11.22mT has been considered.

Figure II .
Figure II.Variation of the fractional singlet yield percentage with magnetic field strength ranging from 0.15 mT to 100 mT for the Tyr Pathway.The figure is plotted for the choice of k = 2 × 10 6 s −1 and r = 3 × 10 5 s −1 : (a) System containing only 24 Mg and 26 Mg isotopes is considered; (b) System containing 25 Mg isotope with the HFCC, a B = −11.22mT has been considered.

Figure III .
Figure III.Ratio of the fractional singlet yield percentages, S plotted against k(s −1 ) and r(s −1 ) for Tyr radical case: (a) System containing only 24 Mg and 26 Mg isotopes is considered; (b) System containing 25 Mg isotope with the HFCC, a B = −11.22mT has been considered.

Figure V .
Figure V. Variation of the FSY percentage with magnetic field strength ranging from 0.15 mT to 100 mT plotted with k = 7.5 × 10 6 s −1 and r = 2 × 10 5 s −1 .(a) and (b) are plotted considering a system with only 24 Mg and 26 Mg isotopes; (c) and (d) are plotted considering system with only 25 Mg isotope at its natural abundance of 10%.The HFI contribution from the same is a B = −11.22mT.

Table 1 .
Approximate quantification of the experimental observations reported by Hirai et al.

Table 2 .
contain data from the experiments conducted by Hirai et.al.The data is approximate and has been extracted from plots in the two papers that visualize the variation of Ca 2+ ion concentration in the presence of the SMF, as NMDA concentration is varied.Approximate quantification of the experimental observations reported by Hirai et al.