Equilibrium Dynamics of β-N-Methylamino-L-Alanine (BMAA) and Its Carbamate Adducts at Physiological Conditions

Elevated incidences of Amyotrophic Lateral Sclerosis/Parkinsonism Dementia complex (ALS/PDC) is associated with β-methylamino-L-alanine (BMAA), a non-protein amino acid. In particular, the native Chamorro people living in the island of Guam were exposed to BMAA by consuming a diet based on the cycad seeds. Carbamylated forms of BMAA are glutamate analogues. The mechanism of neurotoxicity of the BMAA is not completely understood, and BMAA acting as a glutamate receptor agonist may lead to excitotoxicity that interferes with glutamate transport systems. Though the interaction of BMAA with bicarbonate is known to produce carbamate adducts, here we demonstrate that BMAA and its primary and secondary adducts coexist in solution and undergoes a chemical exchange among them. Furthermore, we determined the rates of formation/cleavage of the carbamate adducts under equilibrium conditions using two-dimensional proton exchange NMR spectroscopy (EXSY). The coexistence of the multiple forms of BMAA at physiological conditions adds to the complexity of the mechanisms by which BMAA functions as a neurotoxin.


Introduction
The neurotoxicity of β-N-methylamino-L-alanine (BMAA) from the seeds of cycad plants was first demonstrated by Vega and Bell [1].The higher than normal incidence of the Amyotrophic Lateral Sclerosis/Parkinsonism Dementia complex (ALS/PDC), particularly within the local Chamorro people of Guam, is associated with the presence of BMAA [2,3].The role of BMAA in the onset as well as progression of neurodegenerative diseases including Alzheimer's disease (AD) is strongly supported by many studies [4][5][6][7][8][9].Based on the numerous in vitro investigations, the detrimental effects of BMAA are implicated due to the mechanism of excitotoxicity via the activation of excitatory amino acid (EAA) receptors [4,[10][11][12][13][14]. Excitotoxicity could be a significant contributor to neurodegenerative diseases as evidenced by the presence of increased levels of glutamate in the cerebrospinal fluid of ALS patients [15,16].
The interaction of BMAA with bicarbonate ions is critical in the modality of BMAA's role in excitotoxicity.Pioneering work by Weiss and Choi [17] demonstrated the requirement of bicarbonate ions as a cofactor for the activity of BMAA.All the assays performed since the discovery by Weiss and Choi [17], included bicarbonate ions (HCO 3 -) in the media at concentrations ranging from 10 mM to 25 mM [13].The presence of bicarbonate ions at close to physiological pH therefore created ideal conditions for the formation of carbamate adducts [18].Myers and Nelson [18] characterized the interaction between BMAA and the bicarbonate ions and identified the formation of β-carbamate of BMAA.These findings led to a high similarity between the chemical structures of β-carbamate of the BMAA and glutamic acid.This resemblance might hold the potential clue to the mechanism of how BMAA affects glutamate receptors [13,[17][18][19][20].
Molecular interactions with bicarbonate have been studied due to the importance of spontaneous carbamate formation in biological reactions.Carbamate reactions are reversible and occur as the result of a nucleophilic attack by uncharged amines on the carbon dioxide produced from bicarbonate, in vivo [21].The exogenic reaction between the uncharged amine and carbon dioxide to form carbamates has considerable free energy contribution from the dissociation reaction involving carbamic acid (pK a < 4.8) [22].Under physiological conditions (pH ~7), the amines are mostly protonated (pK a ~8) and the CO 2 exists in the hydrated carbonic acid, H 2 CO 3 (pK a1 ~6.3).These conditions tend to be unfavorable for carbamate formation at physiological conditions.Nonetheless the chemical environment can lower the pK a of the protonated amines and thus benefit the carbamate formation which may be further stabilized by other intermolecular non-covalent interactions [23].Carbamate formation plays an important role in many biological functions such as carboxylation of the active sites (lysine side chain εamino group), biosynthesis of biotin and purine biosynthesis [23][24][25], as well as methanogenic archaea mediated reduction of CO 2 to methane [26].The carbamate formation in most cases are non-enzymatic reactions as demonstrated in the case of the carbamate (N-carboxymethanofuran) in methanogenic archaea [24].
Nunn and co-workers first studied the changes in the NMR spectra of BMAA in the presence of bicarbonate ions and demonstrated the formation of the primary carbamate adduct [25,26].They also first hypothesized the formation of a second adduct of BMAA, but the presence was not demonstrated by NMR spectroscopy, as two-dimensional methods then were not routinely employed.In a later study by Myers and Nelson [18], 13 C labeled bicarbonate was used to examine the formation of the BMAA/bicarbonate adducts.Myers and Nelson observed the formation of two carbamate adducts and found these products were formed when the amine, from the BMAA, binds to carbon dioxide produced from bicarbonate.
In this manuscript, we present a study of BMAA: HCO 3 -interaction using high-resolution NMR spectroscopy.In addition to confirming the earlier observations that BMAA: HCO 3 -interactions lead to the formation of both α and β carbamates, we observed that these adducts coexist in the solution state at physiological conditions.Proton chemical shifts of the carbamate adducts are distinctly different from the free BMAA and thus enable the characterization of the exchange kinetics.Two-dimensional saturation transfer exchange spectroscopy (EXSY) [27,28] is used to determine the pseudo first order rate constants of the αand βcarbamate formation/cleavage.The conformational preference of the carbamate adducts and the intermolecular dynamics may provide an important insight to the understanding of the neurotoxicity of BMAA.

