In [1]:
from bussilab import lohman
import matplotlib.pyplot as plt
import numpy as np
In [2]:
print(lohman.lohman.__doc__)
Lohman model for helicases.
Compute the fraction of unwound helices as a function of time.
See [Lucius et al, Biophys J 2003](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1303449/).
Parameters
----------
t: float or sequence or np.ndarray
Time. If a sequence or np.ndarray is provided, the function is computed
for all values and an array is returned.
ku: float
Unwinding rate.
kd: float
Dissociation rate.
n: int, optional
Step size
boundaries: tuple with 2 elements
Result is mapped to this range.
Returns
-------
float or np.ndarray
The fraction of unfolded helices at a given time `t`. If an array is provided for `t`,
an array is returned containing the fractions at all the times.
If `boundaries` is provided, the fraction is linearly mapped into the
`boundaries[0], boundaries[1]` range.
In [3]:
t=np.linspace(0,4)
plt.plot(t,lohman.lohman(t,ku=1,kd=0.1),label="k_U=1, k_D=0.1, N=1")
plt.plot(t,lohman.lohman(t,ku=2,kd=0.1),label="k_U=2, k_D=0.1, N=1")
plt.plot(t,lohman.lohman(t,ku=1,kd=0.2),label="k_U=1, k_D=0.2, N=1")
plt.plot(t,lohman.lohman(t,ku=1,kd=0.1,n=2),label="k_U=1, k_D=0.1, N=2")
plt.plot(t,lohman.lohman(t,ku=1,kd=0.1,boundaries=(0.3,0.7)),
label="k_U=1, k_D=0.1, N=1, [0.3,0.7]")
plt.legend()
plt.show()
In [ ]: