Fragments beat lead-like compounds in a screen against OGG1

The twin rise of make-on-demand libraries and speedy in silico docking has supercharged fragment screening and optimization: we’ve written previously about V-SYNTHES, Crystal Structure First and a related method. Another advance is described by Jens Carlsson (Uppsala University) and a large group of multinational collaborators in an (open access) Nat. Commun. paper.

 

The researchers were interested in 8-oxoguanine DNA glycosylase (OGG1), a DNA-repair enzyme and potential anti-inflammatory and anticancer target. They started with a crystal structure into which they docked 14 million fragments (MW < 250 Da) or 235 million lead-like molecules (250-350 Da) from ZINC15. Multiple conformations and thousands of orientations were sampled for each molecule. In all, 13 trillion fragment complexes and 149 trillion lead-like complexes were evaluated using DOCK3.7, a process that took just 2 hours and 11 hours on a 3500 core cluster.

 

After removing PAINS and molecules similar to previously reported OGG1 inhibitors, the top-scoring 0.05-0.07% molecules from each screen were clustered and, after manual evaluation, 29 fragments and 36 lead-like compounds were purchased from make-on-demand catalogs. These were tested at 495 µM (for fragments) or 99 µM (for larger molecules) in a DSF screen. None of the lead-like compounds significantly stabilized the protein, while several fragments did. Four of the fragments were successfully crystallized with OGG1, and in all cases the key interactions predicted in the computational screens were confirmed in the actual crystal structures.

 

Compound 1 showed the greatest stabilization of OGG1 (2.8 ºC) and some inhibition in an enzymatic assay, but not enough to calculate an IC₅₀. Searching for analogs that contained compound 1 as a substructure in the Enamine REAL database of 11 billion compounds produced few hits, but, as before, thinking in fragments proved fruitful. Searching for molecules containing just the core heterocycle and amide (colored blue below) yielded nearly 43,000 possibilities. Docking these and making and testing a few dozen led to compound 5, with mid-micromolar inhibition. Further iterations led to low micromolar compound 7.

At this point the researchers turned from make-on-demand libraries to synthetically accessible virtual libraries to fine-tune the molecule. After docking 6720 virtual molecules, they synthesized and tested 16, of which 12 were more potent than compound 7, with five of them being submicromolar. Compound 23 showed low micromolar activity in two different cell assays and was selective against four other DNA repair enzymes.

 

The same high-throughput docking approach was applied to three other protein targets: SMYD3, NUDT5, and PHIP. In each case crystal structures of bound fragments were available to use as starting points. Multiple compounds with improved docking scores compared to the initial fragments were identified, though no compounds were actually synthesized and tested.

 

The success in finding compound 1 demonstrates experimentally the advantage fragments have in efficiently searching chemical space. The researchers note that 97% of the >30 billion currently available make-on-demand compounds have molecular weights >350 Da, while only 50 million are < 250 Da. Screening all of these fragments in silico is possible; screening everything, less so. Although the fragment hits for OGG1 were weak, this isn’t always the case, as noted here. The fact that fragment 1 could be advanced to a sub-micromolar inhibitor after synthesizing just a few dozen molecules also testifies to the efficiency of in silico approaches.

 

The paper contains lots of useful details and suggestions for streamlining the process and is well worth perusing if you are trying to find hits against a structurally-enabled protein.