Sample preparation
β-N-methylamino-L-alanine (BMAA): bicarbonate samples were prepared with a BMAA concentration of 5 mM with varying bicarbonate concentrations.HCO 3 -concentrations were varied from 0.5 mM to 300 mM by adjusting the concentration of the added sodium bicarbonate in the solution.The BMAA:HCO 3 -samples were prepared in 90% H 2 O and 10% D 2 O or in 100% D 2 O and pH (pD) adjusted to 7.6.A final volume of 600 μL is used for the NMR studies.In all the calculations the pH to pD conversion was not performed because of the constant term that relates pH to pD (Pd = pH+0.41)[29].This is a valid approximation considering the decreased acidities of acids in D 2 O [30,31].

NMR spectroscopy
All NMR experiments were performed in a 400-MHz ( 1 H resonance frequency) VNMRS spectrometer (Varian-Agilent).2D EXSY spectra [27] were collected using standard NOESY pulse sequence in the phase-sensitive mode (States-TPPI).The probe temperature was 30°C.The mixing time for the EXSY experiment (τ m ) was 400 ms to reduce the interference from zeroquantum induced coherence transfer between the J-coupled spins and other intramolecular NOE transfers [32].Two-dimensional TOCSY experiment was performed with a DIPSI multiple pulse mixing sequence (mixing time 80 ms) [33].All the 2D data were collected with 256 complex points along the indirect dimension (t 1 ), 2048 complex points in t 2 over a spectral width of 10 ppm and signal averaged over 32 transients with a recycle delay of 5 s.Timedomain data were zero-filled to 1024 points in ω 1 and 4096 points in ω 2 .A 90°-shifted squared sine bell window function was applied prior to a Fourier transformation in each dimension.All the spectra were referenced to TSP  propionic-2, 2, 3, 3-d 4 acid sodium salt) signal at 0.0 ppm.One-dimensional NMR spectra of freshly prepared samples were performed before and after running the 2D experiments to confirm that the sample conditions were unaltered during the course of the data collection.The spectra were processed and analyzed using Mestrenova (Mestrelab Research, Santiago de Compostela, Spain).

