ZheLi-Lab/Weighted_cc
Folders and files
| Name | Name | 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 dG_wcc1 dG_wcc2 path_dependent_error path_independent_error 3K -9.6490 -9.8057 -9.9030 0.6800 0.6300 4M -9.2534 -9.1101 -9.1422 0.6800 0.6300 3N -7.2975 -7.3467 -7.4097 0.4600 0.5000 3A -8.8300 -8.8300 -8.8300 0.0000 0.4600 3M -9.3731 -9.3372 -9.3444 0.4600 0.5000 4P -9.5708 -9.6788 -9.7445 0.9300 0.6300 4L -8.2269 -8.2082 -8.3992 0.9300 0.6300 4N -8.4780 -8.5057 -8.5536 0.8500 0.5100 4I -7.7850 -7.8019 -7.9035 0.6800 0.5000 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.6300 4M -9.2534 -9.1101 -9.1422 0.6800 0.6300 3N -7.2975 -7.3467 -7.4097 0.4600 0.5000 3A -8.8300 -8.8300 -8.8300 0.0000 0.4600 3M -9.3731 -9.3372 -9.3444 0.4600 0.5000 4P -9.5708 -9.6788 -9.7445 0.9300 0.6300 4L -8.2269 -8.2082 -8.3992 0.9300 0.6300 4N -8.4780 -8.5057 -8.5536 0.8500 0.5100 4I -7.7850 -7.8019 -7.9035 0.6800 0.5000 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.6300 4M -9.2534 -9.1101 -9.1422 0.6800 0.6300 3N -7.2975 -7.3467 -7.4097 0.4600 0.5000 3A -8.8300 -8.8300 -8.8300 0.0000 0.4600 3M -9.3731 -9.3372 -9.3444 0.4600 0.5000 4P -9.5708 -9.6788 -9.7445 0.9300 0.6300 4L -8.2269 -8.2082 -8.3992 0.9300 0.6300 4N -8.4780 -8.5057 -8.5536 0.8500 0.5100 4I -7.7850 -7.8019 -7.9035 0.6800 0.5000 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.