Quiz 4: Computing in Carbon Help Center

Learn more

Warning: The hard deadline has passed. You can attempt it, but you will not get credit for it. You are welcome to try it as a learning exercise.

Question 1

When we talk about "spikes", we are referring to the change in some property of the neuron over time. When we typically plot a spike, the x-axis represents time. What does the y-axis represent?

Question 2



In this circuit diagram representing a piece of neuronal cell membrane, the battery, resistor, and capacitor are roughly analogous to the ___, ___, and ___ respectively.

Question 3

Let's imagine there is another ion that is relevant to determining a neuron's membrane potential in addition to those discussed in the lecture. We'll call the ion Im+ (for Imaginary). The equilibrium potential of Im+ (EIm) is -100 mV. Assume the resting potential of the neuron is -65 mV. When specialized Im+ channels open, the Im+ conductance will increase. This will ___ the cell, thus ___ its membrane potential.

Question 4

Suppose Im+ channels are composed of 5 subunits that open and close independently, as well as an additional "ball-in-socket" gating mechanism. Each of the 5 subunits has a voltage-dependent open probability u and closed probability (1−u), while the ball-in-socket gating mechanism has a voltage-dependent open probability z and closed probability (1−z). Which expression could most likely be used to express the Im+ current across the membrane?

Question 5



Let's come back to the real-life Na+ (sodium) channel, whose voltage-dependent dynamics are shown above. m∞ and h∞ (the steady states of m and h) are shown for different V, along with the associated voltage-dependent time constants τm and τh. Remember that m is the probability that any individual channel subunit is open and h is the probability that the additional ball-in-socket gating mechanism is open. If the neuron is at its resting potential (around -65 mV), and we deliver a current injection that decreases the membrane voltage to -75 mV, which of the following will be the first to reach its new steady state?

Question 6

Refer again to the figure shown above the previous question. Remember that Na+ current depolarizes the cell and is the principal driver for the upward portion of a spike. Both m and h must be high for there to be a lot of Na+ current. h∞ becomes 0 when voltage is close to or greater than -30 mV. How, then, is it possible for the membrane to depolarize beyond V = -30 mV during a spike (spikes peak closer to V = 40 mV)?

Question 7

True or false: All neural coding can essentially be reduced to variations in firing rate. Thus, precise spike timing is usually unimportant: what matters is a "rate code."

Question 8



Recall the exponential integrate and fire neuron model, schematized above. How many stable fixed points does the system have?

Question 9

QUESTIONS 9 AND 10 ARE OPTIONAL AND ARE WORTH ZERO POINTS
The FitzHugh-Nagumo model is a 2-dimensional dynamical neuron model. It is defined by the following two differential equations:

dVdt=V(a−V)(V−1)−w+I

dwdt=bV−cw

where V is voltage, w represents an outward hyperpolarizing current, I is injected current, and a, b, and c are constants. Which of the following is an expression for the "w-nullcline?"

Question 10

QUESTIONS 9 AND 10 ARE OPTIONAL AND ARE WORTH ZERO POINTS


The above figure is a phase plane based on the FitzHugh-Nagumo model from the previous question. A vector field is shown that gives a sense of the flow of the system. If we observe a spike in this system, the trajectory will travel counterclockwise around the phase plane. The first part, or “upstroke,” of the spike occurs in which of the following regions of the phase plane (regions are labeled in the figure)?

Question 11

The next five questions utilize the following code to model a passive neuronal membrane as an RC-circuit. (Remember that in the membrane model, the resistor and capacitor are in parallel.)

MATLAB: membrane.m
Python (all versions): membrane.py

This code demonstrates how a membrane responds to a constant current input that is turned on for a fixed time interval and then turned off. What if the current were not turned off? What would the steady state voltage of the membrane be?

Use the values given in the script to compute your answer (C = 0.1 nF, R = 100 MΩ, I = 10 nA). You should give your answer in mV. Do not include units in your answer.

Question 12

Change the values for the membrane's resistance and capacitance (R and C), and find out how this influences the response of the membrane. Does it reach a stable value more quickly or more slowly after multiplying R by 5?

Question 13

Does it reach a stable value more quickly or more slowly after dividing C by 10?

Question 14

Does it reach a stable value more quickly or more slowly after multiplying R by 10 AND dividing C by 10?

Question 15

An experimental method for calculating a membrane’s time constant τ (when R and/or C are not known) is to start at zero and record the time at which the membrane potential V reaches a value approximately equal to 0.6321∗Vpeak=0.6321∗IR, where I is the constant injected current. Check if this method works by injecting different amounts of current I and changing the values for R and C. Once you’ve convinced yourself that the experimental τ appears to be identical to the theoretical τ(=RC) in all these cases, provide a theoretical justification for why this method works.

To do this, find the solution to the differential equation for the passive membrane:
dVdt=−VRC+IC

V(0)=0


After solving the differential equation you should be able to use the fact that Vpeak=IR and e−1e=.6321 to complete the derivation and show that V(τ)=0.6321∗IR.

Which of these equations is the solution to the given differential equation?

Question 16

In the next three questions, you will explore the integrate-and-fire neuron model. To do this, you should start by downloading the following code and tweaking it to run "experiments" on the neuron:

MATLAB: intfire.m
Python: intfire.py

What is the largest current that will fail to cause the neuron to spike? Give your answer in pA and round down to the nearest 10 pA. Do not include units in your answer.

You should vary the input current gradually from very low to high values to find this value and then validate your answer with an analytical solution.

Question 17

What is the maximum firing rate (spike count/trial duration) of this neuron?

Give your answer in Hz and round your answer to the nearest integer value. Do not include the units in your answer.

Question 18

Now let's consider the case that your neuron is not receiving simply a constant input, but a barrage of asynchronous inputs from many other neurons. Model this with the following code by adding a white noise component to the input current (the constant part of the input current is reset to one nA):

MATLAB: intfireNoise.m
Python: intfire_noise.py

Plot the interspike interval distribution, that is, the distribution of the time intervals between consecutive spikes, for a range of different noise amplitudes. Hints: this noise component is already implemented in the code: just switch it on. You can make use of the MATLAB functions diff and hist (in Python np.diff and plt.hist). You will probably want to increase tstop, the length of time you are integrating for, by quite a bit to get a well-sampled histogram. What best describes how this distribution changes as you increase the amplitude of the noise input? (Stay within a range between 0 and 5 nA.)
    
You cannot submit your work until you agree to the Honor Code. Thanks!