ChromaWalker is a Python package for exploring chromatin architecture via the Markov State Model analysis of Hi-C data as delineated in
![image](https://private-user-images.githubusercontent.com/8078280/321658417-1384a3e1-2079-41d6-baae-c9fb7d781e15.png?jwt=eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJpc3MiOiJnaXRodWIuY29tIiwiYXVkIjoicmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbSIsImtleSI6ImtleTUiLCJleHAiOjE3MjQ4MDk3NDksIm5iZiI6MTcyNDgwOTQ0OSwicGF0aCI6Ii84MDc4MjgwLzMyMTY1ODQxNy0xMzg0YTNlMS0yMDc5LTQxZDYtYmFhZS1jOWZiN2Q3ODFlMTUucG5nP1gtQW16LUFsZ29yaXRobT1BV1M0LUhNQUMtU0hBMjU2JlgtQW16LUNyZWRlbnRpYWw9QUtJQVZDT0RZTFNBNTNQUUs0WkElMkYyMDI0MDgyOCUyRnVzLWVhc3QtMSUyRnMzJTJGYXdzNF9yZXF1ZXN0JlgtQW16LURhdGU9MjAyNDA4MjhUMDE0NDA5WiZYLUFtei1FeHBpcmVzPTMwMCZYLUFtei1TaWduYXR1cmU9ZDcwMjhkOGMxYTFmOWFmMWJmMmNiMzk1NjFlZjZkOTY0NjM0MDVhNWEwODM4NjMyOGI4OWIzMmQ0YmJkYTc3NCZYLUFtei1TaWduZWRIZWFkZXJzPWhvc3QmYWN0b3JfaWQ9MCZrZXlfaWQ9MCZyZXBvX2lkPTAifQ.7pXCQ5IlRfj1AnUCYIKwpQhCPjPD5AJ-EMm1rlSVmFQ)
Tan, Z. W.; Guarnera, E.; Berezovsky, I. N. Exploring Chromatin Hierarchical Organization via Markov State Modelling. PLoS Computational Biology 2018, 14 (12). https://doi.org/10.1371/journal.pcbi.1006686.
In this work, we proposed a new approach for extracting robust genomic partitions from HiCdata, seeking to capture the hallmarks of chromatin structure and organization by considering the entire interaction landscape of this complex system. The objectives are to identify and study structural features of chromatin from Hi-C interaction data and to find a connection between these features and data on epigenetic regulation. We introduced a Markov State Model (MSM) approach with minimal assumptions and parameters on the chromatin interaction network, aiming at identifing structural partitions and their interactions. By analogy with a biomolecule moving and interacting in condensed chromatin, the MSM allows one to explore chromatin structure using a “probe” randomly walking in the contact energy landscape derived from Hi-C data. Given the multiscale nature of the data-derived contact energy landscape and the metastability of the corresponding MSM, we can identify regions of dense intra- and inter-chromosomal interactions, linkers between these regions, as well as the overall topology of individual chromosomes and the complex structures that chromosomes form by interacting with each other. The MSM approach and the metastability analisys of the contact energy landscape derived from Hi-C data are a development of a previous work in the context of protein dynamics (see Guarnera, E.; Vanden-Eijnden, E. Optimized Markov State Models for Metastable Systems. Journal of Chemical Physics 2016, 145 (2). https://doi.org/10.1063/1.4954769).
These instructions will get you a copy of the project up and running on your local machine for development and testing purposes. Deployment on a server should require just the same steps.
This library is built for Python 2.7.
List of required Python libraries (versions):
- Numpy (1.14)
- Scipy (1.0)
- Pandas (0.22)
- Matplotlib (2.0)
- Install the SciPy stack
The most convenient way is to use pip:
python -m pip install --user numpy scipy matplotlib pandas
- Copy package source to Python path
- Basic usage of the package can be inferred by looking at the code that's run when you execute ChromaWalker.py as a main program. In other words, look at how the program is run in ChromaWalker.py, under the conditional statement:
if __name__ == “__main__”:
MIT license