Skip to content

Commit 2045983

Browse files
IAlibayatravitz
andauthored
Apply suggestions from code review
Co-authored-by: Alyssa Travitz <31974495+atravitz@users.noreply.github.com>
1 parent 8ad344c commit 2045983

1 file changed

Lines changed: 10 additions & 11 deletions

File tree

cookbook/openfold3_to_openfe.ipynb

Lines changed: 10 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -17,7 +17,7 @@
1717
"source": [
1818
"## Introduction\n",
1919
"\n",
20-
"This notebook offers a brief preview of some of the steps one would need to go through to calculate binding affinities using [OpenFE](https://docs.openfree.energy/en/latest/) with [OpenFold3](https://openfold-3.readthedocs.io/en/latest/) generated structures.\n",
20+
"This notebook offers a brief preview of some of the steps you may need to go through to calculate binding affinities using [OpenFE](https://docs.openfree.energy/en/latest/) with [OpenFold3](https://openfold-3.readthedocs.io/en/latest/) generated structures.\n",
2121
"\n",
2222
"Notably, we will demonstrate how to generate a set of predictions, clean them up for use in molecular dynamics simulations, and then apply them to both relative binding free energy (RBFE) and absolute binding free energy (ABFE) methods."
2323
]
@@ -29,7 +29,7 @@
2929
"source": [
3030
"<div class=\"alert alert-block alert-danger\"> \n",
3131
"⚠️ <b>Important Note:</b>\n",
32-
" Please note that this notebook does not offer an exhaustive look into all everyone you may need to do when parsing OpenFold3 inputs into OpenFE. Where possible, we offer insights into some of the other hurdles you may face, but you are likely to encounter more in production settings.\n",
32+
" Please note that this notebook does not offer an exhaustive look into everything you may need to do when parsing OpenFold3 inputs into OpenFE. Where possible, we offer insights into some of the other hurdles you may face, but you are likely to encounter more in production settings.\n",
3333
"</div>"
3434
]
3535
},
@@ -47,7 +47,7 @@
4747
"\n",
4848
"**The tutorial is broken into the following steps:**\n",
4949
"\n",
50-
"1. Run cofolding with OpenFold3 on a set of TYK2 ligands\n",
50+
"1. Running cofolding with OpenFold3 on a set of TYK2 ligands\n",
5151
"\n",
5252
"2. Processing the OpenFold3 output for use with OpenFE\n",
5353
"\n",
@@ -138,7 +138,7 @@
138138
"id": "80da7967-8046-44e4-aef7-c81ac6712f98",
139139
"metadata": {},
140140
"source": [
141-
"We then create the json file that will be the OpenFold input."
141+
"We then create the JSON file that will be the OpenFold input."
142142
]
143143
},
144144
{
@@ -231,13 +231,11 @@
231231
" - This is necessary for RBFE simulations, where a single structure must be kept constant throughout the transformations. Note that for ABFE simulations, this step can be avoided.\n",
232232
" - For the sake of this demonstration, we will be choosing the predicted structure that is bound to the ligand occupying the largest volume. The idea here being that hopefully the binding site of that structure can accomodate all the other ligands. We note that in practice this is not always true.\n",
233233
"\n",
234-
"2. Aligning all the structures to the chosen structure\n",
234+
"2. Aligning the structures to the reference & extracting the necessary structures\n",
235235
"\n",
236-
" - This is so that we can extract the ligands and use them with the chosen structure.\n",
236+
" - This is so that we can extract the ligands and use them with the chosen structure.\n",.\n",
237237
"\n",
238-
"3. Extracting the ligands to SDF files and the chosen protein structure to a PDB file.\n",
239-
"\n",
240-
"4. Fixing protonation & capping issues with the PDB file.\n",
238+
"3. Fixing & protonating the reference protein structure.\n",
241239
"\n",
242240
"\n",
243241
"#### Choosing a reference protein structure\n",
@@ -324,7 +322,7 @@
324322
"import pathlib\n",
325323
"from rdkit.Chem.Descriptors3D import Asphericity\n",
326324
"\n",
327-
"# gather all the complexes and align to a single structure which has the largest ligand by volume?\n",
325+
"# gather all the complexes and align to a single structure which has the largest ligand by volume \n",
328326
"ligand_templates = {\n",
329327
" \"lig_ejm_31\": \"CC(=O)Nc1cc(NC(=O)c2c(Cl)cccc2Cl)ccn1\",\n",
330328
" \"lig_ejm_42\": \"CCC(=O)Nc1cc(NC(=O)c2c(Cl)cccc2Cl)ccn1\",\n",
@@ -500,7 +498,8 @@
500498
],
501499
"source": [
502500
"# now run the reference protein through pdbfixer\n",
503-
"! pdbfixer assets/of3_output/p_aligned/protein.pdb --add-atoms=all --add-residues --output=assets/of3_output/p_aligned/protein_fixed.pdb"
501+
"%% bash\n"
502+
"pdbfixer assets/of3_output/p_aligned/protein.pdb --add-atoms=all --add-residues --output=assets/of3_output/p_aligned/protein_fixed.pdb"
504503
]
505504
},
506505
{

0 commit comments

Comments
 (0)