Computational prediction models for assessing endocrine disrupting potential of chemicals

Structure-based modeling to assess binding and endocrine disrupting potential of polycyclic aromatic hydrocarbons in Daniorerio

Environmental contaminants, such as polycyclic aromatic hydrocarbons (PAHs), have raised concerns regarding their potential endocrine-disrupting effects on aquatic organisms, including fish. In this study, molecular docking and molecular dynamics techniques were employed to evaluate the endocrine-disrupting potential of PAHs in zebrafish, as a model organism. A virtual screening with 72 PAHs revealed a correlation between the number of PAH aromatic rings and their binding affinity to proteins involved in endocrine regulation. Furthermore, PAHs with the highest binding affinities for each protein were identified: cyclopenta[cd]pyrene for AR (-9.7 kcal/mol), benzo(g)chrysene for ERα (-11.5 kcal/mol), dibenzo(a,e)pyrene for SHBG (-8.7 kcal/mol), dibenz(a,h)anthracene for StAR (-11.2 kcal/mol), and 2,3-benzofluorene for TRα (-9.8 kcal/mol). Molecular dynamics simulations confirmed the stability of the protein-ligand complexes formed by the PAHs with the highest binding affinities throughout the simulations. Additionally, the effectiveness of the protocol used in this study was demonstrated by the receiver operating characteristic curve (ROC) analysis, which effectively distinguished decoys from true ligands. Therefore, this research provides valuable insights into the endocrine-disrupting potential of PAHs in fish, highlighting the importance of assessing their impact on aquatic ecosystems.

Recognizing high-priority disinfection byproducts based on experimental and predicted endocrine disrupting data: Virtual screening and in vitro study

So far, about 130 disinfection by-products (DBPs) and several DBPs-groups have had their potential endocrine-disrupting effects tested on some endocrine endpoints. However, it is still not clear which specific DBPs, DBPs-groups/subgroups may be the most toxic substances or groups/subgroups for any given endocrine endpoint. In this study, we attempt to address this issue. First, a list of relevant DBPs was updated, and 1187 DBPs belonging to 4 main-groups (aliphatic, aromatic, alicyclic, heterocyclic) and 84 subgroups were described. Then, the high-priority endocrine endpoints, DBPs-groups/subgroups, and specific DBPs were determined from 18 endpoints, 4 main-groups, 84 subgroups, and 1187 specific DBPs by a virtual-screening method. The results demonstrate that most of DBPs could not disturb the endocrine endpoints in question because the proportion of active compounds associated with the endocrine endpoints ranged from 0 (human thyroid receptor beta) to 32% (human transthyretin (hTTR)). All the endpoints with a proportion of active compounds greater than 10% belonged to the thyroid system, highlighting that the potential disrupting effects of DBPs on the thyroid system should be given more attention. The aromatic and alicyclic DBPs may have higher priority than that of aliphatic and heterocyclic DBPs by considering the activity rate and potential for disrupting effects. There were 2 (halophenols and estrogen DBPs), 12, and 24 subgroups that belonged to high, moderate, and low priority classes, respectively. For individual DBPs, there were 23 (2%), 193 (16%), and 971 (82%) DBPs belonging to the high, moderate, and low priority groups, respectively. Lastly, the hTTR binding affinity of 4 DBPs was determined by an in vitro assay and all the tested DBPs exhibited dose-dependent binding potency with hTTR, which was consistent with the predicted result. Thus, more efforts should be performed to reveal the potential endocrine disruption of those high research-priority main-groups, subgroups, and individual DBPs.

Comparison of photoinduced and electrochemically induced degradation of venlafaxine

