|
1 | 1 | # Generating Partial Charges with the OpenFE CLI
|
2 | 2 |
|
3 |
| -This tutorial will show you how to use the OpenFE CLI to generate and store partial charges for a series of ligands |
4 |
| -which can be used with OpenFE protocols. It is recommended to use a single set of charges for each ligand to ensure reproducibility |
5 |
| -between repeats or consistent charges between different legs of a calculation involving the same ligand, like a relative binding |
6 |
| -affinity calculation for example. As such both the `plan-rbfe-network` and `plan-rhfe-network` commands will calculate |
7 |
| -partial charges for ligands making it expensive to run multiple network mappings while finding the optimal one for the resources |
8 |
| -available. Hence, the `charge-molecules` command offers a way to generate and store the charges into an SDF file which |
9 |
| -can be used with the rest of the OpenFE CLI and python API. |
| 3 | +It is recommended to use a single set of charges for each ligand to ensure reproducibility between repeats or consistent |
| 4 | +charges between different legs of a calculation involving the same ligand, like a relative binding affinity calculation for example. |
| 5 | +As such both the `plan-rbfe-network` and `plan-rhfe-network` commands will calculate partial charges for ligands making it expensive |
| 6 | +to run multiple network mappings while finding the optimal one for the resources available. |
| 7 | + |
| 8 | +Here we present a CLI tool to allow you to do this ahead of time, reducing overheads and further improving reproducibility. |
| 9 | +This tutorial will show you how to use the OpenFE CLI command `charge-molecules` to generate and store partial charges for a series of ligands |
| 10 | +into an SDF file which can be used with OpenFE protocols. |
10 | 11 |
|
11 | 12 | ## Charging Molecules
|
12 | 13 |
|
@@ -42,7 +43,8 @@ lig_ejm_42
|
42 | 43 |
|
43 | 44 | Generating partial charges with the `am1bcc` method can be slow as they require a semi-empirical quantum chemical calculation,
|
44 | 45 | we can however take advantage of multiprocessing to calculate the charges in parallel for each ligand which offers a
|
45 |
| -significant speed-up. The number of processors available for the workflow can be specified using the `-n` flag: |
| 46 | +significant speed-up. The number of processors available for the workflow can be specified using the `-n` flag. For |
| 47 | +example to spread out the calculation over 4 cores: |
46 | 48 |
|
47 | 49 | ```bash
|
48 | 50 | openfe charge-molecules -M tyk2_ligands.sdf -o charged_tyk2_ligands.sdf -n 4
|
@@ -82,7 +84,7 @@ Using Options:
|
82 | 84 | Partial Charge Generation: am1bcc
|
83 | 85 | ```
|
84 | 86 |
|
85 |
| -The full range of partial charge settings can be found in the snipit bellow, note that some may require installing extra packages. |
| 87 | +The full range of partial charge settings can be found in the snippet bellow, note that some may require installing extra packages. |
86 | 88 |
|
87 | 89 | ```yaml
|
88 | 90 | partial_charge:
|
|
0 commit comments