Evariste Technologies compounds

@mattodd

Hi all,

We’ve already caught up with @mattodd about this but will just give a quick intro for everybody else. Evariste Technologies is a start-up focusing on a probabilistic approach to drug discovery. Briefly, the platform we’ve built, Frobenius, takes an existing dataset and identifies the most promising starting point/s, then designs a bunch of new compounds and scores them according to the likelihood of achieving a set of pre-specified endpoints.

We were really interested in the recent publication detailing the open competition run by the OSM team and thought we’d have a crack at the problem ourselves. The compounds attached are the output generated by Frobenius when it’s presented with the series 4 data. More specifically, we’ve taken the two most promising starting points, applied the various compound designers and selected a subset of the highest scoring compounds (filtered by a medicinal chemist for synthetic feasibility etc). The number associated with each compound is the probability of it achieving a pIC50 of 8.

As we mentioned to Mat, we’re keen to get some of these synthesised and are able to contribute towards cost of synthesis. We’re also more than happy for anyone interested in these compounds to use them as inspiration for similar structures, if this is the case, we’d really appreciate being kept in the loop as we can very readily score the idea in Frobenius.

If anyone’s interested in knowing more about the modelling than the (very) brief overview I’ve given here we’re happy to discuss in detail.

Best wishes,
Alfie

https://www.evaristetechnologies.com
https://www.linkedin.com/in/alfie-brennan-746ba6b1

Series4_EVT_Suggestions.xlsx

Malaria suggestions pIC50 8 .pdf

@abrennan5

Hi @abrennan5 - this is very interesting, thanks for posting. It'd be really great to make some of these, yes.

Could someone please take a look and generate a pic of the molecules we can quickly post into an issue, just to help with digesting the ideas simply (I know there's a PDF, I'd just like a cdx. Alfie if you have one, you can drag it here, but you may have to zip it up for Github to accept it).

I'm intrigued by the meta-subst rings in the northeast, which is something we played with a bit in @maratsydney 's work, but have not explored a lot. The accessibility of those will depend on availability of materials, I suspect.

For EVT-004 there's a CF3 there. I thought we'd done that?

Also of interest is the EVT-007, EVT-008, EVT-009. Didn't we take that N away (in EVT-009) already? Replacement with C-F and C-Cl is something we've not done.

Hi Mat,

Compressed cdx below. Each compound designer should drop suggestions that are already in the dataset but if any of these have cropped up before let me know and we can fix that.

Malaria suggestions pIC50 8 .zip

I've updated the original post with the structures of the suggestions. Just some quick comments:

For EVT-004 there's a CF3 there. I thought we'd done that?
We have 1 compound with an OCF3 group as in EVT-004 but it has an amide on the LHS (MMV675963/OSM-S-271/TM 55-1).

