add ./For_Submission/Koch_reviewer_comments to address reviewer comments
This commit is contained in:
nkoch1
parent
78bd33da4e
commit
a49a1b8034
BIN
Frontiers_manuscript/For_Submission/AUC_correlation.jpg
Normal file
BIN
Frontiers_manuscript/For_Submission/AUC_correlation.jpg
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 441 KiB |
BIN
Frontiers_manuscript/For_Submission/AUC_correlation.pdf
Normal file
BIN
Frontiers_manuscript/For_Submission/AUC_correlation.pdf
Normal file
Binary file not shown.
956
Frontiers_manuscript/For_Submission/FrontiersinHarvardReview.cls
Normal file
956
Frontiers_manuscript/For_Submission/FrontiersinHarvardReview.cls
Normal file
@ -0,0 +1,956 @@
|
||||
\newcommand\classname{FrontiersinHarvard} %this file is based in bioninfo.cls
|
||||
\newcommand\lastmodifieddate{2016/11/10}
|
||||
\newcommand\versionnumber{2.6}
|
||||
|
||||
% For printing crop marks
|
||||
\newif\if@cropmarkson \@cropmarksontrue
|
||||
|
||||
\NeedsTeXFormat{LaTeX2e}[2013/26/03]
|
||||
\ProvidesClass{\classname}[\lastmodifieddate\space\versionnumber]
|
||||
|
||||
\setlength{\paperheight}{11truein}
|
||||
\setlength{\paperwidth}{8.5truein}
|
||||
|
||||
\newif\if@final
|
||||
|
||||
\DeclareOption{draft}{\PassOptionsToPackage{draft}{graphicx}}
|
||||
\DeclareOption{paper}{\PassOptionsToPackage{paper}{graphicx}}
|
||||
\DeclareOption{a4paper}{\PassOptionsToPackage{a4}{crop}}
|
||||
\DeclareOption{centre}{\PassOptionsToPackage{center}{crop}}
|
||||
\DeclareOption{crop}{\PassOptionsToPackage{cam}{crop}\global\@cropmarksontrue}
|
||||
\DeclareOption{nocrop}{\PassOptionsToPackage{off}{crop}\global\@cropmarksonfalse}
|
||||
\DeclareOption{info}{\PassOptionsToPackage{info}{crop}}
|
||||
\DeclareOption{noinfo}{\PassOptionsToPackage{noinfo}{crop}}
|
||||
\DeclareOption{final}{\global\@finaltrue}
|
||||
\DeclareOption{utf8}{\PassOptionsToPackage{utf8}{inputenc}}
|
||||
|
||||
\ExecuteOptions{a4paper,nocrop,centre,info,paper}
|
||||
|
||||
\ProcessOptions
|
||||
|
||||
% Load all necessary packages
|
||||
\RequirePackage{inputenc,crop,graphicx,amsmath,array,color,amssymb,flushend,stfloats,amsthm,chngpage,times,datetime,epstopdf}
|
||||
\usepackage{parskip}[=v1]
|
||||
\def\helvetica{\fontfamily{phv}\selectfont}
|
||||
\def\helveticaitalic{\fontfamily{phv}\itshape\selectfont}
|
||||
\def\helveticabold{\fontfamily{phv}\bfseries\selectfont}
|
||||
\def\helveticabolditalic{\fontfamily{phv}\bfseries\itshape\selectfont}
|
||||
\usepackage[usenames,dvipsnames]{xcolor}
|
||||
\usepackage[labelsep=quad,indention=10pt]{caption}
|
||||
\usepackage[labelfont=bf,list=true]{subcaption}
|
||||
|
||||
\newcommand\@ptsize{0}
|
||||
|
||||
% Set twoside printing
|
||||
\@twosidetrue
|
||||
|
||||
% Marginal notes are on the outside edge
|
||||
\@mparswitchfalse
|
||||
|
||||
\reversemarginpar
|
||||
|
||||
\renewcommand\normalsize{%
|
||||
\@setfontsize\normalsize{12}{12}%
|
||||
\abovedisplayskip 10\p@ \@plus2\p@ \@minus5\p@
|
||||
\abovedisplayshortskip \z@ \@plus3\p@
|
||||
\belowdisplayshortskip 6\p@ \@plus3\p@ \@minus3\p@
|
||||
\belowdisplayskip \abovedisplayskip
|
||||
\let\@listi\@listI}
|
||||
\normalsize
|
||||
\let\@bls\baselineskip
|
||||
|
||||
|
||||
\newcommand\small{%
|
||||
\@setfontsize\small{11}{12}%
|
||||
\abovedisplayskip 11\p@ minus 3\p@
|
||||
\belowdisplayskip \abovedisplayskip
|
||||
\abovedisplayshortskip \z@ plus 2\p@
|
||||
\belowdisplayshortskip 4\p@ plus 2\p@ minus2\p@
|
||||
\def\@listi{\topsep 4.5\p@ plus 2\p@ minus 1\p@
|
||||
\itemsep \parsep
|
||||
\topsep 4\p@ plus 2\p@ minus 2\p@}}
|
||||
|
||||
\newcommand\footnotesize{%
|
||||
\@setfontsize\footnotesize{8}{10}%
|
||||
\abovedisplayskip 6\p@ minus 3\p@
|
||||
\belowdisplayskip\abovedisplayskip
|
||||
\abovedisplayshortskip \z@ plus 3\p@
|
||||
\belowdisplayshortskip 6\p@ plus 3\p@ minus 3\p@
|
||||
\def\@listi{\topsep 3\p@ plus 1\p@ minus 1\p@
|
||||
\parsep 2\p@ plus 1\p@ minus 1\p@\itemsep \parsep}}
|
||||
|
||||
\def\scriptsize{\@setfontsize\scriptsize{11pt}{11pt}}
|
||||
\def\tiny{\@setfontsize\tiny{10pt}{11pt}}
|
||||
\def\large{\@setfontsize\large{13pt}{14pt}}
|
||||
\def\Large{\@setfontsize\Large{14pt}{16}}
|
||||
\def\LARGE{\@setfontsize\LARGE{15pt}{17pt}}
|
||||
\def\huge{\@setfontsize\huge{22pt}{22pt}}
|
||||
\def\Huge{\@setfontsize\Huge{30pt}{30pt}}
|
||||
|
||||
|
||||
\DeclareOldFontCommand{\rm}{\normalfont\rmfamily}{\mathrm}
|
||||
\DeclareOldFontCommand{\sf}{\normalfont\sffamily}{\mathsf}
|
||||
\DeclareOldFontCommand{\tt}{\normalfont\ttfamily}{\mathtt}
|
||||
\DeclareOldFontCommand{\bf}{\normalfont\bfseries}{\mathbf}
|
||||
\DeclareOldFontCommand{\it}{\normalfont\itshape}{\mathit}
|
||||
\DeclareOldFontCommand{\sl}{\normalfont\slshape}{\@nomath\sl}
|
||||
\DeclareOldFontCommand{\sc}{\normalfont\scshape}{\@nomath\sc}
|
||||
|
||||
% Line spacing
|
||||
\setlength\lineskip{1\p@}
|
||||
\setlength\normallineskip{1\p@}
|
||||
\renewcommand\baselinestretch{}
|
||||
|
||||
|
||||
% Set inter-para skips
|
||||
\setlength\smallskipamount{3\p@ \@plus 1\p@ \@minus 1\p@}
|
||||
\setlength\medskipamount{6\p@ \@plus 2\p@}
|
||||
\setlength\bigskipamount{12\p@ \@plus 4\p@ \@minus 4\p@}
|
||||
|
||||
% Paragraph dimensions and inter-para spacing
|
||||
\setlength\parindent{3mm}
|
||||
\setlength{\parskip}{\medskipamount}
|
||||
|
||||
% Page break penalties
|
||||
\@lowpenalty 51
|
||||
\@medpenalty 151
|
||||
\@highpenalty 301
|
||||
|
||||
% Disallow widows and orphans
|
||||
\clubpenalty 10000
|
||||
\widowpenalty 10000
|
||||
|
||||
% Disable page breaks before equations, allow pagebreaks after
|
||||
% equations and discourage widow lines before equations.
|
||||
\displaywidowpenalty 100
|
||||
\predisplaypenalty 10000
|
||||
\postdisplaypenalty 2500
|
||||
|
||||
% Allow breaking the page in the middle of a paragraph
|
||||
\interlinepenalty 0
|
||||
|
||||
% Disallow breaking the page after a hyphenated line
|
||||
\brokenpenalty 10000
|
||||
|
||||
% Hyphenation; don't split words into less than three characters
|
||||
\lefthyphenmin=3
|
||||
\righthyphenmin=3
|
||||
|
||||
%
|
||||
% Set page layout dimensions
|
||||
%
|
||||
\setlength\headheight{16\p@} % height of running head
|
||||
\setlength\topmargin{1.3pc} % head margin
|
||||
\addtolength\topmargin{-1in} % subtract out the 1 inch driver margin
|
||||
|
||||
\setlength\topskip{10\p@} % height of first line of text
|
||||
\setlength\headsep{19\p@} % space below running head --
|
||||
|
||||
\setlength\footskip{34\p@} % space above footer line
|
||||
\setlength\maxdepth{1\topskip} % pages can be short or deep by half a line?
|
||||
|
||||
\setlength\textwidth{42pc} % text measure excluding margins
|
||||
|
||||
\setlength\textheight{56\baselineskip} % 54 lines on a full page,
|
||||
\addtolength\textheight{\topskip} % including the first
|
||||
% line on the page
|
||||
|
||||
% Set the margins
|
||||
\setlength\marginparsep{3\p@}
|
||||
\setlength\marginparpush{3\p@}
|
||||
\setlength\marginparwidth{35\p@}
|
||||
|
||||
\setlength\oddsidemargin{4.5pc}
|
||||
\addtolength\oddsidemargin{-1in} % subtract out the 1 inch driver margin
|
||||
\setlength\@tempdima{\paperwidth}
|
||||
\addtolength\@tempdima{-\textwidth}
|
||||
\addtolength\@tempdima{-4.5pc}
|
||||
\setlength\evensidemargin{\@tempdima}
|
||||
\addtolength\evensidemargin{-1in}
|
||||
|
||||
\setlength\columnsep{1.5pc} % space between columns for double-column text
|
||||
\setlength\columnseprule{0\p@} % width of rule between two columns
|
||||
|
||||
% Footnotes
|
||||
\setlength\footnotesep{9\p@} % space between footnotes
|
||||
% space between text and footnote
|
||||
\setlength{\skip\footins}{12\p@ \@plus 6\p@ \@minus 1\p@}
|
||||
|
||||
% Float placement parameters
|
||||
|
||||
% The total number of floats that can be allowed on a page.
|
||||
\setcounter{totalnumber}{10}
|
||||
% The maximum number of floats at the top and bottom of a page.
|
||||
\setcounter{topnumber}{5}
|
||||
\setcounter{bottomnumber}{5}
|
||||
% The maximum part of the top or bottom of a text page that can be
|
||||
% occupied by floats. This is set so that at least four lines of text
|
||||
% fit on the page.
|
||||
\renewcommand\topfraction{.9}
|
||||
\renewcommand\bottomfraction{.9}
|
||||
% The minimum amount of a text page that must be occupied by text.
|
||||
% This should accomodate four lines of text.
|
||||
\renewcommand\textfraction{.06}
|
||||
% The minimum amount of a float page that must be occupied by floats.
|
||||
\renewcommand\floatpagefraction{.94}
|
||||
|
||||
% The same parameters repeated for double column output
|
||||
\renewcommand\dbltopfraction{.9}
|
||||
\renewcommand\dblfloatpagefraction{.9}
|
||||
|
||||
% Space between floats
|
||||
\setlength\floatsep {12\p@ \@plus 2\p@ \@minus 2\p@}
|
||||
% Space between floats and text
|
||||
\setlength\textfloatsep{20\p@ \@plus 2\p@ \@minus 4\p@}
|
||||
% Space above and below an inline figure
|
||||
\setlength\intextsep {18\p@ \@plus 2\p@ \@minus 2\p@}
|
||||
|
||||
% For double column floats
|
||||
\setlength\dblfloatsep {12\p@ \@plus 2\p@ \@minus 2\p@}
|
||||
\setlength\dbltextfloatsep{20\p@ \@plus 2\p@ \@minus 4\p@}
|
||||
|
||||
% Space left at top, bottom and inbetween floats on a float page.
|
||||
\setlength\@fptop{0\p@} % no space above float page figures
|
||||
\setlength\@fpsep{12\p@ \@plus 1fil}
|
||||
\setlength\@fpbot{0\p@}
|
||||
|
||||
% The same for double column
|
||||
\setlength\@dblfptop{0\p@}
|
||||
\setlength\@dblfpsep{12\p@ \@plus 1fil}
|
||||
\setlength\@dblfpbot{0\p@}
|
||||
|
||||
% Override settings in mathtime back to TeX defaults
|
||||
\DeclareMathSizes{5} {5} {5} {5}
|
||||
\DeclareMathSizes{6} {6} {5} {5}
|
||||
\DeclareMathSizes{7} {7} {5} {5}
|
||||
\DeclareMathSizes{8} {8} {6} {5}
|
||||
\DeclareMathSizes{9} {9} {6.5} {5}
|
||||
\DeclareMathSizes{10} {10} {7.5} {5}
|
||||
\DeclareMathSizes{12} {12} {9} {7}
|
||||
|
||||
% Page styles
|
||||
\def\ps@headings
|
||||
{%
|
||||
\def\@oddfoot{\vbox to 12.5\p@{\rule{\textwidth}{0.5\p@}\vss
|
||||
\hbox to \textwidth{\helveticabold\small {Frontiers}\hfill \thepage}%
|
||||
}}%
|
||||
\def\@evenfoot{\vbox to 12.5\p@{\rule{\textwidth}{0.5\p@}\vss
|
||||
\hbox to \textwidth{\helveticabold\small {Frontiers}\hfill \thepage}%
|
||||
}}%
|
||||
\def\@evenhead{\vbox{\hbox to \textwidth{\fontsize{11}{10}\selectfont
|
||||
\helveticabold{\fontshape{it}\selectfont
|
||||
\strut\firstAuthorLast}\hfill\strut\rightmark}\vspace{6.5\p@}\rule{\textwidth}{0.5\p@}}}%
|
||||
\def\@oddhead{\vbox{\hbox to \textwidth{\fontsize{11}{10}\selectfont
|
||||
\helveticabold{\fontshape{it}\selectfont
|
||||
\strut\firstAuthorLast}\hfill\strut\rightmark}\vspace{6.5\p@}\rule{\textwidth}{0.5\p@}}}%
|
||||
\def\titlemark##1{\markboth{##1}{##1}}%
|
||||
\def\authormark##1{\gdef\leftmark{##1}}%
|
||||
}
|
||||
|
||||
\def\ps@opening
|
||||
{%
|
||||
\def\@oddfoot{\vbox to 13\p@{\hbox{\rule{\textwidth}{1\p@}}\vss
|
||||
\hbox to \textwidth{\helvetica
|
||||
\fontsize{7}{9}\fontshape{n}\selectfont \space %
|
||||
\hfill\small\helveticabold\thepage}%
|
||||
}}%
|
||||
\def\@evenfoot{\vbox to 13\p@{\rule{\textwidth}{1\p@}\vss
|
||||
\hbox to \textwidth{\helvetica\thepage\hfill
|
||||
\fontsize{7}{9}\fontshape{n}\selectfont \space }%
|
||||
}}%
|
||||
\let\@evenhead\relax
|
||||
\let\@oddhead\relax}
|
||||
|
||||
% Page range
|
||||
\newif\iflastpagegiven \lastpagegivenfalse
|
||||
\newcommand\firstpage[1]{%
|
||||
\gdef\@firstpage{#1}%
|
||||
\ifnum\@firstpage>\c@page
|
||||
\setcounter{page}{#1}%
|
||||
\ClassWarning{BIO}{Increasing pagenumber to \@firstpage}%
|
||||
\else \ifnum\@firstpage<\c@page
|
||||
\ClassWarning{BIO}{Firstpage lower than pagenumber}\fi\fi
|
||||
\xdef\@firstpage{\the\c@page}%
|
||||
}
|
||||
\def\@firstpage{1}
|
||||
\def\pagenumbering#1{%
|
||||
\global\c@page \@ne
|
||||
\gdef\thepage{\csname @#1\endcsname \c@page}%
|
||||
\gdef\thefirstpage{%
|
||||
\csname @#1\endcsname \@firstpage}%
|
||||
\gdef\thelastpage{%
|
||||
\csname @#1\endcsname \@lastpage}%
|
||||
}
|
||||
|
||||
\newcommand\lastpage[1]{\xdef\@lastpage{#1}%
|
||||
\global\lastpagegiventrue}
|
||||
\def\@lastpage{0}
|
||||
\def\setlastpage{\iflastpagegiven\else
|
||||
\edef\@tempa{@lastpage@}%
|
||||
\expandafter
|
||||
\ifx \csname \@tempa \endcsname \relax
|
||||
\gdef\@lastpage{0}%
|
||||
\else
|
||||
\xdef\@lastpage{\@nameuse{@lastpage@}}%
|
||||
\fi
|
||||
\fi }
|
||||
\def\writelastpage{%
|
||||
\iflastpagegiven \else
|
||||
\immediate\write\@auxout%
|
||||
{\string\global\string\@namedef{@lastpage@}{\the\c@page}}%
|
||||
\fi
|
||||
}
|
||||
\def\thepagerange{%
|
||||
\ifnum\@lastpage =0 {\ \bf ???} \else
|
||||
\ifnum\@lastpage = \@firstpage \ \thefirstpage\else
|
||||
\thefirstpage--\thelastpage \fi\fi}
|
||||
|
||||
\AtBeginDocument{\setlastpage
|
||||
\pagenumbering{arabic}%
|
||||
}
|
||||
\AtEndDocument{%
|
||||
\writelastpage
|
||||
\if@final
|
||||
\clearemptydoublepage
|
||||
\else
|
||||
\clearpage
|
||||
\fi}
|
||||
|
||||
|
||||
%
|
||||
% Sectional units
|
||||
%
|
||||
|
||||
% Counters
|
||||
\newcounter{section}
|
||||
\newcounter{subsection}[section]
|
||||
\newcounter{subsubsection}[subsection]
|
||||
\newcounter{paragraph}[subsubsection]
|
||||
\newcounter{subparagraph}[paragraph]
|
||||
\newcounter{figure}
|
||||
\newcounter{subfigure}
|
||||
\newcounter{table}
|
||||
|
||||
% Form of the numbers
|
||||
\newcommand\thepage{\arabic{page}}
|
||||
\renewcommand\thesection{\arabic{section}}
|
||||
\renewcommand\thesubsection{{\thesection.\arabic{subsection}}}
|
||||
\renewcommand\thesubsubsection{{\thesubsection.\arabic{subsubsection}}}
|
||||
\renewcommand\theparagraph{\thesubsubsection.\arabic{paragraph}}
|
||||
\renewcommand\thesubparagraph{\theparagraph.\arabic{subparagraph}}
|
||||
\renewcommand\theequation{\arabic{equation}}
|
||||
|
||||
% Form of the words
|
||||
\newcommand\contentsname{Contents}
|
||||
\newcommand\listfigurename{List of Figures}
|
||||
\newcommand\listtablename{List of Tables}
|
||||
\newcommand\partname{Part}
|
||||
\newcommand\appendixname{Appendix}
|
||||
\newcommand\abstractname{Abstract}
|
||||
\newcommand\refname{References}
|
||||
\newcommand\bibname{References}
|
||||
\newcommand\indexname{Index}
|
||||
\newcommand\figurename{Figure}
|
||||
\newcommand\tablename{Table}
|
||||
|
||||
% Clearemptydoublepage should really clear the running heads too
|
||||
\newcommand{\clearemptydoublepage}{\newpage{\pagestyle{empty}\cleardoublepage}}
|
||||
|
||||
% Frontmatter, mainmatter and backmatter
|
||||
|
||||
\newif\if@mainmatter \@mainmattertrue
|
||||
|
||||
\newcommand\frontmatter{%
|
||||
\clearpage
|
||||
\@mainmatterfalse
|
||||
\pagenumbering{roman}}
|
||||
|
||||
\newcommand\mainmatter{%
|
||||
\clearpage
|
||||
\@mainmattertrue
|
||||
\pagenumbering{arabic}}
|
||||
|
||||
\newcommand\backmatter{%
|
||||
\clearpage
|
||||
\@mainmatterfalse}
|
||||
|
||||
\newdateformat{mydate}{\THEDAY \ \monthname[\THEMONTH] \THEYEAR}
|
||||
\newdateformat{mydate2}{\monthname[\THEMONTH] \THEYEAR}
|
||||
|
||||
%********
|
||||
% TITLE *
|
||||
%********
|
||||
\newlength{\dropfromtop}
|
||||
\setlength{\dropfromtop}{\z@}
|
||||
|
||||
% Application Notes
|
||||
\newif\if@appnotes
|
||||
\newcommand{\application}{%
|
||||
% \setlength{\dropfromtop}{-2.25pc}%
|
||||
\global\@appnotestrue}
|
||||
|
||||
% Short Title is for the Title of the Article
|
||||
% Title is for the Title of the Journal
|
||||
|
||||
\long\def\title{\@ifnextchar[{\short@title}{\@@title}}
|
||||
\def\short@title[#1]{\titlemark{#1}\@@@title}
|
||||
\def\@@title#1{\authormark{#1}\@@@title{#1}}
|
||||
\long\def\@@@title#1{\gdef\@title{#1}}
|
||||
|
||||
\long\def\author{\@ifnextchar[{\short@uthor}{\@uthor}}
|
||||
\def\short@uthor[#1]{\authormark{#1}\@@author}
|
||||
\def\@uthor#1{\authormark{#1}\@@author{#1}}
|
||||
\long\def\@@author#1{\gdef\@author{#1}}
|
||||
|
||||
\def\vol#1{\global\def\@vol{#1}}
|
||||
\def\issue#1{\global\def\@issue{#1}}
|
||||
%\def\address#1{\global\def\@issue{#1}}
|
||||
\def\history#1{\global\def\@history{#1}}
|
||||
%\def\correspondance#1{\global\def\@correspondance{#1}}
|
||||
%\def\extraAuth#1{\global\def\@extraAuth{#1}}
|
||||
%\def\address#1{\global\def\@address{#1}}
|
||||
|
||||
|
||||
\def\date#1{\gdef\@date{#1}}
|
||||
|
||||
%***********
|
||||
% Colors *
|
||||
%***********
|
||||
|
||||
% Other colors
|
||||
\definecolor{black}{cmyk}{0, 0, 0, 1}
|
||||
\definecolor{gray}{cmyk}{0, 0, 0, 0.25}
|
||||
\definecolor{darkgray}{cmyk}{0, 0, 0, 0.5}
|
||||
\newlength{\extraspace}
|
||||
\setlength{\extraspace}{\z@}
|
||||
|
||||
\newcommand\maketitle{\par
|
||||
\begingroup
|
||||
\renewcommand\thefootnote{\@fnsymbol\c@footnote}%
|
||||
\def\@makefnmark{\rlap{\@textsuperscript{\normalfont\@thefnmark}}}%
|
||||
\long\def\@makefntext##1{\parindent 3mm\noindent
|
||||
% \@textsuperscript{\normalfont\@thefnmark}\raggedright##1}%
|
||||
\@textsuperscript{\normalfont\@thefnmark}##1}%
|
||||
\if@twocolumn
|
||||
\ifnum \col@number=\@ne
|
||||
\@maketitle
|
||||
\else
|
||||
\twocolumn[\@maketitle]%
|
||||
\fi
|
||||
\else
|
||||
\newpage
|
||||
\global\@topnum\z@ % Prevents figures from going at top of page.
|
||||
\@maketitle
|
||||
\fi
|
||||
\thispagestyle{opening}\@thanks
|
||||
\endgroup
|
||||
\setcounter{footnote}{0}%
|
||||
\global\let\thanks\relax
|
||||
\global\let\maketitle\relax
|
||||
\global\let\@maketitle\relax
|
||||
% \global\let\@address\@empty
|
||||
\global\let\@history\@empty
|
||||
% \global\let\@extraAuth\@empty
|
||||
\global\let\@topic\@empty
|
||||
\global\let\@thanks\@empty
|
||||
\global\let\@author\@empty
|
||||
\global\let\@date\@empty
|
||||
\global\let\@title\@empty
|
||||
\global\let\@pubyear\@empty
|
||||
% \global\let\address\relax
|
||||
\global\let\topic\relax
|
||||
% \global\let\extraAuth\relax
|
||||
\global\let\title\relax
|
||||
\global\let\author\relax
|
||||
\global\let\date\relax
|
||||
\global\let\pubyear\relax
|
||||
\global\let\@copyrightline\@empty
|
||||
\global\let\and\relax
|
||||
\@afterindentfalse\@afterheading
|
||||
}
|
||||
|
||||
\newlength{\aboveskipchk}%for checking oddpage or evenpage top skip
|
||||
\setlength{\aboveskipchk}{\z@}%
|
||||
|
||||
\def\@maketitle{%
|
||||
\let\footnote\thanks
|
||||
\clearemptydoublepage
|
||||
\checkoddpage\ifcpoddpage\setlength{\aboveskipchk}{-7pc}\else\setlength{\aboveskipchk}{-3pc}\fi%for checking oddpage or evenpage top skip%%
|
||||
\vspace*{\aboveskipchk}%
|
||||
\vspace{\dropfromtop}%
|
||||
\hbox to \textwidth{%
|
||||
|
||||
\parbox[t]{15pc}{%
|
||||
\helvetica
|
||||
\hfil
|
||||
\flushleft \includegraphics[width=12pc,angle=0]{./logo1.eps}
|
||||
|
||||
|
||||
}}
|
||||
|
||||
\rule{\textwidth}{1\p@}\par%
|
||||
\helvetica
|
||||
\hbox to \textwidth{%
|
||||
\parbox[t]{36.5pc}{%
|
||||
\vspace*{1sp}
|
||||
{\helveticabold\fontsize{20}{21}\color{black}\selectfont\raggedright \@title \par}%
|
||||
\vspace{4.5\p@}
|
||||
{\helveticabold\fontsize{12}{15}\selectfont\raggedright \@author \par}%
|
||||
\vspace{4\p@}
|
||||
% {\helvetica\fontsize{12}{12}\selectfont\raggedright\slshape\@address \Address \\ \par}%
|
||||
% \vspace{6\p@}
|
||||
% {\helvetica\fontsize{12}{10}\selectfont\raggedright {Correspondence*:\\ }\@correspondance \corrAuthor \\ \corrEmail \par}
|
||||
% \vspace{4\p@}
|
||||
% {\helvetica\fontsize{12}{12}\selectfont\raggedright\@extraAuth \par}%
|
||||
% \vspace{8\p@}
|
||||
}%
|
||||
}
|
||||
\vspace{14.5\p@}%
|
||||
}
|
||||
|
||||
%***********
|
||||
% Abstract *
|
||||
%***********
|
||||
\newcommand{\absection}[1]{%
|
||||
\par\noindent{\bfseries #1}\space\ignorespaces}
|
||||
|
||||
\newenvironment{abstract}{%
|
||||
\begingroup
|
||||
\let\section\absection
|
||||
\fontfamily{\sfdefault}\fontsize{12}{12}\sffamily\selectfont
|
||||
{\fontseries{b}\selectfont ABSTRACT}\par}
|
||||
{\endgroup\smallskip\@afterheading\@afterindentfalse\vskip 12pt plus -2pt minus 2pt}
|
||||
|
||||
% Section macros
|
||||
|
||||
% Lowest level heading that takes a number by default
|
||||
\setcounter{secnumdepth}{5}
|
||||
|
||||
\renewcommand{\@seccntformat}[1]{\csname the#1\endcsname\quad}
|
||||
|
||||
\def\section{%
|
||||
\@startsection{section}{1}{\z@}
|
||||
{-10\p@ plus -3\p@}{3\p@}
|
||||
{\reset@font\raggedright\helveticabold\fontsize{13}{13}\color{black}\selectfont\MakeUppercase}}
|
||||
|
||||
\def\subsection{%
|
||||
\@startsection{subsection}{2}{\z@}
|
||||
{4\p@ plus 2\p@}{2\p@}
|
||||
{\reset@font\raggedright\mathversion{bold}\helveticabold\fontsize{11.5}{12}}}
|
||||
|
||||
\def\subsubsection{%
|
||||
\@startsection{subsubsection}{3}{\z@}
|
||||
{2\p@ plus 2\p@}{1\p@}
|
||||
{\reset@font\raggedright\mathversion{bold}\helvetica\fontsize{11.5}{12}}}
|
||||
|
||||
\def\paragraph{%
|
||||
\@startsection{paragraph}{4}{\z@}
|
||||
{2\p@ plus 1\p@}{1\p@}
|
||||
{\reset@font\itshape\helveticabold\fontsize{11.5}{12}}}
|
||||
|
||||
\def\subparagraph{%
|
||||
\@startsection{subparagraph}{5}{\z@}
|
||||
{2\p@ plus 1\p@}{1\p@}
|
||||
{\reset@font\itshape\helvetica\fontsize{11.5}{12}}}
|
||||
|
||||
\def\textcolon{\text{\rm :}}
|
||||
|
||||
% ********************
|
||||
% Figures and tables *
|
||||
% ********************
|
||||
|
||||
% Table and array parameters
|
||||
\setlength\arraycolsep{.5em}
|
||||
\setlength\tabcolsep{.5em}
|
||||
\setlength\arrayrulewidth{.5pt}
|
||||
\setlength\doublerulesep{2.5pt}
|
||||
\setlength\extrarowheight{\z@}
|
||||
\renewcommand\arraystretch{1}
|
||||
|
||||
%\newlength{\abovecaptionskip}
|
||||
%\newlength{\belowcaptionskip}
|
||||
%\setlength{\abovecaptionskip}{13pt}
|
||||
%\setlength{\belowcaptionskip}{10.5pt}
|
||||
|
||||
\long\def\@makecaption#1#2{\vspace{\abovecaptionskip}%
|
||||
\begingroup
|
||||
\normalsize
|
||||
\textbf{#1.}\enskip{#2}\par
|
||||
\endgroup}
|
||||
|
||||
|
||||
% Table rules
|
||||
\def\toprule{\noalign{\ifnum0=`}\fi\hrule \@height 0.5pt \hrule \@height 6pt \@width 0pt \futurelet
|
||||
\@tempa\@xhline}
|
||||
\def\midrule{\noalign{\ifnum0=`}\fi \hrule \@height 6.75pt \@width 0pt \hrule \@height 0.5pt
|
||||
\hrule \@height 6pt \@width 0pt \futurelet \@tempa\@xhline}
|
||||
\def\botrule{\noalign{\ifnum0=`}\fi \hrule \@height 5.75pt \@width 0pt \hrule \@height 0.5pt \futurelet
|
||||
\@tempa\@xhline}
|
||||
\def\hrulefill{\leavevmode\leaders\hrule height .5pt\hfill\kern\z@}
|
||||
|
||||
\def\thefigure{\@arabic\c@figure}
|
||||
\def\fps@figure{tbp}
|
||||
\def\ftype@figure{1}
|
||||
\def\ext@figure{lof}
|
||||
\def\fnum@figure{\figurename~\thefigure}
|
||||
\def\figure{\@float{figure}}
|
||||
\let\endfigure\end@float
|
||||
\@namedef{figure*}{\@dblfloat{figure}}
|
||||
\@namedef{endfigure*}{\end@dblfloat}
|
||||
\def\thetable{\@arabic\c@table}
|
||||
\def\fps@table{tbp}
|
||||
\def\ftype@table{2}
|
||||
\def\ext@table{lot}
|
||||
\def\fnum@table{Table~\thetable}
|
||||
\def\table{\let\source\tablesource\@float{table}}
|
||||
\def\endtable{\end@float}
|
||||
\@namedef{table*}{\@dblfloat{table}}
|
||||
\@namedef{endtable*}{\end@dblfloat}
|
||||
|
||||
\newif\if@rotate \@rotatefalse
|
||||
\newif\if@rotatecenter \@rotatecenterfalse
|
||||
\def\rotatecenter{\global\@rotatecentertrue}
|
||||
\def\rotateendcenter{\global\@rotatecenterfalse}
|
||||
\def\rotate{\global\@rotatetrue}
|
||||
\def\endrotate{\global\@rotatefalse}
|
||||
\newdimen\rotdimen
|
||||
\def\rotstart#1{\special{ps: gsave currentpoint currentpoint translate
|
||||
#1 neg exch neg exch translate}}
|
||||
\def\rotfinish{\special{ps: currentpoint grestore moveto}}
|
||||
\def\rotl#1{\rotdimen=\ht#1\advance\rotdimen by \dp#1
|
||||
\hbox to \rotdimen{\vbox to\wd#1{\vskip \wd#1
|
||||
\rotstart{270 rotate}\box #1\vss}\hss}\rotfinish}
|
||||
\def\rotr#1{\rotdimen=\ht #1\advance\rotdimen by \dp#1
|
||||
\hbox to \rotdimen{\vbox to \wd#1{\vskip \wd#1
|
||||
\rotstart{90 rotate}\box #1\vss}\hss}\rotfinish}
|
||||
|
||||
\newif\ifsub@@figure \sub@@figuretrue
|
||||
\def\thesubfigure{\@arabic\c@figure\alph{subfigure}}
|
||||
\def\fps@subfigure{tbp}
|
||||
\def\ftype@subfigure{1}
|
||||
\def\ext@subfigure{lof}
|
||||
\def\fnum@subfigure{Figure \thesubfigure}
|
||||
\def\subfigure{\ifsub@@figure\global\sub@@figurefalse\stepcounter{figure}\fi%
|
||||
\@float{subfigure}}
|
||||
\let\endsubfigure\end@float
|
||||
\@namedef{subfigure*}{\ifsub@@figure\global\sub@@figurefalse%
|
||||
\stepcounter{figure}\fi \@dblfloat{subfigure}}
|
||||
\@namedef{endsubfigure*}{\end@dblfloat}
|
||||
|
||||
\def\endfigure{\global\sub@@figuretrue\end@float}
|
||||
\@namedef{endfigure*}{\global\sub@@figuretrue\end@dblfloat}
|
||||
|
||||
\newdimen\tempdime
|
||||
\newbox\temptbox
|
||||
|
||||
% From ifmtarg.sty
|
||||
% Copyright Peter Wilson and Donald Arseneau, 2000
|
||||
\begingroup
|
||||
\catcode`\Q=3
|
||||
\long\gdef\@ifmtarg#1{\@xifmtarg#1QQ\@secondoftwo\@firstoftwo\@nil}
|
||||
\long\gdef\@xifmtarg#1#2Q#3#4#5\@nil{#4}
|
||||
\long\gdef\@ifnotmtarg#1{\@xifmtarg#1QQ\@firstofone\@gobble\@nil}
|
||||
\endgroup
|
||||
|
||||
\def\tablesize{\@setfontsize\tablesize{10\p@}{10\p@}}
|
||||
|
||||
\newenvironment{processtable}[3]{\setbox\temptbox=\hbox{{\tablesize #2}}%
|
||||
\tempdime\wd\temptbox\@processtable{#1}{#2}{#3}{\tempdime}}
|
||||
{\relax}
|
||||
|
||||
\newcommand{\@processtable}[4]{%
|
||||
\if@rotate
|
||||
\setbox4=\vbox to \hsize{\vss\hbox to \textheight{%
|
||||
\begin{minipage}{#4}%
|
||||
\@ifmtarg{#1}{}{\caption{#1}}{\tablesize #2}%
|
||||
\vskip7\p@\noindent
|
||||
\parbox{#4}{\fontsize{10}{10}\selectfont #3\par}%
|
||||
\end{minipage}}\vss}%
|
||||
\rotr{4}
|
||||
\else
|
||||
\hbox to \hsize{\hss\begin{minipage}[t]{#4}%
|
||||
\vskip2.9pt
|
||||
\@ifmtarg{#1}{}{\caption{#1}}{\tablesize #2}%
|
||||
\vskip6\p@\noindent
|
||||
\parbox{#4}{\fontsize{8}{9}\selectfont #3\par}%
|
||||
\end{minipage}\hss}\fi}%
|
||||
|
||||
\newcolumntype{P}[1]{>{\raggedright\let\\\@arraycr\hangindent1em}p{#1}}
|
||||
|
||||
% ******************************
|
||||
% List numbering and lettering *
|
||||
% ******************************
|
||||
\def\labelenumi{{\rm\arabic{enumi}.}}
|
||||
\def\theenumi{\arabic{enumi}}
|
||||
\def\labelenumii{{\rm\alph{enumii}.}}
|
||||
\def\theenumii{\alph{enumii}}
|
||||
\def\p@enumii{\theenumi}
|
||||
\def\labelenumiii{{\rm(\arabic{enumiii})}}
|
||||
\def\theenumiii{\roman{enumiii}}
|
||||
\def\p@enumiii{\theenumi(\theenumii)}
|
||||
\def\labelenumiv{{\rm(\arabic{enumiv})}}
|
||||
\def\theenumiv{\Alph{enumiv}}
|
||||
\def\p@enumiv{\p@enumiii\theenumiii}
|
||||
\def\labelitemi{{\small$\bullet$}}
|
||||
\def\labelitemii{{\small$\bullet$}}
|
||||
\def\labelitemiii{{\small$\bullet$}}
|
||||
\def\labelitemiv{{\small$\bullet$}}
|
||||
|
||||
\def\@listI{\leftmargin\leftmargini \topsep\medskipamount}
|
||||
\let\@listi\@listI
|
||||
\@listi
|
||||
\def\@listii{\topsep\z@\leftmargin\leftmarginii}
|
||||
\def\@listiii{\leftmargin\leftmarginiii \topsep\z@}
|
||||
\def\@listiv{\leftmargin\leftmarginiv \topsep\z@}
|
||||
\def\@listv{\leftmargin\leftmarginv \topsep\z@}
|
||||
\def\@listvi{\leftmargin\leftmarginvi \topsep\z@}
|
||||
|
||||
\setlength{\leftmargini}{3mm}
|
||||
\setlength{\leftmarginii}{\z@}
|
||||
\setlength{\leftmarginiii}{\z@}
|
||||
\setlength{\leftmarginiv}{\z@}
|
||||
|
||||
% Changes to the list parameters for enumerate
|
||||
\def\enumargs{%
|
||||
\partopsep \z@
|
||||
\itemsep 3\p@
|
||||
\parsep \z@
|
||||
\labelsep 0.5em
|
||||
\listparindent \parindent
|
||||
\itemindent \z@
|
||||
\topsep 2\p@
|
||||
}
|
||||
|
||||
\def\enumerate{%
|
||||
\@ifnextchar[{\@numerate}{\@numerate[0]}}
|
||||
|
||||
\def\@numerate[#1]{%
|
||||
\ifnum \@enumdepth >3 \@toodeep\else
|
||||
\advance\@enumdepth \@ne
|
||||
\edef\@enumctr{enum\romannumeral\the\@enumdepth}
|
||||
\list{\csname label\@enumctr\endcsname}{%
|
||||
\enumargs
|
||||
\setlength{\leftmargin}{\csname leftmargin\romannumeral\the\@enumdepth\endcsname}
|
||||
\usecounter{\@enumctr}
|
||||
\settowidth\labelwidth{#1}
|
||||
\addtolength{\leftmargin}{\labelwidth}
|
||||
\addtolength{\leftmargin}{\labelsep}
|
||||
\def\makelabel##1{\hss \llap{##1}}}%
|
||||
\fi
|
||||
}
|
||||
\let\endenumerate\endlist
|
||||
|
||||
% Changes to the list parameters for itemize
|
||||
\def\itemargs{%
|
||||
\partopsep \z@
|
||||
\itemsep 3\p@
|
||||
\parsep \z@
|
||||
\labelsep 0.5em
|
||||
\rightmargin \z@
|
||||
\listparindent \parindent
|
||||
\itemindent \z@
|
||||
\topsep 2\p@
|
||||
}
|
||||
|
||||
\def\itemize{%
|
||||
\@ifnextchar[{\@itemize}{\@itemize[$\bullet$]}}
|
||||
|
||||
\def\@itemize[#1]{%
|
||||
\ifnum \@itemdepth >3 \@toodeep\else
|
||||
\advance\@itemdepth \@ne
|
||||
\edef\@itemctr{item\romannumeral\the\@itemdepth}
|
||||
\list{\csname label\@itemctr\endcsname}{%
|
||||
\itemargs
|
||||
\setlength{\leftmargin}{\csname leftmargin\romannumeral\the\@itemdepth\endcsname}
|
||||
\settowidth\labelwidth{#1}
|
||||
\addtolength{\leftmargin}{\labelwidth}
|
||||
\addtolength{\leftmargin}{\labelsep}
|
||||
\def\makelabel##1{\hss \llap{##1}}}%
|
||||
\fi
|
||||
}
|
||||
\let\enditemize\endlist
|
||||
|
||||
\newenvironment{unlist}{%
|
||||
\begin{list}{}%
|
||||
{\setlength{\labelwidth}{\z@}%
|
||||
\setlength{\labelsep}{\z@}%
|
||||
\setlength{\topsep}{\medskipamount}%
|
||||
\setlength{\itemsep}{3\p@}%
|
||||
\setlength{\leftmargin}{2em}%
|
||||
\setlength{\itemindent}{-2em}}}
|
||||
{\end{list}}
|
||||
|
||||
|
||||
% ***********************
|
||||
% Quotes and Quotations *
|
||||
% ***********************
|
||||
\def\quotation{\par\begin{list}{}{
|
||||
\setlength{\topsep}{\medskipamount}
|
||||
\setlength{\leftmargin}{2em}%
|
||||
\setlength{\rightmargin}{\z@}%
|
||||
\setlength\labelwidth{0pt}%
|
||||
\setlength\labelsep{0pt}%
|
||||
\listparindent\parindent}%
|
||||
\item[]}
|
||||
\def\endquotation{\end{list}}
|
||||
\let\quote\quotation
|
||||
\let\endquote\endquotation
|
||||
|
||||
\skip\@mpfootins = \skip\footins
|
||||
\fboxsep=6\p@
|
||||
\fboxrule=1\p@
|
||||
|
||||
% *******************
|
||||
% Table of contents *
|
||||
% *******************
|
||||
\newcommand\@pnumwidth{4em}
|
||||
\newcommand\@tocrmarg{2.55em plus 1fil}
|
||||
\newcommand\@dotsep{1000}
|
||||
\setcounter{tocdepth}{4}
|
||||
|
||||
\def\numberline#1{\hbox to \@tempdima{{#1}}}
|
||||
|
||||
\def\@authortocline#1#2#3#4#5{%
|
||||
\vskip 1.5\p@
|
||||
\ifnum #1>\c@tocdepth \else
|
||||
{\leftskip #2\relax \rightskip \@tocrmarg \parfillskip -\rightskip
|
||||
\parindent #2\relax\@afterindenttrue
|
||||
\interlinepenalty\@M
|
||||
\leavevmode
|
||||
\@tempdima #3\relax
|
||||
\advance\leftskip \@tempdima \null\nobreak\hskip -\leftskip
|
||||
{\itshape #4}\nobreak
|
||||
\leaders\hbox{$\m@th
|
||||
\mkern \@dotsep mu\hbox{.}\mkern \@dotsep
|
||||
mu$}\hfill
|
||||
\nobreak
|
||||
\hb@xt@\@pnumwidth{\hfil}%
|
||||
\par}%
|
||||
\fi}
|
||||
|
||||
\newcommand*\l@author{\@authortocline{2}{0pt}{30pt}}
|
||||
\newcommand*\l@section{\@dottedtocline{3}{11pt}{20pt}}
|
||||
\newcommand*\l@subsection{\@dottedtocline{4}{31pt}{29pt}}
|
||||
\newcommand*\l@subsubsection[2]{}
|
||||
|
||||
|
||||
|
||||
% ***********
|
||||
% Footnotes *
|
||||
% ***********
|
||||
|
||||
\def\footnoterule{\noindent\rule{\columnwidth}{0.5pt}}
|
||||
\def\@makefnmark{\@textsuperscript{\normalfont\@thefnmark}}%
|
||||
\newcommand\@makefntext[1]{\noindent{\@makefnmark}\enskip#1}
|
||||
|
||||
% ***********
|
||||
% References *
|
||||
% ***********
|
||||
|
||||
\providecommand{\newblock}{}
|
||||
\newenvironment{thebibliography}{%
|
||||
\section{\bibname}%
|
||||
\begingroup
|
||||
\small
|
||||
\begin{list}{}{%
|
||||
\setlength{\topsep}{\z@}%
|
||||
\setlength{\labelsep}{\z@}%
|
||||
\settowidth{\labelwidth}{\z@}%
|
||||
\setlength{\leftmargin}{4mm}%
|
||||
\setlength{\itemindent}{-4mm}}\small}
|
||||
{\end{list}\endgroup}
|
||||
|
||||
\RequirePackage{natbib}
|
||||
|
||||
|
||||
% **********
|
||||
% Appendix *
|
||||
% **********
|
||||
\newif\ifappend
|
||||
\def\appendix{\par
|
||||
\setcounter{section}{0}
|
||||
\setcounter{subsection}{0}
|
||||
\appendtrue
|
||||
}
|
||||
|
||||
%Math parameters
|
||||
|
||||
\setlength{\jot}{5\p@}
|
||||
\mathchardef\@m=1500 % adapted value
|
||||
|
||||
\def\frenchspacing{\sfcode`\.\@m \sfcode`\?\@m \sfcode`\!\@m
|
||||
\sfcode`\:\@m \sfcode`\;\@m \sfcode`\,\@m}
|
||||
|
||||
% Theorems
|
||||
\def\th@plain{%
|
||||
%% \let\thm@indent\noindent % no indent
|
||||
\thm@headfont{\quad\scshape}% heading font is bold
|
||||
\thm@notefont{\upshape\mdseries}% same as heading font
|
||||
\thm@headpunct{.}% no period after heading
|
||||
\thm@headsep 5\p@ plus\p@ minus\p@\relax
|
||||
\itshape % body font
|
||||
}
|
||||
|
||||
\vbadness=9999
|
||||
\tolerance=9999
|
||||
\doublehyphendemerits=10000
|
||||
\doublehyphendemerits 640000 % corresponds to badness 800
|
||||
\finalhyphendemerits 1000000 % corresponds to badness 1000
|
||||
|
||||
\flushbottom
|
||||
\frenchspacing
|
||||
\ps@headings
|
||||
\twocolumn
|
||||
|
||||
% Screen PDF compatability
|
||||
\newcommand{\medline}[1]{%
|
||||
\unskip\unskip\ignorespaces}
|
||||
|
||||
|
||||
%%%%for smaller size text
|
||||
\newenvironment{methods}{%
|
||||
\begingroup
|
||||
\def\section{%
|
||||
\@startsection{section}{1}{\z@}
|
||||
{-24\p@ plus -3\p@}{4\p@}
|
||||
{\reset@font\raggedright\helveticabold\fontsize{10}{12}\selectfont\MakeUppercase}}
|
||||
\def\subsection{%
|
||||
\@startsection{subsection}{2}{\z@}
|
||||
{-5\p@ plus -2\p@}{4\p@}
|
||||
{\reset@font\raggedright\mathversion{bold}\fontseries{b}\fontsize{10}{12}\selectfont}}
|
||||
\def\subsubsection{%
|
||||
\@startsection{subsubsection}{3}{\z@}
|
||||
{-6\p@ plus -1\p@}{-1em}
|
||||
{\reset@font\normalfont\normalsize\itshape}}
|
||||
\footnotesize
|
||||
\par}
|
||||
{\par\endgroup\bigskip\@afterheading\@afterindentfalse}
|
||||
|
||||
|
||||
|
||||
\graphicspath{{g:/artwork/oup/bioinfo/}}
|
||||
|
||||
\language=2
|
||||
|
||||
\hyphenation{Figure Table Figures Tables}
|
||||
|
||||
\newcommand{\href}[2]{#2}
|
||||
|
||||
\renewenvironment{proof}[1][\proofname]{\par
|
||||
\normalfont \topsep6\p@\@plus6\p@\relax
|
||||
\labelsep 0.5em
|
||||
\trivlist
|
||||
\item[\hskip\labelsep\hskip1em\textsc{#1}.]\ignorespaces
|
||||
}{\endtrivlist\@endpefalse}
|
||||
|
||||
%%Different Bonds
|
||||
|
||||
\def\sbond{\ensuremath{\raise.25ex\hbox{${-}\!\!\!\!{-}$}}\kern -.9pt}
|
||||
\def\dbond{\ensuremath{\raise.25ex\hbox{=$\!$=}}}
|
||||
\def\tbond{\ensuremath{\raise.20ex\hbox{${\equiv}\!\!\!{\equiv}$}}}
|
||||
|
||||
\newcommand{\query}[2][0pt]{}%
|
||||
|
||||
\renewcommand{\dag}{{\mathversion{normal}$^{\dagger}$}}
|
||||
|
||||
\endinput
|
||||
BIN
Frontiers_manuscript/For_Submission/Koch_frontiers.pdf
Normal file
BIN
Frontiers_manuscript/For_Submission/Koch_frontiers.pdf
Normal file
Binary file not shown.
BIN
Frontiers_manuscript/For_Submission/Koch_frontiers.synctex.gz
Normal file
BIN
Frontiers_manuscript/For_Submission/Koch_frontiers.synctex.gz
Normal file
Binary file not shown.
319
Frontiers_manuscript/For_Submission/Koch_frontiers.tex
Normal file
319
Frontiers_manuscript/For_Submission/Koch_frontiers.tex
Normal file
@ -0,0 +1,319 @@
|
||||
\documentclass[utf8]{FrontiersinHarvard}
|
||||
\DeclareUnicodeCharacter{03B2}{\(\beta\)}
|
||||
\DeclareUnicodeCharacter{03B1}{\(\alpha\)}
|
||||
\DeclareUnicodeCharacter{00C5}{\AA}
|
||||
\usepackage[english]{babel}
|
||||
\usepackage{url,hyperref,lineno,microtype,subcaption}
|
||||
\usepackage[onehalfspacing]{setspace}
|
||||
\usepackage{amsmath}
|
||||
|
||||
\usepackage{xspace}
|
||||
\usepackage{longtable}
|
||||
\usepackage{boldline}
|
||||
\usepackage{geometry}
|
||||
\usepackage{makecell}
|
||||
\usepackage{float}
|
||||
\usepackage{placeins}
|
||||
|
||||
\usepackage[tableposition=top]{caption}
|
||||
|
||||
|
||||
\newcommand{\Kv}{\(\text{K}_{\text{V}}\text{1.1}\)\xspace}
|
||||
\newcommand{\IKv}{\(\text{I}_{\text{K}_{\text{V}}\text{1.1}}\)\xspace}
|
||||
\newcommand{\drheo}{\(\Delta\)rheobase\xspace}
|
||||
\newcommand{\ndAUC}{normalized \(\Delta\)AUC\xspace}
|
||||
|
||||
\newcommand{\note}[2][]{\textbf{[#1: #2]}}
|
||||
\newcommand{\notenk}[1]{\note[NK]{#1}}
|
||||
\newcommand{\notels}[1]{\note[LS]{#1}}
|
||||
\newcommand{\notejb}[1]{\note[JB]{#1}}
|
||||
|
||||
\linenumbers
|
||||
|
||||
\def\keyFont{\fontsize{8}{11}\helveticabold }
|
||||
\def\firstAuthorLast{Koch {et~al.}} %use et al only if is more than 1 author
|
||||
\def\Authors{Nils A. Koch\,$^{1,2}$, Lukas Sonnenberg\,$^{1,2}$, Ulrike B.S. Hedrich\,$^{3}$, Stephan Lauxmann\,$^{1,3}$ and Jan Benda\,$^{1,2,*}$}
|
||||
\def\Address{$^{1}$Institute for Neurobiology, University of T{\"u}bingen, 72072 T{\"u}bingen, Germany \\
|
||||
$^{2}$Bernstein Center for Computational Neuroscience T{\"u}bingen, 72076 T{\"u}bingen, Germany \\
|
||||
$^{3}$Department of Neurology and Epileptology, Hertie Institute for Clinical Brain Research, University of T{\"u}bingen, 72076 T{\"u}bingen, Germany}
|
||||
|
||||
% The Corresponding Author should be marked with an asterisk
|
||||
% Provide the exact contact address (this time including street name and city zip code) and email of the corresponding author
|
||||
\def\corrAuthor{Jan Benda}
|
||||
|
||||
\def\corrEmail{ jan.benda@uni-tuebingen.de}
|
||||
|
||||
|
||||
|
||||
\begin{document}
|
||||
\onecolumn
|
||||
\firstpage{1}
|
||||
|
||||
\title {Loss or Gain of Function? Effects of Ion Channel Mutations on Neuronal Firing Depend on the Neuron Type}
|
||||
|
||||
\author[\firstAuthorLast ]{\Authors} %This field will be automatically populated
|
||||
\address{} %This field will be automatically populated
|
||||
\correspondance{} %This field will be automatically populated
|
||||
\extraAuth{}
|
||||
|
||||
|
||||
\maketitle
|
||||
|
||||
|
||||
\begin{abstract}
|
||||
\section{}
|
||||
Clinically relevant mutations to voltage-gated ion channels, called channelopathies, alter ion channel function, properties of ionic currents and neuronal firing. The effects of ion channel mutations are routinely assessed and characterized as loss of function (LOF) or gain of function (GOF) at the level of ionic currents. However, emerging personalized medicine approaches based on LOF/GOF characterization have limited therapeutic success. Potential reasons are among others that the translation from this binary characterization to neuronal firing is currently not well understood --- especially when considering different neuronal cell types. Here we investigate the impact of neuronal cell type on the firing outcome of ion channel mutations with simulations of a diverse collection of conductance-based neuron models. We systematically analyzed the effects of changes in ion current properties on firing in different neuronal types. Additionally, we simulated the effects of known mutations in the \textit{KCNA1} gene encoding the \Kv potassium channel subtype associated with episodic ataxia type~1 (EA1). These simulations revealed that the outcome of a given change in ion channel properties on neuronal excitability depends on neuron type, i.e. the properties and expression levels of the unaffected ionic currents. Consequently, neuron-type specific effects are vital to a full understanding of the effects of channelopathies on neuronal excitability and are an important step towards improving the efficacy and precision of personalized medicine approaches.
|
||||
|
||||
|
||||
|
||||
|
||||
\tiny
|
||||
\keyFont{ \section{Keywords:} Channelopathy, Epilepsy, Ataxia, Potassium Current, Neuronal Simulations, Conductance-based Models, Neuronal heterogeneity }
|
||||
\end{abstract}
|
||||
\section{Introduction}
|
||||
The properties and combinations of voltage-gated ion channels are vital in determining neuronal excitability \citep{bernard_channelopathies_2008, carbone_ion_2020, rutecki_neuronal_1992, pospischil_minimal_2008}. However, ion channel function can be disturbed, for instance through genetic alterations, resulting in altered neuronal firing behavior \citep{carbone_ion_2020}. In recent years, next generation sequencing has led to an increase in the discovery of clinically relevant ion channel mutations and has provided the basis for pathophysiological studies of genetic epilepsies, pain disorders, dyskinesias, intellectual disabilities, myotonias, and periodic paralyses \citep{bernard_channelopathies_2008, carbone_ion_2020}.
|
||||
Ongoing efforts of many research groups have contributed to the current understanding of underlying disease mechanism in channelopathies. However, a complex pathophysiological landscape has emerged for many channelopathies and is likely a reason for limited therapeutic success with standard care.
|
||||
|
||||
Ion channel variants are frequently classified in heterologous expression systems as either a loss of function (LOF) or a gain of function (GOF) in the respective ionic current \citep{Musto2020, Kullmann2002, Waxman2011, Kim2021}. This LOF/GOF classification is often directly used to predict the effects on neuronal firing \citep{Niday2018, Wei2017, Wolff2017,Masnada2017}, which in turn is important for understanding the pathophysiology of these disorders and for identification of potential therapeutic targets \citep{Orsini2018, Yang2018, Colasante2020, Yu2006}. Experimentally, the effects of channelopathies on neuronal firing are assessed using primary neuronal cultures \citep{Scalmani2006, Smith2018, Liu2019} or \textit{in vitro} recordings from slices of transgenic mouse lines \citep{Mantegazza2019, Xie2010,Lory2020, Habib2015, Hedrich14874}, but are restricted to limited number of different neuron types. Neuron types differ in many aspects. They may differ in their composition of ionic currents \citep{yao2021taxonomy, Cadwell2016, BICCN2021, Scala2021} and therefore likely respond differently to changes in the properties of a single ionic current. Expression level of an affected gene \citep{Layer2021} and relative amplitudes of ionic currents \citep{rutecki_neuronal_1992, pospischil_minimal_2008,Kispersky2012, golowasch_failure_2002, barreiro_-current_2012} indeed dramatically influence the firing behavior and dynamics of neurons. Mutations in different sodium channel genes have been experimentally shown to affect firing in a neuron-type specific manner based on differences in expression levels of the affected gene \citep{Layer2021}, but also on other neuron-type specific mechanisms \citep{Hedrich14874, makinson_scn1a_2016}.
|
||||
|
||||
Neuron-type specificity is likely vital for successful precision medicine treatment approaches. For example, Dravet syndrome was identified as the consquence of LOF mutations in \textit{SCN1A} \citep{Claes2001,Fujiwara2003,Ohmori2002}, however limited success in the treatment of Dravet syndrome persisted \citep{Claes2001,Oguni2001} in part due to lack of understanding that inhibitory interneurons and not pyramidal neurons had altered excitability as a result of LOF \textit{SCN1A} mutations \citep{Yu2006, Colasante2020}.
|
||||
|
||||
Taken together, these examples demonstrate the need to study the effects of ion channel mutations in many different neuron types --- a daunting if not impossible experimental challenge. In the context of this diversity, simulations of conductance-based neuronal models are a powerful tool bridging the gap between altered ionic currents and firing in a systematic and efficient way. Furthermore, simlutions allow to predict the potential effects of drugs needed to alleviate the pathophysiology of the respective mutation \citep{johannesen_genotype-phenotype_2021, lauxmann_therapeutic_2021, Bayraktar}.
|
||||
|
||||
In this study, we therefore investigated how the outcome of ionic current kinetic changes on firing depend on neuronal cell type by (1) characterizing firing responses with two measures, (2) simulating the response of a repertoire of different neuronal models to changes in single current parameters as well as (3) to more complex changes in this case as they were observed for specific \textit{KCNA1} mutations that are associated with episodic ataxia type~1 \citep{Browne1994, Browne1995, lauxmann_therapeutic_2021}.
|
||||
|
||||
\section{Material and Methods}
|
||||
All modelling and simulation was done in parallel with custom written Python 3.8 (Python Programming Language; RRID:SCR\_008394) software, run on a Cent-OS 7 server with an Intel(R) Xeon (R) E5-2630 v2 CPU.
|
||||
|
||||
\subsection{Different Neuron Models}
|
||||
A group of neuronal models representing the major classes of cortical and thalamic neurons including regular spiking pyramidal (RS pyramidal; model D), regular spiking inhibitory (RS inhibitory; model B), and fast spiking (FS; model C) neurons were used \citep{pospischil_minimal_2008}. Additionally, a \Kv current (\IKv; \citealt{ranjan_kinetic_2019}) was added to each of these models (RS pyramidal +\Kv; model H, RS inhibitory +\Kv; model E, and FS +\Kv; model G respectively). A cerebellar stellate cell model from \citet{alexander_cerebellar_2019} is used (Cb stellate; model A) in this study. This neuron model was also extended by a \Kv current \citep{ranjan_kinetic_2019}, either in addition to the A-type potassium current (Cb stellate +\Kv; model F) or by replacing the A-type potassium current (Cb stellate \(\Delta\)\Kv; model J). A subthalamic nucleus (STN; model L) neuron model as described by \citet{otsuka_conductance-based_2004} was also used. The STN neuron model (model L) was additionally extended by a \Kv current \citep{ranjan_kinetic_2019}, either in addition to the A-type potassium current (STN +\Kv; model I) or by replacing the A-type potassium current (STN \(\Delta\)\Kv; model K). Model letter naming corresponds to panel lettering in Figure \ref{fig:diversity_in_firing}. The properties and maximal conductances of each model are detailed in Table \ref{tab:g} and the gating properties are unaltered from the original Cb stellate (model A) and STN (model L) models \citep{alexander_cerebellar_2019, otsuka_conductance-based_2004}. For enabling the comparison of models with the typically reported electrophysiological data fitting reported and for ease of further gating curve manipulations, a modified Boltzmann function
|
||||
|
||||
\begin{equation}\label{eqn:Boltz}
|
||||
x_\infty = {\left(\frac{1-a}{1+{\exp\left[{\frac{V-V_{1/2}}{k}}\right]}} +a\right)^j}
|
||||
\end{equation}
|
||||
|
||||
with slope \(k\), voltage for half-maximal activation or inactivation (\(V_{1/2}\)), exponent \(j\), and persistent current \(0 \leq a \leq 1\) were fitted to the original formulism for RS pyramidal (model D), RS inhibitory (model B) and FS (model C) models from \citet{pospischil_minimal_2008}. The properties of \IKv were fitted to the mean wild type biophysical parameters of \Kv described in \citet{lauxmann_therapeutic_2021}. The fitted gating parameters are detailed in Table \ref{tab:gating}. Each of the original single-compartment models used here can reproduce physiological firing behavior of the neurons they represent (Figure \ref{fig:diversity_in_firing}; \citealt{pospischil_minimal_2008, alexander_cerebellar_2019, otsuka_conductance-based_2004}) and capture key aspects of the dynamics of these neuron types.
|
||||
|
||||
\subsection{Firing Frequency Analysis}
|
||||
The membrane responses to 200 equidistant two second long current steps were simulated using the forward-Euler method with a \(\Delta \textrm{t} = 0.01\)\,ms from steady state. Current steps ranged from 0 to 1\,nA (step size 5\,pA) for all models except for the RS inhibitory neuron models where a range of 0 to 0.35\,nA (step size 1.75\,pA) was used to ensure repetitive firing across the range of input currents. For each current step, action potentials were detected as peaks with at least 50\,mV prominence, or the relative height above the lowest contour line encircling it, and a minimum interspike interval of 1\,ms. The interspike interval was computed and used to determine the instantaneous firing frequencies elicited by the current step.
|
||||
|
||||
To ensure accurate firing frequencies at low firing rates and reduced spike sampling bias, steady-state firing was defined as the mean firing frequency in a 500\,ms window in the last second of the current steps starting at the inital action potential in this last second.
|
||||
Firing characterization was performed in the last second of current steps to ensure steady-state firing is captured and adaptation processes are neglected in our analysis. Alteration in current magnitudes can have different effects on rheobase and the initial slope of the fI curve \citep{Kispersky2012}.
|
||||
For this reason, we quantified neuronal firing using the rheobase as well as the area under the curve (AUC) of the initial portion of the fI curve as a measure of the initial slope of the fI curve Figure \ref{fig:firing_characterization}~A.
|
||||
|
||||
The smallest current at which steady state firing occured was identified and the current step interval preceding the occurrence of steady state firing was simulated at higher resolution (100 current steps) to determine the current at which steady state firing began. Firing was simulated with 100 current steps from this current upwards for 1/5 of the overall current range. Over this range a fI curve was constructed and the integral, or area under the curve (AUC), of the fI curve over this interval was computed with the composite trapezoidal rule and used as a measure of firing rate independent from rheobase.
|
||||
|
||||
To obtain the rheobase at a higher current resolution than the fI curve, the current step interval preceding the occurrence of action potentials was explored at higher resolution with 100 current steps spanning the interval (step sizes of 0.05 pA and 0.0175 pA, respectively). Membrane responses to these current steps were then analyzed for action potentials and the rheobase was considered the lowest current step for which an action potential was elicited.
|
||||
|
||||
All models exhibited tonic steady-state firing with default parameters. In limited instances, variations of parameters elicited periodic bursting, however these instances were excluded from further analysis.
|
||||
|
||||
|
||||
\subsection{Sensitivity Analysis and Comparison of Models}
|
||||
Properties of ionic currents common to all models (\(\text{I}_{\text{Na}}\), \(\text{I}_{\text{Kd}}\), \(\text{I}_{\text{A}}\)/\IKv, and \(\text{I}_{\text{Leak}}\)) were systematically altered in a one-factor-at-a-time sensitivity analysis for all models. The gating curves for each current were shifted (\(\Delta V_{1/2}\)) from -10 to 10\,mV in increments of 1\,mV. The voltage dependence of the time constant associated with the shifted gating curve was correspondingly shifted. The slope (\(k\)) of the gating curves were altered from half to twice the initial slope. Similarly, the maximal current conductance (\(g\)) was also scaled from half to twice the initial value. For both slope and conductance alterations, alterations consisted of 21 steps spaced equally on a \(\textrm{log}_2\) scale. We neglected the variation of time constants for the practical reason that estimation and assessment of time constants and changes to them is not straightforward \citep{Clerx2019, Whittaker2020}.
|
||||
|
||||
\subsection{Model Comparison}
|
||||
Changes in rheobase (\drheo) were calculated in relation to the original model rheobase. The contrast of each AUC value (\(AUC_i\)) was computed in comparison to the AUC of the unaltered wild type model (\(AUC_{wt}\))
|
||||
\begin{equation}\label{eqn:AUC_contrast}
|
||||
\textrm{normalized } \Delta \text{AUC} = \frac{AUC_i - AUC_{wt}}{AUC_{wt}}
|
||||
\end{equation}
|
||||
|
||||
To assess whether the effects of a given alteration on \ndAUC or \drheo were robust across models, the correlation between \ndAUC or \drheo and the magnitude of the alteration of a current property was computed for each alteration in each model and compared across alteration types.
|
||||
The Kendall's \(\tau\) coefficient, a non-parametric rank correlation, is used to describe the relationship between the magnitude of the alteration and AUC or rheobase values. A Kendall \(\tau\) value of -1 or 1 is indicative of monotonically decreasing and increasing relationships respectively.
|
||||
|
||||
\subsection{\textit{KCNA1} Mutations}\label{subsec:mut}
|
||||
Known episodic ataxia type~1 associated \textit{KCNA1} mutations and their electrophysiological characterization have been reviewed in \citet{lauxmann_therapeutic_2021}. The mutation-induced changes in \IKv amplitude and activation slope (\(k\)) were normalized to wild type measurements and changes in activation \(V_{1/2}\) were used relative to wild type measurements. Although initially described to lack fast activation, \Kv displays prominent inactivation at physiologically relevant temperatures \citep{ranjan_kinetic_2019}. The effects of a mutation were also applied to \(\text{I}_{\text{A}}\) when present as both potassium currents display inactivation. In all cases, the mutation effects were applied to half of the \Kv or \(\text{I}_{\text{A}}\) under the assumption that the heterozygous mutation results in 50\% of channels carrying the mutation. Frequency-current curves for each mutation in each model were obtained through simulation and used to characterize firing behavior as described above. For each model the differences in mutation AUC to wild type AUC were normalized by wild type AUC (\ndAUC) and mutation rheobases were compared to wild type rheobase values (\drheo). Pairwise Kendall rank correlations (Kendall \(\tau\)) were used to compare the correlation in the effects of \Kv mutations on AUC and rheobase between models.
|
||||
|
||||
\subsection{Code Accessibility}
|
||||
The simulation and analysis code including full specification of the models is freely available online at \href{https://github.com/nkoch1/LOFGOF2023}{https://github.com/nkoch1/LOFGOF2023}. %The code is available as \ref{code_zip}.
|
||||
|
||||
\section{Results}
|
||||
|
||||
To examine the role of neuron-type specific ionic current environments on the impact of altered ion currents properties on firing behavior:
|
||||
(1) firing responses were characterized with rheobase and \(\Delta\)AUC, (2) a set of neuronal models was used and properties of channels common across models were altered systematically one at a time, and (3) the effects of a set of episodic ataxia type~1 associated \textit{KCNA1} mutations on firing was then examined across different neuronal models with different ionic current environments.
|
||||
|
||||
\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 (Figure \ref{fig:diversity_in_firing}). The models chosen for this study all fire tonically and do not exhibit bursting (see methods for details and naming of the models). Models are qualitatively sorted based on their firing curves and labeled model A through L accordingly. Some models, such as models A and B, display type I firing, whereas others such as models J and L exhibit 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; \citealt{ermentrout_type_1996, Rinzel_1998}). 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 models I, J, and K have large hysteresis (Figure \ref{fig:diversity_in_firing} and Supplementary Figure S1). Different types of underlying current dynamics are known to generate these different firing types and hysteresis \cite{ERMENTROUT2002, ermentrout_type_1996, Izhikevich2006}. This broad range of single-compartmental models represents the distinct dynamics of various neuron types across diverse brain regions.
|
||||
|
||||
\subsection{Characterization of Neuronal Firing Properties}
|
||||
Neuronal firing is a complex phenomenon, and a quantification of firing properties is required for comparisons across neuron types and between different conditions. Here we focus on two aspects of firing that are routinely measured in clinical settings \citep{Bryson_2020}: rheobase, the smallest injected current at which the neuron 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 (Figure \ref{fig:firing_characterization}~A). The characterization of the firing properties of a neuron by using 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) by two independent measures. Note that AUC is essentially quantifying the slope of a neuron's fI curve.
|
||||
|
||||
Using these two measures we quantified 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 of 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 (Figure \ref{fig:firing_characterization}~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 gain of function (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, loss of function (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.
|
||||
|
||||
|
||||
\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 delayed rectifier potassium current in the model G to more depolarized values, then the rheobase of the resulting fI curves shifted to lower currents \(-\)\drheo, making the neuron more sensitive to weak inputs, but at the same time the slope of the fI curves was reduced (\(-\)\ndAUC), which resulted in a reduced firing rate (Figure \ref{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 model G is in the bottom left quadrant of Figure \ref{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 Figure \ref{fig:AUC_correlation}~B). For each of the many models we got a different relation between the changes in AUC and the shifts in half maximal potential \(V_{1/2}\) (thin lines in Figure \ref{fig:AUC_correlation}~B). To further summarize these different dependencies of the various models we characterized each of these curves by a single number, the \( \text{Kendall} \ \tau \) correlation coefficient. A monotonically increasing curve resulted in a \( \text{Kendall} \ \tau \) close to \(+1\) a monotonously decreasing curve in \( \text{Kendall} \ \tau \approx -1 \), and a non-monotonous, non-linear relation in \( \text{Kendall} \ \tau \) close to zero (compare lines in Figure \ref{fig:AUC_correlation}~B with dots in black box in panel C).
|
||||
|
||||
Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affected the AUC (Figure \ref{fig:AUC_correlation}), but how exactly the AUC was affected usually depended on the specific neuronal model. Increasing the slope factor of the \Kv activation curve for example increased the AUC in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (Figure \ref{fig:AUC_correlation}~D,E,F). Similar consistent positive correlations could 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 resulted in negative correlations with the AUC in all models (\( \text{Kendall} \ \tau \approx -1\)).
|
||||
|
||||
Qualitative differences could be found, for example, when increasing the maximal conductance of the delayed rectifier (Figure \ref{fig:AUC_correlation}~G,H,I). In some model neurons this increased AUC (\( \text{Kendall} \ \tau \approx +1\)), whereas in others AUC was decreased (\( \text{Kendall} \ \tau \approx -1\)). In model I, AUC depended in a non-linear way on the maximal conductance of the delayed rectifier, resulting in a \( \text{Kendall} \ \tau \) close to zero. Even more dramatic qualitative differences between models resulted from shifts of the activation curve of the delayed rectifier, as discussed already above (Figure \ref{fig:AUC_correlation}~A,B,C). Some model neurons did almost not depend on changes in Kd-current half activation \(V_{1/2}\) or showed strong non-linear dependencies, both resulting in \( \text{Kendall} \ \tau\) close to zero. Many model neurons showed strongly negative correlations, and a few displayed positive correlations with shifting the activation curve of the delayed rectifier.
|
||||
|
||||
Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affected rheobase (Figure \ref{fig:rheobase_correlation}). However, in contrast to AUC, qualitatively consistent effects on rheobase across models could be observed. An increasing of the maximal conductance of the leak current in the model A increased the rheobase (Figure \ref{fig:rheobase_correlation}~G). When these changes were plotted against the change in maximal conductance a monotonically increasing relationship was evident (thick teal line in Figure \ref{fig:rheobase_correlation}~H). This monotonically increasing relationship was evident in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (thin lines in Figure \ref{fig:rheobase_correlation}~H). Similarly, positive correlations were consistently found across models for maximal conductances of delayed rectifier Kd, \Kv, and A type currents, whereas the maximal conductance of the sodium current was consistently associated with negative correlations (\( \text{Kendall} \ \tau \approx -1\); Figure \ref{fig:rheobase_correlation}~I), i.e. rheobase decreased with increasing maximum conductance in all models.
|
||||
|
||||
Although changes in half maximal potential \(V_{1/2}\) and slope factor \(k\) generally correlated with rheobase similarly across models there were some exceptions. Rheobase was affected with both with positive and negative correlations in different models as a result of changing slope factor of Na\(^+\)-current inactivation (positive: models A--H and J; negative: models I, K and L), \Kv-current inactivation (positive: models I and K; negative: models E--G, J, H), and A-current activation (positive: models A, F and L; negative: model I; Figure \ref{fig:rheobase_correlation}~F). Departures from monotonic relationships also occurred in some models as a result of Kd-current activation \(V_{1/2}\) (e.g. model J) and slope factor \(k\) (models F and G), \Kv-current inactivation slope factor \(k\) (model K), and A-current activation slope factor \(k\) (model L). 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.
|
||||
|
||||
\subsection{\textit{KCNA1} Mutations}
|
||||
Mutations in \textit{KCNA1} are associated with episodic ataxia type~1 (EA1) and have been characterized biophysically (as reviewed by \citet{lauxmann_therapeutic_2021}). Here they were used as a test case in the effects of various ionic current environments on neuronal firing and on the outcomes of channelopathies. The changes in AUC and rheobase from wild type values for reported EA1 associated \textit{KCNA1} mutations were heterogeneous across models containing \Kv, but generally showed decreases in rheobase (Figure \ref{fig:simulation_model_comparision}A--I). Pairwise non-parametric Kendall \(\tau\) rank correlations between the simulated effects of these \Kv mutations on rheobase were highly correlated across models (Figure \ref{fig:simulation_model_comparision}J) indicating that EA1 associated \textit{KCNA1} mutations generally decrease rheobase across diverse neuron-types. However, the effects of the \Kv mutations on AUC were more heterogenous as reflected by both weak and strong positive and negative pairwise correlations between models (Figure \ref{fig:simulation_model_comparision}K), suggesting that the effects of ion-channel variant on super-threshold neuronal firing depend both quantitatively and qualitatively on the specific composition of ionic currents in a given neuron.
|
||||
|
||||
|
||||
\section{Discussion}
|
||||
|
||||
To compare the effects of ion channel mutations on neuronal firing of different neuron types, we used a diverse set of conductance-based models to systematically characterize the effects of changes in individual channel properties. Additionally, we simulated the effects of specific episodic ataxia type~1 associated (EA1) \textit{KCNA1} mutations. Changes to single ionic current properties, as well as known EA1 associated \textit{KCNA1} mutations showed consistent effects on the rheobase across neuron types, whereas the effects on the slope of the steady-state fI-curve depended on the neuron type. Our results demonstrate that loss of function (LOF) and gain of function (GOF) on the biophysical level cannot be uniquely transferred to the level of neuronal firing. Thus, the effects caused by different mutations depend on the properties of the other ion channels expressed in a neuron and are therefore depend on the channel ensemble of a specific neuron type.
|
||||
|
||||
\subsection{Firing Frequency Analysis}
|
||||
Although differences in neuronal firing can be characterized by an area under the curve of the fI curve for a fixed current range, this approach characterizes firing as a mixture of key features: rheobase and the initial slope of the fI curve. By probing rheobase directly and using an AUC relative to rheobase, we disambiguate these features and enable insights into the effects on rheobase and initial fI curve steepness. This increases the specificity of our understanding of how ion channel mutations alter firing across neuron types and enable classification as described in Figure \ref{fig:firing_characterization}. Importantly, in cases when ion channel mutations alter rheobase and initial fI curve steepness in ways that oppose effects on firing (upper left and bottom right quadrants of Figure \ref{fig:firing_characterization}~B), this disamgibuation is important for understanding the outcome of the mutation. In these cases, the regime the neuron is operating in is vital in determining the neuron's firing outcome. If it is in its excitable regime and only occasionally generates an action potential, then the effect on the rheobase is more important. If it is firing periodically with high rates, then the change in AUC might be more relevant.
|
||||
|
||||
\subsection{Modelling Limitations}
|
||||
The models used here are simple and they all capture key aspects of the firing dynamics for their respective neuron. The simple models fall short of capturing the complex physiology, biophysics and heterogeneity of real neurons, nor do they take into account subunit stoichiometry, auxillary subunits, membrane composition which influence the biophysics of ionic currents \citep{Al-Sabi_2013, Oliver_2004, Pongs_2009, Rettig1994}. However, for the purpose of understanding how different neuron-types, or current environments, contribute to the diversity in firing outcomes of ion channel mutations, the fidelity of the models to the physiological neurons they represent is of a minor concern. For exploring possible neuron-type specific effects, variety in currents and dynamics across models is of utmost importance. With this context in mind, the collection of models used here are labelled as models A-L to highlight that the physiological neurons they represent is not of chief concern, but rather that the collection of models with different attributes respond heterogeneously to the same perturbation. Additionally, the development of more realistic models is a high priority and will enable neuron-type specific predictions that may aid precision medicine approaches. Thus, weight should not be put on any single predicted firing outcome here in a specific model, but rather on the differences in outcomes that occur across the neuron-type spectrum the models used here represent. % Menchaca_2012,
|
||||
|
||||
\subsection{Neuronal Diversity}
|
||||
The nervous system consists of a vastly diverse and heterogenous collection of neurons with variable properties and characteristics including diverse combinations and expression levels of ion channels which are vital for neuronal firing dynamics.
|
||||
|
||||
Advances in high-throughput techniques have enabled large-scale investigation into single-cell properties across the CNS \citep{Poulin2016} that have revealed large diversity in neuronal gene expression, morphology and neuronal types in the motor cortex \citep{Scala2021}, neocortex \cite{Cadwell2016, Cadwell2020}, GABAergic neurons in the cortex and retina \citep{Huang2019, Laturnus2020}, cerebellum \citep{Kozareva2021}, spinal cord \citep{Alkaslasi2021}, visual cortex \citep{Gouwens2019} as well as the retina \citep{Baden2016, Voigt2019, Berens2017, Yan2020a, Yan2020b}.
|
||||
|
||||
Diversity across neurons is not limited to gene expression and can also be seen electrophysiologically \citep{Tripathy2017, Gouwens2018, Tripathy2015, Scala2021, Cadwell2020, Gouwens2019, Baden2016, Berens2017} with correlations existing between gene expression and electrophysiological properties \citep{Tripathy2017}. At the ion channel level, diversity exists not only between the specific ion channels the different neuron types express but heterogeneity also exists in ion channel expression levels within neuron types \citep{marder_multiple_2011, goaillard_ion_2021,barreiro_-current_2012}. As ion channel properties and expression levels are key determinents of neuronal dynamics and firing \citep{Balachandar2018, Gu2014, Zeberg2015, Aarhem2007, Qi2013, Gu2014a, Zeberg2010, Zhou2020, Kispersky2012} neurons with different ion channel properties and expression levels display different firing properties.
|
||||
|
||||
To capture the diversity in neuronal ion channel expression and its relevance in the outcome of ion channel mutations, we used multiple neuronal models with different ionic currents and underlying firing dynamics here.
|
||||
|
||||
\subsection{Ionic Current Environments Determine the Effect of Ion Channel Mutations}
|
||||
|
||||
To our knowledge, no comprehensive evaluation of how ionic current environment and neuron type affect the outcome of ion channel mutations have been reported. However, comparisons between the effects of such mutations between certain neuron types were described. For instance, the R1648H mutation in \textit{SCN1A} does not alter the excitability of cortical pyramidal neurons, but causes hypoexcitability of adjacent inhibitory GABAergic neurons \citep{Hedrich14874}. In the CA3 region of the hippocampus, the equivalent mutation in \textit{SCN8A}, R1627H, increases the excitability of pyramidal neurons and decreases the excitability of parvalbumin positive interneurons \cite{makinson_scn1a_2016}. Additionally, the L858H mutation in \(\text{Na}_\text{V}\text{1.7}\), associated with erythermyalgia, has been shown to cause hypoexcitability in sympathetic ganglion neurons and hyperexcitability in dorsal root ganglion neurons \citep{Waxman2007, Rush2006}. The differential effects of L858H \(\text{Na}_\text{V}\text{1.7}\) on firing is dependent on the presence or absence of another sodium channel, namely the \(\text{Na}_\text{V}\text{1.8}\) subunit \citep{Waxman2007, Rush2006}. These findings, in concert with our findings emphasize that the ionic current environment in which a channelopathy occurs is vital in determining the outcomes of the channelopathy on firing.
|
||||
|
||||
Neuron type specific differences in ionic current properties are important in the effects of ion channel mutations. However, within a neuron type heterogeneity in channel expression levels exists and it is often desirable to generate a population of neuronal models and to screen them for plausibility to biological data in order to capture neuronal population diversity \citep{marder_multiple_2011,OLeary2016}. The models we used here are originally generated by characterization of current gating properties and by fitting of maximal conductances to experimental data \citep{pospischil_minimal_2008, ranjan_kinetic_2019, alexander_cerebellar_2019, otsuka_conductance-based_2004}. This practice of fixing maximal conductances based on experimental data is limiting as it does not reproduce the variability in channel expression and neuronal firing behavior of a heterogeneous neuron population \citep{verma_computational_2020}. For example, a model derived from the mean conductances in a neuronal sub-population within the stomatogastric ganglion, the so-called "one-spike bursting" neurons fire three spikes instead of one per burst due to an L-shaped distribution of sodium and potassium conductances \citep{golowasch_failure_2002}. Multiple sets of conductances can give rise to the same patterns of activity also termed degeneracy and differences in neuronal dynamics may only be evident with perturbations \citep{marder_multiple_2011, goaillard_ion_2021}. The variability in ion channel expression often correlates with the expression of other ion channels \citep{goaillard_ion_2021} and neurons whose behavior is similar may possess correlated variability across different ion channels resulting in stability in the neuronal phenotype \citep{lamb_correlated_2013, soofi_co-variation_2012, taylor_how_2009}. The variability of ionic currents and degeneracy of neurons may account, at least in part, for the observation that the effect of toxins within a neuronal type is frequently not constant \citep{khaliq_relative_2006, puopolo_roles_2007, ransdell_neurons_2013}.
|
||||
|
||||
\subsection{Effects of \textit{KCNA1} Mutations}
|
||||
Changes in delayed rectifier potassium currents, analogous to those seen in LOF \textit{KCNA1} mutations, change the underlying firing dynamics of the Hodgkin Huxley model result in reduced thresholds for repetitive firing and thus contribute to increased excitability \citep{hafez_altered_2020}. Although the Hodgkin Huxley delayed rectifier lacks inactivation, the increases in excitability observed by \citet{hafez_altered_2020} are in line with our simulation-based predictions of the outcomes of \textit{KCNA1} mutations. LOF \textit{KCNA1} mutations generally increase neuronal excitability, however the varying susceptibility on rheobase and different effects on AUC of the fI-curve of \textit{KCNA1} mutations across models are indicative that a certain neuron type specific complexity exists. Increased excitability is seen experimentally with \Kv null mice \citep{smart_deletion_1998, zhou_temperature-sensitive_1998}, with pharmacological \Kv block \citep{chi_manipulation_2007, morales-villagran_protection_1996} and by \citet{hafez_altered_2020} with simulation-based predictions of \textit{KCNA1} mutations. Contrary to these results, \citet{zhao_common_2020} predicted \textit{in silico} that the depolarizing shifts seen as a result of \textit{KCNA1} mutations broaden action potentials and interfere negatively with high frequency action potential firing. However, they varied stimulus duration between different models and therefore comparability of firing rates is lacking in this study.
|
||||
|
||||
In our simulations, different current properties alter the impact of \textit{KCNA1} mutations on firing as evident in the differences seen in the impact of \(\text{I}_\text{A}\) and \IKv in the Cb stellate and STN model families on \textit{KCNA1} mutation firing. This highlights that not only knowledge of the biophysical properties of a channel but also its neuronal expression and other neuronal channels present is vital for the holistic understanding of the effects of a given ion channel mutation both at the single cell and network level.
|
||||
|
||||
\subsection{Loss or Gain of Function Characterizations Do Not Fully Capture Ion Channel Mutation Effects on Firing}
|
||||
The effects of changes in channel properties depend in part on the neuronal model in which they occur and can be seen in the variance of correlations (especially in AUC of the fI-curve) across models for a given current property change. Therefore, relative conductances and gating properties of currents in the ionic current environment in which an alteration in current properties occurs plays an important role in determining the outcome on firing. The use of LOF and GOF is useful at the level of ion channels to indicate whether a mutation results in more or less ionic current. However, the extension of this thinking onto whether mutations induce LOF or GOF at the level of neuronal firing based on the ionic current LOF/GOF is problematic because of this dependency of neuronal firing changes on the ionic channel environment. Thus, the direct leap from current level LOF/GOF characterizations to effects on firing without experimental or modelling-based evidence, although tempting, should be refrained from and viewed with caution when reported. This is especially relevant in the recent development of personalized medicine for channelopathies, where a patient's specific channelopathy is identified and used to tailor treatments \citep{Weber2017, Ackerman2013, Helbig2020, Gnecchi2021, Musto2020, Brunklaus2022, Hedrich2021}. In these cases the effects of specific ion channel mutations are often characterized based on ionic currents in expression systems and classified as LOF or GOF to aid in treatment decisions \citep{johannesen_genotype-phenotype_2021, Brunklaus2022, Musto2020}. Although positive treatment outcomes occur with sodium channel blockers in patients with GOF \(\text{Na}_{\text{V}}\text{1.6}\) mutations, patients with both LOF and GOF \(\text{Na}_{\text{V}}\text{1.6}\) mutations can benefit from treatment with sodium channel blockers \citep{johannesen_genotype-phenotype_2021}. This example suggests that the relationship between effects at the level of ion channels and effects at the level of firing and therapeutics is not linear or evident without further contextual information.
|
||||
|
||||
Therefore, the transfer of LOF or GOF from the current to the firing level should be used with caution; the neuron type in which the mutant ion channel is expressed may provide valuable insight into the functional consequences of an ion channel mutation. Experimental assessment of the effects of a patient's specific ion channel mutation \textit{in vivo} is not generally feasible at a large scale. Therefore, modelling approaches investigating the effects of patient specific channelopathies provide a viable method bridging between characterization of changes in biophysical properties of ionic currents and the firing consequences of these effects. In both experimental and modelling studies on the effects of ion channel mutations on neuronal firing the specific dependency on neuron type should be considered.
|
||||
|
||||
The effects of altered ion channel properties on firing is generally influenced by the other ionic currents present in the neuron. In channelopathies the effect of a given ion channel mutation on neuronal firing therefore depends on the neuron type in which those changes occur \citep{Hedrich14874, makinson_scn1a_2016, Waxman2007, Rush2006}. Although certain complexities of neurons such as differences in neuron-type sensitivities to current property changes, interactions between ionic currents, cell morphology and subcellular ion channel distribution are neglected here, it is likely that this increased complexity \textit{in vivo} would contribute to the neuron-type dependent effects on neuronal firing. The complexity and nuances of the nervous system, including neuron-type dependent firing effects of channelopathies explored here, likely underlie shortcomings in treatment approaches in patients with channelopathies. Accounting for neuron-type dependent firing effects provides an opportunity to improve the efficacy and precision in personalized medicine approaches.
|
||||
|
||||
With this study we suggest that neuron-type specific effects are vital to a full understanding of the effects of channelopathies at the level of neuronal firing. Furthermore, we highlight the use of modelling approaches to enable relatively fast and efficient insight into channelopathies.
|
||||
|
||||
|
||||
\bibliographystyle{Frontiers-Harvard}
|
||||
\bibliography{Koch_ref}
|
||||
\newpage
|
||||
\section*{Figures}
|
||||
%%% NB logo1.eps is required in the path in order to correctly compile front page header %%%
|
||||
|
||||
|
||||
\begin{figure}[H]
|
||||
\centering
|
||||
\includegraphics[width=0.875\linewidth]{diversity_in_firing.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Diversity in Neuronal Model Firing. Spike trains (left), frequency-current (fI) curves (right) for Cb stellate \textbf{(A)}, RS inhibitory \textbf{(B)}, FS \textbf{(C)}, RS pyramidal \textbf{(D)}, RS inhibitory +\Kv \textbf{(E)}, Cb stellate +\Kv \textbf{(F)}, FS +\Kv \textbf{(G)}, RS pyramidal +\Kv \textbf{(H)}, STN +\Kv \textbf{(I)}, Cb stellate \(\Delta\)\Kv \textbf{(J)}, STN \(\Delta\)\Kv \textbf{(K)}, and STN \textbf{(L)} neuron models. Models are sorted qualitatively based on their fI curves. Black markers on the fI curves indicate the current step at which the spike train occurs. The green marker indicates the current at which firing begins in response to an ascending current ramp, whereas the red marker indicates the current at which firing ceases in response to a descending current ramp (see Supplementary Figure S1).} % \textcolor{red}{(see Supplementary Figure \ref{ramp_firing}).}}
|
||||
\label{fig:diversity_in_firing}
|
||||
\end{figure}
|
||||
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=0.5\linewidth]{firing_characterization_arrows.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Characterization of firing with AUC and rheobase. \textbf{(A)} The area under the curve (AUC) of the repetitive firing frequency-current (fI) curve. \textbf{(B)}
|
||||
Changes in firing as characterized by \(\Delta\)AUC and \(\Delta\)rheobase occupy four quadrants separated by no changes in AUC and rheobase. Representative schematic fI curves in red with respect to a reference (or wild type) fI curve (blue) depict the general changes associated with each quadrant.}
|
||||
\label{fig:firing_characterization}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{AUC_correlation.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of altered channel kinetics on AUC in various neuron models. The fI curves corresponding to shifts in model G delayed rectifier Kd half activation \(V_{1/2}\) \textbf{(A)}, changes \Kv activation slope factor \(k\) in model G \textbf{(D)}, and changes in maximal conductance of delayed rectifier K current in the model I \textbf{(G)} are shown. The fI curves from the smallest (grey) to the largest (black) alterations are seen for \textbf{(A,D,} and \textbf{G}) in accordance to the greyscale of the x-axis in B, E, and H. The \ndAUC of fI curves is plotted against delayed rectifier K half activation potential (\(\Delta V_{1/2}\); \textbf{B}), \Kv activation slope factor \(k\) (k/\(\text{k}_{\text{WT}}\); \textbf{E}) and maximal conductance \(g\) of the delayed rectifier Kd current (g/\(\text{g}_{\text{WT}}\); \textbf{H}) for all models (thin lines) with relationships from the fI curve examples ( \textbf{A, D, G} respectively) highlighted by thick lines with colors corresponding to the box highlighting each set of fI curves. The Kendall rank correlation (Kendall \(\tau\)) coefficients between shifts in half maximal potential \(V_{1/2}\) and \ndAUC \textbf{(C)}, slope factor k and \ndAUC \textbf{(F)} as well as maximal current conductances and \ndAUC \textbf{(I)} for each model and current property is computed. The relationships between \(\Delta V_{1/2}\), k/\(\text{k}_{\text{WT}}\), and g/\(\text{g}_{\text{WT}}\) and \ndAUC for the Kendall rank correlations highlighted in the black boxes are depicted in \textbf{(B)}, \textbf{(E)} and \textbf{(H)} respectively. }
|
||||
\label{fig:AUC_correlation}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{rheobase_correlation.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of altered channel kinetics on rheobase. The fI curves corresponding to shifts in model G \Kv activation \(V_{1/2}\) \textbf{(A)}, changes \Kv inactivation slope factor \(k\) in model F \textbf{(D)}, and changes in maximal conductance of the leak current in model A \textbf{(G)} are shown. The fI curves from the smallest (grey) to the largest (black) alterations are seen for ( \textbf{A,D,} and \textbf{G}) in accordance to the greyscale of the x-axis in \textbf{B}, \textbf{E}, and \textbf{H}. The \drheo of fI curves is plotted against \Kv half activation potential (\(\Delta V_{1/2}\); \textbf{B}), \Kv inactivation slope factor \(k\) (k/\(\text{k}_{\text{WT}}\); \textbf{E}) and maximal conductance \(g\) of the leak current (g/\(\text{g}_{\text{WT}}\); \textbf{H}) for all models (thin lines) with relationships from the fI curve examples ( \textbf{A, D, G} respectively) highlighted by thick lines with colors corresponding to the box highlighting each set of fI curves. The Kendall rank correlation (Kendall \(\tau\)) coefficients between shifts in half maximal potential \(V_{1/2}\) and \drheo \textbf{(C)}, slope factor k and \drheo \textbf{(F)} as well as maximal current conductances and \drheo \textbf{(I)} for each model and current property is computed. The relationships between \(\Delta V_{1/2}\), k/\(\text{k}_{\text{WT}}\), and g/\(\text{g}_{\text{WT}}\) and \drheo for the Kendall rank correlations highlighted in the black boxes are depicted in \textbf{(B)}, \textbf{(E)} and \textbf{(H)} respectively.}
|
||||
\label{fig:rheobase_correlation}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{simulation_model_comparison.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of episodic ataxia type~1 associated \textit{KCNA1} mutations on firing. Effects of \textit{KCNA1} mutations on AUC (percent change in normalized \(\Delta\)AUC) and rheobase (\(\Delta\)Rheobase) compared to wild type for model H \textbf{(A)}, model E \textbf{(B)}, model G \textbf{(C)}, model A \textbf{(D)}, model F \textbf{(E)}, model J \textbf{(F)}, model L \textbf{(G)}, model I \textbf{(H)} and model K \textbf{(I)}. All \textit{KCNA1} Mutations are marked in grey with the V174F, F414C, E283K, and V404I \textit{KCNA1} mutations highlighted in color for each model. Pairwise Kendall rank correlation coefficients (Kendall \(\tau\)) between the effects of \textit{KCNA1} mutations on rheobase and on AUC are shown in J and K respectively. Marker shape is indicative of model/firing type, and grey dashed lines denote the quadrants of firing characterization (see Figure \ref{fig:firing_characterization}).}
|
||||
\label{fig:simulation_model_comparision}
|
||||
\end{figure}
|
||||
|
||||
|
||||
|
||||
\FloatBarrier
|
||||
\section*{Tables}
|
||||
|
||||
%\input{g_table}
|
||||
\begin{table}[ht!]
|
||||
% \fontsize{10pt}{10pt}
|
||||
\caption[Cell properties and conductances of neuronal models]{Cell properties and conductances of regular spiking pyramidal neuron (RS Pyramidal; model D), regular spiking inhibitory neuron (RS Inhibitory; model B), fast spiking neuron (FS; model C) each with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (RS Pyramidal +\Kv; model H, RS Inhibitory +\Kv; model E, FS +\Kv; model G respectively), cerebellar stellate cell (Cb Stellate; model A), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (Cb Stellate +\Kv; model F) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_\textrm{A}\) (Cb Stellate \(\Delta\)\Kv; model J), and subthalamic nucleus neuron (STN; model L), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (STN +\Kv; model I) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_{\textrm{A}}\) (STN \Kv; model K) models. All conductances are given in \(\textrm{mS}/\textrm{cm}^2\). Capacitances (\(C_m\)) and \(\tau_{max, M}\) are given in pF and ms respectively.}
|
||||
\centering
|
||||
% \begin{tabular}[x]{@{}l@{}} \\ \end{tabular}
|
||||
\linespread{1.5}\selectfont
|
||||
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\resizebox{\textwidth}{!}{\begin{tabular}{cccccccccc}
|
||||
% in mS/cm^2
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& \begin{tabular}[x]{@{}c@{}} RS\\Pyra-\\midal\\(+\Kv) \end{tabular} & \begin{tabular}[x]{@{}c@{}} RS\\Inhib-\\itory\\(+\Kv)\end{tabular} & \begin{tabular}[x]{@{}c@{}}FS\\(+\Kv) \end{tabular}& \begin{tabular}[x]{@{}c@{}} Cb\\Stellate \end{tabular}& \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\+\Kv \end{tabular} & \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\\(\Delta\)\Kv \end{tabular} & STN &\begin{tabular}[x]{@{}c@{}} STN\\+\Kv \end{tabular} &\begin{tabular}[x]{@{}c@{}} STN\\\(\Delta\)\Kv \end{tabular} \\
|
||||
Model & D (H) & B (E) & C (G) & A & F & J & L & I & K \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(g_{Na}\) & \(56\) & \(10\) & \(58\) & \(3.4\) & \(3.4\) & \(3.4\) & \(49\) & \(49\) & \(49\) \\
|
||||
\(g_{Kd}\) & \(6\) (\(5.4\)) & 2.1 (\(1.89\)) & 3.9 (\(3.51\)) & \(9.0556\) & \(8.15\) &\(9.0556\) & \(57\) & \(56.43\) & \(57\) \\
|
||||
\(g_{K_V1.1}\) & --- (\(0.6\)) & --- (\(0.21\)) & --- (\(0.39\)) & --- & \(0.90556\) & \(1.50159\) & --- & \(0.57\) & \(0.5\) \\
|
||||
\(g_{A}\) & --- & --- & --- & \(15.0159\) & \(15.0159\) & --- & \(5\) & \(5\) & --- \\
|
||||
\(g_{M}\) & \(0.075\) & \(0.0098\) &\(0.075\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(g_{L}\) & --- & --- & --- & --- & --- & --- & \(5\) & \(5 \) & \(5\) \\
|
||||
\(g_{T}\) & --- & --- & --- & \(0.45045\) & \(0.45045\) & \(0.45045\) & \(5\) & \(5\) & \(5\) \\
|
||||
\(g_{Ca,K}\) & --- & --- & --- & --- & --- & --- & \(1\) & \(1\) & \(1\) \\
|
||||
\(g_{Leak}\) &\(0.0205\) & \(0.0205\) & \(0.038\) & \(0.07407\) & \(0.07407\) & \(0.07407\) & \(0.035\) & \(0.035\) & \(0.035\) \\%\Xhline{1\arrayrulewidth}
|
||||
\(\tau_{max, M}\)& \(608\) & \(934\) & \(502\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(C_m\) & \(118.44\) & \(119.99\) & \(101.71\)& \(177.83\) & \(177.83\) & \(177.83\) & \(118.44\)& \(118.44\)& \(118.44\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\end{tabular}}}
|
||||
|
||||
|
||||
\label{tab:g}
|
||||
\end{table}
|
||||
|
||||
|
||||
|
||||
%\input{gating_table}
|
||||
\begin{table}[ht!]
|
||||
\caption[Gating Properties]{ For comparability to typical electrophysiological data fitting reported and for ease of further gating curve manipulations, a sigmoid function (Eqn.\ref{eqn:Boltz}) %Boltzmann \(x_\infty = {\left(\frac{1-a}{1+{exp[{\frac{V-V_{1/2}}{k}}]}} +a\right)^j}\)
|
||||
with slope \(k\), voltage for half-maximal activation or inactivation (\(V_{1/2}\)), exponent \(j\), and persistent current \(0 \leq a \leq 1\) were fitted for the models originating from \citet{pospischil_minimal_2008} (models B, C, D, E, G, H) where \(\alpha_x\) and \(\beta_x\) are used. Gating parameters for \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) are taken from \citet{ranjan_kinetic_2019} and fit to mean wild type parameters in \citet{lauxmann_therapeutic_2021}. Model gating parameters not listed are taken directly from source publication.}
|
||||
\centering
|
||||
\linespread{1.5}\selectfont
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\begin{tabular}{c c c c c c }
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& Gating & \(V_{1/2}\) [mV]& \(k\) & \(j\) & \(a\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
%Pospischil
|
||||
& \(\textrm{I}_{\textrm{Na}}\) activation &\(-34.33054521\) & \(-8.21450277\) & \(1.42295686\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{Na}}\) inactivation &\(-34.51951036\) & \(4.04059373\) & \(1\) & \(0.05\) \\
|
||||
Models & \(\textrm{I}_{\textrm{Kd}}\) activation &\(-63.76096946\) & \(-13.83488194\) & \(7.35347425\) & --- \\
|
||||
B, C, D, E, G, H & \(\textrm{I}_{\textrm{L}}\) activation &\(-39.03684525\) & \(-5.57756176\) & \(2.25190197\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{L}}\) inactivation &\(-57.37\) & \(20.98\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{M}}\) activation &\(-45\) & \(-9.9998807337\) & \(1\) & --- \\ %-45 with 10 mV shift to contributes to resting potential
|
||||
% & & & & &\\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) & \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) activation &\(-30.01851852\) & \(-7.73333333\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) Inactivation &\(-46.85851852\) & \(7.67266667\) & \(1\) & \(0.245\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
|
||||
|
||||
\end{tabular}}
|
||||
|
||||
\label{tab:gating}
|
||||
\end{table}
|
||||
|
||||
\end{document}
|
||||
319
Frontiers_manuscript/For_Submission/Koch_frontiers_revised.tex
Normal file
319
Frontiers_manuscript/For_Submission/Koch_frontiers_revised.tex
Normal file
@ -0,0 +1,319 @@
|
||||
\documentclass[utf8]{FrontiersinHarvard}
|
||||
\DeclareUnicodeCharacter{03B2}{\(\beta\)}
|
||||
\DeclareUnicodeCharacter{03B1}{\(\alpha\)}
|
||||
\DeclareUnicodeCharacter{00C5}{\AA}
|
||||
\usepackage[english]{babel}
|
||||
\usepackage{url,hyperref,lineno,microtype,subcaption}
|
||||
\usepackage[onehalfspacing]{setspace}
|
||||
\usepackage{amsmath}
|
||||
|
||||
\usepackage{xspace}
|
||||
\usepackage{longtable}
|
||||
\usepackage{boldline}
|
||||
\usepackage{geometry}
|
||||
\usepackage{makecell}
|
||||
\usepackage{float}
|
||||
\usepackage{placeins}
|
||||
|
||||
\usepackage[tableposition=top]{caption}
|
||||
|
||||
|
||||
\newcommand{\Kv}{\(\text{K}_{\text{V}}\text{1.1}\)\xspace}
|
||||
\newcommand{\IKv}{\(\text{I}_{\text{K}_{\text{V}}\text{1.1}}\)\xspace}
|
||||
\newcommand{\drheo}{\(\Delta\)rheobase\xspace}
|
||||
\newcommand{\ndAUC}{normalized \(\Delta\)AUC\xspace}
|
||||
|
||||
\newcommand{\note}[2][]{\textbf{[#1: #2]}}
|
||||
\newcommand{\notenk}[1]{\note[NK]{#1}}
|
||||
\newcommand{\notels}[1]{\note[LS]{#1}}
|
||||
\newcommand{\notejb}[1]{\note[JB]{#1}}
|
||||
|
||||
\linenumbers
|
||||
|
||||
\def\keyFont{\fontsize{8}{11}\helveticabold }
|
||||
\def\firstAuthorLast{Koch {et~al.}} %use et al only if is more than 1 author
|
||||
\def\Authors{Nils A. Koch\,$^{1,2}$, Lukas Sonnenberg\,$^{1,2}$, Ulrike B.S. Hedrich\,$^{3}$, Stephan Lauxmann\,$^{1,3}$ and Jan Benda\,$^{1,2,*}$}
|
||||
\def\Address{$^{1}$Institute for Neurobiology, University of T{\"u}bingen, 72072 T{\"u}bingen, Germany \\
|
||||
$^{2}$Bernstein Center for Computational Neuroscience T{\"u}bingen, 72076 T{\"u}bingen, Germany \\
|
||||
$^{3}$Department of Neurology and Epileptology, Hertie Institute for Clinical Brain Research, University of T{\"u}bingen, 72076 T{\"u}bingen, Germany}
|
||||
|
||||
% The Corresponding Author should be marked with an asterisk
|
||||
% Provide the exact contact address (this time including street name and city zip code) and email of the corresponding author
|
||||
\def\corrAuthor{Jan Benda}
|
||||
|
||||
\def\corrEmail{ jan.benda@uni-tuebingen.de}
|
||||
|
||||
|
||||
|
||||
\begin{document}
|
||||
\onecolumn
|
||||
\firstpage{1}
|
||||
|
||||
\title {Loss or Gain of Function? Effects of Ion Channel Mutations on Neuronal Firing Depend on the Neuron Type}
|
||||
|
||||
\author[\firstAuthorLast ]{\Authors} %This field will be automatically populated
|
||||
\address{} %This field will be automatically populated
|
||||
\correspondance{} %This field will be automatically populated
|
||||
\extraAuth{}
|
||||
|
||||
|
||||
\maketitle
|
||||
|
||||
|
||||
\begin{abstract}
|
||||
\section{}
|
||||
Clinically relevant mutations to voltage-gated ion channels, called channelopathies, alter ion channel function, properties of ionic currents and neuronal firing. The effects of ion channel mutations are routinely assessed and characterized as loss of function (LOF) or gain of function (GOF) at the level of ionic currents. However, emerging personalized medicine approaches based on LOF/GOF characterization have limited therapeutic success. Potential reasons are among others that the translation from this binary characterization to neuronal firing is currently not well understood --- especially when considering different neuronal cell types. Here we investigate the impact of neuronal cell type on the firing outcome of ion channel mutations with simulations of a diverse collection of conductance-based neuron models. We systematically analyzed the effects of changes in ion current properties on firing in different neuronal types. Additionally, we simulated the effects of known mutations in the \textit{KCNA1} gene encoding the \Kv potassium channel subtype associated with episodic ataxia type~1 (EA1). These simulations revealed that the outcome of a given change in ion channel properties on neuronal excitability depends on neuron type, i.e. the properties and expression levels of the unaffected ionic currents. Consequently, neuron-type specific effects are vital to a full understanding of the effects of channelopathies on neuronal excitability and are an important step towards improving the efficacy and precision of personalized medicine approaches.
|
||||
|
||||
|
||||
|
||||
|
||||
\tiny
|
||||
\keyFont{ \section{Keywords:} Channelopathy, Epilepsy, Ataxia, Potassium Current, Neuronal Simulations, Conductance-based Models, Neuronal heterogeneity }
|
||||
\end{abstract}
|
||||
\section{Introduction}
|
||||
The properties and combinations of voltage-gated ion channels are vital in determining neuronal excitability \citep{bernard_channelopathies_2008, carbone_ion_2020, rutecki_neuronal_1992, pospischil_minimal_2008}. However, ion channel function can be disturbed, for instance through genetic alterations, resulting in altered neuronal firing behavior \citep{carbone_ion_2020}. In recent years, next generation sequencing has led to an increase in the discovery of clinically relevant ion channel mutations and has provided the basis for pathophysiological studies of genetic epilepsies, pain disorders, dyskinesias, intellectual disabilities, myotonias, and periodic paralyses \citep{bernard_channelopathies_2008, carbone_ion_2020}.
|
||||
Ongoing efforts of many research groups have contributed to the current understanding of underlying disease mechanism in channelopathies. However, a complex pathophysiological landscape has emerged for many channelopathies and is likely a reason for limited therapeutic success with standard care.
|
||||
|
||||
Ion channel variants are frequently classified in heterologous expression systems as either a loss of function (LOF) or a gain of function (GOF) in the respective ionic current \citep{Musto2020, Kullmann2002, Waxman2011, Kim2021}. This LOF/GOF classification is often directly used to predict the effects on neuronal firing \citep{Niday2018, Wei2017, Wolff2017,Masnada2017}, which in turn is important for understanding the pathophysiology of these disorders and for identification of potential therapeutic targets \citep{Orsini2018, Yang2018, Colasante2020, Yu2006}. Experimentally, the effects of channelopathies on neuronal firing are assessed using primary neuronal cultures \citep{Scalmani2006, Smith2018, Liu2019} or \textit{in vitro} recordings from slices of transgenic mouse lines \citep{Mantegazza2019, Xie2010,Lory2020, Habib2015, Hedrich14874}, but are restricted to limited number of different neuron types. Neuron types differ in many aspects. They may differ in their composition of ionic currents \citep{yao2021taxonomy, Cadwell2016, BICCN2021, Scala2021} and therefore likely respond differently to changes in the properties of a single ionic current. Expression level of an affected gene \citep{Layer2021} and relative amplitudes of ionic currents \citep{rutecki_neuronal_1992, pospischil_minimal_2008,Kispersky2012, golowasch_failure_2002, barreiro_-current_2012} indeed dramatically influence the firing behavior and dynamics of neurons. Mutations in different sodium channel genes have been experimentally shown to affect firing in a neuron-type specific manner based on differences in expression levels of the affected gene \citep{Layer2021}, but also on other neuron-type specific mechanisms \citep{Hedrich14874, makinson_scn1a_2016}.
|
||||
|
||||
Neuron-type specificity is likely vital for successful precision medicine treatment approaches. For example, Dravet syndrome was identified as the consquence of LOF mutations in \textit{SCN1A} \citep{Claes2001,Fujiwara2003,Ohmori2002}, however limited success in the treatment of Dravet syndrome persisted \citep{Claes2001,Oguni2001} in part due to lack of understanding that inhibitory interneurons and not pyramidal neurons had altered excitability as a result of LOF \textit{SCN1A} mutations \citep{Yu2006, Colasante2020}.
|
||||
|
||||
Taken together, these examples demonstrate the need to study the effects of ion channel mutations in many different neuron types --- a daunting if not impossible experimental challenge. In the context of this diversity, simulations of conductance-based neuronal models are a powerful tool bridging the gap between altered ionic currents and firing in a systematic and efficient way. Furthermore, simlutions allow to predict the potential effects of drugs needed to alleviate the pathophysiology of the respective mutation \citep{johannesen_genotype-phenotype_2021, lauxmann_therapeutic_2021, Bayraktar}.
|
||||
|
||||
In this study, we therefore investigated how the outcome of ionic current kinetic changes on firing depend on neuronal cell type by (1) characterizing firing responses with two measures, (2) simulating the response of a repertoire of different neuronal models to changes in single current parameters as well as (3) to more complex changes in this case as they were observed for specific \textit{KCNA1} mutations that are associated with episodic ataxia type~1 \citep{Browne1994, Browne1995, lauxmann_therapeutic_2021}.
|
||||
|
||||
\section{Material and Methods}
|
||||
All modelling and simulation was done in parallel with custom written Python 3.8 (Python Programming Language; RRID:SCR\_008394) software, run on a Cent-OS 7 server with an Intel(R) Xeon (R) E5-2630 v2 CPU.
|
||||
|
||||
\subsection{Different Neuron Models}
|
||||
A group of neuronal models representing the major classes of cortical and thalamic neurons including regular spiking pyramidal (RS pyramidal; model D), regular spiking inhibitory (RS inhibitory; model B), and fast spiking (FS; model C) neurons were used \citep{pospischil_minimal_2008}. Additionally, a \Kv current (\IKv; \citealt{ranjan_kinetic_2019}) was added to each of these models (RS pyramidal +\Kv; model H, RS inhibitory +\Kv; model E, and FS +\Kv; model G respectively). A cerebellar stellate cell model from \citet{alexander_cerebellar_2019} is used (Cb stellate; model A) in this study. This neuron model was also extended by a \Kv current \citep{ranjan_kinetic_2019}, either in addition to the A-type potassium current (Cb stellate +\Kv; model F) or by replacing the A-type potassium current (Cb stellate \(\Delta\)\Kv; model J). A subthalamic nucleus (STN; model L) neuron model as described by \citet{otsuka_conductance-based_2004} was also used. The STN neuron model (model L) was additionally extended by a \Kv current \citep{ranjan_kinetic_2019}, either in addition to the A-type potassium current (STN +\Kv; model I) or by replacing the A-type potassium current (STN \(\Delta\)\Kv; model K). Model letter naming corresponds to panel lettering in Figure \ref{fig:diversity_in_firing}. The properties and maximal conductances of each model are detailed in Table \ref{tab:g} and the gating properties are unaltered from the original Cb stellate (model A) and STN (model L) models \citep{alexander_cerebellar_2019, otsuka_conductance-based_2004}. For enabling the comparison of models with the typically reported electrophysiological data fitting reported and for ease of further gating curve manipulations, a modified Boltzmann function
|
||||
|
||||
\begin{equation}\label{eqn:Boltz}
|
||||
x_\infty = {\left(\frac{1-a}{1+{\exp\left[{\frac{V-V_{1/2}}{k}}\right]}} +a\right)^j}
|
||||
\end{equation}
|
||||
|
||||
with slope \(k\), voltage for half-maximal activation or inactivation (\(V_{1/2}\)), exponent \(j\), and persistent current \(0 \leq a \leq 1\) were fitted to the original formulism for RS pyramidal (model D), RS inhibitory (model B) and FS (model C) models from \citet{pospischil_minimal_2008}. The properties of \IKv were fitted to the mean wild type biophysical parameters of \Kv described in \citet{lauxmann_therapeutic_2021}. The fitted gating parameters are detailed in Table \ref{tab:gating}. Each of the original single-compartment models used here can reproduce physiological firing behavior of the neurons they represent (Figure \ref{fig:diversity_in_firing}; \citealt{pospischil_minimal_2008, alexander_cerebellar_2019, otsuka_conductance-based_2004}) and capture key aspects of the dynamics of these neuron types.
|
||||
|
||||
\subsection{Firing Frequency Analysis}
|
||||
The membrane responses to 200 equidistant two second long current steps were simulated using the forward-Euler method with a \(\Delta \textrm{t} = 0.01\)\,ms from steady state. Current steps ranged from 0 to 1\,nA (step size 5\,pA) for all models except for the RS inhibitory neuron models where a range of 0 to 0.35\,nA (step size 1.75\,pA) was used to ensure repetitive firing across the range of input currents. For each current step, action potentials were detected as peaks with at least 50\,mV prominence, or the relative height above the lowest contour line encircling it, and a minimum interspike interval of 1\,ms. The interspike interval was computed and used to determine the instantaneous firing frequencies elicited by the current step.
|
||||
|
||||
To ensure accurate firing frequencies at low firing rates and reduced spike sampling bias, steady-state firing was defined as the mean firing frequency in a 500\,ms window in the last second of the current steps starting at the inital action potential in this last second.
|
||||
Firing characterization was performed in the last second of current steps to ensure steady-state firing is captured and adaptation processes are neglected in our analysis. Alteration in current magnitudes can have different effects on rheobase and the initial slope of the fI curve \citep{Kispersky2012}.
|
||||
For this reason, we quantified neuronal firing using the rheobase as well as the area under the curve (AUC) of the initial portion of the fI curve as a measure of the initial slope of the fI curve Figure \ref{fig:firing_characterization}~A.
|
||||
|
||||
The smallest current at which steady state firing occured was identified and the current step interval preceding the occurrence of steady state firing was simulated at higher resolution (100 current steps) to determine the current at which steady state firing began. Firing was simulated with 100 current steps from this current upwards for 1/5 of the overall current range. Over this range a fI curve was constructed and the integral, or area under the curve (AUC), of the fI curve over this interval was computed with the composite trapezoidal rule and used as a measure of firing rate independent from rheobase.
|
||||
|
||||
To obtain the rheobase at a higher current resolution than the fI curve, the current step interval preceding the occurrence of action potentials was explored at higher resolution with 100 current steps spanning the interval (step sizes of 0.05 pA and 0.0175 pA, respectively). Membrane responses to these current steps were then analyzed for action potentials and the rheobase was considered the lowest current step for which an action potential was elicited.
|
||||
|
||||
All models exhibited tonic steady-state firing with default parameters. In limited instances, variations of parameters elicited periodic bursting, however these instances were excluded from further analysis.
|
||||
|
||||
|
||||
\subsection{Sensitivity Analysis and Comparison of Models}
|
||||
Properties of ionic currents common to all models (\(\text{I}_{\text{Na}}\), \(\text{I}_{\text{Kd}}\), \(\text{I}_{\text{A}}\)/\IKv, and \(\text{I}_{\text{Leak}}\)) were systematically altered in a one-factor-at-a-time sensitivity analysis for all models. The gating curves for each current were shifted (\(\Delta V_{1/2}\)) from -10 to 10\,mV in increments of 1\,mV. The voltage dependence of the time constant associated with the shifted gating curve was correspondingly shifted. The slope (\(k\)) of the gating curves were altered from half to twice the initial slope. Similarly, the maximal current conductance (\(g\)) was also scaled from half to twice the initial value. For both slope and conductance alterations, alterations consisted of 21 steps spaced equally on a \(\textrm{log}_2\) scale. We neglected the variation of time constants for the practical reason that estimation and assessment of time constants and changes to them is not straightforward \citep{Clerx2019, Whittaker2020}.
|
||||
|
||||
\subsection{Model Comparison}
|
||||
Changes in rheobase (\drheo) were calculated in relation to the original model rheobase. The contrast of each AUC value (\(AUC_i\)) was computed in comparison to the AUC of the unaltered wild type model (\(AUC_{wt}\))
|
||||
\begin{equation}\label{eqn:AUC_contrast}
|
||||
\textrm{normalized } \Delta \text{AUC} = \frac{AUC_i - AUC_{wt}}{AUC_{wt}}
|
||||
\end{equation}
|
||||
|
||||
To assess whether the effects of a given alteration on \ndAUC or \drheo were robust across models, the correlation between \ndAUC or \drheo and the magnitude of the alteration of a current property was computed for each alteration in each model and compared across alteration types.
|
||||
The Kendall's \(\tau\) coefficient, a non-parametric rank correlation, is used to describe the relationship between the magnitude of the alteration and AUC or rheobase values. A Kendall \(\tau\) value of -1 or 1 is indicative of monotonically decreasing and increasing relationships respectively.
|
||||
|
||||
\subsection{\textit{KCNA1} Mutations}\label{subsec:mut}
|
||||
Known episodic ataxia type~1 associated \textit{KCNA1} mutations and their electrophysiological characterization have been reviewed in \citet{lauxmann_therapeutic_2021}. The mutation-induced changes in \IKv amplitude and activation slope (\(k\)) were normalized to wild type measurements and changes in activation \(V_{1/2}\) were used relative to wild type measurements. Although initially described to lack fast activation, \Kv displays prominent inactivation at physiologically relevant temperatures \citep{ranjan_kinetic_2019}. The effects of a mutation were also applied to \(\text{I}_{\text{A}}\) when present as both potassium currents display inactivation. In all cases, the mutation effects were applied to half of the \Kv or \(\text{I}_{\text{A}}\) under the assumption that the heterozygous mutation results in 50\% of channels carrying the mutation. Frequency-current curves for each mutation in each model were obtained through simulation and used to characterize firing behavior as described above. For each model the differences in mutation AUC to wild type AUC were normalized by wild type AUC (\ndAUC) and mutation rheobases were compared to wild type rheobase values (\drheo). Pairwise Kendall rank correlations (Kendall \(\tau\)) were used to compare the correlation in the effects of \Kv mutations on AUC and rheobase between models.
|
||||
|
||||
\subsection{Code Accessibility}
|
||||
The simulation and analysis code including full specification of the models is freely available online at \href{https://github.com/nkoch1/LOFGOF2023}{https://github.com/nkoch1/LOFGOF2023}. %The code is available as \ref{code_zip}.
|
||||
|
||||
\section{Results}
|
||||
|
||||
To examine the role of neuron-type specific ionic current environments on the impact of altered ion currents properties on firing behavior:
|
||||
(1) firing responses were characterized with rheobase and \(\Delta\)AUC, (2) a set of neuronal models was used and properties of channels common across models were altered systematically one at a time, and (3) the effects of a set of episodic ataxia type~1 associated \textit{KCNA1} mutations on firing was then examined across different neuronal models with different ionic current environments.
|
||||
|
||||
\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 (Figure \ref{fig:diversity_in_firing}). The models chosen for this study all fire tonically and do not exhibit bursting (see methods for details and naming of the models). Models are qualitatively sorted based on their firing curves and labeled model A through L accordingly. Some models, such as models A and B, display type I firing, whereas others such as models J and L exhibit 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; \citealt{ermentrout_type_1996, Rinzel_1998}). 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 models I, J, and K have large hysteresis (Figure \ref{fig:diversity_in_firing} and Supplementary Figure S1). Different types of underlying current dynamics are known to generate these different firing types and hysteresis \cite{ERMENTROUT2002, ermentrout_type_1996, Izhikevich2006}. This broad range of single-compartmental models represents the distinct dynamics of various neuron types across diverse brain regions.
|
||||
|
||||
\subsection{Characterization of Neuronal Firing Properties}
|
||||
Neuronal firing is a complex phenomenon, and a quantification of firing properties is required for comparisons across neuron types and between different conditions. Here we focus on two aspects of firing that are routinely measured in clinical settings \citep{Bryson_2020}: rheobase, the smallest injected current at which the neuron 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 (Figure \ref{fig:firing_characterization}~A). The characterization of the firing properties of a neuron by using 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) by two independent measures. Note that AUC is essentially quantifying the slope of a neuron's fI curve.
|
||||
|
||||
Using these two measures we quantified 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 of 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 (Figure \ref{fig:firing_characterization}~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 gain of function (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, loss of function (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.
|
||||
|
||||
|
||||
\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 delayed rectifier potassium current in the model G to more depolarized values, then the rheobase of the resulting fI curves shifted to lower currents \(-\)\drheo, making the neuron more sensitive to weak inputs, but at the same time the slope of the fI curves was reduced (\(-\)\ndAUC), which resulted in a reduced firing rate (Figure \ref{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 model G is in the bottom left quadrant of Figure \ref{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 Figure \ref{fig:AUC_correlation}~B). For each of the many models we got a different relation between the changes in AUC and the shifts in half maximal potential \(V_{1/2}\) (thin lines in Figure \ref{fig:AUC_correlation}~B). To further summarize these different dependencies of the various models we characterized each of these curves by a single number, the \( \text{Kendall} \ \tau \) correlation coefficient. A monotonically increasing curve resulted in a \( \text{Kendall} \ \tau \) close to \(+1\) a monotonously decreasing curve in \( \text{Kendall} \ \tau \approx -1 \), and a non-monotonous, non-linear relation in \( \text{Kendall} \ \tau \) close to zero (compare lines in Figure \ref{fig:AUC_correlation}~B with dots in black box in panel C).
|
||||
|
||||
Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affected the AUC (Figure \ref{fig:AUC_correlation}), but how exactly the AUC was affected usually depended on the specific neuronal model. Increasing the slope factor of the \Kv activation curve for example increased the AUC in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (Figure \ref{fig:AUC_correlation}~D,E,F). Similar consistent positive correlations could 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 resulted in negative correlations with the AUC in all models (\( \text{Kendall} \ \tau \approx -1\)).
|
||||
|
||||
Qualitative differences could be found, for example, when increasing the maximal conductance of the delayed rectifier (Figure \ref{fig:AUC_correlation}~G,H,I). In some model neurons this increased AUC (\( \text{Kendall} \ \tau \approx +1\)), whereas in others AUC was decreased (\( \text{Kendall} \ \tau \approx -1\)). In model I, AUC depended in a non-linear way on the maximal conductance of the delayed rectifier, resulting in a \( \text{Kendall} \ \tau \) close to zero. Even more dramatic qualitative differences between models resulted from shifts of the activation curve of the delayed rectifier, as discussed already above (Figure \ref{fig:AUC_correlation}~A,B,C). Some model neurons did almost not depend on changes in Kd-current half activation \(V_{1/2}\) or showed strong non-linear dependencies, both resulting in \( \text{Kendall} \ \tau\) close to zero. Many model neurons showed strongly negative correlations, and a few displayed positive correlations with shifting the activation curve of the delayed rectifier.
|
||||
|
||||
Changes in gating half activation potential \(V_{1/2}\) and slope factor \(k\) as well as the maximum conductance \(g\) affected rheobase (Figure \ref{fig:rheobase_correlation}). However, in contrast to AUC, qualitatively consistent effects on rheobase across models could be observed. An increasing of the maximal conductance of the leak current in the model A increased the rheobase (Figure \ref{fig:rheobase_correlation}~G). When these changes were plotted against the change in maximal conductance a monotonically increasing relationship was evident (thick teal line in Figure \ref{fig:rheobase_correlation}~H). This monotonically increasing relationship was evident in all models (\( \text{Kendall} \ \tau \approx +1\)), but with different slopes (thin lines in Figure \ref{fig:rheobase_correlation}~H). Similarly, positive correlations were consistently found across models for maximal conductances of delayed rectifier Kd, \Kv, and A type currents, whereas the maximal conductance of the sodium current was consistently associated with negative correlations (\( \text{Kendall} \ \tau \approx -1\); Figure \ref{fig:rheobase_correlation}~I), i.e. rheobase decreased with increasing maximum conductance in all models.
|
||||
|
||||
Although changes in half maximal potential \(V_{1/2}\) and slope factor \(k\) generally correlated with rheobase similarly across models there were some exceptions. Rheobase was affected with both with positive and negative correlations in different models as a result of changing slope factor of Na\(^+\)-current inactivation (positive: models A--H and J; negative: models I, K and L), \Kv-current inactivation (positive: models I and K; negative: models E--G, J, H), and A-current activation (positive: models A, F and L; negative: model I; Figure \ref{fig:rheobase_correlation}~F). Departures from monotonic relationships also occurred in some models as a result of Kd-current activation \(V_{1/2}\) (e.g. model J) and slope factor \(k\) (models F and G), \Kv-current inactivation slope factor \(k\) (model K), and A-current activation slope factor \(k\) (model L). 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.
|
||||
|
||||
\subsection{\textit{KCNA1} Mutations}
|
||||
Mutations in \textit{KCNA1} are associated with episodic ataxia type~1 (EA1) and have been characterized biophysically (as reviewed by \citet{lauxmann_therapeutic_2021}). Here they were used as a test case in the effects of various ionic current environments on neuronal firing and on the outcomes of channelopathies. The changes in AUC and rheobase from wild type values for reported EA1 associated \textit{KCNA1} mutations were heterogeneous across models containing \Kv, but generally showed decreases in rheobase (Figure \ref{fig:simulation_model_comparision}A--I). Pairwise non-parametric Kendall \(\tau\) rank correlations between the simulated effects of these \Kv mutations on rheobase were highly correlated across models (Figure \ref{fig:simulation_model_comparision}J) indicating that EA1 associated \textit{KCNA1} mutations generally decrease rheobase across diverse neuron-types. However, the effects of the \Kv mutations on AUC were more heterogenous as reflected by both weak and strong positive and negative pairwise correlations between models (Figure \ref{fig:simulation_model_comparision}K), suggesting that the effects of ion-channel variant on super-threshold neuronal firing depend both quantitatively and qualitatively on the specific composition of ionic currents in a given neuron.
|
||||
|
||||
|
||||
\section{Discussion}
|
||||
|
||||
To compare the effects of ion channel mutations on neuronal firing of different neuron types, we used a diverse set of conductance-based models to systematically characterize the effects of changes in individual channel properties. Additionally, we simulated the effects of specific episodic ataxia type~1 associated (EA1) \textit{KCNA1} mutations. Changes to single ionic current properties, as well as known EA1 associated \textit{KCNA1} mutations showed consistent effects on the rheobase across neuron types, whereas the effects on the slope of the steady-state fI-curve depended on the neuron type. Our results demonstrate that loss of function (LOF) and gain of function (GOF) on the biophysical level cannot be uniquely transferred to the level of neuronal firing. Thus, the effects caused by different mutations depend on the properties of the other ion channels expressed in a neuron and are therefore depend on the channel ensemble of a specific neuron type.
|
||||
|
||||
\subsection{Firing Frequency Analysis}
|
||||
Although differences in neuronal firing can be characterized by an area under the curve of the fI curve for a fixed current range, this approach characterizes firing as a mixture of key features: rheobase and the initial slope of the fI curve. By probing rheobase directly and using an AUC relative to rheobase, we disambiguate these features and enable insights into the effects on rheobase and initial fI curve steepness. This increases the specificity of our understanding of how ion channel mutations alter firing across neuron types and enable classification as described in Figure \ref{fig:firing_characterization}. Importantly, in cases when ion channel mutations alter rheobase and initial fI curve steepness in ways that oppose effects on firing (upper left and bottom right quadrants of Figure \ref{fig:firing_characterization}~B), this disamgibuation is important for understanding the outcome of the mutation. In these cases, the regime the neuron is operating in is vital in determining the neuron's firing outcome. If it is in its excitable regime and only occasionally generates an action potential, then the effect on the rheobase is more important. If it is firing periodically with high rates, then the change in AUC might be more relevant.
|
||||
|
||||
\subsection{Modelling Limitations}
|
||||
The models used here are simple and they all capture key aspects of the firing dynamics for their respective neuron. The simple models fall short of capturing the complex physiology, biophysics and heterogeneity of real neurons, nor do they take into account subunit stoichiometry, auxillary subunits, membrane composition which influence the biophysics of ionic currents \citep{Al-Sabi_2013, Oliver_2004, Pongs_2009, Rettig1994}. However, for the purpose of understanding how different neuron-types, or current environments, contribute to the diversity in firing outcomes of ion channel mutations, the fidelity of the models to the physiological neurons they represent is of a minor concern. For exploring possible neuron-type specific effects, variety in currents and dynamics across models is of utmost importance. With this context in mind, the collection of models used here are labelled as models A-L to highlight that the physiological neurons they represent is not of chief concern, but rather that the collection of models with different attributes respond heterogeneously to the same perturbation. Additionally, the development of more realistic models is a high priority and will enable neuron-type specific predictions that may aid precision medicine approaches. Thus, weight should not be put on any single predicted firing outcome here in a specific model, but rather on the differences in outcomes that occur across the neuron-type spectrum the models used here represent. % Menchaca_2012,
|
||||
|
||||
\subsection{Neuronal Diversity}
|
||||
The nervous system consists of a vastly diverse and heterogenous collection of neurons with variable properties and characteristics including diverse combinations and expression levels of ion channels which are vital for neuronal firing dynamics.
|
||||
|
||||
Advances in high-throughput techniques have enabled large-scale investigation into single-cell properties across the CNS \citep{Poulin2016} that have revealed large diversity in neuronal gene expression, morphology and neuronal types in the motor cortex \citep{Scala2021}, neocortex \cite{Cadwell2016, Cadwell2020}, GABAergic neurons in the cortex and retina \citep{Huang2019, Laturnus2020}, cerebellum \citep{Kozareva2021}, spinal cord \citep{Alkaslasi2021}, visual cortex \citep{Gouwens2019} as well as the retina \citep{Baden2016, Voigt2019, Berens2017, Yan2020a, Yan2020b}.
|
||||
|
||||
Diversity across neurons is not limited to gene expression and can also be seen electrophysiologically \citep{Tripathy2017, Gouwens2018, Tripathy2015, Scala2021, Cadwell2020, Gouwens2019, Baden2016, Berens2017} with correlations existing between gene expression and electrophysiological properties \citep{Tripathy2017}. At the ion channel level, diversity exists not only between the specific ion channels the different neuron types express but heterogeneity also exists in ion channel expression levels within neuron types \citep{marder_multiple_2011, goaillard_ion_2021,barreiro_-current_2012}. As ion channel properties and expression levels are key determinents of neuronal dynamics and firing \citep{Balachandar2018, Gu2014, Zeberg2015, Aarhem2007, Qi2013, Gu2014a, Zeberg2010, Zhou2020, Kispersky2012} neurons with different ion channel properties and expression levels display different firing properties.
|
||||
|
||||
To capture the diversity in neuronal ion channel expression and its relevance in the outcome of ion channel mutations, we used multiple neuronal models with different ionic currents and underlying firing dynamics here.
|
||||
|
||||
\subsection{Ionic Current Environments Determine the Effect of Ion Channel Mutations}
|
||||
|
||||
To our knowledge, no comprehensive evaluation of how ionic current environment and neuron type affect the outcome of ion channel mutations have been reported. However, comparisons between the effects of such mutations between certain neuron types were described. For instance, the R1648H mutation in \textit{SCN1A} does not alter the excitability of cortical pyramidal neurons, but causes hypoexcitability of adjacent inhibitory GABAergic neurons \citep{Hedrich14874}. In the CA3 region of the hippocampus, the equivalent mutation in \textit{SCN8A}, R1627H, increases the excitability of pyramidal neurons and decreases the excitability of parvalbumin positive interneurons \cite{makinson_scn1a_2016}. Additionally, the L858H mutation in \(\text{Na}_\text{V}\text{1.7}\), associated with erythermyalgia, has been shown to cause hypoexcitability in sympathetic ganglion neurons and hyperexcitability in dorsal root ganglion neurons \citep{Waxman2007, Rush2006}. The differential effects of L858H \(\text{Na}_\text{V}\text{1.7}\) on firing is dependent on the presence or absence of another sodium channel, namely the \(\text{Na}_\text{V}\text{1.8}\) subunit \citep{Waxman2007, Rush2006}. These findings, in concert with our findings emphasize that the ionic current environment in which a channelopathy occurs is vital in determining the outcomes of the channelopathy on firing.
|
||||
|
||||
Neuron type specific differences in ionic current properties are important in the effects of ion channel mutations. However, within a neuron type heterogeneity in channel expression levels exists and it is often desirable to generate a population of neuronal models and to screen them for plausibility to biological data in order to capture neuronal population diversity \citep{marder_multiple_2011,OLeary2016}. The models we used here are originally generated by characterization of current gating properties and by fitting of maximal conductances to experimental data \citep{pospischil_minimal_2008, ranjan_kinetic_2019, alexander_cerebellar_2019, otsuka_conductance-based_2004}. This practice of fixing maximal conductances based on experimental data is limiting as it does not reproduce the variability in channel expression and neuronal firing behavior of a heterogeneous neuron population \citep{verma_computational_2020}. For example, a model derived from the mean conductances in a neuronal sub-population within the stomatogastric ganglion, the so-called "one-spike bursting" neurons fire three spikes instead of one per burst due to an L-shaped distribution of sodium and potassium conductances \citep{golowasch_failure_2002}. Multiple sets of conductances can give rise to the same patterns of activity also termed degeneracy and differences in neuronal dynamics may only be evident with perturbations \citep{marder_multiple_2011, goaillard_ion_2021}. The variability in ion channel expression often correlates with the expression of other ion channels \citep{goaillard_ion_2021} and neurons whose behavior is similar may possess correlated variability across different ion channels resulting in stability in the neuronal phenotype \citep{lamb_correlated_2013, soofi_co-variation_2012, taylor_how_2009}. The variability of ionic currents and degeneracy of neurons may account, at least in part, for the observation that the effect of toxins within a neuronal type is frequently not constant \citep{khaliq_relative_2006, puopolo_roles_2007, ransdell_neurons_2013}.
|
||||
|
||||
\subsection{Effects of \textit{KCNA1} Mutations}
|
||||
Changes in delayed rectifier potassium currents, analogous to those seen in LOF \textit{KCNA1} mutations, change the underlying firing dynamics of the Hodgkin Huxley model result in reduced thresholds for repetitive firing and thus contribute to increased excitability \citep{hafez_altered_2020}. Although the Hodgkin Huxley delayed rectifier lacks inactivation, the increases in excitability observed by \citet{hafez_altered_2020} are in line with our simulation-based predictions of the outcomes of \textit{KCNA1} mutations. LOF \textit{KCNA1} mutations generally increase neuronal excitability, however the varying susceptibility on rheobase and different effects on AUC of the fI-curve of \textit{KCNA1} mutations across models are indicative that a certain neuron type specific complexity exists. Increased excitability is seen experimentally with \Kv null mice \citep{smart_deletion_1998, zhou_temperature-sensitive_1998}, with pharmacological \Kv block \citep{chi_manipulation_2007, morales-villagran_protection_1996} and by \citet{hafez_altered_2020} with simulation-based predictions of \textit{KCNA1} mutations. Contrary to these results, \citet{zhao_common_2020} predicted \textit{in silico} that the depolarizing shifts seen as a result of \textit{KCNA1} mutations broaden action potentials and interfere negatively with high frequency action potential firing. However, they varied stimulus duration between different models and therefore comparability of firing rates is lacking in this study.
|
||||
|
||||
In our simulations, different current properties alter the impact of \textit{KCNA1} mutations on firing as evident in the differences seen in the impact of \(\text{I}_\text{A}\) and \IKv in the Cb stellate and STN model families on \textit{KCNA1} mutation firing. This highlights that not only knowledge of the biophysical properties of a channel but also its neuronal expression and other neuronal channels present is vital for the holistic understanding of the effects of a given ion channel mutation both at the single cell and network level.
|
||||
|
||||
\subsection{Loss or Gain of Function Characterizations Do Not Fully Capture Ion Channel Mutation Effects on Firing}
|
||||
The effects of changes in channel properties depend in part on the neuronal model in which they occur and can be seen in the variance of correlations (especially in AUC of the fI-curve) across models for a given current property change. Therefore, relative conductances and gating properties of currents in the ionic current environment in which an alteration in current properties occurs plays an important role in determining the outcome on firing. The use of LOF and GOF is useful at the level of ion channels to indicate whether a mutation results in more or less ionic current. However, the extension of this thinking onto whether mutations induce LOF or GOF at the level of neuronal firing based on the ionic current LOF/GOF is problematic because of this dependency of neuronal firing changes on the ionic channel environment. Thus, the direct leap from current level LOF/GOF characterizations to effects on firing without experimental or modelling-based evidence, although tempting, should be refrained from and viewed with caution when reported. This is especially relevant in the recent development of personalized medicine for channelopathies, where a patient's specific channelopathy is identified and used to tailor treatments \citep{Weber2017, Ackerman2013, Helbig2020, Gnecchi2021, Musto2020, Brunklaus2022, Hedrich2021}. In these cases the effects of specific ion channel mutations are often characterized based on ionic currents in expression systems and classified as LOF or GOF to aid in treatment decisions \citep{johannesen_genotype-phenotype_2021, Brunklaus2022, Musto2020}. Although positive treatment outcomes occur with sodium channel blockers in patients with GOF \(\text{Na}_{\text{V}}\text{1.6}\) mutations, patients with both LOF and GOF \(\text{Na}_{\text{V}}\text{1.6}\) mutations can benefit from treatment with sodium channel blockers \citep{johannesen_genotype-phenotype_2021}. This example suggests that the relationship between effects at the level of ion channels and effects at the level of firing and therapeutics is not linear or evident without further contextual information.
|
||||
|
||||
Therefore, the transfer of LOF or GOF from the current to the firing level should be used with caution; the neuron type in which the mutant ion channel is expressed may provide valuable insight into the functional consequences of an ion channel mutation. Experimental assessment of the effects of a patient's specific ion channel mutation \textit{in vivo} is not generally feasible at a large scale. Therefore, modelling approaches investigating the effects of patient specific channelopathies provide a viable method bridging between characterization of changes in biophysical properties of ionic currents and the firing consequences of these effects. In both experimental and modelling studies on the effects of ion channel mutations on neuronal firing the specific dependency on neuron type should be considered.
|
||||
|
||||
The effects of altered ion channel properties on firing is generally influenced by the other ionic currents present in the neuron. In channelopathies the effect of a given ion channel mutation on neuronal firing therefore depends on the neuron type in which those changes occur \citep{Hedrich14874, makinson_scn1a_2016, Waxman2007, Rush2006}. Although certain complexities of neurons such as differences in neuron-type sensitivities to current property changes, interactions between ionic currents, cell morphology and subcellular ion channel distribution are neglected here, it is likely that this increased complexity \textit{in vivo} would contribute to the neuron-type dependent effects on neuronal firing. The complexity and nuances of the nervous system, including neuron-type dependent firing effects of channelopathies explored here, likely underlie shortcomings in treatment approaches in patients with channelopathies. Accounting for neuron-type dependent firing effects provides an opportunity to improve the efficacy and precision in personalized medicine approaches.
|
||||
|
||||
With this study we suggest that neuron-type specific effects are vital to a full understanding of the effects of channelopathies at the level of neuronal firing. Furthermore, we highlight the use of modelling approaches to enable relatively fast and efficient insight into channelopathies.
|
||||
|
||||
|
||||
\bibliographystyle{Frontiers-Harvard}
|
||||
\bibliography{Koch_ref}
|
||||
\newpage
|
||||
\section*{Figures}
|
||||
%%% NB logo1.eps is required in the path in order to correctly compile front page header %%%
|
||||
|
||||
|
||||
\begin{figure}[H]
|
||||
\centering
|
||||
\includegraphics[width=0.875\linewidth]{diversity_in_firing.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Diversity in Neuronal Model Firing. Spike trains (left), frequency-current (fI) curves (right) for Cb stellate \textbf{(A)}, RS inhibitory \textbf{(B)}, FS \textbf{(C)}, RS pyramidal \textbf{(D)}, RS inhibitory +\Kv \textbf{(E)}, Cb stellate +\Kv \textbf{(F)}, FS +\Kv \textbf{(G)}, RS pyramidal +\Kv \textbf{(H)}, STN +\Kv \textbf{(I)}, Cb stellate \(\Delta\)\Kv \textbf{(J)}, STN \(\Delta\)\Kv \textbf{(K)}, and STN \textbf{(L)} neuron models. Models are sorted qualitatively based on their fI curves. Black markers on the fI curves indicate the current step at which the spike train occurs. The green marker indicates the current at which firing begins in response to an ascending current ramp, whereas the red marker indicates the current at which firing ceases in response to a descending current ramp (see Supplementary Figure S1).} % \textcolor{red}{(see Supplementary Figure \ref{ramp_firing}).}}
|
||||
\label{fig:diversity_in_firing}
|
||||
\end{figure}
|
||||
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=0.5\linewidth]{firing_characterization_arrows.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Characterization of firing with AUC and rheobase. \textbf{(A)} The area under the curve (AUC) of the repetitive firing frequency-current (fI) curve. \textbf{(B)}
|
||||
Changes in firing as characterized by \(\Delta\)AUC and \(\Delta\)rheobase occupy four quadrants separated by no changes in AUC and rheobase. Representative schematic fI curves in red with respect to a reference (or wild type) fI curve (blue) depict the general changes associated with each quadrant.}
|
||||
\label{fig:firing_characterization}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{AUC_correlation.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of altered channel kinetics on AUC in various neuron models. The fI curves corresponding to shifts in model G delayed rectifier Kd half activation \(V_{1/2}\) \textbf{(A)}, changes \Kv activation slope factor \(k\) in model G \textbf{(D)}, and changes in maximal conductance of delayed rectifier K current in the model I \textbf{(G)} are shown. The fI curves from the smallest (grey) to the largest (black) alterations are seen for \textbf{(A,D,} and \textbf{G}) in accordance to the greyscale of the x-axis in B, E, and H. The \ndAUC of fI curves is plotted against delayed rectifier K half activation potential (\(\Delta V_{1/2}\); \textbf{B}), \Kv activation slope factor \(k\) (k/\(\text{k}_{\text{WT}}\); \textbf{E}) and maximal conductance \(g\) of the delayed rectifier Kd current (g/\(\text{g}_{\text{WT}}\); \textbf{H}) for all models (thin lines) with relationships from the fI curve examples ( \textbf{A, D, G} respectively) highlighted by thick lines with colors corresponding to the box highlighting each set of fI curves. The Kendall rank correlation (Kendall \(\tau\)) coefficients between shifts in half maximal potential \(V_{1/2}\) and \ndAUC \textbf{(C)}, slope factor k and \ndAUC \textbf{(F)} as well as maximal current conductances and \ndAUC \textbf{(I)} for each model and current property is computed. The relationships between \(\Delta V_{1/2}\), k/\(\text{k}_{\text{WT}}\), and g/\(\text{g}_{\text{WT}}\) and \ndAUC for the Kendall rank correlations highlighted in the black boxes are depicted in \textbf{(B)}, \textbf{(E)} and \textbf{(H)} respectively. }
|
||||
\label{fig:AUC_correlation}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{rheobase_correlation.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of altered channel kinetics on rheobase. The fI curves corresponding to shifts in model G \Kv activation \(V_{1/2}\) \textbf{(A)}, changes \Kv inactivation slope factor \(k\) in model F \textbf{(D)}, and changes in maximal conductance of the leak current in model A \textbf{(G)} are shown. The fI curves from the smallest (grey) to the largest (black) alterations are seen for ( \textbf{A,D,} and \textbf{G}) in accordance to the greyscale of the x-axis in \textbf{B}, \textbf{E}, and \textbf{H}. The \drheo of fI curves is plotted against \Kv half activation potential (\(\Delta V_{1/2}\); \textbf{B}), \Kv inactivation slope factor \(k\) (k/\(\text{k}_{\text{WT}}\); \textbf{E}) and maximal conductance \(g\) of the leak current (g/\(\text{g}_{\text{WT}}\); \textbf{H}) for all models (thin lines) with relationships from the fI curve examples ( \textbf{A, D, G} respectively) highlighted by thick lines with colors corresponding to the box highlighting each set of fI curves. The Kendall rank correlation (Kendall \(\tau\)) coefficients between shifts in half maximal potential \(V_{1/2}\) and \drheo \textbf{(C)}, slope factor k and \drheo \textbf{(F)} as well as maximal current conductances and \drheo \textbf{(I)} for each model and current property is computed. The relationships between \(\Delta V_{1/2}\), k/\(\text{k}_{\text{WT}}\), and g/\(\text{g}_{\text{WT}}\) and \drheo for the Kendall rank correlations highlighted in the black boxes are depicted in \textbf{(B)}, \textbf{(E)} and \textbf{(H)} respectively.}
|
||||
\label{fig:rheobase_correlation}
|
||||
\end{figure}
|
||||
|
||||
\begin{figure}[!ht]
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{simulation_model_comparison.jpg}
|
||||
\linespread{1.}\selectfont
|
||||
\caption[]{Effects of episodic ataxia type~1 associated \textit{KCNA1} mutations on firing. Effects of \textit{KCNA1} mutations on AUC (percent change in normalized \(\Delta\)AUC) and rheobase (\(\Delta\)Rheobase) compared to wild type for model H \textbf{(A)}, model E \textbf{(B)}, model G \textbf{(C)}, model A \textbf{(D)}, model F \textbf{(E)}, model J \textbf{(F)}, model L \textbf{(G)}, model I \textbf{(H)} and model K \textbf{(I)}. All \textit{KCNA1} Mutations are marked in grey with the V174F, F414C, E283K, and V404I \textit{KCNA1} mutations highlighted in color for each model. Pairwise Kendall rank correlation coefficients (Kendall \(\tau\)) between the effects of \textit{KCNA1} mutations on rheobase and on AUC are shown in J and K respectively. Marker shape is indicative of model/firing type, and grey dashed lines denote the quadrants of firing characterization (see Figure \ref{fig:firing_characterization}).}
|
||||
\label{fig:simulation_model_comparision}
|
||||
\end{figure}
|
||||
|
||||
|
||||
|
||||
\FloatBarrier
|
||||
\section*{Tables}
|
||||
|
||||
%\input{g_table}
|
||||
\begin{table}[ht!]
|
||||
% \fontsize{10pt}{10pt}
|
||||
\caption[Cell properties and conductances of neuronal models]{Cell properties and conductances of regular spiking pyramidal neuron (RS Pyramidal; model D), regular spiking inhibitory neuron (RS Inhibitory; model B), fast spiking neuron (FS; model C) each with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (RS Pyramidal +\Kv; model H, RS Inhibitory +\Kv; model E, FS +\Kv; model G respectively), cerebellar stellate cell (Cb Stellate; model A), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (Cb Stellate +\Kv; model F) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_\textrm{A}\) (Cb Stellate \(\Delta\)\Kv; model J), and subthalamic nucleus neuron (STN; model L), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (STN +\Kv; model I) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_{\textrm{A}}\) (STN \Kv; model K) models. All conductances are given in \(\textrm{mS}/\textrm{cm}^2\). Capacitances (\(C_m\)) and \(\tau_{max, M}\) are given in pF and ms respectively.}
|
||||
\centering
|
||||
% \begin{tabular}[x]{@{}l@{}} \\ \end{tabular}
|
||||
\linespread{1.5}\selectfont
|
||||
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\resizebox{\textwidth}{!}{\begin{tabular}{cccccccccc}
|
||||
% in mS/cm^2
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& \begin{tabular}[x]{@{}c@{}} RS\\Pyra-\\midal\\(+\Kv) \end{tabular} & \begin{tabular}[x]{@{}c@{}} RS\\Inhib-\\itory\\(+\Kv)\end{tabular} & \begin{tabular}[x]{@{}c@{}}FS\\(+\Kv) \end{tabular}& \begin{tabular}[x]{@{}c@{}} Cb\\Stellate \end{tabular}& \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\+\Kv \end{tabular} & \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\\(\Delta\)\Kv \end{tabular} & STN &\begin{tabular}[x]{@{}c@{}} STN\\+\Kv \end{tabular} &\begin{tabular}[x]{@{}c@{}} STN\\\(\Delta\)\Kv \end{tabular} \\
|
||||
Model & D (H) & B (E) & C (G) & A & F & J & L & I & K \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(g_{Na}\) & \(56\) & \(10\) & \(58\) & \(3.4\) & \(3.4\) & \(3.4\) & \(49\) & \(49\) & \(49\) \\
|
||||
\(g_{Kd}\) & \(6\) (\(5.4\)) & 2.1 (\(1.89\)) & 3.9 (\(3.51\)) & \(9.0556\) & \(8.15\) &\(9.0556\) & \(57\) & \(56.43\) & \(57\) \\
|
||||
\(g_{K_V1.1}\) & --- (\(0.6\)) & --- (\(0.21\)) & --- (\(0.39\)) & --- & \(0.90556\) & \(1.50159\) & --- & \(0.57\) & \(0.5\) \\
|
||||
\(g_{A}\) & --- & --- & --- & \(15.0159\) & \(15.0159\) & --- & \(5\) & \(5\) & --- \\
|
||||
\(g_{M}\) & \(0.075\) & \(0.0098\) &\(0.075\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(g_{L}\) & --- & --- & --- & --- & --- & --- & \(5\) & \(5 \) & \(5\) \\
|
||||
\(g_{T}\) & --- & --- & --- & \(0.45045\) & \(0.45045\) & \(0.45045\) & \(5\) & \(5\) & \(5\) \\
|
||||
\(g_{Ca,K}\) & --- & --- & --- & --- & --- & --- & \(1\) & \(1\) & \(1\) \\
|
||||
\(g_{Leak}\) &\(0.0205\) & \(0.0205\) & \(0.038\) & \(0.07407\) & \(0.07407\) & \(0.07407\) & \(0.035\) & \(0.035\) & \(0.035\) \\%\Xhline{1\arrayrulewidth}
|
||||
\(\tau_{max, M}\)& \(608\) & \(934\) & \(502\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(C_m\) & \(118.44\) & \(119.99\) & \(101.71\)& \(177.83\) & \(177.83\) & \(177.83\) & \(118.44\)& \(118.44\)& \(118.44\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\end{tabular}}}
|
||||
|
||||
|
||||
\label{tab:g}
|
||||
\end{table}
|
||||
|
||||
|
||||
|
||||
%\input{gating_table}
|
||||
\begin{table}[ht!]
|
||||
\caption[Gating Properties]{ For comparability to typical electrophysiological data fitting reported and for ease of further gating curve manipulations, a sigmoid function (Eqn.\ref{eqn:Boltz}) %Boltzmann \(x_\infty = {\left(\frac{1-a}{1+{exp[{\frac{V-V_{1/2}}{k}}]}} +a\right)^j}\)
|
||||
with slope \(k\), voltage for half-maximal activation or inactivation (\(V_{1/2}\)), exponent \(j\), and persistent current \(0 \leq a \leq 1\) were fitted for the models originating from \citet{pospischil_minimal_2008} (models B, C, D, E, G, H) where \(\alpha_x\) and \(\beta_x\) are used. Gating parameters for \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) are taken from \citet{ranjan_kinetic_2019} and fit to mean wild type parameters in \citet{lauxmann_therapeutic_2021}. Model gating parameters not listed are taken directly from source publication.}
|
||||
\centering
|
||||
\linespread{1.5}\selectfont
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\begin{tabular}{c c c c c c }
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& Gating & \(V_{1/2}\) [mV]& \(k\) & \(j\) & \(a\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
%Pospischil
|
||||
& \(\textrm{I}_{\textrm{Na}}\) activation &\(-34.33054521\) & \(-8.21450277\) & \(1.42295686\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{Na}}\) inactivation &\(-34.51951036\) & \(4.04059373\) & \(1\) & \(0.05\) \\
|
||||
Models & \(\textrm{I}_{\textrm{Kd}}\) activation &\(-63.76096946\) & \(-13.83488194\) & \(7.35347425\) & --- \\
|
||||
B, C, D, E, G, H & \(\textrm{I}_{\textrm{L}}\) activation &\(-39.03684525\) & \(-5.57756176\) & \(2.25190197\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{L}}\) inactivation &\(-57.37\) & \(20.98\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{M}}\) activation &\(-45\) & \(-9.9998807337\) & \(1\) & --- \\ %-45 with 10 mV shift to contributes to resting potential
|
||||
% & & & & &\\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) & \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) activation &\(-30.01851852\) & \(-7.73333333\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) Inactivation &\(-46.85851852\) & \(7.67266667\) & \(1\) & \(0.245\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
|
||||
|
||||
\end{tabular}}
|
||||
|
||||
\label{tab:gating}
|
||||
\end{table}
|
||||
|
||||
\end{document}
|
||||
2233
Frontiers_manuscript/For_Submission/Koch_ref.bib
Normal file
2233
Frontiers_manuscript/For_Submission/Koch_ref.bib
Normal file
File diff suppressed because it is too large
Load Diff
2220
Frontiers_manuscript/For_Submission/Koch_ref.bib.sav.tmp
Normal file
2220
Frontiers_manuscript/For_Submission/Koch_ref.bib.sav.tmp
Normal file
File diff suppressed because it is too large
Load Diff
BIN
Frontiers_manuscript/For_Submission/Koch_reviewer_comments.pdf
Normal file
BIN
Frontiers_manuscript/For_Submission/Koch_reviewer_comments.pdf
Normal file
Binary file not shown.
Binary file not shown.
115
Frontiers_manuscript/For_Submission/Koch_reviewer_comments.tex
Normal file
115
Frontiers_manuscript/For_Submission/Koch_reviewer_comments.tex
Normal file
@ -0,0 +1,115 @@
|
||||
\documentclass[utf8]{FrontiersinHarvardReview}
|
||||
\DeclareUnicodeCharacter{03B2}{\(\beta\)}
|
||||
\DeclareUnicodeCharacter{03B1}{\(\alpha\)}
|
||||
\DeclareUnicodeCharacter{00C5}{\AA}
|
||||
\usepackage[english]{babel}
|
||||
\usepackage{url,hyperref,lineno,microtype,subcaption}
|
||||
\usepackage[onehalfspacing]{setspace}
|
||||
\usepackage{amsmath}
|
||||
|
||||
\usepackage{xspace}
|
||||
\usepackage{longtable}
|
||||
\usepackage{boldline}
|
||||
\usepackage{geometry}
|
||||
\usepackage{makecell}
|
||||
\usepackage{float}
|
||||
\usepackage{placeins}
|
||||
|
||||
\usepackage[tableposition=top]{caption}
|
||||
|
||||
|
||||
\newcommand{\Kv}{\(\text{K}_{\text{V}}\text{1.1}\)\xspace}
|
||||
\newcommand{\IKv}{\(\text{I}_{\text{K}_{\text{V}}\text{1.1}}\)\xspace}
|
||||
\newcommand{\drheo}{\(\Delta\)rheobase\xspace}
|
||||
\newcommand{\ndAUC}{normalized \(\Delta\)AUC\xspace}
|
||||
|
||||
\newcommand{\note}[2][]{\textbf{[#1: #2]}}
|
||||
\newcommand{\notenk}[1]{\note[NK]{#1}}
|
||||
\newcommand{\notels}[1]{\note[LS]{#1}}
|
||||
\newcommand{\notejb}[1]{\note[JB]{#1}}
|
||||
|
||||
%\linenumbers
|
||||
|
||||
\def\keyFont{\fontsize{8}{11}\helveticabold }
|
||||
\def\firstAuthorLast{Koch {et~al.}} %use et al only if is more than 1 author
|
||||
\def\Authors{Nils A. Koch, Lukas Sonnenberg, Ulrike B.S. Hedrich, Stephan Lauxmann and Jan Benda}
|
||||
%\def\Authors{Nils A. Koch\,$^{1,2}$, Lukas Sonnenberg\,$^{1,2}$, Ulrike B.S. Hedrich\,$^{3}$, Stephan Lauxmann\,$^{1,3}$ and Jan Benda\,$^{1,2,*}$}
|
||||
%\def\Address{$^{1}$Institute for Neurobiology, University of T{\"u}bingen, 72072 T{\"u}bingen, Germany \\
|
||||
%$^{2}$Bernstein Center for Computational Neuroscience T{\"u}bingen, 72076 T{\"u}bingen, Germany \\
|
||||
%$^{3}$Department of Neurology and Epileptology, Hertie Institute for Clinical Brain Research, University of T{\"u}bingen, 72076 T{\"u}bingen, Germany}
|
||||
|
||||
%% The Corresponding Author should be marked with an asterisk
|
||||
%% Provide the exact contact address (this time including street name and city zip code) and email of the corresponding author
|
||||
%\def\corrAuthor{Jan Benda}
|
||||
%
|
||||
%\def\corrEmail{ jan.benda@uni-tuebingen.de}
|
||||
|
||||
|
||||
|
||||
\begin{document}
|
||||
\onecolumn
|
||||
\firstpage{1}
|
||||
|
||||
\title{REVIEW COMMENTS: Loss or Gain of Function? Effects of Ion Channel Mutations on Neuronal Firing Depend on the Neuron Type}
|
||||
|
||||
\author[\firstAuthorLast ]{\Authors} %This field will be automatically populated
|
||||
%\address{} %This field will be automatically populated
|
||||
%\correspondance{} %This field will be automatically populated
|
||||
%\extraAuth{}
|
||||
|
||||
|
||||
\maketitle
|
||||
|
||||
|
||||
\section{Reviewer 1}
|
||||
|
||||
|
||||
\textit{Although the proposed study for publication presents an up-to-date approach, it does not bring any innovation beyond the previous studies (Hedrich et al., 2014; Makinson et al., 2016; Waxman, 2007; Rush et al., 2006).}
|
||||
|
||||
We thank Reviewer 1 for their comments. Although previous studies (Hedrich et al., 2014; Makinson et al., 2016; Waxman, 2007; Rush et al., 2006) referred to by the reviewer and discussed in our manuscript investigate differences in outcomes of ion channel mutations in different neuronal types, these studies have addressed specific instances in which this is the case. For instance, Waxman, 2007 and Rush et al., 2006 demonstrate the dependence of Nav1.7 mutations on the presence/absence of Nav1.8. However whether the outcome of the mutation is graded with different levels of Nav1.8 or whether this effect is modulated by other channels is not investigated nor known. In this context, our work furthers the current knowledge that certain ion channel mutations can have different effects on firing in specific contexts (i.e. cell types) by generalizing this viewpoint to demonstrate that the effects of any ion channel mutation will depend in the neuron-type and its ion channel composition. This suggests that holistic understanding of the effects of ion channel mutations likely requires a more specific neuron-type dependent level of understanding than is generally used currently.
|
||||
|
||||
|
||||
%Previous studies (Hedrich et al., 2014; Makinson et al., 2016; Waxman, 2007; Rush et al., 2006) referred to by the reviewer and discussed in our manuscript investigate differences in outcomes of ion channel mutations in different neuronal types.
|
||||
%This is a vital step in the understanding of the effects of these channelopathies, however the modelling work in our manuscript extends this view that not only the presence or absence of a channel can alter the effects of ion channel mutations, but also the relative expression levels of the mutated channel as well as the context of channels (i.e. the other channels expressed and their relative abundances). \textcolor{red}{The manuscript has been updated to describe this elaboration of the role of other ion channels in the outcomes of ion channel mutations more explicitly}
|
||||
%
|
||||
\textit{The authors say that taking into account neuron type-dependent firing effects offers an opportunity to increase efficacy and precision in personalized medicine approaches, but they do not elaborate on how this applies to the clinic.}
|
||||
|
||||
Currently personalized medicine approaches, although promising, have limited efficacy. Increased efficacy in these approaches likely will result with more therapeutic specificity, however different avenues exist to increase specificity. Although we do not offer a directly actionable clinical outcome, we propose that understanding of the neuron-type level firing effects of ion channel mutations may provide an opportunity for therapeutic strategies that provide more specificity and thus greater efficacy.
|
||||
|
||||
\textit{Considering that even a single mutation in the organism will affect many neuron types (cells) and the signaling pathways to which it is associated, how can all other possible effects be eliminated?}
|
||||
|
||||
We thank the reviewer for their question. The complexity noted in the reviewers question is, in our opinion, vital to the outcome of any ion channel mutation. With this study we demonstrate this to be the case at a neuron-type level and suggest that distillation of this complexity into a LOF or GOF characterization in terms of firing is not meaningful unless tied to a neuron type. Although we make no claim as to at which point certain aspects of the complexity become less or not relevant to the outcome of an ion channel mutation, we use the simulation presented within this study to encourage investigation into neuron type specific effects and to use LOF/GOF firing characterizations with caution.
|
||||
|
||||
\textit{Failure to analyze the effects of the obtained data on living organisms casts a shadow over the reliability of the results. The study would be much more valuable if the authors could at least support their results with animal experiments. Otherwise, it will not go beyond being an evaluation study based on already known mutations. This limits the originality of the study.}
|
||||
|
||||
Although we agree that data obtained from animal experiments is invaluable in the investigation of channelopathies and in neuroscience and neurology generally, the modelling and simulation in this study serves to provide theoretical framework to explain why neuron-type level investigation is important in understanding of ion channel mutations. Furthermore, efforts to reduce animal experimentation are ongoing and modelling such as that presented in this study enables the reduction of animal experimentation. Although we recognize the importance and value in collection of new mutations, further investigation and analysis of known mutations and their effects is also valuable and essential for understanding the effects of ion channel mutations at multiple levels of scale.
|
||||
|
||||
|
||||
\section{Reviewer 2}
|
||||
\textit{The work presented by the authors "Loss or Gain of Function? Effects of Ion Channel Mutations on Neuronal Firing Depend on the Neuron Type" is an approach to elucidate how the loss or gain of function mutation in ionic currents can alter the firing rate of action potentials in different neuronal types. They analyze changes in the rheobase and the area under the curve of the fI curves. Although the approach is adequate to distinguish LOF or GOF, intermediate cases are not easy to analyze. I emphasize that the authors ponder a large part of the limitations of the method used in a broad discussion. However, I believe that the authors should make some clarifications that allow the reader to better understand the information from the results.}
|
||||
|
||||
|
||||
|
||||
\textit{Authors should include a schematic representing the general biological characteristics of each model tested. Even this schematic could be appended to the model description section}
|
||||
|
||||
Todo?
|
||||
|
||||
%We thank Reviewer 2 for their comments.
|
||||
|
||||
\textit{In Figure 1 the authors should comment in more detail on the behavior of the fI curve (B)}
|
||||
|
||||
This model, the RS inhibitory model ceases to fire with large current steps as evidences in the fI curve in Figure 1B. Conceptually, this can be seen in neurons as different processes including depolarization block. \textcolor{red}{This has been update in the manuscript.}
|
||||
|
||||
\textit{Bearing in mind that the authors mention that each model may have hysteresis depending on the stimulation ramp that was applied in figure 1, under what conditions was the rheobase determined in figure 2-5? It is also not clear which model they use to build figure 2A.}
|
||||
|
||||
Although hysteresis is important both from a dynamical systems perspective and in a physiological context, we exclude further analysis of hysteresis in our simulations past Figure 1 to be as applicable as possible to commonly used protocols in the assessment of ion channel mutations. For Figures 2-5, rheobase is calculated by a 2 step process as outlined in the methods. Briefly described this methods uses a sequence of coarse current step amplitudes to identify the step interval in which firing begins. This interval is then explored at high resolution with a second sequence of current step inputs. \textcolor{red}{The manuscript has been updated to explicitly state this.}
|
||||
|
||||
\textit{A final outline summarizing the main results of the main models and the general limitations of this approach could be included.}
|
||||
|
||||
Todo?
|
||||
|
||||
\section{Reviewer 3}
|
||||
|
||||
|
||||
|
||||
\end{document}
|
||||
BIN
Frontiers_manuscript/For_Submission/diversity_in_firing.jpg
Normal file
BIN
Frontiers_manuscript/For_Submission/diversity_in_firing.jpg
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 110 KiB |
BIN
Frontiers_manuscript/For_Submission/diversity_in_firing.pdf
Normal file
BIN
Frontiers_manuscript/For_Submission/diversity_in_firing.pdf
Normal file
Binary file not shown.
Binary file not shown.
|
After Width: | Height: | Size: 125 KiB |
Binary file not shown.
33
Frontiers_manuscript/For_Submission/g_table.tex
Normal file
33
Frontiers_manuscript/For_Submission/g_table.tex
Normal file
@ -0,0 +1,33 @@
|
||||
|
||||
\begin{table}[ht!]
|
||||
% \fontsize{10pt}{10pt}
|
||||
\caption[Cell properties and conductances of neuronal models]{Cell properties and conductances of regular spiking pyramidal neuron (RS Pyramidal; model D), regular spiking inhibitory neuron (RS Inhibitory; model B), fast spiking neuron (FS; model C) each with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (RS Pyramidal +\Kv; model H, RS Inhibitory +\Kv; model E, FS +\Kv; model G respectively), cerebellar stellate cell (Cb Stellate; model A), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (Cb Stellate +\Kv; model F) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_\textrm{A}\) (Cb Stellate \(\Delta\)\Kv; model J), and subthalamic nucleus neuron (STN; model L), with additional \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) (STN +\Kv; model I) and with \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\) replacement of \(\textrm{I}_{\textrm{A}}\) (STN \Kv; model K) models. All conductances are given in \(\textrm{mS}/\textrm{cm}^2\). Capacitances (\(C_m\)) and \(\tau_{max, M}\) are given in pF and ms respectively.}
|
||||
\centering
|
||||
% \begin{tabular}[x]{@{}l@{}} \\ \end{tabular}
|
||||
\linespread{1.5}\selectfont
|
||||
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\resizebox{\textwidth}{!}{\begin{tabular}{cccccccccc}
|
||||
% in mS/cm^2
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& \begin{tabular}[x]{@{}c@{}} RS\\Pyra-\\midal\\(+\Kv) \end{tabular} & \begin{tabular}[x]{@{}c@{}} RS\\Inhib-\\itory\\(+\Kv)\end{tabular} & \begin{tabular}[x]{@{}c@{}}FS\\(+\Kv) \end{tabular}& \begin{tabular}[x]{@{}c@{}} Cb\\Stellate \end{tabular}& \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\+\Kv \end{tabular} & \begin{tabular}[x]{@{}c@{}}Cb\\Stellate\\\(\Delta\)\Kv \end{tabular} & STN &\begin{tabular}[x]{@{}c@{}} STN\\+\Kv \end{tabular} &\begin{tabular}[x]{@{}c@{}} STN\\\(\Delta\)\Kv \end{tabular} \\
|
||||
Model & D (H) & B (E) & C (G) & A & F & J & L & I & K \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(g_{Na}\) & \(56\) & \(10\) & \(58\) & \(3.4\) & \(3.4\) & \(3.4\) & \(49\) & \(49\) & \(49\) \\
|
||||
\(g_{Kd}\) & \(6\) (\(5.4\)) & 2.1 (\(1.89\)) & 3.9 (\(3.51\)) & \(9.0556\) & \(8.15\) &\(9.0556\) & \(57\) & \(56.43\) & \(57\) \\
|
||||
\(g_{K_V1.1}\) & --- (\(0.6\)) & --- (\(0.21\)) & --- (\(0.39\)) & --- & \(0.90556\) & \(1.50159\) & --- & \(0.57\) & \(0.5\) \\
|
||||
\(g_{A}\) & --- & --- & --- & \(15.0159\) & \(15.0159\) & --- & \(5\) & \(5\) & --- \\
|
||||
\(g_{M}\) & \(0.075\) & \(0.0098\) &\(0.075\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(g_{L}\) & --- & --- & --- & --- & --- & --- & \(5\) & \(5 \) & \(5\) \\
|
||||
\(g_{T}\) & --- & --- & --- & \(0.45045\) & \(0.45045\) & \(0.45045\) & \(5\) & \(5\) & \(5\) \\
|
||||
\(g_{Ca,K}\) & --- & --- & --- & --- & --- & --- & \(1\) & \(1\) & \(1\) \\
|
||||
\(g_{Leak}\) &\(0.0205\) & \(0.0205\) & \(0.038\) & \(0.07407\) & \(0.07407\) & \(0.07407\) & \(0.035\) & \(0.035\) & \(0.035\) \\%\Xhline{1\arrayrulewidth}
|
||||
\(\tau_{max, M}\)& \(608\) & \(934\) & \(502\) & --- & --- & --- & --- & --- & --- \\
|
||||
\(C_m\) & \(118.44\) & \(119.99\) & \(101.71\)& \(177.83\) & \(177.83\) & \(177.83\) & \(118.44\)& \(118.44\)& \(118.44\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\end{tabular}}}
|
||||
|
||||
|
||||
\label{tab:g}
|
||||
\end{table}
|
||||
|
||||
28
Frontiers_manuscript/For_Submission/gating_table.tex
Normal file
28
Frontiers_manuscript/For_Submission/gating_table.tex
Normal file
@ -0,0 +1,28 @@
|
||||
\begin{table}[ht!]
|
||||
\caption[Gating Properties]{ For comparability to typical electrophysiological data fitting reported and for ease of further gating curve manipulations, a sigmoid function (Eqn.\ref{eqn:Boltz}) %Boltzmann \(x_\infty = {\left(\frac{1-a}{1+{exp[{\frac{V-V_{1/2}}{k}}]}} +a\right)^j}\)
|
||||
with slope \(k\), voltage for half-maximal activation or inactivation (\(V_{1/2}\)), exponent \(j\), and persistent current \(0 \leq a \leq 1\) were fitted for the models originating from \citet{pospischil_minimal_2008} (models B, C, D, E, G, H) where \(\alpha_x\) and \(\beta_x\) are used. Gating parameters for \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) are taken from \citet{ranjan_kinetic_2019} and fit to mean wild type parameters in \citet{lauxmann_therapeutic_2021}. Model gating parameters not listed are taken directly from source publication.}
|
||||
\centering
|
||||
\linespread{1.5}\selectfont
|
||||
\fontsize{10pt}{12pt}\selectfont{
|
||||
\begin{tabular}{c c c c c c }
|
||||
\Xhline{1\arrayrulewidth}
|
||||
& Gating & \(V_{1/2}\) [mV]& \(k\) & \(j\) & \(a\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
%Pospischil
|
||||
& \(\textrm{I}_{\textrm{Na}}\) activation &\(-34.33054521\) & \(-8.21450277\) & \(1.42295686\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{Na}}\) inactivation &\(-34.51951036\) & \(4.04059373\) & \(1\) & \(0.05\) \\
|
||||
Models & \(\textrm{I}_{\textrm{Kd}}\) activation &\(-63.76096946\) & \(-13.83488194\) & \(7.35347425\) & --- \\
|
||||
B, C, D, E, G, H & \(\textrm{I}_{\textrm{L}}\) activation &\(-39.03684525\) & \(-5.57756176\) & \(2.25190197\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{L}}\) inactivation &\(-57.37\) & \(20.98\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{M}}\) activation &\(-45\) & \(-9.9998807337\) & \(1\) & --- \\ %-45 with 10 mV shift to contributes to resting potential
|
||||
% & & & & &\\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
\(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) & \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) activation &\(-30.01851852\) & \(-7.73333333\) & \(1\) & --- \\
|
||||
& \(\textrm{I}_{\textrm{K}_{\textrm{V}}\textrm{1.1}}\ \) Inactivation &\(-46.85851852\) & \(7.67266667\) & \(1\) & \(0.245\) \\
|
||||
\Xhline{1\arrayrulewidth}
|
||||
|
||||
|
||||
\end{tabular}}
|
||||
|
||||
\label{tab:gating}
|
||||
\end{table}
|
||||
BIN
Frontiers_manuscript/For_Submission/ramp_firing.jpg
Normal file
BIN
Frontiers_manuscript/For_Submission/ramp_firing.jpg
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 82 KiB |
BIN
Frontiers_manuscript/For_Submission/rheobase_correlation.jpg
Normal file
BIN
Frontiers_manuscript/For_Submission/rheobase_correlation.jpg
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 424 KiB |
BIN
Frontiers_manuscript/For_Submission/rheobase_correlation.pdf
Normal file
BIN
Frontiers_manuscript/For_Submission/rheobase_correlation.pdf
Normal file
Binary file not shown.
Binary file not shown.
|
After Width: | Height: | Size: 436 KiB |
Binary file not shown.
Loading…
Reference in New Issue
Block a user