Drug Design Update: Tool Talk

Started by free_radical

free_radical Lv 1

Hello everyone!

We are now ready to deploy a series of tools for small molecule drug discovery! Our goal is to release new tools on a rolling basis for the next couple of weeks. Because the tools are still being tested and not guaranteed to be bug free, we have created an “experimental” user group. This user group is open to everyone who wants to test the drug discovery tools; however, we do ask, that you report any bugs that you find or any suggestions that you have for the interface. After all the tools have been tested, we will release the drug discovery tools to Foldit’s main client.

You should expect this build to be experimental. This means that you have a chance of losing your designs. Because of this, it is highly recommend that you keep a separate install outside of your main client. This will ensure that you do not lose your main game puzzle progress.

In order to get access for the experimental group, follow the steps for access to the devprev build. Instead of replacing “main” with “devprev” in the options.txt file, type “experimental”.

"update_group" : "experimental"

The first tool is a very simple tool and is an expansion of the very first puzzle’s tool released for small molecule drug discovery. It is simply named, “Ligand Design Tool”. This tool provides the ability to change the identity of elements, create bonds between atoms, add a predefined set of fragments, delete atoms, and delete bonds. The Ligand Design Tool is available in both the Selection Interface and in the original Foldit interface.

You may want to change your view to the ligand design view. This will let you see the hydrogens on the small molecule, which can be used to extend the ligand by a single atom. You can access this by changing the view options in the advanced options menu.

In addition to releasing the first tool, we have three new tutorials. The tutorials are a work in progress and will be updated when we have some down time. The tutorials are particularly exciting as they follow the progression of a scientific team as they design an inhibitor for the FKB binding protein.

The puzzle that we are using for testing is Dihydrofulate Reductase (DHFR). Read on to get some more information on DHFR!

Finally, here is what to expect in the upcoming weeks. We will have a new blog post for each topic explaining in detail the science behind each tool and the tool's purpose:

Ligand minimization with a new force field, MMFF, and a set of new filters to guide you in designing small molecules, the rule of five filter and a similarity filter. We will have an explanation of all these features later on, but here is a video of the MMFF minimization.

There is also the ligand queueing interface. This tool allows us to give you small molecules that have been pre-identified (through experiments or virtual high throughput screening) that might bind in the target protein. This tool will also allow you to share your small molecule designs between your teams.

Finally, we are also creating a tool that lets you design small molecules like medicinal chemists. This is called reaction based drug design and provides a synthetic pathway for organic chemists to design the small molecule you create.

Additionally, I should be around in chat on Friday (April 29, 2016) around 2PM Eastern to help answer questions! Think of it like an "office hour", where you can drop in and get things answered versus an actual scientist chat.

free_radical Lv 1

Dihydrofolate reductase is an enzyme that converts dihydrofolate into tetrayhydrofolate (THF). THF is the active form of folic acid and a compound required for making thymine, DNA nucleotides, and certain amino acids. Important for cancer treatment, THF is also responsible for cell growth and division. The DHFR enzyme is the sole source of THF for your body; thus, inhibiting DHFR can reduce the proliferating cells found in cancer tissue.

Knowing the mechanism by which an enzyme works can be used as a starting point when designing small molecule drugs for protein targets. For DHFR, THF is created through catalysis with a cofactor NADPH. Small molecules that are similar to THF can “trick” the enzyme into trying to catalyze the drug instead of the original substrate. In fact, one of the inhibitors designed by chemists (Methotraxate) looks very similar to THF. Methotrexate is very similar to THF and only differs by a single methyl (CH3) and amino group (NH2).

For this series of puzzles, NADPH will be present, but you will not be able to design it. Instead, your goal is to design a small molecule inhibitor that would replace THF. A starting fragment is given to you to help in the design. Keep in mind that the closer you can mimic the original substrate (THF), the more successful your designs will be. If you find bugs, please report them here so that we can fix them!

altejoh Lv 1

And not just for drug testing, but as a potential modelling system for chemicals in general. Imagine adding in other first principles calculations to predict things like reactivity or spectra plots, while still being able to directly interact with your structure.

Bruno Kestemont Lv 1