The European Union requires environmental monitoring of the antidepressant drug venlafaxine. Advanced oxidation processes provide a remedy against the spread of micropollutants. In this study, the photoinduced and electrochemical decompositions of venlafaxine were investigated in terms of mechanism and efficacy using high-performance liquid chromatography coupled to high-resolution multifragmentation mass spectrometry. Kinetic analysis, structure elucidation, matrix variation, and radical scavenging indicated the dominance of a hydroxyl-mediated indirect mechanism during photodegradation and hydroxyl and direct electrochemical oxidation for electrochemical degradation. Oxidants, sulfate, and chloride ions acted as accelerants, which reduced venlafaxine half-lives from 62 to 25 min. Humic acid decelerated degradation during ultra-violet irradiation up to 50%, but accelerated during electrochemical oxidation up to 56%. In silico quantitative structure activity relationship analysis predicted decreased environmental hazard after advanced oxidation process treatment. In general, photoirradiation proved more efficient due to faster decomposition and slightly less toxic transformation products. Yet, matrix effects would have to be carefully evaluated when potential applications as a fourth purification stage were to be considered.

Effects of endocrine active contaminating pesticides on RACK1 expression and immunological consequences in THP-1 cells

We have previously demonstrated that RACK1, which expression is under steroid hormone control, plays an important role in the activation of immune cells and its expression can be useful to evaluate the immunotoxic profile of endocrine disrupting chemicals (EDCs). Hence, we investigated the effects of three contaminating and persistent pesticides: the fungicide vinclozolin (VIN), the herbicide atrazine (ATR) and the insecticide cypermethrin (CYP) on RACK1 expression and on innate immune response. VIN resulted in modest alteration of RACK1 while ATR and CYP reduced in a dose dependent manner RACK1 expression, ultimately leading to the decrease in lipopolysaccharide-induced IL-8 and TNF-α release and CD86 and CD54 surface marker expression. Moreover, our data indicate that, after exposure to EDCs, alterations of RACK1 expression can also occur with mechanisms not directly mediated by an interaction with a nuclear or membrane steroid receptors. Therefore, RACK1 could represent a useful EDCs screening tool to evaluate their immunotoxic potential and to dissect their mechanisms of action.

In Silico Investigation of Some Compounds from the N-Butanol Extract of Centaurea tougourensis Boiss. & Reut

Bioinformatics as a newly emerging discipline is considered nowadays a reference to characterize the physicochemical and pharmacological properties of the actual biocompounds contained in plants, which has helped the pharmaceutical industry a lot in the drug development process. In this study, a bioinformatics approach known as in silico was performed to predict, for the first time, the physicochemical properties, ADMET profile, pharmacological capacities, cytotoxicity, and nervous system macromolecular targets, as well as the gene expression profiles, of four compounds recently identified from Centaurea tougourensis via the gas chromatography–mass spectrometry (GC–MS) approach. Thus, four compounds were tested from the n-butanol (n-BuOH) extract of this plant, named, respectively, Acridin-9-amine, 1,2,3,4-tetrahydro-5,7-dimethyl- (compound 1), 3-[2,3-Dihydro-2,2-dimethylbenzofuran-7-yl]-5-methoxy-1,3,4-oxadiazol-2(3H)-one (compound 2), 9,9-Dimethoxybicyclo[3.3.1]nona-2,4-dione (compound 3), and 3-[3-Bromophenyl]-7-chloro-3,4-dihydro-10-hydroxy-1,9(2H,10H)-acridinedione (compound 4). The insilico investigation revealed that the four tested compounds could be a good candidate to regulate the expression of key genes and may also exert significant cytotoxic effects against several tumor celllines. In addition, these compounds could also be effective in the treatment of some diseases related to diabetes, skin pathologies, cardiovascular, and central nervous system disorders. The bioactive compounds of plant remain the best alternative in the context of the drug discovery and development process.

Machine Learning Models for Predicting Cytotoxicity of Nanomaterials

The wide application of nanomaterials in consumer and medical products has raised concerns about their potential adverse effects on human health. Thus, more and more biological assessments regarding the toxicity of nanomaterials have been performed. However, the different ways the evaluations were performed, such as the utilized assays, cell lines, and the differences of the produced nanoparticles, make it difficult for scientists to analyze and effectively compare toxicities of nanomaterials. Fortunately, machine learning has emerged as a powerful tool for the prediction of nanotoxicity based on the available data. Among different types of toxicity assessments, nanomaterial cytotoxicity was the focus here because of the high sensitivity of cytotoxicity assessment to different treatments without the need for complicated and time-consuming procedures. In this review, we summarized recent studies that focused on the development of machine learning models for prediction of cytotoxicity of nanomaterials. The goal was to provide insight into predicting potential nanomaterial toxicity and promoting the development of safe nanomaterials.