Chemical equilibrium dynamics of carbamate adduct formations
The formation of α and β carbamate adducts of BMAA are shown in Fig 1, with an equilibrium constant of K α and K β , respectively.The homogeneous interactions of CO 2 in aqueous BMAA solutions include a complex series of reactions to form carbamates, hydration reactions involving hydroxide, water and CO 2 , as well as amines [34,35].Thus, the collection of reaction mechanisms involving the various ions and water must be accounted for and quantified to describe the chemical equilibrium of αand βcarbamate formation (see Fig 1).In the presence of H + and CO 2 , the formation of the αand βcarbamates is assumed to follow the equilibrium constants defined as follows: Where [α-carbamate] and [β-carbamate] are the concentrations of the αand βcarbamates in equilibrium with the BMAA ([BMAA]) and carbon dioxide ([CO 2 ]).Based on the NMR data (Fig 2 ) and the previous studies [18,25] there is no evidence of the formation of αβ-carbamate (doubly carbamylated) form of BMAA and is therefore not considered as part of the reaction mechanism.
In order to determine the equilibrium constants (Eq (1)), all the concentration values are required.However, the NMR spectrum cannot differentiate between the protonated vs. deprotonated forms of the carbamate adducts (Fig 1) and therefore all the BMAA concentrations must be related to the starting total BMAA concentration ([BMAA] T ).During the chemical exchange process, the total BMAA concentration is Where [BMAA-α + ] and [BMAA-β + ] are the concentrations of BMAA protonated at the primary and secondary amines, respectively (Fig 1).The equilibrium process of the adduct formation will depend on the total amount of carbon dioxide present in the experimental vessel.The total concentration of CO 2 ([CO 2 ] T ) can also be defined in terms of the various reactions (Fig 1, center) as It is important to note the integral role of bicarbonate ion in Eq (3) (Fig 1) in the formation of CO2 [25,26,36].Using the definitions of total concentrations of BMAA and CO 2 from Eqs (2) and (3), the equilibrium process for the formation of carbamate adducts can be written as Where K α Ã and K β Ã are the equilibrium constants involving the total concentrations of BMAA and CO 2 .Eq (4) is the representation of Eq (1) in terms of the total concentration of BMAA and CO 2 .In Eq (4), k 1i Ã and k 2i Ã with i = α or β are defined as , the equilibrium constants of the adduct formation (Eq (4)) can be written as The constants k 1i ÃÃ and k 2i Ã (i = α or β) (i = α or β) are the pseudo first order exchange constants of carbamate adduct formation and pseudo first order exchange constant of carbamate cleavage, respectively.The total concentration of BMAA (Eq (2)) can be estimated using the pK a of protonated and non-protonated amines.Similarly, the total concentration of the carbon dioxide is estimated using the equilibration process of carbamate and carbamic acid.Two-dimensional EXSY spectrum is used to determine the rate constants k 1α ÃÃ , k 2α Ã (and k 1β ÃÃ , k 2β Ã ) (Eqs ( 4) and ( 6)), while using the equilibrium concentrations of [CO 2 ] T , [CO 2 ] and [H + ].It must be noted that these calculations are valid at any given sample conditions (BMAA:HCO 3 -concentration).In the event of concentration dependent change in the equilibrium, this approach measures the quasi-equilibrium kinetics.

Calculation of exchange rates from EXSY data
The equilibrium process between the BMAA and its αand βcarbamates affects the nuclear spin relaxation process via the mechanism of chemical exchange.For a system of three spins (I = ½), the nuclear spin relaxation process in terms of the rate constants of chemical kinetics can be written in terms of the relaxation matrix as [27,28,37] where τ m is the mixing time used in the EXSY experiment and the elements of Ã (matrix) are defined as [27,28] with 'i' and 'j' are the three chemically exchanging sites from three molecules: free BMAA, αcarbamate and β-carbamate adducts.V ij are the cross peak volumes between chemical shifts of 'i' and 'j' and p j 's are the relative population of the i th species.The relative populations of the three molecules can be defined as Using the definitions from Eqs ( 8)- (11), the matrices Ã and Ř in Eq (7) can be written as The relaxation matrix Ř contains the kinetic parameters of the chemical exchange and longitudinal relaxation rates.In matrix Ã, the quantities V ij (i, j = [BMAA] T , or [α-carbamate] or [β-carbamate]) are the two-dimensional cross-peak volumes in the EXSY experiment and p i ([BMAA] T , or [α-carbamate] or [β-carbamate]) are the relative populations determined from the 1D-NMR spectrum obtained at the same conditions (or from the diagonal peak volumes of the EXSY spectrum (V ij )).The rate constants k 1βα Ã and k 2αβ Ã are the magnetization transfer resulting from the coupled linear differential equation relating non-equilibrium nuclear spin magnetization.These second order (relayed) transfer of magnetization rates are included in the calculations to account that the total magnetization within the exchanging spins is conserved.
The Ř matrix can be obtained by diagonalizing Ã as Where Û is the square matrix of eigenvectors of Ã, such that Û -1 Ã Û = λ = diag (λ j ) and ln(Λ) = diag (ln λ j ), with λ j the j th eigenvalue of Ã [37][38][39].Eq (13) can be solved by several standard routines and herein we have used "R-Statistical Programming" for the process [40].The concentrations of [BMAA] T , [α-carbamate] and [β-carbamate]) in the 600 μL reaction vessels (NMR tubes) were obtained from the integration of the 1D spectra (total BMAA concentration = [BMAA] T , or [α-carbamate] or [β-carbamate]).Concentration of BMAA (unprotonated) were calculated using the known pK a of amines [41].We further assume that the pK a for primary and secondary amines are the same and temperature independent.Total carbon dioxide concentration was estimated using the total NaHCO 3 in the solution (pK a1 of H 2 CO 3 = 6.34) [42].At equilibrium conditions close to the neutral pH, the ratios [H 2 CO 3 ]/[CO 2 ] and [H 2 CO 3 ]/[HCO 3 -] are small (< 0.005 or less) [43].Furthermore, contributions from CO 3 2may also be neglected, as the pK a2 of H 2 CO 3 is 10.3.Upon substituting these values and the corresponding rate constants determined from the EXSY spectrum, the equilibrium constant can be determined in a straightforward manner.
Uncertainty in the measured values were determined by error propagation method (using R).The standard deviation in the spectral data was measured by estimating the noise in the 1D or 2D data by randomly selecting five different regions of the spectra.The standard deviation in each of the measured values was simulated randomly using a Monte Carlo method (5000 samples) with respect to the mean and standard deviation of the of the individual variables in the respective measurements.Typically, the accuracy of the measured peak intensity in the 1D spectra is ~1% and the variation in the diagonal and cross peak volumes of the 2D data are in the range of 3-7%.The details of the calculations, in spreadsheet format along with the R-code, used for matrix diagonalization are provided in S1 Table.