Didn't we take that N away (in EVT-009) already?
Yes, we have 4 compounds with X = H instead of that N atom and a 4-CN on the RHS (from Patrick Thomson's work). Examples include both ether and amide linkers on the LHS but all are inactive.

I'm not sure if the synthetic accessibility of the benzylic OCHF2 group has changed but that was one of the things that we weren't able to remake from the inherited compounds. If we can't find a way to make that side-chain then all the compounds based on Starting Point 2 may not be accessible.

I'll have a look at whether the aldehydes needed for the different cores are available in the meantime and update later.

Great to meet you Edwin, thanks for uploading the image. Re the compounds containing the benzylic OCF2H, the OCH3 analogue was also one of the top ranked hits, see below for a set of suggestions from that starting point. Perhaps unsurprisingly, the model scores modifications of the methoxy group relatively highly as it views changes here as likely to increase potency. There were a significant number of designs in this region, all with around the same 4 - 7% chance of having pIC50 > 8. I selected the smaller changes to avoid having too much of an impact on solubility and logD. I've also (hopefully) fixed the filtering such that there shouldn't be any compounds already present in the dataset included here but please let me know if that isn't the case.

Let me know your thoughts around this set of molecules as well as just replacing the OCF2H with OCH3 in EVT-011/021.

OMe_suggestions.zip

@mattodd

Hi @mattodd and @edwintse,

I hope you're both well. Just wanted to follow up with our most recent work here, we're actually planning to publish the process of generating these ideas in a blog post at some point soon but I wanted to get them uploaded here first as it builds on so much of your hard work. You'll note the compounds suggested here are slightly different to those above. There are a couple of reasons for this; most importantly, instead of just optimising for improved potency, we asked the model to also aim for molecules with improved solubility and logD in the range of 0 - 4 (essentially 1 - 3 + prediction error). We have also used an updated version of our platform which predicts our error more accurately.

Finally, when selecting compounds we used a 'pessimistic design' algorithm which picks molecules for synthesis predicated on the failure of all the preceding molecules. This filters our list of 30 designs (based on three different starting points) down to 10. You'll note that in the PDF there are two predicted potency values, black is independent of the result for the earlier molecules, red is presuming that the preceding molecules failed to hit the desired endpoint (pIC50 > 8).

None of these feature the troublesome benzylic OCF2 due to the prioritisation of molecules which were predicted to be more soluble.

Happy to field any questions and share a wider selection of designs if that would be helpful.

Best wishes,
Alfie
compound_picker_suggestions.pdf
compound_picker_suggestions.zip

@abrennan5

Hi @abrennan5, I've finally had a look into potential synthesis of these compounds and it doesn't look like we'll be able to make any from your last post. I think the main issue is with the difluorophenyl ring - these substrates are typically more expensive to get a hold of compared to phenyl ring analogues, and this is especially the case when there are additional substituents on the ring. Just based on gut feeling, the compounds in your original post from Starting Point 1 are also less desirable as changes to the NE portion tend to decrease potency. I'm starting to think that the Starting Point 1 compound isn't the most ideal in terms of synthesising analogues. Any chance you'd be able to redesign some new compounds with this in mind?

@edwintse

Hi @edwintse, thanks for taking the time to review the molecules and potential routes, we really appreciate your input here. We can absolutely go from a different starting point, I'll remove both fluorines from that ring and see where we end up - most of the suggestions will design something back in in those positions, but the overall complexity should be lower.

I'll also look at building a virtual library using the published route, this sort of forward synthesis is something we're introducing more and more often on our projects. We can also take into account cost when doing this (very much a beta version of this feature at the minute and it is highly dependent on where the building block comes from).

I'll get something back to you later this week.

@edwintse

Hi @edwintse,

I ran the modelling again and selected some starting points which should (hopefully) be somewhat easier/cheaper to access. The two most promising analogues without any of the complex/expensive SMs discussed previously were:

indole in the NE position and the phenyl ring with the pendant ethyl alcohol in the NW
4-chlorophenyl in the NE position and the simple benzylethyl ether in the NW.

I applied our design algorithms but I also built virtual libraries using in stock alcohols and aldehydes available from Enamine so at least some of the compounds attached should be easy to plug into the existing route. The building block price filter is still work in progress so I wasn't able to apply it here.

The modelling was the same as the previous post so the designs have been biased for improved solubility and then selected based on predicted potency. Changes to the NE portion still crop up and I suspect this is the case because the modelling views regions of steep SAR as relatively promising for finding highly potent molecules. The reason for this is that they might well be rubbish, but they could also be great if they pick up the right interactions close to the surface of the protein. There are a few anilines present which I would normally filter out, but in this case I chose to leave them in as they are also likely to be easy to make.

Happy to hear your thoughts on these, have a great weekend!
Alfie
new_starting_points.zip

@abrennan5

Hi @abrennan5, thanks for the new suggestions! We might have a go at making the higher scoring ones like the cyclic urea (0.07). Interesting that the 0.14 compound is predicted to be active. We've made the same compound but with the NH2CH2 at the 3-position and it ended up being inactive.

@edwintse

@edwintse Glad to hear a few of these might make a target list! I initially couldn't find the molecule you mention so I've had a look and realised that there were a small number (about 25) that had a potency value in the Dundee column but not in the collated 'PfaI EC50 uMol (Mean)' column we'd been using.

I've updated the modelling and the indole predictions are broadly similar with slightly (2-fold) lower success estimates. The drop off is broadly consistent with the other suggestions, with the exception of the 0.14 compound, which is now 0.03 and probably better reflects your understanding of the SAR. No other starting points jump up the list, most of the missing data was for inactive molecules.

Hope that helps clear things up!

Just to update on where we're at with these compounds:

When trying to make the respective cores via the typical route, the condensation works fine but no reaction occurs under the oxidative cyclisation conditions.
A quick and dirty method to couple the diol to the core (for Suzuki coupling) was more dirty than quick.
I've made the mono-THP protected diol via the long route. Coupling to the core works fine but bromination with NBS gave a mix of products where the THP group had fallen off.
The Suzuki on the monobrominated free alcohol looks promising but the scale was too small.
I'm currently scaling up the coupled core and will retry the bromination using reXST'd NBS.

@edwintse

Hi @edwintse, thanks for the update and for organising the testing! Some interesting if not super exciting results. I think it's fair to say that the SAR really does not track in this series. I'll update the modelling and get back to you on Monday.

Thanks again!

@mattodd

@mattodd @edwintse Hope you're both well, here are some of my thoughts on the above data:

We would have loved to find something highly potent, but we're fairly pleased overall with the accuracy of the predictions. Our mean values were slightly high but we had large standard deviations for everything as well. Our modelling has also moved on slightly since these compounds were made, for example our prediction for EGT-555-1 would have been 4.8 +- 1.5.
Very interesting to see the comparison of 552-1 and 552-3. The fact the SAR is so clearly not additive suggests it's worth exploring further combinations in these positions.
A key question is; did we select the "best" (most potent/most information rich) compounds from the limited library we designed from this core or are there still really good compounds hiding further down the list? Adding this set of data (I've taken the average of the two results for each molecule as a starting point) and re-running the modelling on the same library picks out some 3,5 substituted molecules, often with the CF3 in one of the two positions. The top compound is 3-CF3-5-CN (5.4 +- 1.6).
Do you want to carry on with a few more cycles? If so would that be a few more compounds based on the existing library or some changes to the ether portion that use similar chemistry.

Really appreciate all your help on this project, it's great to have some data we can discuss openly.

@abrennan5

That's great @abrennan5. Interesting to get your thoughts here. I'd love to engage in another round and will check in with @edwintse. Are you able to post structures or SMILES for the new ones so that we can consider cost and synthetic difficulty?

@mattodd

@mattodd @edwintse

The structures here are the 10 things we would make next and the csv file contains the full library scored using the updated models. The values below the structures are: predicted pIC50, the error bar around that prediction, and the conditioned probability of achieving pIC50 > 7.

By conditioned probability I mean that, having selected the 'best' molecule, we assume that it fails to hit the endpoints we want, and rebuild the model before selecting the second one. This helps remove very similar molecules that have clustered at the top of the list. The full library contains the raw probability of success scores so you can look further down the list if you want.

Our predicted sigma is much higher than we normally see on datasets with 400ish compounds. This isn't necessarily a bad thing as it hopefully reflects the relatively steep SAR you seem to encounter in this part of the pocket. It also takes into account that the SAR between the NE and NW groups is not at all additive.

If we were to look into other chemistry, it would be to take a few of the more potent groups we have already found or those that show up in the next round and try them with different ethers. Perhaps those that include the benzylic substitution found in some of the more potent compounds.

Let me know your thoughts!

boronic_library_rescore.csv.zip

Minor edit: Added compound numbers. Also, we didn't filter out the aryl chlorides in the list below as they are semi-compatible with this chemistry but obviously it might lead to some useful but potentially annoying side products.

@abrennan5

Hi @abrennan5, thanks for generating some new compounds! I've had a look for the commercial availability of the corresponding boronic esters and they are shown below.

I was looking at the csv file and there's a compound in the list that isn't in the figure. The compound numbered 620 (2-OMe;3-OH) in the figure isn't in the list, and the compound numbered 250 in the list (3-Me; 5-NO2) isn't in the figure but is in the list?

The pyrrole might be interesting (closest we have already are imidazole and pyrazole, both were inactive).
The triazolopyridine might also be interesting
The 3,5-Cl;4-OCHF2 one might be good too

@edwintse

Hi @edwintse,

Sorry that's my mistake! The list is supposed to be the full library but I've accidentally taken a subset instead. The full library is now attached to this message, although the compounds are numbered differently. See if you're interested in any of the others in there.

You can also sort the compounds by "potency_mu" which is our mean prediction. You'll see that the top scorers here are more similar to the last round - number 10 and number 7 in original picture (your post, July 21st) are still predicted to be about 1 uM but unlikely to be a substantial more active than that. It might be interesting to include a couple of the higher confidence/exploitation compounds and a few of the riskier bets.

Maybe a list that looks something like 7 + 10 from the first list, as well as 1, 2, 5, and 8 (all in stock and fairly cheap) from the second round of designs. Does that sound reasonable to you? If we're lucky, 7 and 10 will also provide the de-chlorinated analogues which are predicted to be modestly potent as well.

boronic_library_rescore.csv.zip

@abrennan5

Hi @abrennan5, here's the list of the 6 compounds and their prices and codes from Enamine. I'm making more of the brominated core so it'll be ready to Suzuki when these arrive. Thanks again!

@mattodd

@mattodd I've uploaded the data curves from Mark for the latest compounds we got back. MMV1903416 and MMV1903417 were run 4 times each. The last repeats for those were inconsistent so Mark said to take the values that are highlighted in red (2 each which would be averaged).

Ed Report.xlsx

@abrennan5

Hi @abrennan5, sorry for the delay. I've just finished purifying all the compounds and just need to NMR check them but hopefully they're all ok. I managed to get a number of bonus compounds which were the result of dechlorination or double Suzuki coupling.

@mattodd In total there should be 12 compounds for testing (11 of these + 1 positive control)

@edwintse

Thanks @edwintse, this is brilliant! Great job getting everything separated - if those double suzuki bonus compounds are any good we're going to have a job getting the logP down.

I'll run all the compounds through the model again, although the double suzukis will fall outside the domain of applicability I'd imagine. Looking forward to seeing the results!

@edwintse

@edwintse @abrennan5 Now i'm curious as to what a naked or mono-chlorinated biphenyl would do

CLogP goes straight into the skip but curious nonetheless

Para-biphenyl-boronic acid is not that costly:

@mattodd

@mattodd @abrennan5 @miquelduranfrigola @GemmaTuron Results are in!! This round looks particularly promising. Quite a few of them are <1 uM.

The double Suzuki products are (I suppose) expectedly less potent
OSM-LO-69 is a compound that Jo made early on and it had an activity of 0.281 uM

@edwintse

Thanks @edwintse and @mattodd, this looks promising and there seems to be some interesting SAR to discuss. A lot of us are away on holiday at the minute but I'll follow up with an assessment of the predictions next week.

So far this project has been a real test of the modelling so it's great to see some fairly potent molecules. Thanks again!

@edwintse

@edwintse http://www.fluorochem.co.uk/Products/Product?Code=209754 (4‐cyclopropoxyphenyl)boronic acid is commercially available - might be worth trying to make something with it since it seems like para position is important for potency?

@edwintse

@edwintse @mattodd Hi both,

Image attached shows our predicted 50% confidence intervals for the new, non-biphenyl molecules. A comparison of 11 and 12 is illustrative for why we think predicting model error is so important. 11 has the highest mean prediction but a low confidence interval (exploiting well understood chemical space) whereas 12 has a modest mean prediction but significant uncertainty (more explorative) and as it turns out it's the best compound in the whole set. Overall really pleased with the results and very grateful for your help so far.

We're more than happy to do another round but I appreciate you might well need some time to focus on the suggestions from @miquelduranfrigola and @GemmaTuron, which is absolutely fair enough given the amount of compounds that have already been made based on our modelling. Happy to discuss this whenever you think would be suitable.

Link below is to a blog that we've just put on our website writing up a summary of the work so far.

https://www.evaristetechnologies.com/post/automated-design-of-potent-antimalarials

Navigation Menu

Repository navigation

Evariste Technologies compounds #29

27 remaining items

Metadata

Assignees

Labels

Type

Projects

Milestone

Relationships

Development

Evariste Technologies compounds #29

Description

Activity

mattodd commented on Jan 12, 2021

abrennan5 commented on Jan 12, 2021

edwintse commented on Jan 12, 2021

abrennan5 commented on Jan 12, 2021

abrennan5 commented on Apr 1, 2021

edwintse commented on Apr 27, 2021

abrennan5 commented on Apr 27, 2021

abrennan5 commented on Apr 30, 2021

edwintse commented on May 6, 2021

abrennan5 commented on May 7, 2021

edwintse commented on Jun 10, 2021

27 remaining items

abrennan5 commented on Nov 19, 2021

abrennan5 commented on Nov 29, 2021

mattodd commented on Nov 29, 2021

abrennan5 commented on Nov 30, 2021

edwintse commented on Dec 1, 2021

abrennan5 commented on Dec 1, 2021

edwintse commented on Dec 10, 2021

edwintse commented on Dec 13, 2021

edwintse commented on Jan 31, 2022

abrennan5 commented on Jan 31, 2022

MFernflower commented on Feb 1, 2022

edwintse commented on Mar 7, 2022

abrennan5 commented on Mar 7, 2022

MFernflower commented on Mar 7, 2022

abrennan5 commented on Mar 18, 2022

Metadata

Metadata

Assignees

Labels

Type

Projects

Milestone

Relationships

Development

Issue actions