-
Notifications
You must be signed in to change notification settings - Fork 7
License
zlisysu/Weighted_cc
Folders and files
Name | Name | Last commit message | Last commit date | |
---|---|---|---|---|
Repository files navigation
Please ensure that your input file contains at least three columns: The first two columns are the names of the mutation ligands, and the third colunm are their pair-wise energy. Example: 1. Example with only three data columns: ligand1, ligand2 and pair-wise energy(Input file: Example/example_without_w). Usage: python wcc_main.py -f Example/example_without_w -r 3A -e -8.83 output: Node dG_cc path_dependent_error path_independent_error 3K -9.6490 0.6800 0.3917 4M -9.2534 0.6800 0.3917 3N -7.2975 0.4600 0.2542 3A -8.8300 0.0000 0.2126 3M -9.3731 0.4600 0.2542 4P -9.5708 0.9300 0.3917 4L -8.2269 0.9300 0.3917 4N -8.4780 0.8500 0.2580 4I -7.7850 0.6800 0.2542 2. Example with five data columns: ligand1, ligand2, pair-wise energy, bar_std, slide_std(Input file: Example/example_with_w) Usage: python wcc_main.py -f Example/example_with_w -r 3A -e -8.83 output: Node dG_cc dG_wcc1 dG_wcc2 path_dependent_error path_independent_error 3K -9.6490 -9.8057 -9.9030 0.6800 0.3917 4M -9.2534 -9.1101 -9.1422 0.6800 0.3917 3N -7.2975 -7.3467 -7.4097 0.4600 0.2542 3A -8.8300 -8.8300 -8.8300 0.0000 0.2126 3M -9.3731 -9.3372 -9.3444 0.4600 0.2542 4P -9.5708 -9.6788 -9.7445 0.9300 0.3917 4L -8.2269 -8.2082 -8.3992 0.9300 0.3917 4N -8.4780 -8.5057 -8.5536 0.8500 0.2580 4I -7.7850 -7.8019 -7.9035 0.6800 0.2542 3. If you want obtain pair-wise results after calculation, add "-p yes" option Usage: python wcc_main.py -f Example/example_with_w -r 3A -e -8.83 -p yes output: Printing Pairwise Energies: Pair ddG_cc ddG_wcc1 ddG_wcc2 pair_error 3K-4M 0.3956 0.6956 0.7607 0.5600 3N-3A -1.5325 -1.4833 -1.4203 0.4600 3M-3A 0.5431 0.5072 0.5144 0.4600 3K-3N 2.3515 2.4590 2.4932 0.5000 4P-3K -0.0782 -0.1269 -0.1585 0.6300 4M-4P -0.3174 -0.5688 -0.6022 0.6300 3K-4L 1.4221 1.5975 1.5038 0.6300 4M-3M -0.1197 -0.2271 -0.2022 0.5000 4N-3K -1.1710 -1.3001 -1.3494 0.5100 4I-3N 0.4875 0.4552 0.4938 0.5000 4I-3M -1.5881 -1.5354 -1.4409 0.5000 4M-4L 1.0264 0.9018 0.7432 0.6300 4M-4N 0.7753 0.6043 0.5888 0.5100 **************************************************************************************************** Node dG_cc dG_wcc1 dG_wcc2 path_dependent_error path_independent_error 3K -9.6490 -9.8057 -9.9030 0.6800 0.3917 4M -9.2534 -9.1101 -9.1422 0.6800 0.3917 3N -7.2975 -7.3467 -7.4097 0.4600 0.2542 3A -8.8300 -8.8300 -8.8300 0.0000 0.2126 3M -9.3731 -9.3372 -9.3444 0.4600 0.2542 4P -9.5708 -9.6788 -9.7445 0.9300 0.3917 4L -8.2269 -8.2082 -8.3992 0.9300 0.3917 4N -8.4780 -8.5057 -8.5536 0.8500 0.2580 4I -7.7850 -7.8019 -7.9035 0.6800 0.2542 References ---------- Li, Yishui & Liu, Runduo & Liu, Jie & Luo, Haibin & Wu, Chengkun & Li, Zhe. An Open Source Graph-Based Weighted Cycle Closure Method for Relative Binding Free Energy Calculations. Journal of Chemical Information and Modeling. 2023, 63, 2, 561–570 Refer to the publication above for a detailed description of the wcc method and the parameters and please cite it to support our work if you use this software in your research.
About
No description, website, or topics provided.
Resources
License
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published