NMR characterization of the BMAA-carbamate adducts
Proton NMR spectra of BMAA in the presence of bicarbonate ions show spectral features from the free BMAA and the αand βcarbamate adducts (Fig 2).The NMR spectra of BMAA and carbamates are similar to the NMR results obtained by Nunn and co-workers [25,26] except for the increased spectral resolution ( 1 H frequency at 400 MHz vs. 300 MHz).We did not find any evidence of the presence of doubly carbamylated form of BMAA in the NMR spectra, which is consistent with the previous observations [18,25].The spectral patterns of all the protons are distinctly different between the three molecules in terms of chemical shift and coupling constants.At 400 MHz, the two 1 Hβ protons and 1 Hα proton show a typical ABX spin-system with the 1 Hα proton chemical shifts of the α-carbamate, β-carbamate and the free BMAA at 4.24 ppm (dd), 3.90 ppm (dd) and Two-dimensional total correlation spectroscopy (TOCSY) experiment was used to identify the sub-set of the intramolecular J-coupled correlations with the BMAA, α-carbamate and βcarbamate molecules (S1 Fig) .Free BMAA in equilibrium with the two BMAA carbamate adducts was established by the two-dimensional EXSY experiment (S2 Fig) .The EXSY experiment shows saturation transfer of magnetization between the protons that undergo chemical exchange between the free-BMAA, and its αand β-carbamate adducts.The free-BMAA peaks are identified in comparison with the NMR spectrum of the BMAA in absence of bicarbonate ions (data not shown).The β-carbamate formation brings the carbonyl oxygen in closer vicinity to the methyl group (Fig 1) leading to a downfield shift.Therefore, the methyl resonance at 2.86 ppm is assigned to the β-carbamate, confirming the earlier observation by Nunn et al. [25,26].The chemical shift assignments for rest of the protons are consistent with the cross peak patterns observed in the TOCSY experiment (S1 Fig).

