From 58a435a0c95517ffd10326b38dc121954cc965f0 Mon Sep 17 00:00:00 2001 From: nkoch1 Date: Sun, 22 May 2022 22:50:41 -0400 Subject: [PATCH] Rheobase results section rework --- manuscript.tex | 28 ++++++++++++++++++---------- 1 file changed, 18 insertions(+), 10 deletions(-) diff --git a/manuscript.tex b/manuscript.tex index 020cceb..18f885b 100644 --- a/manuscript.tex +++ b/manuscript.tex @@ -285,7 +285,7 @@ To examine the role of cell-type specific ionic current environments on the impa \end{figure} \subsection*{Variety of model neurons} -Neuronal firing is heterogenous across the CNS and a set of neuronal models with heterogenous firing due to different ionic currents is desirable to reflect this heterogeneity. The set of single-compartment, conductance-based neuronal models used here has considerable diversity as evident in the variability seen across neuronal models both in spike trains and their fI curves (\Cref{fig:diversity_in_firing}). The models chosen all fire tonically and do not exhibit bursting. See methods for details and naming of the models. Some models, such as Cb stellate and RS inhibitory models, display type I firing whereas others such as Cb stellate \(\Delta\)\Kv and STN models have type II firing. Type I firing is characterized by continuous fI curves (i.e. firing rate increases from 0 in a continuous fashion) whereas type II firing is characterized by a discontinuity in the fI curve (i.e. a jump occurs from no firing to firing at a certain frequency). The other models used here lie on a continuum between these prototypical firing classifications. Most neuronal models exhibit hysteresis with ascending and descending ramps eliciting spikes at different current thresholds, however the STN +\Kv, STN \(\Delta\)\Kv, and Cb stellate \(\Delta\)\Kv models have large hysteresis (\Cref{fig:diversity_in_firing}, \Cref{fig:ramp_firing}). \notejb{No bifurctions, but mention different dynamics!} \notenk{Like this?} Different underlying dynamics generate these different firing types and hysteresis \cite{ERMENTROUT2002, ermentrout_type_1996, Izhikevich2006}. +Neuronal firing is heterogenous across the CNS and a set of neuronal models with heterogenous firing due to different ionic currents is desirable to reflect this heterogeneity. The set of single-compartment, conductance-based neuronal models used here has considerable diversity as evident in the variability seen across neuronal models both in spike trains and their fI curves (\Cref{fig:diversity_in_firing}). The models chosen all fire tonically and do not exhibit bursting. See methods for details and naming of the models. Some models, such as Cb stellate and RS inhibitory models, display type I firing whereas others such as Cb stellate \(\Delta\)\Kv and STN models have type II firing. Type I firing is characterized by continuous fI curves (i.e. firing rate increases from 0 in a continuous fashion) whereas type II firing is characterized by a discontinuity in the fI curve (i.e. a jump occurs from no firing to firing at a certain frequency). The other models used here lie on a continuum between these prototypical firing classifications. Most neuronal models exhibit hysteresis with ascending and descending ramps eliciting spikes at different current thresholds, however the STN +\Kv, STN \(\Delta\)\Kv, and Cb stellate \(\Delta\)\Kv models have large hysteresis (\Cref{fig:diversity_in_firing}, \Cref{fig:ramp_firing}). \notejb{No bifurctions, but mention different dynamics!} \notenk{Bifurcations removed and new sentence like this?} Different underlying dynamics generate these different firing types and hysteresis \cite{ERMENTROUT2002, ermentrout_type_1996, Izhikevich2006}. \subsection*{Characterization of Neuronal Firing Properties} \begin{figure}[tp] @@ -297,23 +297,23 @@ Neuronal firing is heterogenous across the CNS and a set of neuronal models with \linespread{1.}\selectfont \caption[]{Characterization of firing with AUC and rheobase. (A) The area under the curve (AUC) of the repetitive firing frequency-current (fI) curve. (B) Changes in firing as characterized by \(\Delta\)AUC and \(\Delta\)rheobase occupy 4 quadrants separated by no changes in AUC and rheobase. Representative schematic fI curves in blue with respect to a reference fI curve (black) depict the general changes associated with each quadrant.} - \label{fig:firing_characterizaton} + \label{fig:firing_characterization} \end{figure} -Neuronal firing is a complex phenomenon and a quantification of firing properties is required for comparisons across cell types and between conditions. Here we focus on two aspects of firing: rheobase, the smallest injected current at which the cell fires an action potential, and the shape of the frequency-current (fI) curve as quantified by the area under the curve (AUC) for a fixed range of input currents above rheobase (\Cref{fig:firing_characterizaton}A). \notejb{this way we make the AUC independent from changes in rheobase} \notenk{How about ``This enables a AUC measurement independent from rheobase.''}The characterization of firing by rheobase and AUC allows to characterize both a neuron's excitability in the sub-threshold regime (rheobase) and periodic firing in the super-threshold regime (AUC). Note that AUC is essentially quantifying the slope of a neuron's fI curve. +Neuronal firing is a complex phenomenon and a quantification of firing properties is required for comparisons across cell types and between conditions. Here we focus on two aspects of firing: rheobase, the smallest injected current at which the cell fires an action potential, and the shape of the frequency-current (fI) curve as quantified by the area under the curve (AUC) for a fixed range of input currents above rheobase (\Cref{fig:firing_characterization}A). \notejb{this way we make the AUC independent from changes in rheobase} \notenk{How about ``This enables a AUC measurement independent from rheobase.''}The characterization of firing by rheobase and AUC allows to characterize both a neuron's excitability in the sub-threshold regime (rheobase) and periodic firing in the super-threshold regime (AUC). Note that AUC is essentially quantifying the slope of a neuron's fI curve. -Using these two measures we quantify the effects a changed property of an ionic current has on neural firing by the differences in both rheobase, \drheo, and in AUC, \(\Delta\)AUC, relative to the wild type neuron. \(\Delta\)AUC is in addition normalized to the AUC if the wild type neuron, see Eq.~\eqref{eqn:AUC_contrast}. Each fI curve resulting from an altered ionic current is a point in a two-dimensional coordinate system spanned by \drheo and \notejb{normalized \(\Delta\)AUC} (\Cref{fig:firing_characterizaton}B). An fI curve similar to the one of the wild type neuron is marked by a point close to the origin. In the upper left quadrant, fI curves become steeper (positive difference of AUC values: \(+\Delta\)AUC) and are shifted to lower rheobases (negative difference of rheobases: \(-\)\drheo), unambigously indicating an increased firing that clearly might be classified as a GOF of neuronal firing. The opposite happens in the bottom right quadrant where the slope of fI curves decreases (\(-\Delta\)AUC) and the rheobase is shifted to higher currents (\(+\)\drheo), indicating a decreased, LOF firing. In the lower left (\(-\Delta\)AUC and \(-\)\drheo) and upper right (\(+\Delta\)AUC and \(+\)\drheo) quadrants, the effects on firing are less clear-cut, because the changes in rheobase and AUC have opposite effects on neuronal firing. Changes in a neuron's fI curves in these two quadrants cannot uniquely be described as a gain or loss of excitability. In these cases it depends on the regime the neuron is operating in. If it is in its excitable regime and only occasionaly generates an action potential, then the effect on the rheobase matters more. If it is firing periodically with high rates, then the change in AUC might be more relevant. +Using these two measures we quantify the effects a changed property of an ionic current has on neural firing by the differences in both rheobase, \drheo, and in AUC, \(\Delta\)AUC, relative to the wild type neuron. \(\Delta\)AUC is in addition normalized to the AUC if the wild type neuron, see Eq.~\eqref{eqn:AUC_contrast}. Each fI curve resulting from an altered ionic current is a point in a two-dimensional coordinate system spanned by \drheo and \ndAUC (\Cref{fig:firing_characterizaton}B). An fI curve similar to the one of the wild type neuron is marked by a point close to the origin. In the upper left quadrant, fI curves become steeper (positive difference of AUC values: \(+\Delta\)AUC) and are shifted to lower rheobases (negative difference of rheobases: \(-\)\drheo), unambigously indicating an increased firing that clearly might be classified as a GOF of neuronal firing. The opposite happens in the bottom right quadrant where the slope of fI curves decreases (\(-\Delta\)AUC) and the rheobase is shifted to higher currents (\(+\)\drheo), indicating a decreased, LOF firing. In the lower left (\(-\Delta\)AUC and \(-\)\drheo) and upper right (\(+\Delta\)AUC and \(+\)\drheo) quadrants, the effects on firing are less clear-cut, because the changes in rheobase and AUC have opposite effects on neuronal firing. Changes in a neuron's fI curves in these two quadrants cannot uniquely be described as a gain or loss of excitability. In these cases it depends on the regime the neuron is operating in. If it is in its excitable regime and only occasionaly generates an action potential, then the effect on the rheobase matters more. If it is firing periodically with high rates, then the change in AUC might be more relevant. \subsection*{Sensitivity Analysis} Sensitivity analyses are used to understand how input model parameters contribute to determining the output of a model \citep{Saltelli2002}. In other words, sensitivity analyses are used to understand how sensitive the output of a model is to a change in input or model parameters. One-factor-a-time sensitivity analyses involve altering one parameter at a time and assessing the impact of this parameter on the output. This approach enables the comparison of given alterations in parameters of ionic currents across models. -For example, when shifting the half-activation voltage \(V_{1/2}\) of the potassium current in the FS +\Kv model to more depolarized values, then the rheobase of the resulting fI curves shifts to lower currents \(-\)\drheo), making the neuron more sensitive to weak inputs, but at the same time the slope of the fI curves is reduced (\(-\)\ndAUC), reducing firing rate (\Cref{fig:AUC_correlation}~A). \notejb{Refer to Fig 2B} Plotting the corresponding changes in AUC against the change in \(V_{1/2}\) results in a monotonously falling curve (thick orange line in \Cref{fig:AUC_correlation}~B). For each of the many models we get a different relation between the changes in AUC and the shifts in \(V_{1/2}\) (thin lines in \Cref{fig:AUC_correlation}~B). To further summarize these different dependencies of the various models we characterize each of these curve by a single number, the \( \text{Kendall} \ \tau \) correlation coefficient. A monotonously increasing curve results in a \( \text{Kendall} \ \tau \) close to \(+1\), a monotounsly decreasing curve in \( \text{Kendall} \ \tau \approx -1 \), and a non-monotonous, non-linear relation relation in \( \text{Kendall} \ \tau \) close to zero (compare lines in \Cref{fig:AUC_correlation}~B with dots in black box in panel C). +For example, when shifting the half activation voltage \(V_{1/2}\) of the delayed rectifier potassium current in the FS +\Kv model to more depolarized values, then the rheobase of the resulting fI curves shifts to lower currents \(-\)\drheo, making the neuron more sensitive to weak inputs, but at the same time the slope of the fI curves is reduced (\(-\)\ndAUC), reducing firing rate (\Cref{fig:AUC_correlation}~A). As a result the effect of a depolarizing shift in the delayed rectifier potassium current half activation \(V_{1/2}\) in FS neurons is in the bottom left quadrant of \Cref{fig:firing_characterization}~B and characterization as LOF or GOF in excitability is not possible. Plotting the corresponding changes in AUC against the change in half activation potential \(V_{1/2}\) results in a monotonically falling curve (thick orange line in \Cref{fig:AUC_correlation}~B). For each of the many models we get a different relation between the changes in AUC and the shifts in half maximal potential \(V_{1/2}\) (thin lines in \Cref{fig:AUC_correlation}~B). To further summarize these different dependencies of the various models we characterize each of these curve by a single number, the \( \text{Kendall} \ \tau \) correlation coefficient. A monotonically increasing curve results in a \( \text{Kendall} \ \tau \) close to \(+1\), a monotounsly decreasing curve in \( \text{Kendall} \ \tau \approx -1 \), and a non-monotonous, non-linear relation relation in \( \text{Kendall} \ \tau \) close to zero (compare lines in \Cref{fig:AUC_correlation}~B with dots in black box in panel C). -Changes in gating \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affect AUC (\Cref{fig:AUC_correlation}), but how exactly AUC is affected usually depends on the specific models. Increasing, for example, the slope factor of the \Kv activation curve, increases the AUC in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (\Cref{fig:AUC_correlation}~B). Similar consistent positive correlations can be found for shifts in A-current activation \(V_{1/2}\). Changes in \Kv activation \(V_{1/2}\) and in maximal A-current conductance result in negative correlations with the AUC in all models (\( \text{Kendall} \ \tau \approx -1\)). +Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affect AUC (\Cref{fig:AUC_correlation}), but how exactly AUC is affected usually depends on the specific models. Increasing, for example, the slope factor of the \Kv activation curve, increases the AUC in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (\Cref{fig:AUC_correlation}~D,E,F). Similar consistent positive correlations can be found for shifts in A-current activation \(V_{1/2}\). Changes in \Kv half activation \(V_{1/2}\) and in maximal A-current conductance result in negative correlations with the AUC in all models (\( \text{Kendall} \ \tau \approx -1\)). -Qualitative differences can be found, for example, when increasing the maximal conductance of the delayed rectifier (\Cref{fig:AUC_correlation}~C). In some model neurons this increased AUC (\( \text{Kendall} \ \tau \approx +1\)), whereas in others AUC is decreased (\( \text{Kendall} \ \tau \approx -1\)). In the STN +\Kv model, AUC depends in a non-linear way on the maximal conductance of the delayed rectifier, resulting in an \( \text{Kendall} \ \tau \) close to zero. Even more dramatic qualitative differences between models result from shifts of the activation curve of the delayed rectifier, as discussed already above (\Cref{fig:AUC_correlation}~A). Some model neurons do almost not depend on changes in K-current activation \(V_{1/2}\) or show strongly non-linear dependencies, both resulting in \( \text{Kendall} \ \tau\) close to zero. Many model neurons show strongly negative correlations, and a few show positive correlations with shifting the activation curve of the delayed rectifier. +Qualitative differences can be found, for example, when increasing the maximal conductance of the delayed rectifier (\Cref{fig:AUC_correlation}~G,H,I). In some model neurons this increased AUC (\( \text{Kendall} \ \tau \approx +1\)), whereas in others AUC is decreased (\( \text{Kendall} \ \tau \approx -1\)). In the STN +\Kv model, AUC depends in a non-linear way on the maximal conductance of the delayed rectifier, resulting in an \( \text{Kendall} \ \tau \) close to zero. Even more dramatic qualitative differences between models result from shifts of the activation curve of the delayed rectifier, as discussed already above (\Cref{fig:AUC_correlation}~A,B,C). Some model neurons do almost not depend on changes in K-current half activation \(V_{1/2}\) or show strongly non-linear dependencies, both resulting in \( \text{Kendall} \ \tau\) close to zero. Many model neurons show strongly negative correlations, and a few show positive correlations with shifting the activation curve of the delayed rectifier. -\notejb{the slope factor has a name (``slope factor''), but \(V_{1/2}\) not. How is this called, ``midpoint potential''? ``half activation potential/voltage''}\notenk{How about ``half-maximal potential''?} +\notejb{the slope factor has a name (``slope factor''), but \(V_{1/2}\) not. How is this called, ``midpoint potential''? ``half activation potential/voltage''}\notenk{Changed this to half activation or half inactivation when only talking about activation or inactivation and to half maximal potential when referring to both activation and inactivation.} \notejb{``maximum conductance''} \notenk{fixed this in results section} @@ -333,8 +333,16 @@ Qualitative differences can be found, for example, when increasing the maximal c \label{fig:AUC_correlation} \end{figure} -\notejb{Add more detailed step-by-step description what we did here, like for AUC.} -The rheobase is also affected by changes in channel kinetics (\Cref{fig:rheobase_correlation}). In contrast to AUC, most alterations result in similar changes of rheobase in all models, but there are some noteable exceptions. Changing the slope factor of Na-current inactivation, \Kv-current inactivation, and A-current activation affect rheobase both with positive and negative correlations in different models (\Cref{fig:rheobase_correlation}B). Departures from monotonic relationships occur in some models as a result of K-current activation, \Kv-current inactivation, and A-current activation in some models \notejb{which parameters are changed?}. Maximum conductance affects rheobase similarly across models (\Cref{fig:rheobase_correlation}C). However, identical changes in current gating properties such as activation or inactivation \(V_{1/2}\) or slope factor \(k\) can have differing effects on firing depending on the model in which they occur. \notejb{This we just have said...} + +\notejb{Add more detailed step-by-step description what we did here, like for AUC.}\notenk{started doing this} + +Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affect rheobase (\Cref{fig:rheobase_correlation}), however in contrast to AUC qualitatively consistent effects on rheobase across models are observed. Increasing, for example, the maximal conductance of the leak current in the Cb stellate model increases the rheobase (\Cref{fig:rheobase_correlation}~G). When these changes are plotted against the change in maximal conductance a monontonically increasing relationship is evident (thick teal line in \Cref{fig:AUC_correlation}~H). This monotonically increasing relationship is evident in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (thin lines in \Cref{fig:rheobase_correlation}~H). Similarily, positive correlations are consistently found across models for maximal conductances of delayed rectifier K, \Kv, and A type currents, whereas the maximal conductance of the sodium current consistently is associated with negative correlations (\( \text{Kendall} \ \tau \approx +1\); \Cref{fig:rheobase_correlation}~I). + +Although changes in half maximal potential \(V_{1/2}\) and slope factor \(k\) generally correlate with rheobase similarly across model there are some exceptions. Changing the slope factor of Na-current inactivation, \Kv-current inactivation, and A-current activation affect rheobase both with positive and negative correlations in different models (\Cref{fig:rheobase_correlation}~F). Departures from monotonic relationships also occur in some models as a result of K-current activation \(V_{1/2}\) and slope factor \(k\), \Kv-current inactivation slope factor \(k\), and A-current activation slope factor \(k\) in some models. Thus, identical changes in current gating properties such as the half maximal potential \(V_{1/2}\) or slope factor \(k\) can have differing effects on firing depending on the model in which they occur. \notejb{This we just have said...} + +%The rheobase is also affected by changes in channel kinetics (\Cref{fig:rheobase_correlation}). In contrast to AUC, most alterations result in similar changes of rheobase in all models, but there are some noteable exceptions. Changing the slope factor of Na-current inactivation, \Kv-current inactivation, and A-current activation affect rheobase both with positive and negative correlations in different models (\Cref{fig:rheobase_correlation}~B). Departures from monotonic relationships occur in some models as a result of K-current activation, \Kv-current inactivation, and A-current activation in some models \notejb{which parameters are changed?}. Maximum conductance affects rheobase similarly across models (\Cref{fig:rheobase_correlation}~C). However, identical changes in current gating properties such as the half maximal potential \(V_{1/2}\) or slope factor \(k\) can have differing effects on firing depending on the model in which they occur. \notejb{This we just have said...} + + \begin{figure}[tp] \centering