In silico Modeling method voor computationele aquatische toxicologie van hormoonontregelende stoffen: een op software gebaseerde aanpak met QSAR toolbox

Samenvatting

De kwantitatieve structuur-Activiteitsrelatie (QSAR) modellering is een representatieve bioinformatica-geassisteerde methode in de toxicologische screening. Dit protocol laat zien hoe de Risico's van hormoonontregelende stoffen (EDs) in aquatische omgevingen op een computationeel wijze kunnen worden beoordeeld. Door gebruik te maken van de OESO QSAR toolbox implementeert het protocol een in silico-test voor het analyseren van de toxiciteit van hormoonontregelende stoffen in vissen.

Samenvatting

Computationele analyses van toxicologische processen maakt screening van chemische stoffen met een hoge doorvoer mogelijk en voorspelling van hun eindpunten in biologische systemen. Met name de modellen van de kwantitatieve structuur-activiteitsverhouding (QSAR) zijn in toenemende mate toegepast om de milieueffecten van een overvloed aan giftige stoffen te beoordelen. In de afgelopen jaren, sommige meer gemarkeerde soorten toxicantia zijn hormoonontregelende stoffen (EDs, die chemicaliën die kunnen interfereren met een hormoon-gerelateerde metabolisme). Omdat hormoonontregelende stoffen een aanzienlijke invloed kunnen hebben op de ontwikkeling en reproductie van dieren, is een snelle voorspellen van de nadelige effecten van EDs die in silico-technieken werken vereist Deze studie presenteert een in silico-methode om voorspellingsgegevens te genereren over de effecten van representatieve hormoonontregelende stoffen in aquatische gewervelde dieren, in het bijzonder vissoorten. Het protocol beschrijft een voorbeeld met behulp van de geautomatiseerde workflow van de QSAR toolbox-software die is ontwikkeld door de organisatie voor economische samenwerking en ontwikkeling (OESO) om acute ecotoxiciteits voorspellingen van EDs mogelijk te maken. Als gevolg hiervan worden de volgende gegevens bepaald: (1) berekening van de numerieke correlaties tussen de concentratie voor 50% van de letaliteit (LC₅₀) en de verdelingscoëfficiënt octanol-water (K_ow), (2) uitvoerprestaties waarbij de LC₅₀ waarden bepaald in experimenten worden vergeleken met die gegenereerd door berekeningen, en (3) de afhankelijkheid van oestrogeen receptor bindende affiniteit op de relatie tussen K_ow en LC₅₀.

Inleiding

Nieuwe ontwikkelingen in informatica en computationele technologie hebben de biologische wetenschappen gemachtigd met kwantitatieve methodologieën die hoge precisie en betrouwbaarheid bieden¹. In het bijzonder hebben algoritmen die worden gebruikt in de moleculaire taxonomie en classificatie van eigenschappen geleid tot kwantitatieve structuur-activiteitsrelaties (QSAR) modellen². Deze modellen correleren automatisch de chemische structuren en biologische activiteiten van een bepaalde chemische database en implementeren snel in silico-screening van een breed scala aan chemische substraten volgens hun medicinale of toxicologische acties³. QSAR-hulpmiddelen kunnen voorspellende toxiciteits profielen produceren als functie van vectoren van moleculaire descriptoren (d.w.z. fysisch-chemische parameters) van chemicaliën die van belang zijn voor het numeriek maken van categorische eindpunten⁴. Gewoonlijk wordt elk kwantitatief eindpunt weergegeven als een 2D teststappenvisualization versus wijzigingen in descriptor waarden. Vervolgens wordt een QSAR-model gegenereerd met behulp van (meervoudige) lineaire regressieanalyses. Zodra een gegevensset volledig is benut om een QSAR-model (de Trainingsset) te maken, wordt het model statistisch gevalideerd door de eindpunten te voorspellen van een groep chemicaliën die niet in de Trainingsset is opgenomen (de testset). Het model kan vervolgens worden gebruikt om de biologische activiteiten van niet-geteste verbindingen³te voorspellen.

Onder veel schadelijke chemicaliën, hormoonontregelende stoffen (EDs) zijn gemarkeerd als een groep van toxicantia die kunnen interfereren met talrijke hormoon-gerelateerde metabolismen in zoogdieren, amfibieën, en vis^5,6. EDs zijn bekend voor het opwekken van een verscheidenheid van nadelige effecten, zoals kanker en misvormingen, door het blokkeren of wijzigen van normale hormonale trajecten of het activeren van abnormale hormoon synthese/afbraak signalen. Dientengevolge, deze hormoon-nabootsen chemicaliën kunnen perturb endocriene systemen zodanig dat biologische ontwikkeling en reproductie van Wildlife dierenpopulaties worden belemmerd. In het bijzonder zijn de ecotoxicologische effecten van hormoonontregelende stoffen uitgebreid onderzocht in aquatische gewervelde dieren, die bijna identieke hormoonreceptor structuren hebben met die van zoogdieren, waaronder mensen. Omdat alle hormonale acties optreden bij lage doses in vivo, is het voorspellen van de potentiële toxiciteiten van ED-kandidaten met behulp van Rapid in silico-screening van cruciaal belang voor de gezondheid van het publiek en het milieu.

QSAR-modellen op basis van de toxicologie van EDs zijn uitgevoerd met behulp van zowel 2D-als 3D-descriptoren (bekend als 2D en 3D QSAR, respectievelijk), die de ED ligand bindende affiniteiten van oestrogeen, androgeen en progesteron-receptoren⁷onthullen. Ondanks de zeer nauwkeurige voordelen van 3D QSAR, waarin conformationele en elektrostatische interacties worden overwogen, behoudt 2D QSAR zijn eigen robuustheid in directe wiskundige algoritmen, snelle berekeningen en extreem lage reken krachten. Bovendien zijn 2D-QSAR-modellen flexibel voor gebruik in een breed scala aan toepassingen, terwijl relatief nauwkeurige Voorspellings prestaties worden bereikt.

De OESO QSAR toolbox is momenteel een van de meest gebruikte computersoftware tools die gratis beschikbare en vooraf gebouwde QSAR-modellen^8,9bieden. De Profiler maakt gebruik van 2D descriptor databases. Sinds de release van de eerste versie in 2008 is de software toegepast op het gebied van chemische en biologische industrieën, volksgezondheid en milieu veiligheid voor volledige of gedeeltelijke analyse van de potentiële Risico's van natuurlijke en synthetische verbindingen, met bijzondere belangen in carcinogenese^10,11,12, mutageniteit^13,14,15en ontwikkelingstoxiciteit¹⁶. De toepassing op aquatische toxicologie is ook aangetoond, met focus op bioaccumulatie en biotransformatie¹⁷.