Exciting !
9-1 is enough with only wiggling (that should not be possible for a tutorial on molecule design)
9-2 is enough with only mutating (that should not be possible for a tutorial on molecule design)
9-4 if we wiggle just after the balloon "try wiggling …", we escape the remaining text "tell me more" (however, it's interesting to read). Suggestion "try wiggling" only at the end

Bug for the tool in 9-2: crashed in 2 different situations (unable to reproduce any of them):
-selecting and unselecting around
-selecting and unselecting, selecting 4-5 atoms, changed to O => crash
-pull and adding atoms rounds, undo: at a certain moment, a spot remains at the place where an atom was. Clicking on this small spot crashes

Suggestion for the tool:
1- when selecting and changing an atom, the selection should still remain (in order to be able to quick try all atoms without being obliged to re-select manually)

free_radical Lv 1

Thanks for this bug report…I have a feeling that selecting 4-5 atoms and then hitting an oxygen caused a molecule that isnt chemically feasible to be produced. I had put some checks in for this, but I guess it didnt catch all of them. I am going to try and reproduce that.

Could you expand on this a little:

-pull and adding atoms rounds, undo: at a certain moment, a spot remains at the place where an atom was. Clicking on this small spot crashes

Did you hit undo and an atom was still highlighted?

bertro Lv 1

just clicking for next one.

==» debug.txt:

version: 1
{
"build_id" : "20160426-f037e02514-win_x86-experimental"
"crash_id" : "2107b6f0-02ab-4e0f-a071-bfdf79ac8944"
"current_actions" : ""
"current_tool" : ""
"file" : "core/conformation/Conformation.cc:850"
"machine_id" : "09edfe6b-95dc-4576-888d-aeeae5d74622"
"macro_id" : "-1"
"message" : "Can't create a polymer bond after residue 15 due to incompatible type: ARG:CtermProteinFull"
"player_id" : "0"
"puzzle_id" : "0"
"puzzle_running_time" : "503286"
"running_time" : "503287"
"timestamp" : "1461787297"
}
verify: 4bc937ed75f39b317f3b24f11ebe9d8b

==»log.txt:

args passed:
'F:\FolditNEW\Main C\Foldit\Foldit.exe'
starting the init thread!..
boinc base url: https://fold.it
checking updates…
binary
local: 'db836967ac190e672dc51d1003ec2823'
remote: 'db836967ac190e672dc51d1003ec2823'
database
local: 'f7518d932a1190b04c69a56e71985152'
remote: 'f7518d932a1190b04c69a56e71985152'
resources
local: '1810343af1fd6b549773b27e4ac00063'
remote: '1810343af1fd6b549773b27e4ac00063'
cleaning up old components:
binary 00000000000000000000000000000000
binary db836967ac190e672dc51d1003ec2823
database f7518d932a1190b04c69a56e71985152
resources 00000000000000000000000000000000
resources 1810343af1fd6b549773b27e4ac00063
CRASH: 253386
SoundTheme::load: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_00.ogg
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_01.ogg
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_02.ogg
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_03.ogg
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_04.ogg
loading: cmp-resources-1810343af1fd6b549773b27e4ac00063\resources\sounds/organic_01/rotamer_land_05.ogg
SVM classifier successfully loaded
Feature list successfully loaded
Valid classifier feature list
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue LEU:CtermProteinFull 21
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue GLU:CtermProteinFull 21
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue LYS:CtermProteinFull 21
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue ARG:CtermProteinFull 21
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue ASP:CtermProteinFull 21
core.conformation.Conformation: [ WARNING ] missing heavyatom: OXT on residue LEU:CtermProteinFull 17

ERROR: Can't create a polymer bond after residue 15 due to incompatible type: ARG:CtermProteinFull
ERROR:: Exit from: ......\source\src\core\conformation\Conformation.cc line: 850

UNHANDLED EXCEPTION
1: RaiseException +72 bytes (no line)
2: alIsAuxiliaryEffectSlot +283336 bytes (no line)
3: library_main +36421084 bytes (no line)
4: library_main +34771455 bytes (no line)
5: library_main +26630459 bytes (no line)
6: library_main +27504977 bytes (no line)
7: library_main +27510399 bytes (no line)
8: library_main +27510999 bytes (no line)
9: library_main +27509901 bytes (no line)
10: library_main +27509699 bytes (no line)
11: library_main +2950223 bytes (no line)
12: library_main +2929128 bytes (no line)
13: library_main +403430 bytes (no line)
14: library_main +347656 bytes (no line)
15: library_main +333803 bytes (no line)
16: library_main +223752 bytes (no line)
17: library_main +3256533 bytes (no line)
18: library_main +3256987 bytes (no line)
19: BaseThreadInitThunk +36 bytes (no line)
20: RtlUnicodeStringToInteger +595 bytes (no line)
21: RtlUnicodeStringToInteger +542 bytes (no line)

bertro Lv 1

The puzzle that we are using for testing is Dihydrofulate Reductase (DHFR).

Is there a puzzle because I only see the 3 tutorials? No puzzle in the puzzle menu.

bertro Lv 1

win10 fully updated not on preview. I will try on a new client (after the tutorials it did not show and I closed that client). Still no puzzle showing…