Human transthyretin binding affinity of halogenated thiophenols and halogenated phenols: An in vitro and in silico study

Serious harmful effects have been reported for thiophenols, which are widely used industrial materials. To date, little information is available on whether such chemicals can elicit endocrine-related detrimental effects. Herein the potential binding affinity and underlying mechanism of action between human transthyretin (hTTR) and seven halogenated-thiophenols were examined experimentally and computationally. Experimental results indicated that the halogenated-thiophenols, except for pentafluorothiophenol, were powerful hTTR binders. The differentiated hTTR binding affinity of halogenated-thiophenols and halogenated-phenols were observed. The hTTR binding affinity of mono- and di-halo-thiophenols was higher than that of corresponding phenols; while the opposite relationship was observed for tri- and penta-halo-thiophenols and phenols. Our results also confirmed that the binding interactions were influenced by the degree of ligand dissociation. Molecular modeling results implied that the dominant noncovalent interactions in the molecular recognition processes between hTTR and halogenated-thiophenols were ionic pair, hydrogen bonds and hydrophobic interactions. Finally, a model with acceptable predictive ability was developed, which can be used to computationally predict the potential hTTR binding affinity of other halogenated-thiophenols and phenols. Taken together, our results highlighted that more research is needed to determine their potential endocrine-related harmful effects and appropriate management actions should be taken to promote their sustainable use.

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

An approach to classifying occupational exposures to endocrine disrupting chemicals by sex hormone function using an expert judgment process

Endocrine disrupting chemicals (EDCs) are exogenous substances that interfere with the endocrine system and cause adverse effects. We aimed to classify the effects of 24 known EDCs, prevalent in certain occupations, according to four modes of action (estrogenic, antiestrogenic, androgenic, and/or antiandrogenic). A literature search, stratified into four types of literature was conducted (namely: national and international agency reports; review articles; primary studies; ToxCast^TM). The state of the evidence of each EDC on sex hormone function was summarized and reviewed by an expert panel. For each mode of action, the experts evaluated the likelihood of endocrine disruption in five categories: "No", "Unlikely", "Possibly", "Probably", and "Yes". Seven agents were categorized as "Yes," or having strong evidence for their effects on sex hormone function (antiandrogenic: lead, arsenic, butylbenzyl phthalate, dibutyl phthalate, dicyclohexyl phthalate; estrogenic: nonylphenol, bisphenol A). Nine agents were categorized as "Probable," or having probable evidence (antiandrogenic: bis(2-ethylhexyl)phthalate, nonylphenol, toluene, bisphenol A, diisononyl phthalate; androgenic: cadmium; estrogenic: copper, cadmium and; anti-estrogenic: lead). Two agents (arsenic, polychlorinated biphenyls) had opposing conclusions supporting both "probably" estrogenic and antiestrogenic effects. This synthesis will allow researchers to evaluate the health effects of selected EDCs with an added level of precision related to the mode of action.

Development of a Nicotinic Acetylcholine Receptor nAChR α7 Binding Activity Prediction Model

Despite the well-known adverse health effects associated with tobacco use, addiction to nicotine found in tobacco products causes difficulty in quitting among users. Nicotinic acetylcholine receptors (nAChRs) are the physiological targets of nicotine and facilitate addiction to tobacco products. The nAChR-α7 subtype plays an important role in addiction; therefore, predicting the binding activity of tobacco constituents to nAChR-α7 is an important component for assessing addictive potential of tobacco constituents. We developed an α7 binding activity prediction model based on a large training data set of 843 chemicals with human α7 binding activity data extracted from PubChem and ChEMBL. The model was tested using 1215 chemicals with rat α7 binding activity data from the same databases. Based on the competitive docking results, the docking scores were partitioned to the key residues that play important roles in the receptor-ligand binding. A decision forest was used to train the human α7 binding activity prediction model based on the partition of docking scores. Five-fold cross validations were conducted to estimate the performance of the decision forest models. The developed model was used to predict the potential human α7 binding activity for 5275 tobacco constituents. The human α7 binding activity data for 84 of the 5275 tobacco constituents were experimentally measured to confirm and empirically validate the prediction results. The prediction accuracy, sensitivity, and specificity were 64.3, 40.0, and 81.6%, respectively. The developed prediction model of human α7 may be a useful tool for high-throughput screening of potential addictive tobacco constituents.

