Open probability and dose response curves

A collection of functions for open probability, Popen, or dose-response curve calculations.

popen.EC50(mec, tres, eff='c')

Estimate numerically the equilibrium EC50 for a specified mechanism. If monotonic this is unambiguous. If not monotonic then returned is a concentration for 50% of the peak response to the left of the peak.

Parameters :

mec : instance of type Mechanism

tres : float

Time resolution (dead time).

Returns :

EC50 : float

Concentration at which Popen is 50% of its maximal value.

popen.Popen(mec, tres, conc, eff='c')

Calculate equilibrium open probability (Popen) and correct for unresolved blockages in case of presence of fast pore blocker.

Parameters :

mec : dcpyps.Mechanism

The mechanism to be analysed.

tres : float

Time resolution (dead time).

conc : float

Concentration.

Returns :

Popen : float

Open probability value at a given concentration.

popen.Popen0(mec, tres, eff='c')

Find Popen at concentration = 0.

Parameters :

mec : instance of type Mechanism

tres : float

Time resolution (dead time).

Returns :

P0 : float

Open probability in absence of effector.

popen.decline(mec, tres, eff='c')

Find whether open probability curve increases or decreases with ligand concentration. Popen may decrease if ligand is inhibitor.

Parameters :

mec : instance of type Mechanism

tres : float

Time resolution (dead time).

Returns :

decline : bool

True if Popen curve dectreases with concentration.

popen.maxPopen(mec, tres, eff='c')

Estimate numerically maximum equilibrium open probability. In case Popen curve goes through a maximum, the peak open probability is returned.

Parameters :

mec : instance of type Mechanism

tres : float

Time resolution (dead time).

Returns :

maxPopen : float

Maximum equilibrium open probability.

conc : float

Concentration at which Popen curve reaches maximal value.

popen.nH(mec, tres, eff='c')

Calculate Hill slope, nH, at EC50 of a calculated Popen curve. This is Python implementation of DCPROGS HJC_HILL.FOR subroutine.

Parameters :

mec : instance of type Mechanism

tres : float

Time resolution (dead time).

Returns :

nH : float

Hill slope.

popen.printout(mec, tres, output=<open file '<stdout>', mode 'w' at 0x2b971dd501e0>, eff='c')