This is the official repository of Propensity Optimised Paths (POPs) , a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. The method was presented in the paper:
Wu, N., Yaliraki, S. N., & Barahona, M. (2022). Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity. Journal of Molecular Biology, 167749. https://doi.org/10.1016/j.jmb.2022.167749
The calculation of POPs is based on results from bond-to-bond propensity analysis, which was introduced in Amor et al. Nature Communications, 7, 1–13 (2016) and made available through the easy-to-use webserver ProteinLens.
Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths.
This code is only available in python. We apologise that it is not available in other languages for users who are not familiar with python and would like to use this code.
Bond-to-bond propensity analysis quantifies the non-local effect of instantaneous bond fluctuations propagating through the protein. It is available as a webserver - ProteinLens. Please follow the Tutorial provided in ProteinLens to complete bond-to-bond propensity analysis on your chosen protein and download the results. The results are required for POPs computation and an example of the result folder, sessionIZN1I the bond-to-bond propensity analysis result of h-Ras (PDB ID: 3K8Y), is provided here for illustration.
After completing bond-to-bond propensity analysis, POP.py is the only script required for POPs computation and analysis. Please follow the step by step instruction in POP_example.ipynb in the examples/ directory to complete POPs computation and analysis. Running everything in POP_example.ipynb would give you two folders as the results and they should be the same as the folders given in examples_completed in the examples/ directory. This will make sure everything is set up correctly and you can then proceed with your own protein(s) of interest.
- Nan Wu, GitHub:
nw97nan <https://github.com/nw97nan>
If you use this code in your own work, please cite our paper :
Wu, N., Yaliraki, S. N., & Barahona, M. (2022). Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity. Journal of Molecular Biology, 167749. https://doi.org/10.1016/j.jmb.2022.167749
and the paper on bond-to-bond propensity analysis:
Amor, B. R. C., Schaub, M. T., Yaliraki, S. N., & Barahona, M. (2016). Prediction of allosteric sites and mediating interactions through bond-to-bond propensities. Nature Communications, 7, 1–13. https://doi.org/10.1038/ncomms12477.
The paper that introduced the ProteinLens webserver is:
Mersmann, S. F., Strömich, L., Song, F. J., Wu, N., Vianello, F., Barahona, M., & Yaliraki, S. N. (2021). ProteinLens: a web-based application for the analysis of allosteric signalling on atomistic graphs of biomolecules. Nucleic Acids Research, 49(W1), W551–W558. https://doi.org/10.1093/nar/gkab350