Binding interactions of halo-benzoic acids, halo-benzenesulfonic acids and halo-phenylboronic acids with human transthyretin

The anionic form-dependent binding interaction of halo-phenolic substances with human transthyretin (hTTR) has been observed previously. This indicates that ionizable compounds should be the primary focus in screening potential hTTR disruptors. Here, the potential binding potency of halo-benzoic acids, halo-benzenesulfonic acids/sulfates and halo-phenylboronic acids with hTTR was determined and analyzed by competitive fluorescence displacement assay integrated with computational methods. The laboratorial results indicated that the three test groups of model compounds exhibited a distinct binding affinity to hTTR. All the tested halo-phenylboronic acids, some of the tested halo-benzoic acids and halo-benzenesulfonic acids/sulfates were shown to be inactive with hTTR. Other halo-benzoic acids and halo-benzenesulfonic acids/sulfates were moderate and/or weak hTTR binders. The binding affinity of halo-benzoic acids and halo-benzenesulfonic acids/sulfates with hTTR was similar. The low distribution ability of the model compounds from water to hTTR may be the reason why they exhibited the binding potency observed with hTTR. By introducing other highly hydrophobic compounds, we observed that the binding affinity between compounds and hTTR increased with increasing molecular hydrophobicity. Those results indicated that the highly hydrophobic halo-benzoic acids and halo-benzenesulfonic acids/sulfates may be high-priority hTTR disruptors. Finally, a binary classification model was constructed employing three predictive variables. The sensitivity (S_n), specificity (S_p), predictive accuracy (Q) values of the training set and validation set were >0.83, indicating that the model had good classification performance. Thus, the binary classification model developed here could be used to distinguish whether a given ionizable compound is a potential hTTR binder or not.

Emerging Polar Phenolic Disinfection Byproducts Are High-Affinity Human Transthyretin Disruptors: An in Vitro and in Silico Study

Phenolic disinfection byproducts (phenolic-DBPs) have been identified in recent years. However, the toxicity data for phenolic-DBPs are scarce, hampering their risk assessment and the development of regulations on the acceptable concentration of phenolic-DBPs in water. In this study, the binding potency and underlying interaction mechanism between human transthyretin (hTTR) and five groups of representative phenolic-DBPs (2,4,6-trihalo-phenols, 2,6-dihalo-4-nitrophenols, 3,5-dihalo-4-hydroxybenzaldehydes, 3,5-dihalo-4-hydroxybenzoic acids, halo-salicylic acids) were determined and probed by competitive fluorescence displacement assay integrated with in silico methods. Experimental results implied that 2,4,6-trihalo-phenols, 2,6-dihalo-4-nitrophenols, and 3,5-dihalo-4-hydroxybenzaldehydes have a high binding affinity with hTTR. The hTTR binding potency of the chemicals with electron-withdrawing groups on their molecular structures were higher than that with electron-donor groups. Molecular modeling methods were used to decipher the binding mechanism between model compounds and hTTR. The results documented that ionic pair, hydrogen bonding and hydrophobic interactions were dominant interactions. Finally, a mechanism-based model for predicting the hTTR binding affinity was developed. The determination coefficient ( R²), leave-one-out cross validation Q² ( Q_LOO²), bootstrapping coefficient ( Q_BOOT²), external validation coefficient ( Q_EXT²) and concordance correlation coefficient ( CCC) of the developed model met the acceptable criteria ( Q² > 0.600, R² > 0.700, CCC > 0.850), implying that the model had good goodness-of-fit, robustness, and external prediction performances. All the results indicated that the phenolic-DBPs have the hTTR disrupting effects, and further studies are needed to investigate their other mechanism of endocrine disruption.