Skip to content

Commit

Permalink
chore(ray2020): tweak steps
Browse files Browse the repository at this point in the history
  • Loading branch information
sanjayankur31 committed Apr 16, 2024
1 parent a7f22ed commit 4187553
Showing 1 changed file with 9 additions and 9 deletions.
18 changes: 9 additions & 9 deletions source/Userdocs/Walkthroughs/RayEtAl2020/Conversion.md
Original file line number Diff line number Diff line change
Expand Up @@ -5,8 +5,8 @@ On inspection of the model, we see that it has two biophysically detailed cell m
- the GGN (Giant GABAergic Neuron)
- the KC (Kenyon cell)

## Step 1) Converting the Giant GABAergic Neuron
### 1a) Exporting morphology of the GGN
## Converting the Giant GABAergic Neuron
### Step 1) Exporting morphology of the GGN

Let's start with the GGN first.
It's morphology is defined as an SWC file in [this file](https://github.com/OpenSourceBrain/262670/blob/master/NEURON/mb/cell_templates/GGN_20170309_sc.swc).
Expand Down Expand Up @@ -51,7 +51,7 @@ pynml-plotmorph -i GGN.morph.cell.nml
Visualisation of the GGN using `pynml-plotmorph`
```

### 1b) Adding biophysics to the GGN
### Step 2) Adding biophysics to the GGN

Now that we have the morphology of the GGN exported, we can add the biophysics.
We need to inspect the original model code to learn about the biophysics.
Expand Down Expand Up @@ -160,7 +160,7 @@ Once this is set up, we write the cell to a new file.
The GGN cell has now been converted.
Since the GGN is a simple passive cell, we won't test its biophysics just yet.

## Step 2) Converting the Kenyon Cell
## Converting the Kenyon Cell

The KC cell is defined in the `kc_1_comp.hoc` file.
Whereas the GGN cell had a complex morphology but passive biophysics, the KC cell has very simple morphology---a single compartment---but does contain active channels:
Expand Down Expand Up @@ -206,7 +206,7 @@ proc biophys() {
```

### 2a) Converting ion channels
### Step 1) Converting ion channels

For this cell, we will first convert the various ion channel models.
These are included in the `mod` folder.
Expand Down Expand Up @@ -433,7 +433,7 @@ We have two cases here:
- if `v` is greater than `table_max` (which is 40mV), the rate is the value at `v=40mV`
- otherwise, the rate is calculated from the value of `v`

#### i) Kv, Naf, Kst
#### 1a) Kv, Naf, Kst

The nas channel is now done.
If we look at the other channels---naf, kv, and kst---they follow similar formalisms.
Expand Down Expand Up @@ -485,7 +485,7 @@ Time course of activation variables of kv channel, generated with pynml-channela

In these graphs, the effect of the conditional at 40mV becomes more apparent.

#### ii) Ka
#### 1b) Ka

The last remaining channel is the ka channel.
It's dynamics are defined in the mod file as:
Expand Down Expand Up @@ -565,7 +565,7 @@ Even though the [NMODL language](https://nrn.readthedocs.io/en/8.2.2/python/mode
This means that different people can write the same dynamics in different ways.
On the other hand, NeuroML and LEMS are more formal with more strict structures, and once channels are converted to these formats, they are much easier to understand.

### 2b) Creating the morphology
### Step 2) Creating the morphology

Since the morphology of the KC is a single compartment, we don't need to export it from NEURON.
We can create it ourselves.
Expand Down Expand Up @@ -672,7 +672,7 @@ Now that we have the morphology and the ion channels for the KC, we can add the
This completes the cell.
We export it to a NeuroML file.

### 2c) Testing the model
### Step 3) Testing the model

Finally, we want to test our NeuroML conversion against the original cell to see that it exhibits the same dynamics.
The `test_kc.py` script runs a simple step current simulation with a single KC cell and shows its membrane potentials:
Expand Down

0 comments on commit 4187553

Please sign in to comment.