Nineteen percent (19%) of structures in the training set are duplicates

Out of 4316 cases in the training set, only 3504 compound structures are unique. That means the training set contains 812 duplicate structures (18.81%). While that could be intentionally introduced noise, or even a decoy set, it is also quite common SMILES error and could be simply miss-annotations.

The implications are however severe, some of those contain also different flavor annotations.

Compound #288 and #1664 are exactly the same. This could be a problem of unique SMILES codes, which was not used by the data providers, leading to multiple annotations. Or as mentioned before its noise or decoy intentionally introduced by the competition providers.
image

Either way, if the 812 duplicate structures have different flavor annotations, they have to be removed (garbage-in --> garbage out). Otherwise the trained model is false. If the annotations are the same, one compounds can be removed, in order not to overfit the data with the same component twice.

See details here:
Trust, but verify: On the importance of chemical structure curation in cheminformatics and QSAR modeling research There are many ways to deal with that, one example is given here (duplicate removal in Python):

I have to include a picture below, because the web editor here actively destroys SMILES codes:
Training set structural duplicates (excerpt):

@Tobi, Chiral molecules like in pharmaceutical domain may have different behaviours. This is why you have those β€œchiral” SMILES. The chirality is very important (3D) and this is not duplicated data.

Cheers,

Guillaume

1 Like

@guillaumegodin Thanks, the examples above (I corrected them) clearly show that no double bond or tetrahedral isomers are in the original SMILES code (#288 and #1664). According to Daylight, these are designated by β€œβ€ and β€œ/” for cis/trans and β€œ@” and β€œ@@” for E/Z isomers.

However the (original training) SMILES codes are different, and show a compound that is just mirrored (I corrected that), but has no stereoisomers in the SMILES itself. This is the original data:

image

The correct way to describe the tetrahedral and double bond stereosimers would be by using the up and down bond or the β€œ@” and β€œ@@” signs in the SMILES code. But that is missing. The correct SMILES are in the figure below.

image

Also an additional independent analysis with topology analysis but different stereoisomers kept, came to a similar number 825 structures in the training set are exact duplicates. Maybe providing the correct stereo-isomer form in SD format would be helpful, or the structures are simply not correct. SMILES duplicates were tested in OpenBabel.

3 Likes

Dear all,

You must not take this statement made by @Tobi as correct see my complete response related to that in the other exact duplicate question posted here Training and test set compound overlapping (9%).

First, there is no duplicates. Second we used rdkit Inchi & smiles (we have canonical RDKit SMILES for the isolated stereoisomers, provided information is available and unambiguous) for double check, with multiple teams.

Another very important point canonical smiles is related to toolkit only so if you use another toolkit than RDkit, which is the most used across all AI papers > 60%, we will have other canonical representations due to other canonical procedures. Inchi is unique so you can use them to compare in any toolkit.

I played this game 5/6 years ago for fun: CDK vs RDKit in average have less than 10% identical canonical strings after the canonical procedure taking the same smiles as input for the two toolkits.

Inchi is available in RDkit anyway so double check is very easy!

best regards,

Guillaume

.