Equilibrium dynamics of BMAA: adducts
The equilibrium process of the carbamate formation is slow and depends on the relative concentrations of the BMAA to HCO 3 -ions.The chemical exchange between the three species in the solution is reflected in the two-dimensional EXSY experiment between the methyl protons (Fig 4).The relative intensity of the cross peaks suggests that BMAA to α-carbamate exchange vs. the BMAA to β-carbamate exchange mechanisms are expected to have different rates.The cross peaks between the αand βcarbamates at higher BMAA to HCO 3 -are due to relayed transfer effects and are expected to increase in intensity particularly at longer mixing times of the EXSY experiments.The magnetization must be along the z-axis during the mixing period of the EXSY experiment.This allows for the exchange process between the various non-equilibrium spin states.In addition to the chemical exchange process, both cross-relaxation (NOE effects) and coherence transfer via zero quantum coherences, between the J-coupled spins, are active.The NOE contributions might underestimate the exchange contributions particularly for the small molecules and zero-quantum coherence transfer peaks that have anti-phase dispersive components in the spectra [44].One of the effective methods to reduce the effects of both cross relaxation and zero-quantum induced transfers, particularly in the EXSY experiment, is to perform the EXSY experiments with relatively long mixing times.The EXSY experiments here are performed with mixing times in the range of 300-700 ms (data not shown) and found the mixing a time of 400 ms provided optimum sensitivity.For a given concentration of BMAA and bicarbonate ions, the equilibrium constant between the BMAA to the carbamate adducts can be estimated using the EXSY data (Materials and Methods).Using a BMAA:HCO 3 -ratio of 1:10 (BMAA concentration, 5 mM) from the onedimensional 1 H spectra (methyl groups) and the cross peak intensities of the 2D-EXSY spectra the relative populations of total BMAA (p [BMAA]T , Eq (9)), α-carbamate adduct (p [α-BMAA]T , Eq (10)) and β-carbamate adduct (p [β-BMAA]T , Eq (11)) are estimated to be 47%, 39% and 14%, And Ř (Eq ( 14)) become The noise level in the 1D and 2D EXSY spectra are estimated from the regions of the spectra that do not have any peaks to be < 1% and 3-7%, respectively.The Ř (units of s -1 ) shows the exchange dynamics between the three species in solution.At the ratio of BMAA: HCO 3 -(1:10), the rate of formation of the α-carbamate and β-carbamate adducts are 0.17 ± 0.03 s -1 (k 1α ÃÃ ) and 0.60 ± 0.05 s -1 (k 1β ÃÃ ), while the cleavages are 0.23 ± 0.02 s -1 (k 2α Ã ) and 1.89 ± 0.03 s -1 (k 2β Ã ), respectively (from matrix elements in Eq ( 16)).The half-lives of the formation the α- carbamate and β-carbamate adducts (assuming a first order reaction) are 4.18 ± 0.99 s (t ½,1α ) and 1.17 ± 0.10 s (t ½,1β ), while the cleavages are 3.08 ± 0.34 s (t ½,2α ) and 0.37 ± 0.01 s (t ½,2β ), respectively.The half-life of the α-carbamate formation decreases with increasing bicarbonate formation with increasing carbonate ion concentration while the cleavage increases (Table 1) suggesting that the system is still in quasi-equilibrium condition.Indirect (relayed) magnetization transfer rates between the adducts (k 1βα Ã and k 2αβ Ã ) are much smaller than the direct transfer terms.Upon converting the exchange rates as described previously (Materials and Methods, S1 Table ), the equilibrium constant (K i Ã , Eq (4)) for the formation of α-carbamate and β-carbamate adducts are estimated as 7.70 ± 1.77 and M -1 and 3.11 ± 0.26 M -1 , respectively.A summary of the experimentally measured equilibrium constants of the carbamate adduct formation is given (Table 1) along with the other relevant first and second order rate constants.A spreadsheet describing the calculations along with an R-code for matrix diagonalization is provided in S1 Table .To account for the protonation (K i , Eq (1)), the equilibrium constants of the formation of α-carbamate and β-carbamate are 7.60 ± 1.75 ×10 6 and 3.11 ± 0.26×10 6 , respectively.Using the equilibrium conditions, the reaction Gibbs free energy (ΔG α Ã ) for the αcarbamate formation was estimated to be -5.14 ± 1.18 kJ mol -1 , and -2.89 ± 0.24 kJ mol -1 for the β-carbamate (Table 1).By increasing the concentration of the bicarbonate ions to 200 mM (BMAA:HCO 3 -ratio of 1:20) the concentration of CO 2 increases, which leads to an increase of K 1α Ã to 13.52 ± 1.71 M -1 and K 1β Ã to 3.23 ± 0.09 M -1 and the Gibbs free energy decreases (ΔG α Ã = -6.56 ± 0.83 kJ mol -1 and ΔG β Ã = -2.96± 0.08 kJ mol -1 ).The trend continues further upon increasing the concentration of NaHCO 3 to 300 mM (BMAA: HCO3 -ratio of 1:30) with K 1α Ã = 20.51 ± 0.58 M -1 , ΔG α Ã = -7.61±0.22 kJ.mol -1 , K 1β Ã = 4.41 ± 0.20 M -1 and ΔG β Ã = -3.74± 0.17 kJ.mol -1 (Table 1).A linear trend between the equilibrium constants and the total concentration of NaHCO 3 is observed: K 1α Ã increases approximately 0.07 M -1 for an increase of 1 mM of HCO 3 -ion, while K 1β Ã increases approximately by 0.006 M -1 per 1 mM of HCO 3 -.

