pyLoopSage 0.0.8__tar.gz → 0.1.0__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/PKG-INFO +5 -5
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/README.md +4 -4
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/PKG-INFO +5 -5
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/setup.py +1 -1
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/LICENSE +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/SOURCES.txt +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/dependency_links.txt +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/entry_points.txt +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/requires.txt +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/pyLoopSage.egg-info/top_level.txt +0 -0
- {pyloopsage-0.0.8 → pyloopsage-0.1.0}/setup.cfg +0 -0
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: pyLoopSage
|
|
3
|
-
Version: 0.0
|
|
3
|
+
Version: 0.1.0
|
|
4
4
|
Summary: An energy-based stochastic model of loop extrusion in chromatin.
|
|
5
5
|
Home-page: https://github.com/SFGLab/pyLoopSage
|
|
6
6
|
Author: Sebastian Korsak
|
|
@@ -51,14 +51,14 @@ Please cite the method paper in case that you would like to use this model for y
|
|
|
51
51
|
|
|
52
52
|
## The model
|
|
53
53
|
|
|
54
|
-
We have a polymer chain with
|
|
54
|
+
We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
|
|
55
55
|
|
|
56
56
|
Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
|
|
57
57
|
|
|
58
58
|
* Slide both locations randomly (as 1D random walk) or
|
|
59
59
|
* Rebind somewhere else.
|
|
60
60
|
|
|
61
|
-
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is
|
|
61
|
+
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
|
|
62
62
|
|
|
63
63
|
The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
|
|
64
64
|
|
|
@@ -71,7 +71,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
|
|
|
71
71
|
|
|
72
72
|
$$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
|
|
73
73
|
|
|
74
|
-
The parameters are defined in such a way that when
|
|
74
|
+
The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
|
|
75
75
|
|
|
76
76
|
And the energy difference can be expressed as the energy difference of each term,
|
|
77
77
|
|
|
@@ -96,7 +96,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
|
|
|
96
96
|
|
|
97
97
|
```python
|
|
98
98
|
# Definition of Monte Carlo parameters
|
|
99
|
-
import
|
|
99
|
+
import loopsage.stochastic_simulation as lps
|
|
100
100
|
|
|
101
101
|
N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
|
|
102
102
|
mode = 'Metropolis'
|
|
@@ -17,14 +17,14 @@ Please cite the method paper in case that you would like to use this model for y
|
|
|
17
17
|
|
|
18
18
|
## The model
|
|
19
19
|
|
|
20
|
-
We have a polymer chain with
|
|
20
|
+
We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
|
|
21
21
|
|
|
22
22
|
Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
|
|
23
23
|
|
|
24
24
|
* Slide both locations randomly (as 1D random walk) or
|
|
25
25
|
* Rebind somewhere else.
|
|
26
26
|
|
|
27
|
-
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is
|
|
27
|
+
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
|
|
28
28
|
|
|
29
29
|
The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
|
|
30
30
|
|
|
@@ -37,7 +37,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
|
|
|
37
37
|
|
|
38
38
|
$$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
|
|
39
39
|
|
|
40
|
-
The parameters are defined in such a way that when
|
|
40
|
+
The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
|
|
41
41
|
|
|
42
42
|
And the energy difference can be expressed as the energy difference of each term,
|
|
43
43
|
|
|
@@ -62,7 +62,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
|
|
|
62
62
|
|
|
63
63
|
```python
|
|
64
64
|
# Definition of Monte Carlo parameters
|
|
65
|
-
import
|
|
65
|
+
import loopsage.stochastic_simulation as lps
|
|
66
66
|
|
|
67
67
|
N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
|
|
68
68
|
mode = 'Metropolis'
|
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: pyLoopSage
|
|
3
|
-
Version: 0.0
|
|
3
|
+
Version: 0.1.0
|
|
4
4
|
Summary: An energy-based stochastic model of loop extrusion in chromatin.
|
|
5
5
|
Home-page: https://github.com/SFGLab/pyLoopSage
|
|
6
6
|
Author: Sebastian Korsak
|
|
@@ -51,14 +51,14 @@ Please cite the method paper in case that you would like to use this model for y
|
|
|
51
51
|
|
|
52
52
|
## The model
|
|
53
53
|
|
|
54
|
-
We have a polymer chain with
|
|
54
|
+
We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
|
|
55
55
|
|
|
56
56
|
Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
|
|
57
57
|
|
|
58
58
|
* Slide both locations randomly (as 1D random walk) or
|
|
59
59
|
* Rebind somewhere else.
|
|
60
60
|
|
|
61
|
-
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is
|
|
61
|
+
In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
|
|
62
62
|
|
|
63
63
|
The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
|
|
64
64
|
|
|
@@ -71,7 +71,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
|
|
|
71
71
|
|
|
72
72
|
$$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
|
|
73
73
|
|
|
74
|
-
The parameters are defined in such a way that when
|
|
74
|
+
The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
|
|
75
75
|
|
|
76
76
|
And the energy difference can be expressed as the energy difference of each term,
|
|
77
77
|
|
|
@@ -96,7 +96,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
|
|
|
96
96
|
|
|
97
97
|
```python
|
|
98
98
|
# Definition of Monte Carlo parameters
|
|
99
|
-
import
|
|
99
|
+
import loopsage.stochastic_simulation as lps
|
|
100
100
|
|
|
101
101
|
N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
|
|
102
102
|
mode = 'Metropolis'
|
|
@@ -6,7 +6,7 @@ with open('README.md', 'r', encoding='utf-8') as f:
|
|
|
6
6
|
|
|
7
7
|
setup(
|
|
8
8
|
name='pyLoopSage', # Package name
|
|
9
|
-
version='0.0
|
|
9
|
+
version='0.1.0', # Version of the software
|
|
10
10
|
description='An energy-based stochastic model of loop extrusion in chromatin.',
|
|
11
11
|
long_description=long_description,
|
|
12
12
|
long_description_content_type='text/markdown',
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|