Analyse ToxTwin V1.3 GINEConv OGB d'une anthracycline de référence
La doxorubicine (adriamycine) est l’une des anthracyclines les mieux caractérisées sur le plan toxicologique. Cardiotoxicité dose-dépendante, mutagénicité documentée, perturbation endocrinienne : chaque axe de toxicité a fait l’objet de décennies de recherche clinique et préclinique. Cette molécule constitue donc un terrain de validation exigeant pour un modèle prédictif — non pas parce qu’elle est difficile à classer, mais parce que la concordance mécanistique entre les prédictions et les données expérimentales peut être évaluée endpoint par endpoint.
Le rapport présenté ici soumet la doxorubicine (SMILES canonique : COc1cccc2c1C(=O)c1c(O)c3c(c(O)c1C2=O)CC(O)(C(=O)CO)CC3OC1CC(N)C(O)C(C)O1) au pipeline ToxTwin V1.3, fondé sur une architecture GINEConv OGB entraînée sur les 12 endpoints Tox21, complétée par des modèles dédiés Ames (ICH S2(R1)) et hERG (ICH S7B).
Le résultat le plus significatif n’est pas la détection de la toxicité — attendue pour une molécule cytotoxique — mais la structure du profil prédit. Le triplet critique SR-ARE (1,000) / SR-MMP (0,957) / Ames (0,809) constitue une signature toxicologique cohérente avec la classe des anthracyclines. Le stress oxydant via le cycle rédox de la quinone, la dysfonction mitochondriale par ouverture du mPTP, et la génotoxicité par intercalation et inhibition de la topoisomérase II forment un triptyque que le modèle restitue avec une hiérarchie conforme aux données expérimentales.
La discrimination entre cardiotoxicité mécanistique (mitochondriale/oxydative) et cardiotoxicité électrophysiologique (hERG = 0,227) est cliniquement pertinente. La doxorubicine ne provoque pas de syndrome du QT long par inhibition du canal hERG ; sa cardiotoxicité opère par accumulation de dommages oxydatifs au niveau du cardiomyocyte. Le modèle ne confond pas ces deux voies — ce qui est un indicateur de cohérence architecturale plus informatif qu’un score AUC global.
Les scores élevés sur NR-AhR (0,750) et NR-PPAR-γ (0,750) méritent une lecture attentive. L’activation de l’AhR par la doxorubicine est documentée et implique une induction des CYP1A1/1B1 susceptible de moduler le métabolisme des xénobiotiques co-administrés. Le signal PPAR-γ s’inscrit dans les travaux récents sur la dysrégulation du métabolisme lipidique cardiaque comme facteur contributif de la cardiomyopathie aux anthracyclines — un axe de recherche qui n’est pas encore consolidé mais que le modèle identifie de manière autonome.
Les scores modérés sur les récepteurs endocriniens (NR-AR, NR-ER, aromatase) sont compatibles avec les effets cliniques connus : aménorrhée, infertilité, risque de néoplasies secondaires hormono-dépendantes. Le score NR-AR-LBD élevé (0,667) pose une question pertinente pour les protocoles combinant doxorubicine et thérapie antiandrogénique.
Le score p53 modéré (0,395) constitue la limite la plus visible. La voie p53 est centrale dans le mécanisme d’action antitumoral de la doxorubicine et dans sa génotoxicité. Cette sous-estimation probable reflète une limitation structurelle : les endpoints Tox21 mesurent l’activation à concentration fixe dans des lignées cellulaires spécifiques, ce qui peut ne pas capturer l’ampleur de la réponse p53 observée in vivo.
Plus fondamentalement, le modèle V1.3 opère sur la molécule parente. La doxorubicine est un substrat métabolique complexe — doxorubicinol, aglycone 7-déoxyaglycone, semi-quinone radicalaire. Chacun de ces métabolites possède un profil toxicologique distinct que le modèle ne capture pas. Cette limitation est intrinsèque à l’approche SMILES-in / prédiction-out et ne sera adressée qu’avec l’intégration d’un module de prédiction métabolique en amont (horizon V3.0).
Ce rapport ne constitue pas une découverte toxicologique — le profil de la doxorubicine est connu. Il constitue un cas de validation : la démonstration que le pipeline ToxTwin V1.3, sur une architecture GINEConv OGB avec descripteurs 2D, restitue un profil mécanistiquement cohérent pour une molécule de référence. La valeur du modèle ne réside pas dans sa capacité à prédire le connu, mais dans la structure de ses prédictions — structure qui sera mobilisée pour des molécules dont le profil toxicologique n’est pas encore caractérisé.