Discussion
Carbamate synthesis in biological systems has an integral dependence on the presence of nonprotonated form of the amines, which are the preferred protonation states under alkaline conditions.Aerial carbonation of amides to form carbamates is part of the biological fixation of CO 2 by many organisms.Carbamate formation has been implicated in many biological functions ranging from inflammatory responses to tumor progression [45][46][47][48][49][50].In particular, the formation of α-carbamate by cysteine and its neurotoxicity is bicarbonate-dependent [51].Erving et al. have studied the rate of single carbamate formation in aqueous solution using spectrophotometric methods, including NMR spectroscopy [21].EXSY based determination of both carbamate formation (k 1α and k 1β ) and cleavage (k 2α and k 2β ) are within the range of experimental measurements of aliphatic amines [24].Similar results have been obtained in an EXSY based approach for single carbamate formation in N-carboxymethanofuran (carbamate) from methanofuran [24].The measured experimental equilibrium constant for the β-carbamate  3 and Table 1).The equilibrium constants of the β-carbamate are within the range of previously published values [24,52] while that of the α-carbamate tend to be higher.However, the earlier studies have only focused on the equilibrium process of single carbamate formation while the current approach involves formation/cleavage of two carbamates simultaneously.In this work, in addition to confirming the formation of the carbamate adducts, we demonstrate that BMAA and the αand β-carbamate adducts simultaneously coexist in the solution with an equilibrium dynamic that is dependent on the concentration of CO 2 .
Nunn and coworkers have studied the formation BMAA carbamates at neutral conditions [25,36,53,54].Some of the chemical features that contributes to the carbamate formation in BMAA include, differential degree of ionization of the protons attached to the nitrogen atoms, low pKa of the amino groups [26], and the overall neutral charge at physiological conditions.These features lead to the reaction between the non-protonated α-amino group of BMAA and the bicarbonate which encourages the formation of the carbamate adducts at physiological conditions [36].Therefore the formation of both α-(2-amino) and β-(3-methylamino) carbamates is a probable explanation of why the BMAA activity is measured only in the presence of bicarbonate/CO 2 [25,26,36,53].At equilibrium the relative population of the carbamate adducts depends on the concentration of the bicarbonate ions with the α-carbamate population higher than the β-carbamate (Fig 3).Our data does not show the presence of doubly carbamylated forms of BMAA, which follows a similar observation by Davis et al. [26].Simultaneous presence of the three molecules in solution further suggests that if either one or both the carbamates are removed from the reaction equilibrium (due to binding to EAA receptors), the equilibrium process will continue to replenish the carbamate populations from free BMAA.
Physiological concentration of bicarbonate ions in bodily fluids range from 20-25 mM [55].A BMAA concentration of ~5 mM is required for good signal-to-noise ratio particularly for the 2D EXSY experiments.Thus at a BMAA:HCO 3 -ratio of 1:20, the estimated bicarbonate concentration of 100 mM may be physiologically higher.Biological experiments performed in mixed spinal cord cultures to measure motor neuron (MN) loss suggest the concentration dependence of the bicarbonate ions in the media [12,56].Therefore, the ratio of BMAA to HCO 3 -might be more important than the absolute concentrations of BMAA or bicarbonate ions themselves.
The concentration ratio presented in this work is similar to the experimental conditions of Richter and Mena [13], where the BMAA to bicarbonate ratio was 1 mM to 20 mM.Recently, to investigate the etiopathogenesis of Parkinsonism-dementia (PD), Arif et al. [57] tested the effect of BMAA on protein phosphatase 2A (PP2A) activity and tau hyperphosphorylation in mouse primary neuronal cultures and metabolically active rat brain slices.These experiments are For each i = α (α-carbamate) and i = β (β-carbamate) the rate constants k ÃÃ 1i and k Ã 2i were determined by proton exchange NMR spectroscopy (EXSY).Reaction half-lives were determined by assuming a first order kinetics (t1 = 2 ;i ¼ lnð2Þ = k i ).The equilibrium constants K Ã i were then calculated using k Ã 1i ; k 1i and k 2i along with the estimates obtained by propagating the measured standard error in the NMR spectra (materials and methods).The Gibbs free energy was then estimated using DG Ã i ¼ ÀRTlnðKÞ.The experiments were performed in a total volume of 600 μL containing 5 mM BMAA and varying concentration of NaHCO 3 (100 mM, 200 mM and 300 mM).The sample pD was at 7.6 and the experiments were performed at 303 K. doi:10.1371/journal.pone.0160491.t001performed typically with 1 mM BMAA and to 35mM of NaHCO 3 , leading to a BMAA: HCO 3 -ratio of 1:10 to 1:35 (at pH 7.4).The BMAA: HCO 3 -ratio in these experiments by Arif et al. [57] are in the same range employed in the results presented in this work.
It must be noted that the measured equilibrium constant (Table 1) depends on the concentration HCO 3 -and therefore the system is still in a quasi-equilibrium state.The quasi-equilibrium constant of the α-carbamate formation shows greater than the β-carbamate with increasing the concentration HCO 3 -.Physiological concentrations of BMAA, as low as 10-30 μM, have shown to cause neuronal injury [12,56].Thus with a physiological concentration of bicarbonate in the range of 20-25 mM, the HCO 3 -will be > 1,000 times higher than BMAA (BMAA: HCO 3 -= ~1:1,000-1:2,000).Based on the variation in the relative concentrations of the three species in solution (Fig 3), one may assume the concentration of BMAA and the βcarbamate will be the same (~25% each) while the α-carbamate will be twice that of the BMAA (or β-carbamate) with ~50% of the concentration.These estimates are based on the extrapolation of the experimental data (Fig 3 ) under in vitro conditions.The hypothesized modes of action for BMAA include GluR activation [19] and binding metalloenzymes causing enzyme malfunction [58].In addition to the different modes of action discussed, the required concentration for neurotoxicity is also varied among studies, yielding conflicting results.BMAA can also cause neurodegeneration through misincorporation into proteins.BMAA has been shown to replace serine within proteins [59], and when BMAA is incorporated the proteins can no longer perform their function due to misfolding.It is unknown if BMAA or the carbamate adducts are misincorporated into proteins.There is a possibility BMAA and the carbamate adducts have different probabilities of misincorporation, and the carbamate adducts might misincorporate to replace different amino acids other than serine.If the carbamate adducts can be misincorporated into proteins in addition to free BMAA, there could be a higher percentage of misfolded or inactive proteins.These misfolded proteins might induce the formation of plaques causing cell damage.
BMAA and HCO 3 -interaction produces two carbamate adducts along with free BMAA in equilibrium at physiological conditions.The concentration of the carbamate adducts are determined with varying BMAA:HCO 3 -ratios, and it was determined that the equilibrium of the three molecules depend on the concentration of CO 2 .With an increased amount of HCO 3 more of the carbamate adducts are formed.The equilibrium constants governing the formation of the α-carbamate (K α Ã ) and the β-carbamate (K β Ã ) and the respective reaction Gibbs free energies (ΔG α Ã and ΔG β Ã ) are determined under equilibrium conditions via 2D proton exchange NMR spectroscopy (EXSY).The α-carbamate tends to be a more preferred adduct of BMAA because the primary amine is more likely to be deprotonated over the secondary amine.Equilibrium carbamate formation is important when determining the neurotoxic concentration of BMAA (directly depends on the physiological conditions of [CO 2 ] and pH).This implies a sample with BMAA/HCO 3 -can form a neurotoxic molecule and the neurotoxic concentration is not the initial concentration of BMAA, but a fraction of the initial concentration.The dynamic equilibrium process of BMAA with its carbamate adducts could play a key role in explaining why BMAA can only activate the GluR when in the presence of HCO 3 -and provide insight on BMAA's mode of action for neurotoxicity.S1 Table .Example calculations.Spreadsheet used to determine to equilibrium constants from the EXSY spectrum and along with a code (written in R) used eigenvalue calculation (matrix diagonalization).(PDF)