De QSAR toolbox is nuttig gebleken bij het voorspellen van de toxiciteit op korte termijn van een breed scala van chemicaliën¹⁷, evenals de oestrogeen receptor (er) bindende verwantschap van EDs¹⁸. De acute ecotoxiciteiten van hormoonontregelende stoffen in aquatische gewervelde dieren zijn echter niet geanalyseerd met behulp van de QSAR toolbox. In deze studie, een typische en facile protocol wordt gepresenteerd voor het uitvoeren van QSAR Modeling op de acute nadelige effecten van hormoonontregelende stoffen met een focus in de vissoorten. De studie toont aan dat de QSAR toolbox een zeer toegankelijke software is voor het berekenen en voorspellen van de letaliteit/sterfte van aquatische gewervelde dieren voor sommige representatieve hormoonontregelende stoffen. statistische behandelingsmethoden voor de afgeleide in silico-gegevenssets worden gepresenteerd. Figuur 1 toont het algemene schema voor de algemene werking van de QSAR toolbox. De in Figuur 2 getoonde workflow biedt eenvoudige instructies voor het gebruik van de in silico-assay om acute ecotoxiciteit van doelstoffen zoals hormoonontregelende chemicaliën te voorspellen.

Protocol

1. uitrusting

1. Software: gebruik OESO QSAR toolbox 4,0 of nieuwer (gratis te downloaden van < https://qsartoolbox.org/download/?) en software voor gegevensanalyse.
2. Computer: voor de OESO QSAR toolbox gebruikt u: (i) Systeemtype: 64 bit, Windows 7 of nieuwer; (II) processor: i5 op 2,4 GHz, of een snellere processor of gelijkwaardige AMD CPU; (III) geïnstalleerd geheugen (RAM): 6 GB; (IV) harde schijf (HDD): 20 GB vrije ruimte op de harde schijf (OESO QSAR toolbox 4,3 Release Notes: < https://qsartoolbox.org/file/2019/02/Toolbox-4.3-Release-Notes-1.pdf >).

2. procedure

1. OESO QSAR toolbox
   Opmerking: de QSAR toolbox werkt in zes opeenvolgende stroom modules vanaf input en gevolgd door profilering, Data, categorie definitie, vullen van gegevenskloof, vervolgens rapport, gelegen op de top van de programma-interface.
   1. Verken de bovengenoemde zes fasen via zes werkbalkpictogrammen door met de linkermuisknop te klikken. Kijk eerst naar de stadia van de invoer, de vulling van gegevenskloof en het rapport dat nodig is om de geautomatiseerde workflow "ecotoxicologische endpoint" uit te voeren en om de resultaten ervan te documenteren.
   2. Neem een korte blik op optionele stadia profilering en Data. De profilerings fase biedt een eerste inzicht in het potentieel van de doel stof (eco)-toxiciteit en de milieu-lotkenmerken. Optionele gegevens fase maakt het zoeken naar beschikbare experimentele gegevens met betrekking tot de doel stof mogelijk.
2. Input
   1. Bij het starten van de QSAR toolbox begint de gebruiker standaard bij de input Toolbox stage. De QSAR toolbox maakt automatisch een werkend bestand met de naam "document 1", dat wordt weergegeven in het deelvenster stage Option aan de linkerkant van de programma-interface. Wijzig de naam van het bestand, indien gewenst, door met de rechtermuisknop op het werkbestand te klikken.
   2. Klik op de knop CAS # op de werkbalk acties, voer het CAS-nummer (Chemical Abstract Service) van de doel stof in het beschikbare tekstveld in en klik op zoeken. De tool zal dan zoeken naar de doel stof door CAS-nummer.
   3. Kies indien nodig andere zoekopties die beschikbaar zijn in de actiewerkbalk, zoals zoeken op stoffen naam of vereenvoudigde code voor het regelinvoer systeem (SMILES) voor moleculaire invoer. Glimlachen kan worden ingevoerd als 2D niet-stereochemisch of 3D stereochemische bevattende vormen. Klik respectievelijk op naam of structuur . Gebruik het gereedschap structuur om de doel stof te tekenen.
   4. Met de zoekfunctie worden de zoekresultaten weergegeven via databaserecords in een pop-upvenster. Kies de record rapportage van een ' hoge ' CAS-SMILES-relatie (CS-relatie veld) als er meerdere records voor de doel stof worden opgehaald door het selectievakje aan de linkerkant van de record aan te vinken. Klik op OK.
      Opmerking: vanaf dit punt is alleen mogelijk als de opgehaalde record een GLIMLACH code bevat, omdat de SMILES code (2D niet-stereochemisch bevattende vorm) de basis voor berekening is.
   5. Batch modus: om de in silico-assay voor meerdere doelstoffen uit te voeren, schrijft u een eenvoudige stoffen lijst in een teksteditor waarin elk CAS-nummer in één rij wordt vermeld (aanvullende figuur S3). Sla het tekstbestand met de juiste naam en extensie. txt op de computer op.
   6. Batch modus: Klik op gegevens. Ga vervolgens naar databases in het deelvenster fase optie aan de linkerkant van de programma-interface. Zorg ervoor dat databases die worden vermeld onder ecotoxicologische informatie worden gecontroleerd.
   7. Batch modus: Klik op invoer. Selecteer query in de werkbalk acties. Accepteer de instellingen die zijn ingesteld in stap 2.2.6 door te klikken op Ja in het dialoogvenster.
   8. Batch mode: Kies het tabblad CAS . Upload de stoffen lijst die is opgeslagen als tekstbestand via de ladingslijst van uw computer.
   9. Batch modus: er zijn twee knoppen toevoegen beschikbaar; Klik op de knop toevoegen onder aan het pop-upmenu en klik vervolgens op uitvoeren. In de QSAR toolbox wordt een bericht weergegeven over het aantal stoffen dat voor de zoekopdracht is opgehaald.
      Opmerking: sommige stoffen van de geladen lijst zijn mogelijk niet gevonden door de zoekfunctie of er kunnen meerdere vermeldingen beschikbaar zijn voor één CAS-nummer. Het is niet mogelijk om stoffen uit de opgehaalde verzameling stoffen te verwijderen.
3. Profiling
   Opmerking: de volgende sectie is optioneel. Als dit niet nodig is, gaat u verder met paragraaf 2,5.
   1. Klik op de Toolbox stage knop profilering. Ga naar profilerings methoden in het deelvenster fase optie aan de linkerkant van de programma-interface.
   2. Klik op Alles deselecteren. Controleer alle onder vooraf gedefinieerde profilers en die gerelateerd zijn aan aquatische toxiciteit vermeld onder eindpunt specifiek , zoals "acute aquatische toxiciteit classificatie door Verhaar (gewijzigd)".
   3. Voltooi de selectie. Klik vervolgens op de knop toepassen op de werkbalk acties.
      Opmerking: de QSAR toolbox biedt aanbevelingen voor een set profilers. Deze zijn gemarkeerd in groen (geschikt) en oranje (aannemelijk) bij het kiezen van opties ≫ kleur door: ≫ eindpunt geselecteerd in de gegevensmatrix in de linkerbovenhoek van profilerings methoden. Klik met de linkermuisknop op het gegevensmatrix veld naast het eindpunt van de interesse. Beschikbare eindpunten worden weergegeven in de eindpunt structuur naast het deelvenster fase optie. Het type Profiler- stof geeft aan of de doel stof een ' discrete chemische ' is. De informatie wordt weergegeven in de uitgevouwen eindpunt structuur "profiel", "vooraf gedefinieerd" en "stof type". Alleen als de doel stof een discrete chemische is, kan de geautomatiseerde workflow met succes worden uitgevoerd. "Acute aquatische toxiciteits classificatie door Verhaar (gewijzigd)" geeft een eerste raming van het acute aquatische toxiciteitsmechanisme van de doel stof^19,20. De informatie wordt weergegeven in de uitgebreide eindpunt structuur "profiel", "eindpunt specifiek" en "acute aquatische toxiciteits classificatie door Verhaar (gewijzigd)". Er zijn vijf klassen beschikbaar: (klasse 1) inerte chemicaliën (baseline toxiciteit); (klasse 2) minder inerte chemicaliën; (klasse 3) reactieve chemicaliën; (klasse 4), specifiek handelend chemicaliën; en (klasse 5) voor chemische stoffen die niet kunnen worden ingedeeld.
   4. Klik met de rechtermuisknop op parameter in de eindpunt structuur om, indien gewenst, geïntegreerde 2D-en 3D QSAR-modellen uit te voeren die beschikbaar zijn in de QSAR toolbox. Klik op Bereken/Extraheer alle parameters voor alle chemicaliën in het pop-upmenu.
   5. 2D en 3D QSAR modellen gecompileerd in parameter bieden numerieke waarden. Gebruik "profilerings methoden" voor kwalitatieve informatie (zie stap 2.3.1).
4. Gegevens
   Opmerking: deze sectie is optioneel. Als dit niet nodig is, gaat u verder met paragraaf 2,5.
   1. Klik op de knop gegevens van de werkset fase . Klik vervolgens op verzamelen op de werkbalk acties.
   2. Selecteer alle eindpunten om alle experimentele gegevens te verzamelen en Kies vervolgens voor het verzamelen van eindpunt specifieke experimentele gegevens. Als bijvoorbeeld aquatische toxiciteit de focus van de gebruiker is, klikt u op kies ≫ ecotoxicologische informatie ≫ aquatische toxiciteit > OK.
      Opmerking: kiezen voor het verzamelen van experimentele gegevens voor alle eindpunten kan leiden tot uitgebreide verwerkingstijd. De gebruiker kan de hiërarchie van de eindpunt structuur aanpassen aan het specifieke doel. Hiermee wijzigt u de manier waarop gegevens worden weergegeven.
   3. Indien gewenst, met de rechtermuisknop op het eindpunt van belang in de eindpunt-structuurgebied. Kies structuur hiërarchie instellen in het pop-upmenu. Organiseer de eindpunt structuur op de gewenste manier met de beschikbare voorwaarden en pijlen en klik op OK.
   4. Exporteer indien gewenst de verzamelde gegevens als een Excel-bestand. Klik met de rechtermuisknop op het eindpunt van belang en kies gegevensmatrix exporteren in het pop-upmenu.
   5. Een matrix export wizard wordt geopend en toevoegen van andere eindpunten aan de export lijst. Voltooi de selectie, klik op exporteren en sla het bestand op de computer op.
      Opmerking: het exporteren van gegevens uit alle databases is niet mogelijk. Gegevens die zijn opgehaald uit de database ' ECHA CHEM ' kunnen bijvoorbeeld niet worden opgeslagen.
5. Vullen van gegevenskloof
   1. Klik op de Toolbox stage knop Data gap vulling. Klik vervolgens op geautomatiseerd in de werkbalk acties.
   2. Selecteer ecotoxicologisch eindpunt > vis, LC50 (letale concentratie, 50%) op 96 h voor Pimephales promelas (sterfte). Klik op OK. Een ' Workflowcontroller ' wordt weergegeven en de verwerking duurt enkele minuten, met name in de batchmodus.
      Opmerking: de QSAR toolbox past automatisch een gedefinieerde set profilers toe bij het zoeken naar geschikte stoffen met beschikbare experimentele gegevens voor de voorspelling. De experimentele gegevens [bijv. effect concentraties 96 h LC₅₀ (p. promelas) of 96 h EC₅₀ (p. promelas, mortaliteit)] worden gebruikt om de voorspelling voor de doel stof te genereren door middel van een lineaire benadering of de dichtstbijzijnde Neighbor-methode. Houd er rekening mee dat de methoden van lineaire benadering en dichtstbijzijnde neighbor worden aangeduid als trendanalyse (gelabeld als "T") en read-across (gelabeld als "R"), respectievelijk.
   3. De gebruiker ontvangt een bericht als de voorspelling met succes wordt uitgevoerd. Klik op OK en sluit de workflowcontroller die ' voltooide workflow ' aangeeft door in de rechterbovenhoek op x te klikken.
   4. Batch mode: bij het starten van de geautomatiseerde workflow, wordt de gebruiker gevraagd om het bereik van de stoffen waarover de werkstroom uit te voeren opgeven. Accepteer het volledige bereik van stoffen die standaard in het dialoogvenster zijn geselecteerd door op OK te klikken .
   5. Batch modus: de gebruiker ontvangt geen bericht dat aangeeft of een voorspelling is uitgevoerd of niet is geslaagd. Sluit de workflowcontroller die aangeeft dat de werkstroom is voltooid aan het einde van de batchverwerking door op x te klikken in de rechterbovenhoek.
6. Verslag
   1. Als een voorspelling met succes is uitgevoerd, klikt u op de werkset fase knop rapport.
      Opmerking: er kunnen geen rapporten worden gegenereerd in de batchmodus.
   2. Scrol omlaag en zoek de voorspellingswaarde in het matrixveld in een gele gemarkeerde rij naast eindpunt "96-h." De voorspelde waarde is gelabeld met ' T ' of ' R '. Activeer dit specifieke gegevensmatrix veld door er met de linkermuisknop op te klikken.
   3. Klik op voor Spelling op de werkbalk acties. Pas de inhoud en vormgeving van het rapport aan in de pop-upwizard. Er zijn drie soorten rapporten beschikbaar: (i) voorspelling, (II) categorie en (III) gegevensmatrix.
   4. De wizard stelt de gebruiker in staat de naam en contactgegevens van de auteur in te vullen. Schrijf een korte samenvatting, geef een gedetailleerde uitleg van de mechanistische interpretatie, of geef verantwoording voor de geschiktheid van de voorspelling.
   5. Aanvullende informatie met betrekking tot de uitgevoerde voorspelling opnemen, indien gewenst. De mate van aanvullende informatie is afhankelijk van de gebruiker.
   6. Ga door de wizard door te klikken op Volgende. Klik ten slotte op rapport maken en sla de voorspellings-en categorie rapporten op als PDF-bestanden en de gegevensmatrix als een Excel-spreadsheet op de computer.
   7. Meer informatie over de functionaliteiten van de QSAR toolbox en geautomatiseerde workflows vindt u in de toepassingshandleiding van de OESO QSAR toolbox v. 4 (F1 Help op het toetsenbord). Details over de algoritmen en de logica achter de geautomatiseerde workflow worden beschreven door Dimitrov et al.⁸ en Yordanova et al.⁹.

3. toepassing

1. Gebruik de ondergrens van het betrouwbaarheidsinterval van 95% als u de voorspelde effect concentratie (d.w.z. 96-h LC₅₀ van P. promelas) in de milieurisicobeoordeling gebruikt. Zoek de gegevens op de eerste pagina van het opgeslagen Voorspellings rapport (PDF) in "Voorspellings overzicht", "voorspelde waarde: < gemiddelde > (van < ondergrens > naar < bovengrens >)."
   Opmerking: de hier gegeven opmerkingen zijn gebaseerd op de resultaten van de vergelijking tussen voorspelde en experimentele gegevens voor een reeks doelstoffen die in dit onderzoek worden gerapporteerd. Als u het onderste uiteinde van het 95%-betrouwbaarheids bereik selecteert, wordt de kans groter dat de voorspelde effect concentratie de werkelijke toxiciteit van de stof niet onderschat (Zie de representatieve resultaten). De voorspelde effectieve concentratie van de ondergrens van het 95% betrouwbaarheidsinterval zal daarom een veiligere basis voor risicobeoordeling vormen.

Resultaten

Het in dit onderzoek beschreven voorbeeld werd uitgevoerd voor kwantitatieve analyse en voorspelling van acute toxiciteiten van geselecteerde hormoonontregelende stoffen in vissen. Wanneer de voorspelde gegevenspunten werden uitgezet versus experimentele gegevenspunten als een schaal van log-log, een positieve correlatie tussen beide werd gevonden voor alle vissen en een representatieve soort, namelijk Pimephales promelas (Fathead minnow; Figuur 3). In beide gevallen leek de helling van de lineaire regressie vergelijkbaar te zijn (voorspelde LC₅₀/experimentele LC₅₀ = 0,611 en 0,602 voor alle vissen en P. promelas, respectievelijk). Vanwege de beperkte hoeveelheid experimentele gegevens was het aantal beschikbare waarden van experimentele observatie meestal kleiner dan dat van rekenkundige voorspelling. Toepassing van de tolerantie factor 5-voudig voor de computationele capaciteit²¹ resulteerde respectievelijk in 94% (34/36) en 96% (26/27) van de beschermende voorspelling voor alle vissen en P. promelas. Op basis van deze voorspelling leek 3 ', 5, 7-trihydroxy-4 ', 6-dimethoxyisoflavone en 1,4-benzenediol berekende LC₅₀ -waarden te vertonen die groter zijn dan de tolerantielimiet.

Om veiligheidsbeoordeling op de hoogste betrouwbaarheid mogelijk te maken, werd verdere computationele analyse uitgevoerd door de voorspelde ondergrens van het 95% betrouwbaarheidsinterval van LC₅₀ (in plaats van de gemiddelde waarden in Figuur 3) te plotten versus de experimenteel afgeleide waarden (Figuur 4). Bij deze evaluatie met een verhoogde veiligheids drempel bleek 92% (33/36) van de totale geteste hormoonontregelende stoffen in het beschermingsgebied te vallen in vergelijking met de experimenteel afgeleide waarden, met uitzondering van: 3 ', 5, 7-trihydroxy-4 ', 6- dimethoxyisoflavone; 1,4-benzenediol; en 4-hexylfenol.

Op basis van de beoordelingen van de volledige uit de databank beschikbare soorten, waarden voor de voorspelde en experimentele 96-h Log₁₀LC₅₀ vertoonden lineariteit met de log₁₀K_ow waarden in het domein tussen-1 en 7, wat duidt op een hyperbolische correlatie tussen LC₅₀ en K_ow. Er bestond een algemene tendens waarbij de LC₅₀ daalde voor hogere waarden van K_ow van EDs voor de gegevens verkregen uit zowel computationele voorspellingen als experimenten, wat suggereert dat de acute toxiciteit bij vissoorten voor hormoonontregelende stoffen met een hogere hydrofobiciteit (aanvullend figuur S1).

Door de op regels gebaseerde ER-Profiler die is ingebed in de QSAR toolbox van de OESO, zijn de ER-bindende affiniteiten van de hormoonontregelende stoffen gecategoriseerd als niet-bindende en zwakke, matige, sterke en zeer sterke bindmiddelen, in volgorde van toenemende binding-affiniteit¹⁸. Dienovereenkomstig kan de statistische verdeling van Log₁₀K_ow worden weergegeven als een kwalitatieve classificatie van er binding affiniteit (aanvullende figuur S2). Over het algemeen bleken de veranderingen in de distributie marges van K_ow en hun gemiddelde niveaus niet een gedefinieerde tendens te hebben. Evenzo werden de distributies van voorspelde en experimentele LC₅₀ getoond als de omvang van er binding affiniteit (Figuur 5). In dit geval waren de gemiddelde niveaus van voorspelde LC₅₀ voor er-bindmiddelen hoger dan die van niet-bindmiddelen. Voor de experimentele LC₅₀waren de gemiddelde niveaus van niet-en zwakke bindmiddelen daarentegen hoger dan die van sterkere er bindmiddelen.

Figuur 1: basisschema van de algemene werkstroom van de OESO QSAR toolbox.
Klik hier om een grotere versie van dit cijfer te bekijken.

Figuur 2: workflow.
Getoond is de workflow die de modules en sequenties conceptualiseren om de acute toxiciteiten van hormoonontregelende stoffen in vis te voorspellen met behulp van de OESO QSAR toolbox. Klik hier om een grotere versie van dit cijfer te bekijken.

Figuur 3: voorspelde versus experimentele 96-h LC₅₀ van EDs in tabel 1 voor alle vissen (blauwe diamanten, n = 36) en een geselecteerde soort P. promelas (cyaan diamanten, n = 27).
Voor de voorspelde LC₅₀worden de gemiddelde waarden ("Ave") weergegeven. De onderbroken lijnen vertegenwoordigen lineaire regressies voor de twee groepen: voor alle vissen (licht blauw), voorspelde LC₅₀^Ave = 0,611 x (experimentele LC₅₀) + 0,277 (aangepast r² = 0,408); en voor P. promelas (licht cyaan), voorspelde LC₅₀^Ave = 0,602 x (experimentele LC₅₀) + 0,385 (aangepast r² = 0,441). De vaste diagonale lijn toont eenheid waarin de voorspelde en experimentele waarden gelijk zijn aan²¹. De gestippelde grijze lijn toont de 5-voudige tolerantielimiet van de rekenkundige mogelijkheid¹⁹. Uitschieters: 3 ', 5, 7-trihydroxy-4 ', 6-dimethoxyisoflavone (*) en 1,4-benzenediol (* *). Klik hier om een grotere versie van dit cijfer te bekijken.

Figuur 4: voorspelde (ondergrens van 95% betrouwbaarheidsinterval, "low-95%") vs. experimentele 96-h LC₅₀ van hormoonontregelende stoffen in tabel 1 voor alle vissen (n = 36).
De onderbroken lijn staat voor de lineaire regressie: voorspelde LC₅₀^Low-95% = 0,470 x (experimentele LC₅₀)-0,312, waarbij gecorrigeerd r² = 0,193. De vaste diagonale lijn geeft de eenheid aan waar de voorspelde en experimentele waarden gelijk zijn aan elkaar¹⁹. Uitschieters: 3 ', 5, 7-trihydroxy-4 ', 6-dimethoxyisoflavone (*), 1,4-benzenediol (* *) en 4-hexylphenol (* * *). Klik hier om een grotere versie van dit cijfer te bekijken.

Figuur 5: distributies van voorspelde (vaste dozen, n = 8 – 20 voor elke categorie) en experimentele (onderbroken dozen; n = 3 – 16 voor elke categorie) 96-h LC₅₀ afhankelijk van er bindende affiniteit van EDs in tabel 1 voor alle vissen.
Een vak plot vertegenwoordigt: (A) gemiddelde (klein vierkant met een vette horizontale balk), (B) 1^St en 3^RD kwartiel (onderste en bovenste- uiteinden van de doos, respectievelijk), (C) mediaan (horizontaal segment in de doos), ( D) 5^th en 95^th percentiel (respectievelijk onderste en bovenste foutbalken), (E) 1^St en 99^e percentiel (onderste en bovenste x, respectievelijk) en (F) minimum en maximum (onder-en boven-, respectievelijk). Klik hier om een grotere versie van dit cijfer te bekijken.

No.	CAS-registratienummer		Stofnaam		SMILES formule (2D niet-stereochemische vorm)						Log Kow	Ave voorspelde 96-h LC50 (mg/L)	LAGER 95% BI voorspelde 96-h LC50 (mg/L)	Profiler-oestrogeen receptor binding
1	50-28-2		17-β-Estradiol		CC12CCC3C (CCc4cc (O) ccc34) C1CCC2O						4,01	3,62	1,42	Zeer sterke Binder, OH Group
2	57-63-6		17-α Ethinyl- Estradiol		CC12CCC3C (CCc4cc (O) ccc34) C1CCC2 (O) C # C						3,67	3,00	1,18	Strong Binder, OH Group
3	80-05-7		2, 2-bis (4-hydroxyphe-Nyl) propaan (bisfenol A)		CC (C) (c1ccc (O) CC1) c1ccc (O) CC1						3,32	4,68	1,80	Zeer sterke Binder, OH Group
4	80-46-6		4-tert-Pentylphenol		CCC (C) (C) c1ccc (O) CC1						3,91	2,27	0,87	Zwakke Binder, OH Group
5	140-66-9		4-tert-Octylphenol		CC (c) (C) CC (C) (C) c1ccc (O) CC1						5,28	0,38	0,14	Strong Binder, OH Group
6	446-72-0		Genistein [3 ', 5, 7-trihydroxy-4 ', 6-dimethoxyisoflavone]		Oc1ccc (CC1) C1 = COc2cc (O) CC (O) c2C1 = O						2,84	32,00	10,03	Zeer sterke Binder, OH Group
7	10161-33-8		17β-trenbolon		CC12C = CC3C (CCC4 = CC (= O) CCC = 34) C1CCC2O						2,65	124,72	19,75	Strong Binder, OH Group
8	67747-09-5		Prochloraz (DMI-fungicide)		CCCN (CCOc1c (CL) CC (CL) cc1Cl) C (= O) n1ccnc1						4,1	5,19	1,74	Non Binder, zonder OH-of NH2-groep
9	84852-15-3		4-nonylfenol		CC (C) CCCCCCc1ccc (O) CC1						5,92	0,21	0,07	Strong Binder, OH Group
10	69-72-7		Salicylzuur		OC (= O) c1ccccc1O						2,26	24,07	9,31	Zwakke Binder, OH Group
11	80-09-1		4, 4 '-dihydroxydifenyl methioninesulfon (Bisfenol S)		Oc1ccc (CC1) S (= O) (= O) c1ccc (O) CC1						1,65	48,67	10,67	Zeer sterke Binder, OH Group
12	84-74-2		ftalinezuur, dibutylester		CCCCOC (= O) c1ccccc1C (= O) OCCCC						4,5	0,76	0,06	Non Binder, zonder OH-of NH2-groep
13	92-88-6		4, 4′-dihydroxybifenyl		Oc1ccc (CC1)-c1ccc (O) CC1						2,8	12,05	4,20	Matige Binder, OH grooup
14	94-13-3		4-hydroxybenzoë zuur, PROPYLESTER		CCCOC (= O) c1ccc (O) CC1						3,04	10,32	3,86	Matige Binder, OH grooup
15	98-54-4		4-tert-butylfenol		CC (C) (C) c1ccc (O) CC1						3,31	4,36	1,68	Zwakke Binder, OH Group
16	97-23-4		2, 2 ′-dihydroxy-– 5,5 ′-dichloorodifenyl-methaan		Oc1ccc (CL) cc1Cc1cc (CL) ccc1O						4,26	0,48	0,10	Zeer sterke Binder, OH Group
17	97-53-0		Eugenol		COc1cc (CC = C) ccc1O						2,27	14,70	5,60	Zwakke Binder, OH Group
18	99-76-3		4-hydroxybenzoinezuur, methyl ester		COC (= O) c1ccc (O) CC1						1,96	38,20	14,01	Zwakke Binder, OH Group
19	103-90-2		N-(4-hydroxyphenyl) aceetamide		CC (= O) Nc1ccc (O) CC1						0,46	338,97	43,39	Zwakke Binder, OH Group
20	106-44-5		p-cresol		Cc1ccc (O) CC1						1,94	20,47	7,14	Zwakke Binder, OH Group
21	108-39-4		m-cresol		Cc1cccc (O) C1						1,96	23,45	9,17	Zwakke Binder, OH Group
22	108-45-2		1,3-fenyleendiamine		Nc1cccc (N) C1						-0,33	34,60	0,00	Zwakke Binder, NH2 Group
23	108-46-3		1,3-Dihydroxybenzeen		Oc1cccc (O) C1						0,8	123,03	27,06	Zwakke Binder, OH Group
24	108-91-8		Cyclohexylamine		NC1CCCCC1						1,49	28,08	1,40	Zwakke Binder, NH2 Group
25	119-36-8		salicylzuur, methyl ester		COC (= O) c1ccccc1O						2,55	16,16	5,68	Zwakke Binder, OH Group
26	120-47-8		4-hydroxybenzoëzuur, ethyl ester		CCOC (= O) c1ccc (O) CC1						2,47	19,93	7,40	Zwakke Binder, OH Group
27	120-80-9		1,2-Dihydroxybenzeen		Oc1ccccc1O						0,88	11,14	0,01	Zwakke Binder, OH Group
28	123-31-9		1,4-Dihydroxybenzeen [1,4-benzenediol]		Oc1ccc (O) CC1						0,59	90,75	33,19	Zwakke Binder, OH Group
29	131-53-3		2, 2 ′-dihydroxy-4-methoxybenzophenone		COc1ccc (C (= O) c2ccccc2O) C (O) C1						3,82	3,97	1,46	Zeer sterke Binder, OH Group
30	131-56-6		2, 4-dihydroxybenzophenone		Oc1ccc (c (O) C1) C (= O) c1ccccc1						2,96	12,04	4,73	Strong Binder, OH Group
31	131-57-7		2-hydroxy-4-methoxybenzophenone		COc1ccc (C (= O) c2ccccc2) C (O) C1						3,79	5,96	2,27	Strong Binder, OH Group
32	599-64-4		4-cumylphenol		CC (C) (c1ccccc1) c1ccc (O) CC1						4,12	2,15	0,84	Strong Binder, OH Group
33	2855-13-2		1-Amino-3-aminomethyl-3, 5, 5-trimethyl-cyclohexaan		CC1 (C) CC (N) CC (C) (CN) C1						1,9	30,65	1,53	Moderate Binder, NH2 Group
34	6864-37-5		3, 3 ′-dimethyl-4, 4′-diaminodicyclohexylmethaan		CC1CC (CCC1N) CC1CCC (N) C (C) C1						4,1	1,07	0,05	Sterke Binder, NH2 Group
35	25013-16-5		tert-butyl-4-hydroxyanisol		COc1ccc (O) c (C1) C (C) (C) C						3,5	4,85	1,85	Matige Binder, OH grooup
36	147315-50-2		2-(4, 6-difenyl-1, 3, 5-triazin-2-yl)-5-(hexyloxy) fenol		CCCCCCOc1ccc (c (O) C1)-c1nc (NC (N1)-c1ccccc1)-c1ccccc1						6,24	0,17	0,06	Strong Binder, OH Group
37	88-68-6		2-aminobenzamide		NC (= O) c1ccccc1N						0,35	694,00	84,30	Zwakke Binder, NH2 Group
38	611-99-4		4, 4′-dihydroxybenzophenone		Oc1ccc (CC1) C (= O) c1ccc (O) CC1						2,19	37,74	14,67	Zeer sterke Binder, OH Group
39	27955-94-8		1, 1, 1-tris (4-hydroxyphenol) ethaan		CC (c1ccc (O) CC1) (c1ccc (O) CC1) c1ccc (O) CC1						4,38	2,09	0,82	Zeer sterke Binder, OH Group
40	87-18-3		salicylzuur, 4-tert-Butylphenyl Ester		CC (C) (C) c1ccc (OC (= O) c2ccccc2O) CC1						5,73	0,24	0,09	Strong Binder, OH Group
41	47465-97-4		3, 3-bis (3-methyl-4-hydroxyphenyl) 2-indolinone		Cc1cc (ccc1O) C1 (C (= O) Nc2ccccc12) c1ccc (O) c (C) C1						4,48	2,07	0,77	Zeer sterke Binder, OH Group
42	99-96-7		p-hydroxybenzoinezuur		OC (= O) c1ccc (O) CC1						1,58	8,54	0,00	Zwakke Binder, OH Group
43	80-07-9		1-Chloro-4-(4- chloorfenyl) sulfonylbenz		Clc1ccc (CC1) S (= O) (= O) c1ccc (CL) CC1						3,9	3,92	0,85	Non Binder, zonder OH-of NH2-groep
44	84-65-1		9, 10-anthraquinone		O = C1c2ccccc2C (= O) c2ccccc12						3,39	7,00	3,54	Non Binder, zonder OH-of NH2-groep
45	85-44-9		2-benzofuraan-1, 3-dione		O = C1OC (= O) c2ccccc12						1,6	2,69	0,00	Non Binder, zonder OH-of NH2-groep
46	92-84-2		10H-Fenothiazine		N1c2ccccc2Sc2ccccc12						4,15	1,07	0,08	Non Binder, zonder OH-of NH2-groep
47	2855-13-2		1-Amino-3-aminomethyl-3, 5, 5-trimethyl-cyclohexaan		CC1 (C) CC (N) CC (C) (CN) C1						1,9	30,65	1,53	Moderate Binder, NH2 Group
48	50-27-1		Estriol		CC12CCC3C (CCc4cc (O) ccc34) C1CC (O) C2O						2,45	21,21	8,29	Zeer sterke Binder, OH Group
49	50-50-0		Bèta-Estradiol-3-benzoaat		CC12CCC3C (CCc4cc (OC (= O) c5ccccc5) ccc34) C1CCC2O						5,47	0,36	0,02	Strong Binder, OH Group
50	53-16-7		Estron		CC12CCC3C (CCc4cc (O) ccc34) C1CCC2 = O						3,13	7,78	3,06	Strong Binder, OH Group
51	92-52-4		Bifenyl		c1ccc (CC1)-c1ccccc1						4,01	4,10	0,47	Non Binder, zonder OH-of NH2-groep
52	92-69-3		p-fenylfenol		Oc1ccc (CC1)-c1ccccc1						3,2	5,99	1,82	Matige Binder, OH grooup
53	96-29-7		2-Butanone Oxime		CCC (C) = NEE						0,63	32,67	2,49	Non Binder, niet-cyclische structuur
54	121-75-5		Malathon		CCOC (= O) CC (SP (= S) (OC) OC) C (= O) OCC						2,36	37,73	3,33	Non Binder, niet-cyclische structuur
55	123-07-9		4-ethylfenol		CCc1ccc (O) CC1						2,58	13,63	4,65	Zwakke Binder, OH Group
56	645-56-7		4-n-Propylpehnol		CCCc1ccc (O) CC1						3,2	7,32	2,55	Zwakke Binder, OH Group
57	1638-22-8		p-butyl fenol		CCCCc1ccc (O) CC1						3,65	4,09	1,39	Zwakke Binder, OH Group
58	1912-24-9		Atrazine		CCNc1nc (CL) NC (NC (C) C) N1						2,61	30,87	4,63	Non Binder, zonder OH-of NH2-groep
59	40596-69-8		Methopreen		COC (C) (C) CCCC (C) CC = CC (C) = CC (= O) OC (C) C						5,5	0,08	0,00	Non Binder, niet-cyclische structuur
60	1987-50-4		4-heptylphenol		CCCCCCCc1ccc (O) CC1						5,01	0,66	0,22	Matige Binder, OH grooup
61	92-86-4		p, p'-Dibromobiphenyl		Brc1ccc (CC1)-c1ccc (BR) CC1						5,72	0,11	0,02	Non Binder, zonder OH-of NH2-groep
62	480-41-1		Met		Oc1ccc (CC1) C1CC (= O) C2C (O) CC (O) cc2O1						2,52	27,84	10,87	Zeer sterke Binder, OH Group
63	486-66-8		Daïdzeïne		Oc1ccc (CC1) C1 = COc2cc (O) ccc2C1 = O						2,55	36,47	11,71	Zeer sterke Binder, OH Group
64	491-70-3		Luteoline		Oc1cc (O) c2C (= O) C = C (Oc2c1) c1ccc (O) c (O) C1						2,53	43,75	14,28	Zeer sterke Binder, OH Group
65	491-80-5		Biochanin A		COc1ccc (CC1) C1 = COc2cc (O) CC (O) c2C1 = O						3,41	15,87	3,70	Strong Binder, OH Group
66	520-18-3		Kaempferol		Oc1ccc (CC1) C1Oc2cc (O) CC (O) c2C (= O) C = 1O						1,96	70,98	8,05	Zeer sterke Binder, OH Group
67	2051-60-7		2-chloorobifenyl (PCB 1)		Clc1ccccc1-c1ccccc1						4,53	0,77	0,16	Non Binder, zonder OH-of NH2-groep
68	2051-61-8		3-chloorobifenyl (PCB 2)		Clc1cccc (C1)-c1ccccc1						4,58	0,77	0,16	Non Binder, zonder OH-of NH2-groep
69	2051-62-9		4-Chloro-1, 1 '-bifenyl		Clc1ccc (CC1)-c1ccccc1						4,61	0,77	0,16	Non Binder, zonder OH-of NH2-groep
70	2446-69-7		p-n-Hexylfenol [4-hexylphenol]		CCCCCCc1ccc (O) CC1						4,52	1,22	0,42	Matige Binder, OH grooup
71	14938-35-3		4-n-Amylfenol		CCCCCc1ccc (O) CC1						4,06	2,44	0,89	Zwakke Binder, OH Group
72	17924-92-4		Zearalenon		CC1CCCC (= O) CCCC = Cc2cc (O) CC (O) c2C (= O) O1						3,58	7,22	2,66	Strong Binder, OH Group
73	1743-60-8		Bèta-Estradiol 3-benzoaat 17-nbutyraat		CC (= O) OC1CCC2C3CCc4cc (O) ccc4C3CCC12C						4,95	0,91	0,35	Strong Binder, OH Group
74	479-13-0		Coumestrol		Oc1ccc2c (OC (= O) c3c-2oc2cc (O) ccc32) C1						1,57	52,16	11,44	Zeer sterke Binder, OH Group

Tabel 1: lijst van geëvalueerde hormoonontregelende stoffen. Gemiddeld gemiddelde (AVE) en lagere 95% betrouwbaarheidsinterval (CI) effectieve concentraties (95-h LC₅₀, Pimephales promelas) evenals oestrogeen receptor binding werden voorspeld met de qsar Toolbox versie 4,3 geautomatiseerde workflow. Log₁₀K_ow werd opgehaald via qsar Toolbox versie 4,3 van kowwin v 1.68, 2000, U.S. Environmental Protection Agency. Experimentele Log₁₀K_ow waarden hebben de voorkeur boven voorspelde waarden. De doel stof lijst is samengesteld uit eerder gerapporteerde lijsten van EDs^22,23,24.

Aanvullende informatie. Klik hier om dit bestand te downloaden. 

Discussie

De veelzijdigheid van de OESO QSAR toolbox als analytische software voor Ecotoxicologie wordt hier getoond met specifieke belangstelling voor de nadelige effecten van hormoonontregelende stoffen op aquatische gewervelde dieren. Daarnaast werd een eenvoudig en standaardprotocol aangetoond voor het voorspellen van acute toxiciteit (96-h LC₅₀) van 74 representatieve EDs (tabel 1) voor vissoorten. Dit werd bereikt door het toepassen van categorie bouw, data gap vulling en er profilerings modules ingebed in de QSAR toolbox (Figuur 1, Figuur 2).

De lineaire correlatie tussen log₁₀LC₅₀ en log₁₀K_ow met een negatieve helling (zoals weergegeven in aanvullende figuur S1) is al lange tijd bekend als een standaard kwantitatieve relatie in QSAR-analyses²⁵, waarbij hogere toxiciteit wordt aangetoond hoe meer hydrofobe een bepaalde chemische stof is. Zoals blijkt uit een eenvoudige berekening, is de algemene wiskundige relatie met vergelijking S1 en vergelijking S2 (aanvullende informatie) een geconverteerde uitdrukking uit de volgende Power functie²⁶:

Vanaf het perceel van (vergelijking 2) kan het karakteriseren van een tussen bereik van K_ow²⁶ mogelijk zijn door de parameters a en b aan te passen, waarbij een bepaalde variatie in hydrofobiciteit (of hydrofiele) de eindpunt van acute toxiciteit.

Vergelijkende analyses tussen de computationele voorspellingen en experimentele waarnemingen op de LC₅₀, zoals weergegeven in Figuur 3 en Figuur 4, zijn doorgaans gerapporteerd in studies van QSAR voor verschillende aquatische toxicanten, waaronder technische nonionische oppervlakteactieve stoffen²⁷, triazoolfungiciden²⁸en pesticide metabolieten²¹. Dit type retrospectieve validatie geeft informatie over de mate waarin een bepaalde QSAR-tool kan reiken in termen van vergelijkende prestaties voor experimentele resultaten. In deze studie van acute toxiciteit in vis, de QSAR toolbox bleek te bieden beschermende voorspellingen voor meer dan 90% van geteste EDs in alle vissen en in een enkele soort, Pimephales promelas.

Verdere identificatie van de drie uitschieter chemicaliën in Figuur 3 en Figuur 4, die een hogere voorspelde LC₅₀ gemiddeld en minimaal, respectievelijk, is gebleken. Eerste, de 3 ', 5, 7-trihydroxy-4 ', 6-dimethoxyisoflavone is een soort flavonoïde (meer in het bijzonder, een isoflavone), die wordt beschouwd als over het algemeen veilig en gebruikt in kruidengeneesmiddelen; echter, het heeft nog steeds oestrogeen-gerelateerde bezorgdheid²⁹ en kan leiden tot acute toxiciteit waarschijnlijk door middel van oxidatieve fosforylering ontkoppeling³⁰. Vervolgens is de 1,4-benzenediol, genaamd hydroquinone, een fenol verbinding die een niet-specifieke en cytotoxische immuunrespons in vissen kan veroorzaken³¹. Tot slot, de 4-hexylphenol is gekend om te exposeren voldoende positieve oestrogene activiteit worden geclassificeerd als een ED³². Het is goed bestudeerd dat de belangrijkste reden van de acute toxiciteit van Hydrochinon is de reductie-oxidatie (redox) fietsen. De Hydrochinon wordt geoxideerd tot Benzochinon en teruggebracht tot semi-quinon of Hydrochinon herhaaldelijk, met afputten van cofactoren en het genereren van reactieve zuurstof soorten³³. De andere twee chemicaliën kunnen diepere onderzoeken vereisen om hun werkingsmechanismen in acute ecotoxiciteit te onthullen met behulp van moleculaire docking benaderingen zoals die gebruikt door Panche et al.³⁴, die niet kan worden gedekt door de QSAR toolbox.

EDs interfereren met het endocriene systeem voornamelijk door middel van fysisch-chemische interacties met steroïde receptoren zoals oestrogeen en androgeen receptoren, die van groot belang zijn in QSAR Modeling studies³⁵. Gezien dit, de QSAR toolbox is robuust in termen van facile en snelle classificatie van ER bindende affiniteiten voor een set van chemicaliën op basis van alleen op de 2D descriptoren van moleculaire structuren. Toen dit ER Profiler systeem werd toegepast op onze lijst van hormoonontregelende stoffen, werd er geen duidelijke correlatie gevonden tussen ER bindingsaffiniteit en hydrofobiciteit (aanvullende figuur S2). Dit resultaat kan worden verklaard door het feit dat de vorming van een steroïde-receptor complex is niet een direct gevolg van een hydrofobe hechting bijdrage, maar moet gepaard gaan met een conformationele verandering in de actieve-site receptor structuur³⁶. De receptor binding kan ook worden veroorzaakt door waterstof binding en π-stacking.

Bovendien kan de positie van elke chemische groep op het molecuul de receptor binding beïnvloeden, zelfs als de hydrofobiciteit en het aantal van de waterstof-obligatie acceptoren-donoren hetzelfde blijven. Ten tweede produceerde de ER Profiler tegengestelde trends tussen voorspelde en experimentele LC₅₀ -gemiddelde niveaus met een toenemende er-bindingsaffiniteit (Figuur 5). Dit kan zijn omdat de letaliteit van ouders in een acute toxiciteitstest niet te wijten is aan ER binding, maar eerder aan narcose in de meeste gevallen, of aan redox fietsen in het geval van hydroquinone. Er is bijvoorbeeld een uitgebreidere analyse, inclusief de chronische toxiciteit, vereist voor een grotere set hormoonontregelende stoffen om voorspellende beperkingen van de huidige versie van de QSAR toolbox te definiëren.

Dit voorbereidend onderzoek kan ook gevolgen hebben voor de volksgezondheid, omdat steroïden (androgenen, oestrogenen, progestines, en corticoïden) en hun receptoren vertonen vergelijkbare of zelfs identieke macromoleculaire structuren over gewervelde dieren⁵. Deze soorten analoge endocriene signaleringssystemen kunnen werken met behulp van een gemeenschappelijk mechanisme in belangrijke gebeurtenissen van EDs⁵. Niettemin zijn aanvullende en complementaire methodologieën nodig om dit enorme en complexe aspect te verlichten [bijvoorbeeld door het uitvoeren van computationele modellering van absorptie, distributie, metabolisme en uitscheiding (ADME) en/of ongunstige uitkomst Pathway (AOP)]³⁸. Bovendien, omdat de meeste wetenschappelijke en publieke bezorgdheid over de nadelige effecten van hormoonontregelende stoffen verband houden met hun chronische toxiciteiten, het verbeteren van de databases en algoritmen in de QSAR toolbox en het produceren van betrouwbare lange termijn ecotoxicologie voorspellingen voor EDs zijn beide nodig.

Deze paper demonstreert de toepassing van de QSAR toolbox om ecotoxicologische LC₅₀ waarden voor vissen te vergelijken met Log₁₀K_ow waarden van EDs. in het gehele protocol leidt dit tot zwakke verhoudingen tussen de twee parameters, werd onthuld door eerdere studies (bijv., Kim et al.³⁹) dat Log₁₀K_ow is niet een goede directe voorspeller van aquatische LC₅₀. Ondanks deze beperking biedt dit protocol een algemene beoordeling of "vignet" om te beschrijven hoe het dashboard voor een bepaald doel moet worden gebruikt, omdat het een geldige toepassing is voor het gebruik van de QSAR toolbox voor het onderzoeken van correlaties tussen LC₅₀ (of er-binding -affiniteit) en log₁₀K_ow, of als een hulpmiddel voor snelle acute ecotoxiciteit screening. Niettemin moet worden opgemerkt dat (1) het verlichten van het verband tussen oestrogeen receptor binding en chronische toxiciteit, in plaats van acute toxiciteit (letaliteit), relevanter is zodat duidelijkere correlaties kunnen worden gevonden, en (2) de androgeen receptor, samen met die van oestrogeen, speelt ook een cruciale rol in reproductietoxiciteit. Daarom is het vereist voor de toekomstige versie van de QSAR toolbox om de Voorspellings functies in het licht van die twee punten te verbeteren.

Materialen

Lijst van materialen gebruikt in dit artikel

Naam	Bedrijf	Catalogusnummer	Opmerkingen
Adobe Acrobat Reader DC	Adobe Systems Software Ireland Limited	NA	Vereist om voorspelling en categorierapport te bekijken
Computer	Systeem: Microsoft Corporation	NA	Aanbevolen systeemeenheden: (i) systeemtype: 64 bit, Microsoft Windows 7 of nieuwer, (ii) processor: I5 op 2,4 GHz of snellere processor of gelijkwaardige AMD CPU, (iii) Geïnstalleerd geheugen (RAM): 6 GB RAM, (iv) Harde schijf (HDD): 20 GB vrije schijfruimte
Microsoft Editor	Microsoft Corporation	NA	Vereist voor het uploaden van een stoflijst van CAS-nummers (batchmodus) naar de OECD QSAR Toolbox als .txt-bestand (tekstbestand)
Microsoft Excel 2016	Microsoft Corporation	NA	Vereist om gegevens uit OECD QSAR Toolbox te exporteren als.cvs,.xls of.xlsx-bestanden
OECD QSAR Toolbox versie 4.0 of nieuwer	Organisatie voor Economische Samenwerking en Ontwikkeling	NA	Vereist om OECD QSAR Toolbox Automated Workflows uit te voeren; gratis download: https://qsartoolbox.org/download/
OriginPro 9	OriginLab Corporation	NA	Optioneel programma voor gegevensanalyse; vergelijkbare tools mogelijk