<!-- <picture>
<source media="(prefers-color-scheme: dark)" srcset="./assets/repo_logo_dark.png" width='100%'>
<source media="(prefers-color-scheme: light)" srcset="./assets/repo_logo_light.png" width='100%'>
<img alt="Project logo" src="/assets/" width="100%">
</picture> -->
<br>
[](https://github.com/schwallergroup/CycPepPerm)
<!-- [](https://doi.org/10.48550/arXiv.2304.05376)
[](https://img.shields.io/pypi/v/CycPepPerm)
[](https://img.shields.io/pypi/pyversions/CycPepPerm) -->
[](https://github.com/psf/black)
[](https://github.com/schwallergroup/liac-repo)
[](https://schwallergroup.github.io)
<h1 align="center">
CycPepPerm
</h1>
<br>
Python package to predict membrane permeability of cyclic peptides.
## 👩💻 Installation
We provide the code as a [python package](https://pypi.org/project/cyc-pep-perm/), so the only thing you need is to install it. We recommend creating a new conda environment for that, which allows simple package management for a project. Follow these [instructions](https://docs.anaconda.com/free/anaconda/install/index.html) to install Anaconda. However, the package containing our code can also be installed without creating a project-specific environment. In that case, one just skips the first two lines of the following code:
```bash
$ conda create -n cyc-pep-perm python=3.10
$ conda activate cyc-pep-perm
$ pip install cyc-pep-perm==0.1.0
```
### 🛠️ For Developers
<details>
<summary>See detailed installation instructions</summary>
The repository can be cloned from GitHub and installed with `pip` or `conda`. The code was built with Python 3.10 on Linux but other OS should work as well.
With conda:
```bash
$ git clone git+https://github.com/schwallergroup/CycPepPerm.git
$ cd CycPepPerm
$ conda env create -f environment.yml
$ conda activate cyc_pep_perm
$ pip install -e .
```
or with pip:
```bash
$ git clone git+https://github.com/schwallergroup/CycPepPerm.git
$ cd CycPepPerm
$ conda create -n cyc_pep_perm python=3.10
$ conda activate cyc_pep_perm
$ pip install -r requirements.txt
$ pip install -e .
```
If the options above did not work, please try from scratch:
```bash
$ git clone git+https://github.com/schwallergroup/CycPepPerm.git
$ cd CycPepPerm
$ conda create -c conda-forge -n cyc_pep_perm rdkit=2022.03.5 python=3.10
$ conda activate cyc_pep_perm
$ conda install -c conda-forge scikit-learn=1.0.2
$ conda install -c rdkit -c mordred-descriptor mordred
$ conda install -c conda-forge xgboost
$ conda install -c conda-forge seaborn
$ pip install shap
$ conda install -c conda-forge jupyterlab
$ pip isntall pandas-ods-reader
$ pip install -e .
```
</details>
## 🔥 Usage
> For some more examples on how to process data, train and evaluate the alogrithms, please consult the folder `notebooks/`. This folder also contains a notebook to perform polynomial fits as described in the paper.
All data paths in the following examples are taken from the hard-coded paths that work when one clones this repository. If you use the python package and download the data separately, please change the paths accordingly.
### Data preprocessing
Here we showcase how to handle the data as for our use-case. Some simple reformating is done (see also the notebook `notebooks/01_data_preparation.ipynb`) starting from `.ods` file with DataWarrior output (for data see [Data and Models](#data-and-models)).
```python
import os
from cyc_pep_perm.data.processing import DataProcessing
data_dir = "/path/to/data/folder" # ADAPT TO YOUR PATH!
# this can also be a .csv input
datapath = os.path.join(data_dir, "perm_random80_train_raw.ods")
# instantiate the class and make sure the columns match your inputed file - otherwise change arguments
dp = DataProcessing(datapath=datapath)
# make use of precomputed descriptors from DataWarrior
df = dp.read_data(filename="perm_random80_train_dw.csv")
# calculate Mordred deescripttors
df_mordred = dp.calc_mordred(filename="perm_random80_train_mordred.csv")
```
### Training
Make sure to have the data ready to be used. In order to make the hyperparameter search more extensive, please look into the respective python scripts (e.g. `src/cyc_pep_perm/models/randomforest.py`) and adjust the `PARAMS` dictionary.
```python
import os
from cyc_pep_perm.models.randomforest import RF
data_dir = "/path/to/data/folder" # ADAPT TO YOUR PATH!
train_data = os.path.join(data_dir, "perm_random80_train_dw.csv")
model_dir = "/path/to/model/folder" # ADAPT TO YOUR PATH!
rf_model_trained = os.path.join(model_dir, "rf_random_dw.pkl")
# instantiate class
rf_regressor = RF()
model = rf_regressor.train(
datapath = train_data,
savepath = rf_model_trained,
)
y_pred, rmse, r2 = rf_regressor.evaluate()
# will print training results, e.g.:
>>> RMSE: 8.45
>>> R2: 0.879
```
### Prediction
```python
import os
from cyc_pep_perm.models.randomforest import RF
data_dir = "/path/to/data/folder" # ADAPT TO YOUR PATH!
train_data = os.path.join(data_dir, "perm_random80_train_dw.csv")
model_dir = "/path/to/model/folder" # ADAPT TO YOUR PATH!
rf_model_trained = os.path.join(model_dir, "rf_random_dw.pkl")
# instantiate class
rf_regressor = RF()
# load trained model
rf_regressor.load_model(
modelpath = rf_model_trained,
)
# data to predict on, e.g.:
df = pd.read_csv(train_data)
X = df.drop(columns=["SMILES"])
# predict
y_pred = rf_regressor.predict(X)
```
## Data and Models
All data required for reproducing the results in the paper are provided in the folder `data/`. Beware that due to the random nature of these models, the results might differ from the ones reported in the paper. The files found in `data/` are split into training and test data (randomly split 80/20) and with either the DataWarrior (dw) or the Mordred descriptors. The simple data processing can be found in the notebook `notebooks/01_data_preparation.ipynb`. The DataWarrior descriptors are computed with external software ([DataWarrior](https://openmolecules.org/datawarrior/)). The following files are provided:
- `data/perm_random20_test_dw.csv` - test data with DataWarrior descriptors
- `data/perm_random20_test_mordred.csv` - test data with Mordred descriptors
- `data/perm_random20_test_raw.ods` - test data before processing
- `data/perm_random80_train_dw.csv` - training data with DataWarrior descriptors
- `data/perm_random80_train_mordred.csv` - training data with Mordred descriptors
- `data/perm_random80_train_raw.ods` - training data before processing
The models are provided in the folder `models/` and can be loaded with the `load_model()` method of the respective class. The models provided are:
- `models/rf_random_dw.pkl` - Random Forest trained on DataWarrior descriptors
- `models/rf_random_mordred.pkl` - Random Forest trained on Mordred descriptors
- `models/xgb_random_dw.pkl` - XGBoost trained on DataWarrior descriptors
- `models/xgb_random_mordred.pkl` - XGBoost trained on Mordred descriptors
## ✅ Citation
```bibtex
@Misc{this_repo,
author = { Rebecca M Neeser },
title = { cyc_pep_perm - Python package to predict membrane permeability of cyclic peptides. },
howpublished = {Github},
year = {2023},
url = {https://github.com/schwallergroup/CycPepPerm }
}
```
## 🛠️ For Developers
<details>
<summary>See developer instructions</summary>
### 👐 Contributing
Contributions, whether filing an issue, making a pull request, or forking, are appreciated. See
[CONTRIBUTING.md](https://github.com/schwallergroup/CycPepPerm/blob/master/.github/CONTRIBUTING.md) for more information on getting involved.
### 🥼 Testing
After cloning the repository and installing `tox` with `pip install tox`, the unit tests in the `tests/` folder can be
run reproducibly with:
```shell
$ tox
```
Additionally, these tests are automatically re-run with each commit in a [GitHub Action](https://github.com/schwallergroup/CycPepPerm/actions?query=workflow%3ATests).
### 📖 Building the Documentation
The documentation can be built locally using the following:
```shell
$ git clone git+https://github.com/schwallergroup/CycPepPerm.git
$ cd CycPepPerm
$ tox -e docs
$ open docs/build/html/index.html
```
The documentation automatically installs the package as well as the `docs`
extra specified in the [`setup.cfg`](setup.cfg). `sphinx` plugins
like `texext` can be added there. Additionally, they need to be added to the
`extensions` list in [`docs/source/conf.py`](docs/source/conf.py).
### 📦 Making a Release
After installing the package in development mode and installing
`tox` with `pip install tox`, the commands for making a new release are contained within the `finish` environment
in `tox.ini`. Run the following from the shell:
```shell
$ tox -e finish
```
This script does the following:
1. Uses [Bump2Version](https://github.com/c4urself/bump2version) to switch the version number in the `setup.cfg`,
`src/cyc_pep_perm/version.py`, and [`docs/source/conf.py`](docs/source/conf.py) to not have the `-dev` suffix
2. Packages the code in both a tar archive and a wheel using [`build`](https://github.com/pypa/build)
3. Uploads to PyPI using [`twine`](https://github.com/pypa/twine). Be sure to have a `.pypirc` file configured to avoid the need for manual input at this
step
4. Push to GitHub. You'll need to make a release going with the commit where the version was bumped.
5. Bump the version to the next patch. If you made big changes and want to bump the version by minor, you can
use `tox -e bumpversion -- minor` after.
</details>
Raw data
{
"_id": null,
"home_page": "https://github.com/schwallergroup/CycPepPerm",
"name": "cyc-pep-perm",
"maintainer": "Rebecca M Neeser",
"docs_url": null,
"requires_python": ">=3.8",
"maintainer_email": "rebecca.neeser@epfl.ch",
"keywords": "snekpack, cookiecutter, chemistry, ai, ml, deep-learning",
"author": "Rebecca M Neeser",
"author_email": "rebecca.neeser@epfl.ch",
"download_url": "https://files.pythonhosted.org/packages/2f/90/090b353f5d1ec55f8261e178cc67773d8996ab438a76a31de7ae2942ba70/cyc_pep_perm-0.1.0.tar.gz",
"platform": null,
"description": "\n<!-- <picture>\n <source media=\"(prefers-color-scheme: dark)\" srcset=\"./assets/repo_logo_dark.png\" width='100%'>\n <source media=\"(prefers-color-scheme: light)\" srcset=\"./assets/repo_logo_light.png\" width='100%'>\n <img alt=\"Project logo\" src=\"/assets/\" width=\"100%\">\n</picture> -->\n\n<br>\n\n[](https://github.com/schwallergroup/CycPepPerm)\n<!-- [](https://doi.org/10.48550/arXiv.2304.05376)\n[](https://img.shields.io/pypi/v/CycPepPerm)\n[](https://img.shields.io/pypi/pyversions/CycPepPerm) -->\n[](https://github.com/psf/black)\n[](https://github.com/schwallergroup/liac-repo)\n[](https://schwallergroup.github.io)\n\n\n<h1 align=\"center\">\n CycPepPerm\n</h1>\n\n\n<br>\n\nPython package to predict membrane permeability of cyclic peptides.\n\n## \ud83d\udc69\u200d\ud83d\udcbb Installation\n\nWe provide the code as a [python package](https://pypi.org/project/cyc-pep-perm/), so the only thing you need is to install it. We recommend creating a new conda environment for that, which allows simple package management for a project. Follow these [instructions](https://docs.anaconda.com/free/anaconda/install/index.html) to install Anaconda. However, the package containing our code can also be installed without creating a project-specific environment. In that case, one just skips the first two lines of the following code:\n\n```bash\n$ conda create -n cyc-pep-perm python=3.10\n$ conda activate cyc-pep-perm\n$ pip install cyc-pep-perm==0.1.0\n```\n\n### \ud83d\udee0\ufe0f For Developers\n<details>\n <summary>See detailed installation instructions</summary>\nThe repository can be cloned from GitHub and installed with `pip` or `conda`. The code was built with Python 3.10 on Linux but other OS should work as well.\n\nWith conda:\n\n```bash\n$ git clone git+https://github.com/schwallergroup/CycPepPerm.git\n$ cd CycPepPerm\n$ conda env create -f environment.yml\n$ conda activate cyc_pep_perm\n$ pip install -e .\n```\n\nor with pip:\n\n```bash\n$ git clone git+https://github.com/schwallergroup/CycPepPerm.git\n$ cd CycPepPerm\n$ conda create -n cyc_pep_perm python=3.10\n$ conda activate cyc_pep_perm\n$ pip install -r requirements.txt\n$ pip install -e .\n```\n\nIf the options above did not work, please try from scratch:\n\n```bash\n$ git clone git+https://github.com/schwallergroup/CycPepPerm.git\n$ cd CycPepPerm\n$ conda create -c conda-forge -n cyc_pep_perm rdkit=2022.03.5 python=3.10\n$ conda activate cyc_pep_perm\n$ conda install -c conda-forge scikit-learn=1.0.2\n$ conda install -c rdkit -c mordred-descriptor mordred\n$ conda install -c conda-forge xgboost\n$ conda install -c conda-forge seaborn\n$ pip install shap\n$ conda install -c conda-forge jupyterlab\n$ pip isntall pandas-ods-reader\n$ pip install -e .\n```\n</details>\n\n## \ud83d\udd25 Usage\n\n> For some more examples on how to process data, train and evaluate the alogrithms, please consult the folder `notebooks/`. This folder also contains a notebook to perform polynomial fits as described in the paper.\n\nAll data paths in the following examples are taken from the hard-coded paths that work when one clones this repository. If you use the python package and download the data separately, please change the paths accordingly.\n\n### Data preprocessing\n\nHere we showcase how to handle the data as for our use-case. Some simple reformating is done (see also the notebook `notebooks/01_data_preparation.ipynb`) starting from `.ods` file with DataWarrior output (for data see [Data and Models](#data-and-models)).\n\n```python\nimport os\n\nfrom cyc_pep_perm.data.processing import DataProcessing\n\ndata_dir = \"/path/to/data/folder\" # ADAPT TO YOUR PATH!\n\n# this can also be a .csv input\ndatapath = os.path.join(data_dir, \"perm_random80_train_raw.ods\")\n\n# instantiate the class and make sure the columns match your inputed file - otherwise change arguments\ndp = DataProcessing(datapath=datapath)\n\n# make use of precomputed descriptors from DataWarrior\ndf = dp.read_data(filename=\"perm_random80_train_dw.csv\")\n\n# calculate Mordred deescripttors\ndf_mordred = dp.calc_mordred(filename=\"perm_random80_train_mordred.csv\")\n```\n\n### Training\n\nMake sure to have the data ready to be used. In order to make the hyperparameter search more extensive, please look into the respective python scripts (e.g. `src/cyc_pep_perm/models/randomforest.py`) and adjust the `PARAMS` dictionary.\n\n```python\nimport os\n\nfrom cyc_pep_perm.models.randomforest import RF\n\ndata_dir = \"/path/to/data/folder\" # ADAPT TO YOUR PATH!\ntrain_data = os.path.join(data_dir, \"perm_random80_train_dw.csv\")\nmodel_dir = \"/path/to/model/folder\" # ADAPT TO YOUR PATH!\nrf_model_trained = os.path.join(model_dir, \"rf_random_dw.pkl\")\n\n# instantiate class\nrf_regressor = RF()\n\nmodel = rf_regressor.train(\n datapath = train_data,\n savepath = rf_model_trained,\n)\n\ny_pred, rmse, r2 = rf_regressor.evaluate()\n# will print training results, e.g.:\n>>> RMSE: 8.45\n>>> R2: 0.879\n```\n\n### Prediction\n\n```python\nimport os\n\nfrom cyc_pep_perm.models.randomforest import RF\n\ndata_dir = \"/path/to/data/folder\" # ADAPT TO YOUR PATH!\ntrain_data = os.path.join(data_dir, \"perm_random80_train_dw.csv\")\nmodel_dir = \"/path/to/model/folder\" # ADAPT TO YOUR PATH!\nrf_model_trained = os.path.join(model_dir, \"rf_random_dw.pkl\")\n\n# instantiate class\nrf_regressor = RF()\n\n# load trained model\nrf_regressor.load_model(\n modelpath = rf_model_trained,\n)\n\n# data to predict on, e.g.:\ndf = pd.read_csv(train_data)\nX = df.drop(columns=[\"SMILES\"])\n\n# predict\ny_pred = rf_regressor.predict(X)\n```\n\n## Data and Models\n\nAll data required for reproducing the results in the paper are provided in the folder `data/`. Beware that due to the random nature of these models, the results might differ from the ones reported in the paper. The files found in `data/` are split into training and test data (randomly split 80/20) and with either the DataWarrior (dw) or the Mordred descriptors. The simple data processing can be found in the notebook `notebooks/01_data_preparation.ipynb`. The DataWarrior descriptors are computed with external software ([DataWarrior](https://openmolecules.org/datawarrior/)). The following files are provided:\n\n- `data/perm_random20_test_dw.csv` - test data with DataWarrior descriptors\n- `data/perm_random20_test_mordred.csv` - test data with Mordred descriptors\n- `data/perm_random20_test_raw.ods` - test data before processing\n- `data/perm_random80_train_dw.csv` - training data with DataWarrior descriptors\n- `data/perm_random80_train_mordred.csv` - training data with Mordred descriptors\n- `data/perm_random80_train_raw.ods` - training data before processing\n\nThe models are provided in the folder `models/` and can be loaded with the `load_model()` method of the respective class. The models provided are:\n\n- `models/rf_random_dw.pkl` - Random Forest trained on DataWarrior descriptors\n- `models/rf_random_mordred.pkl` - Random Forest trained on Mordred descriptors\n- `models/xgb_random_dw.pkl` - XGBoost trained on DataWarrior descriptors\n- `models/xgb_random_mordred.pkl` - XGBoost trained on Mordred descriptors\n\n## \u2705 Citation\n\n```bibtex\n@Misc{this_repo,\n author = { Rebecca M Neeser },\n title = { cyc_pep_perm - Python package to predict membrane permeability of cyclic peptides. },\n howpublished = {Github},\n year = {2023},\n url = {https://github.com/schwallergroup/CycPepPerm }\n}\n```\n\n\n## \ud83d\udee0\ufe0f For Developers\n\n\n<details>\n <summary>See developer instructions</summary>\n\n### \ud83d\udc50 Contributing\n\nContributions, whether filing an issue, making a pull request, or forking, are appreciated. See\n[CONTRIBUTING.md](https://github.com/schwallergroup/CycPepPerm/blob/master/.github/CONTRIBUTING.md) for more information on getting involved.\n\n### \ud83e\udd7c Testing\n\nAfter cloning the repository and installing `tox` with `pip install tox`, the unit tests in the `tests/` folder can be\nrun reproducibly with:\n\n```shell\n$ tox\n```\n\nAdditionally, these tests are automatically re-run with each commit in a [GitHub Action](https://github.com/schwallergroup/CycPepPerm/actions?query=workflow%3ATests).\n\n### \ud83d\udcd6 Building the Documentation\n\nThe documentation can be built locally using the following:\n\n```shell\n$ git clone git+https://github.com/schwallergroup/CycPepPerm.git\n$ cd CycPepPerm\n$ tox -e docs\n$ open docs/build/html/index.html\n```\n\nThe documentation automatically installs the package as well as the `docs`\nextra specified in the [`setup.cfg`](setup.cfg). `sphinx` plugins\nlike `texext` can be added there. Additionally, they need to be added to the\n`extensions` list in [`docs/source/conf.py`](docs/source/conf.py).\n\n### \ud83d\udce6 Making a Release\n\nAfter installing the package in development mode and installing\n`tox` with `pip install tox`, the commands for making a new release are contained within the `finish` environment\nin `tox.ini`. Run the following from the shell:\n\n```shell\n$ tox -e finish\n```\n\nThis script does the following:\n\n1. Uses [Bump2Version](https://github.com/c4urself/bump2version) to switch the version number in the `setup.cfg`,\n `src/cyc_pep_perm/version.py`, and [`docs/source/conf.py`](docs/source/conf.py) to not have the `-dev` suffix\n2. Packages the code in both a tar archive and a wheel using [`build`](https://github.com/pypa/build)\n3. Uploads to PyPI using [`twine`](https://github.com/pypa/twine). Be sure to have a `.pypirc` file configured to avoid the need for manual input at this\n step\n4. Push to GitHub. You'll need to make a release going with the commit where the version was bumped.\n5. Bump the version to the next patch. If you made big changes and want to bump the version by minor, you can\n use `tox -e bumpversion -- minor` after.\n</details>\n",
"bugtrack_url": null,
"license": "MIT",
"summary": "Python package to predict membrane permeability of cyclic peptides.",
"version": "0.1.0",
"project_urls": {
"Documentation": "https://cyc_pep_perm.readthedocs.io",
"Download": "https://github.com/schwallergroup/CycPepPerm/releases",
"Homepage": "https://github.com/schwallergroup/CycPepPerm",
"Source": "https://github.com/schwallergroup/CycPepPerm",
"Tracker": "https://github.com/schwallergroup/CycPepPerm/issues"
},
"split_keywords": [
"snekpack",
" cookiecutter",
" chemistry",
" ai",
" ml",
" deep-learning"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "d1286a6736bfe37b83ebf308a4e7c698d1dec0caff55be3632a04f302381841a",
"md5": "27882b9ec6e59d241e87bb681d3ad490",
"sha256": "cc4e692e65aafb43d52e7d73e8a3f5183cde8c58fb194f730d071ff87988e5cd"
},
"downloads": -1,
"filename": "cyc_pep_perm-0.1.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "27882b9ec6e59d241e87bb681d3ad490",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.8",
"size": 21890,
"upload_time": "2024-04-24T12:36:30",
"upload_time_iso_8601": "2024-04-24T12:36:30.003203Z",
"url": "https://files.pythonhosted.org/packages/d1/28/6a6736bfe37b83ebf308a4e7c698d1dec0caff55be3632a04f302381841a/cyc_pep_perm-0.1.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "2f90090b353f5d1ec55f8261e178cc67773d8996ab438a76a31de7ae2942ba70",
"md5": "e1ef9555f78692d1a07a37fd4492eb76",
"sha256": "9a8fa65777bc9b26557e2e52ec1dfc8d3a13339337f993ee517bce2c67bf37ed"
},
"downloads": -1,
"filename": "cyc_pep_perm-0.1.0.tar.gz",
"has_sig": false,
"md5_digest": "e1ef9555f78692d1a07a37fd4492eb76",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.8",
"size": 4419537,
"upload_time": "2024-04-24T12:36:33",
"upload_time_iso_8601": "2024-04-24T12:36:33.035416Z",
"url": "https://files.pythonhosted.org/packages/2f/90/090b353f5d1ec55f8261e178cc67773d8996ab438a76a31de7ae2942ba70/cyc_pep_perm-0.1.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-04-24 12:36:33",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "schwallergroup",
"github_project": "CycPepPerm",
"github_not_found": true,
"lcname": "cyc-pep-perm"
}
JSON