Fig 1 .
Fig 1. Kinetic equilibrium process of the conversion and chemical exchange of β-Methylamino-L-alanine (BMAA) with α-carbamate and βcarbamate upon interaction with sodium bicarbonate.At neutral pH both the αand βcarbamate adducts coexist in solution and undergo detectable chemical exchange with the free BMAA.The exchange rates between the free BMAA and that of the α-carbamate and β-carbamates are defined as k 1α (reverse rate k 2α ) and k 1β (reverse rate k 2β ), with the chemical equilibrium constants of K α and K β , respectively.Two-dimensional saturation transfer exchange spectroscopy (EXSY) is used to characterize the chemical exchange mechanism and equilibrium constants.Figure also shows the other reactions that are responsible for generating the bicarbonate ion and CO 2 under aqueous conditions (middle part) and are shown in block letters.The two-way arrows with continuous lines (i.e., ⇄) to represent equilibria with kinetically measurable rate constants, while arrows with the dashed lines represent equilibrium processes that are faster than the NMR measured chemical exchange parameters.doi:10.1371/journal.pone.0160491.g001
Fig 3 shows the variation in the amount of carbamates (α and β) formed as a function of increasing bicarbonate ions for a fixed concentration of BMAA (5 mM).The continuous curves are polynomial fit to the experimental points and may only be used to suggest an expected trend.As seen in Fig 3, even at a BMAA: HCO 3 -of 1:30, there is a significant equilibrium dynamic process between the free BMAA and the αand βcarbamates.All the resonances ( 1 Hα, 1 Hβ and 1 Hγ) from the BMAA, the α-carbamate, and the β-carbamate show signs of chemical exchange between the three forms (S2 Fig).The methyl proton ( 1 Hγ) resonances of the three molecules are void of overlap from other coherence transfer peaks (Fig 2 Panel B, S2 Fig) and therefore this region can be used to determine the kinetic parameters.

