Skip to main content

Magic methyls and magic carpets


A few days ago, there was this post by Derek Lowe, reviewing a recent paper on magic methyls and their occurrence and impact in medicinal chemistry practice. They're called 'magic' because, although methyls are relatively insignificant in terms of size, polarity or lipophilicity, the addition of one in a compound can sometimes have a dramatic impact in its potency - much more that it would be attributed to any simple desolvation effects.

More generally, the 'magic methyl' phenomenon pops up in discussions about the validity of the molecular similarity principle, descriptors, QSAR - almost everything in the applied Chemoinformatics field - and belongs to the general class of 'activity cliffs'. 

Methylation is a chemical transformation, and transformations along with their impact on a property of choice can be easily mined and studied using the so-called Matched Molecular Pairs analysis (MMPA). We already have a comprehensive database of all the matched pairs and transformations in ChEMBL, so it was relatively straightforward to extract all the methylations (H>>CH3) recorded in ChEMBL_17 and analyse their impact in binding affinity. (b.t.w., MMPs are coming to the ChEMBL interface soon, so look out for this feature if you are interested in this area).

So, in more detail, I extracted all the H>>CH3 pairs and joined them with their pActivities (Ki, IC50, EC50) against human proteins as reported in the literature (our data validity flags were quite useful in this case). The trick here is to only consider molecule pairs tested against the same assay, so that their respective activities are directly comparable and one can safely subtract one from the other.

I ended up with 37,771 data points - much more than another recent publication that looked at this. Here's how the histogram of Delta pActivity (log units) looks like:

As you can see, the scale tilts slightly to the left of zero, meaning that methylation has overall neutral to negative effect on binding affinity. This is not the first time people see this. There are indeed, however, several examples (~2.3K out of 37.8K, to be exact) of magic methyls with more than 10-fold increase in activity. More about this later.

Some of you will ask: 'OK, but what about the context? - methylation of a carbon, nitrogen or oxygen is not the same'. You're right, it's not. So I trellised the above plot by a perception of context - i.e. whether the methylation happens next to an aromatic/aliphatic C or N or next to an oxygen:
The same trend, more or less, is observed with the exception of the aromatic carbon context, whereby methylation seems to have more favourable effect that expected by the overall distribution. Perhaps that could be explained by introducing torsional and planarity changes, etc. For a more thorough explanation of this, see here

Here are some examples of 'magic methyls' in the literature:

The take home message is: Magic methyls, unlike magic carpets, do exist but there are also equally as many, or even more, 'nasty' methyls. However, both of them are just a rather small minority compared to the 'boring' methyls - i.e. methyls with minimal or zero impact on potency.

It's just human nature to remember the few exceptions and outliers and forget the vast evidence to the contrary. However, isolating and understanding such edge cases and black swans is what could make the difference in drug discovery. 

George

Comments

Noel O'Boyle said…
So are "magic methyls" just another way of saying, "there is a part on the RHS of the normal curve which is a long way from the mean"? That is, is it a zero information content phrase.
kott said…
"there is a part on the RHS of the normal curve which is a long way from the mean" - I don't think this is a catchy title...

Popular posts from this blog

New SureChEMBL announcement

(Generated with DALL-E 3 ∙ 30 October 2023 at 1:48 pm) We have some very exciting news to report: the new SureChEMBL is now available! Hooray! What is SureChEMBL, you may ask. Good question! In our portfolio of chemical biology services, alongside our established database of bioactivity data for drug-like molecules ChEMBL , our dictionary of annotated small molecule entities ChEBI , and our compound cross-referencing system UniChem , we also deliver a database of annotated patents! Almost 10 years ago , EMBL-EBI acquired the SureChem system of chemically annotated patents and made this freely accessible in the public domain as SureChEMBL. Since then, our team has continued to maintain and deliver SureChEMBL. However, this has become increasingly challenging due to the complexities of the underlying codebase. We were awarded a Wellcome Trust grant in 2021 to completely overhaul SureChEMBL, with a new UI, backend infrastructure, and new f

A python client for accessing ChEMBL web services

Motivation The CheMBL Web Services provide simple reliable programmatic access to the data stored in ChEMBL database. RESTful API approaches are quite easy to master in most languages but still require writing a few lines of code. Additionally, it can be a challenging task to write a nontrivial application using REST without any examples. These factors were the motivation for us to write a small client library for accessing web services from Python. Why Python? We choose this language because Python has become extremely popular (and still growing in use) in scientific applications; there are several Open Source chemical toolkits available in this language, and so the wealth of ChEMBL resources and functionality of those toolkits can be easily combined. Moreover, Python is a very web-friendly language and we wanted to show how easy complex resource acquisition can be expressed in Python. Reinventing the wheel? There are already some libraries providing access to ChEMBL d

LSH-based similarity search in MongoDB is faster than postgres cartridge.

TL;DR: In his excellent blog post , Matt Swain described the implementation of compound similarity searches in MongoDB . Unfortunately, Matt's approach had suboptimal ( polynomial ) time complexity with respect to decreasing similarity thresholds, which renders unsuitable for production environments. In this article, we improve on the method by enhancing it with Locality Sensitive Hashing algorithm, which significantly reduces query time and outperforms RDKit PostgreSQL cartridge . myChEMBL 21 - NoSQL edition    Given that NoSQL technologies applied to computational chemistry and cheminformatics are gaining traction and popularity, we decided to include a taster in future myChEMBL releases. Two especially appealing technologies are Neo4j and MongoDB . The former is a graph database and the latter is a BSON document storage. We would like to provide IPython notebook -based tutorials explaining how to use this software to deal with common cheminformatics p

Multi-task neural network on ChEMBL with PyTorch 1.0 and RDKit

  Update: KNIME protocol with the model available thanks to Greg Landrum. Update: New code to train the model and ONNX exported trained models available in github . The use and application of multi-task neural networks is growing rapidly in cheminformatics and drug discovery. Examples can be found in the following publications: - Deep Learning as an Opportunity in VirtualScreening - Massively Multitask Networks for Drug Discovery - Beyond the hype: deep neural networks outperform established methods using a ChEMBL bioactivity benchmark set But what is a multi-task neural network? In short, it's a kind of neural network architecture that can optimise multiple classification/regression problems at the same time while taking advantage of their shared description. This blogpost gives a great overview of their architecture. All networks in references above implement the hard parameter sharing approach. So, having a set of activities relating targets and molecules we can tra

ChEMBL 26 Released

We are pleased to announce the release of ChEMBL_26 This version of the database, prepared on 10/01/2020 contains: 2,425,876 compound records 1,950,765 compounds (of which 1,940,733 have mol files) 15,996,368 activities 1,221,311 assays 13,377 targets 76,076 documents You can query the ChEMBL 26 data online via the ChEMBL Interface and you can also download the data from the ChEMBL FTP site . Please see ChEMBL_26 release notes for full details of all changes in this release. Changes since the last release: * Deposited Data Sets: CO-ADD antimicrobial screening data: Two new data sets have been included from the Community for Open Access Drug Discovery (CO-ADD). These data sets are screening of the NIH NCI Natural Product Set III in the CO-ADD assays (src_id = 40, Document ChEMBL_ID = CHEMBL4296183, DOI = 10.6019/CHEMBL4296183) and screening of the NIH NCI Diversity Set V in the CO-ADD assays (src_id = 40, Document ChEMBL_ID = CHEMBL4296182, DOI = 10.601