Fig 3 .
Fig 3. Formation of α-carbamate and β-carbamate in the equilibrium process and its dependence of concentration of bicarbonate ions in the solution.Estimation of the total concentration of different species in solution [BMAA] (black circles), [α-carbamate] (red squares) and [β-carbamate] (blue diamonds) as a function of increasing HCO 3 -([NaHCO 3 ]) concentration.The concentrations of each species were determined using one-dimensional 1 H NMR spectroscopy.The concentration of the initial BMAA was 5 mM and all the experiments were performed at a pD of 7.6 and at 30°C.doi:10.1371/journal.pone.0160491.g003

Fig 4 .
Fig 4. Effect of increasing bicarbonate concentration: Methyl region of the 400 MHz two-dimensional exchange spectroscopy (EXSY) spectra with increasing BMAA: HCO3-ratio: (A) 1:10, (B) 1:20 and (C) 1:30.The BMAA concentration was at 10 mM.EXSY spectrum is recorded at 30°C and with a mixing time of 400 ms.doi:10.1371/journal.pone.0160491.g004 26 M -1 to 4.36 ± 0.09 M -1 while the corresponding values for the α-carbamate (K α Ã = k 1α Ã / k 2α ) are from 7.60 ± 1.75 M -1 to 20.25 ± 0.58 M -1 , at varying ratios of BMAA:HCO 3 -(Figs 2, Fig. TOCSY spectrum.400 MHz two-dimensional total correlation (TOCSY) spectra shows the co-existence of BMAA with the primary and secondary carbamate adducts.The reaction mixture consists of 10 mM BMAA and 200 mM bicarbonate in D 2 O. TOCSY experiments were performed at 30°C with a mixing time of 80 ms.The total spectrum is shown in (A), while the sub-spectra in (B), (C) and (D) shows the J-coupling related to primary carbamate adduct, BMAA (free) and secondary carbamate adduct, respectively.(PDF) S2 Fig. EXSY spectrum.400 MHz two-dimensional exchange spectroscopy (EXSY) spectra shows the chemical exchange between the BMAA with the primary and secondary carbamate adducts.The full spectrum is shown in (a) while the expanded region of the αprotons in (B) and the methyl protons in (C).The sample contains 10 mM BMAA and 200 mM bicarbonate in D 2 O. EXSY spectrum is recorded at 30°C and with a mixing time of 400 ms.The BMAA, primary carbamate adduct (α) and the secondary carbamate adduct (β) are marked panels (B) and (C).(PDF)