Next Article in Journal
Voluntary Physical Exercise Reduces Motor Dysfunction and Hampers Tumor Cell Proliferation in a Mouse Model of Glioma
Next Article in Special Issue
Development of a Human Estrogen Receptor Dimerization Assay for the Estrogenic Endocrine-Disrupting Chemicals Using Bioluminescence Resonance Energy Transfer
Previous Article in Journal
Dementia-Free Life Expectancy among People over 60 Years Old by Sex, Urban and Rural Areas in Jiangxi Province, China
Previous Article in Special Issue
Cross-Sectional and Longitudinal Effects of PCB Exposure on Human Stress Hormones in the German HELPcB Surveillance Program
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 17
Issue 16
10.3390/ijerph17165664
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Basis for Endocrine Disruption by Pesticides Targeting Aromatase and Estrogen Receptor

by
Chao Zhang
1,
Tiziana Schilirò
2,
Marta Gea
2,
Silvia Bianchi
1,
Angelo Spinello
3,
Alessandra Magistrato
3,
Gianfranco Gilardi
1 and
Giovanna Di Nardo
1,*
1
Department of Life Sciences and Systems Biology, University of Torino, 10123 Torino, Italy
2
Department of Public Health and Pediatrics, University of Torino, 10126 Torino, Italy
3
National Research Council-Institute of Materials (CNR-IOM) at International School for Advanced Studies (SISSA), 34165 Trieste, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(16), 5664; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17165664
Submission received: 10 July 2020 / Revised: 1 August 2020 / Accepted: 3 August 2020 / Published: 5 August 2020
(This article belongs to the Special Issue Endocrine Disruptor Chemicals and Estrogenic Activity in the Context of Health Risk Assessment)
Browse Figures
Versions Notes

Abstract

:
The intensive use of pesticides has led to their increasing presence in water, soil, and agricultural products. Mounting evidence indicates that some pesticides may be endocrine disrupting chemicals (EDCs), being therefore harmful for the human health and the environment. In this study, three pesticides, glyphosate, thiacloprid, and imidacloprid, were tested for their ability to interfere with estrogen biosynthesis and/or signaling, to evaluate their potential action as EDCs. Among the tested compounds, only glyphosate inhibited aromatase activity (up to 30%) via a non-competitive inhibition or a mixed inhibition mechanism depending on the concentration applied. Then, the ability of the three pesticides to induce an estrogenic activity was tested in MELN cells. When compared to 17β-estradiol, thiacloprid and imidacloprid induced an estrogenic activity at the highest concentrations tested with a relative potency of 5.4 × 10−10 and 3.7 × 10−9, respectively. Molecular dynamics and docking simulations predicted the potential binding sites and the binding mode of the three pesticides on the structure of the two key targets, providing a rational for their mechanism as EDCs. The results demonstrate that the three pesticides are potential EDCs as glyphosate acts as an aromatase inhibitor, whereas imidacloprid and thiacloprid can interfere with estrogen induced signaling.

Graphical Abstract

1. Introduction

Under the modern lifestyle, humans are exposed to various chemicals such as pesticide residuals in fruits and vegetables, antibiotics in meat and milk, preservatives in cosmetics and personal care products [1,2]. These chemicals are usually in low doses and may not have a short term significant impact on the human body, but they can cause long term damages to health [3]. The effects of low-dose compounds on human health are mainly related to the endocrine system [4,5,6]. These compounds can in fact mimic or influence the action of endogenous hormones through various mechanisms, being therefore referred to as endocrine disrupting chemicals (EDCs) [7,8].
The modern industrial and agricultural system relies heavily on pesticides. The compelling need of high food crop increasingly demands the use of chemicals. This results in the extensive applications of millions of tons of pesticides every year [9,10]. Among the many pesticides available on the market, glyphosate is one of the most widely used herbicides in the world. In plants, it affects the synthesis of essential aromatic amino acids by inhibiting the activity of 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) in the shikimate pathway [11]. As a result, glyphosate is considered harmless to mammals since they do not contain the EPSPS enzyme [12]. Conversely, neonicotinoids are currently the most widely used agricultural insecticides [13]. These selective insecticides specifically bind to the α -subunit of nicotinic acetylcholine receptors (nAChR), which is common in all insects [14,15]. Due to their difficulty to penetrate the blood-brain barrier, they should exert low toxicity in vertebrates. However, different direct, indirect, and suspect toxic effects of these compounds on vertebrate wildlife and human health have been reported [16,17,18]. Neonicotinoids were introduced in the market in the 1990s, therefore their sales volume has enormously proliferated in recent decades, exceeding 25% of the market in 2010 [19]. Imidacloprid has become the world’s best-selling insecticide next to glyphosate [20]. Although the European Union banned the outdoor application of imidacloprid, clothianidin, and thiamethoxam in 2018 due to a high risk for bees that are responsible for pollinating most crops worldwide [21], neonicotinoids remain the most extensively applied insecticides in the world [22].
The intensive use of glyphosate-based herbicides and neonicotinoid insecticides has caused the contamination of soil, water, air, and agricultural products [23,24,25,26,27,28]. The half-life of glyphosate in the field is usually 47 days and it is primarily degraded into aminomethylphosphonic acid (AMPA) and glyoxylic acid by soil microorganisms [29]. Neonicotinoids can exist for a long time in soil and water, and their half-life being as long as three years [22]. Given these observations and the increased use of pesticides, concerns are raising about their potential impact on human health and the environment. In particular, mounting evidence discloses that exposure to pesticides can affect the endocrine system. Toxicological and epidemiological studies indicate possible genotoxic and cytotoxic effects as well as birth defects and neurotoxicity in different cell lines and animal models [18,30].
Aiming at establishing the interference of pollutants with human health, we monitored the effect of glyphosate and two neonicotinoids on two critical targets of the endocrine system: Aromatase, the enzyme responsible for estrogen biosynthesis, and estrogen receptor (ER) alpha, the main protein promoting estrogen signaling.
Aromatase catalyzes the transformation of androgens to estrogens [31]. In vitro experiments showed that glyphosate causes changes of aromatase mRNA levels and activity in placental JEG3 cells and human HEK293 cells, thus interfering with steroid conversion to estradiol [32,33]. In addition, the expression of androgen and ERα was inhibited in glyphosate treated HepG2 cells [34,35] where also the transcription and activity of aromatase were altered [35]. Insecticides such as thiacloprid, thiamethoxam, and imidacloprid affect aromatase expression and activity in a co-culture model of fetoplacental steroidogenesis, increasing estrone and estradiol production, while estriol production is decreased [36]. Thiacloprid, imidacloprid, and thiamethoxam have also been demonstrated to increase aromatase expression and activity in H295R and Hs578t cells [37,38].
Here, a recombinant form of human aromatase is used to test the possible effect of the three pesticides on estrogen biosynthesis, gaining new information at molecular level [39,40,41]. Moreover, all-atom molecular dynamics simulations provided structural insights on the ability of these molecules to target the aromatase enzyme.
At a cellular level, the estrogen-regulated signaling is mainly due to estrogen receptors (ERs). The main ERs are the ERα and ERβ nuclear receptors and the G-protein-coupled estrogen receptor 1, which is a membrane receptor [42]. Estrogen receptor α (ERα) is a nuclear hormone receptor and a ligand-regulated transcription factor, which mediates the activity of estrogens in vital processes (i.e., reproduction, cardiovascular maintenance, bone density/remodeling). ERα is composed of five functional domains, among which the ligand-binding domain, activated upon estrogen binding, stimulates cell growth and proliferation. After the menopause, increased estrogen levels due to a deregulated activity of aromatase bind as an agonist to ERα, exerting a pro-oncogenic effect by either decreasing apoptosis or promoting cell proliferation [43]. Therefore, estrogen selective modulators have been developed and one of them, tamoxifen, is extensively used in adjuvant therapy of breast cancer for its ability to act as an ERα antagonist.
Previous studies have investigated the ability of glyphosate to induce an estrogenic activity mediated by ERs. The results of the study of Thongprakaisang et al., 2013 [44] demonstrated that glyphosate induces an ER-mediated estrogenic activity, mediated by ER activation, similar to 17β-estradiol (E2) on T47D-KBluc cells. However, recently another study, using the same cells, showed that this pesticide induces an estrogenic activity, mediated by ER activation, lower than E2 and that this activation is probably induced by a ligand-independent mechanism [45]. Moreover, additional studies on different transfected cells showed that glyphosate did not induce any ER-mediated estrogenic activity and did not produce any anti-estrogenic effect when tested in combination with E2 [35,46]. Therefore, whether glyphosate can trigger an ER-mediated estrogenic activity remains controversial. To our knowledge, only three studies have been performed on the estrogenic activity of imidacloprid and thiacloprid. In particular, in the study of Mesnage et al., 2018 [47] the proliferative effect of both pesticides was investigated on estrogen-sensitive cells, while Kojima et al., 2004 [46] and Westlund and Yargeau, 2017 [48] assessed the ER-mediated estrogenic activity of imidacloprid and thiacloprid on mammalian or yeast cells, respectively. The results of the three studies demonstrated that the two pesticides induce no proliferative effect and no estrogenic activity was mediated by ERs, while an anti-estrogenic activity was detected testing imidacloprid in combination with E2 on yeast cells.
In this work, in order to increase the knowledge on the estrogenic activity of glyphosate, imidacloprid, and thiacloprid, the gene reporter assay on estrogen-sensitive human breast cancer MCF-7 cells transfected with the ERE-βGlob-Luc-SVNeo plasmid (MELN cells) is applied to test the three pesticides for their possible ER-mediated estrogenic activity.
Possible additive and/or antagonist effects are also investigated. Moreover, docking simulations provide atomic level insights on the potential binding mode of these molecules to the primary ligand (estrogen) binding site as well as to a peripheral allosteric site which may be responsible for the experimental observed additive effect of the pollutant with the endogenous ligand.

2. Materials and Methods

2.1. Materials

All reagents are analytically pure by purchase from Sigma-Aldrich (St. Louis, MO, USA). Stock solutions of chemical compounds were prepared in absolute ethanol or dimethyl sulfoxide (DMSO). Before each experiment, the test sample was diluted into a fresh buffer solution, and the final organic solvent concentration was less than 0.1%. The recombinant human aromatase (Aro) and the human recombinant cytochrome P450 reductase (hCPR) were expressed and purified as previously described [39,49].

2.2. ELISA Assay

An estrone direct competitive ELISA kit (BioVendor, Brno, Czech Republic) was used to evaluate the effect of pesticides on aromatase activity. Different reaction mixtures were set up by mixing 5 nM Aro, 5 nM hCPR, 0.5 mM NADPH, 50 nM androstenedione, and three concentrations of pesticides (500, 1000, and 1500 nM) in a 100 mM potassium buffer (KPi) containing 20% glycerol, 1 mM β-mercaptoethanol at pH 7.0. Reactions were carried out for 10 min at 30 °C, heat-inactivated for 10 min at 90 °C, and centrifuged for 10 min at 11,000 rpm. After centrifugation, the supernatant was diluted 1:8 in the Calibrator A provided by the ELISA kit and the product estrone quantified performing ELISA according to the manufacturer’s instructions. Reactions in the presence of anastrozole or without hCPR were used as negative controls. The concentration of estrone was extrapolated from a calibration curve with known concentrations of estrone.
For the experiment where the catalytic parameters were derived, four substrate concentrations were applied in the reaction mixture (ranging from 25 to 250 nM) in the absence and presence of 1000 and 5000 nM of glyphosate.

2.3. Computational Studies

In order to explain the molecular terms for the action of glyphosate on the aromatase enzyme we docked it into the two possible allosteric sites previously identified [50]. Docking has been performed with the GLIDE software, release 2020-1 (Schrödinger, LLC, New York, NY, USA) using the single-precision protocol [51]. The two neonicotinoids, thiacloprid and imidacloprid, were instead docked into the ERαs active site, using as a starting structure the crystal structure of 17-β-estradiol (EST)-bound ERα dimer (PDB id: 1qku) [52]. In this structure, we have searched for putative allosteric pockets using the SiteMap algorithm [53].
In order to account for its flexibility of the receptor and since is flexibility resulted to be of paramount importance for the identification of novel allosteric inhibitors [54], we performed classical Molecular Dynamics (MD) simulations on the complex with the aromatase enzyme. We employed as a starting structure of our simulation the equilibrated enzyme model which was embedded in a mimic of a membrane bilayer by using the CHARMMGUI webserver [55]. This consisted of POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and 6 wt% of cholesterol (CHL) in order to mimic the endoplasmic reticulum membrane. Physiological protonation states were calculated with the webserver H++ [56]. Asp309 was considered in its neutral form consistently with other literature studies [57]. The glyphosate molecule was considered in the most likely protonation state at physiological pH. According to literature data, the first protonation of the molecule occurs on its phosphate group [58].
The Parm99SB AMBER force field (FF) [59,60] and lipid14 FF [61] were used for the protein and the lipids, respectively.
The Shahrokh et al. parameters were used for the heme moiety and Cys437 [62]. Simulations were done in the presence of the substrate androstenedione (ASD) in the active site and of glyphosate in the allosteric pockets for which the general Amber FF (GAFF) was employed [63]. For the organic ligands the electrostatic potential (ESP) charges [64] were calculated by performing geometry optimization of the substrates at the Hartree-Fock level of theory using a 6-31G* basis set with the Gaussian 09 software (Gaussian Inc., Wallingford, CT, USA) [65]. These were later transformed in RESP charges by using the Antechamber tool [66].
The system was then explicitly solvated using the TIP3P water model, leading to a total of 131,454 atoms. Topology, built with AmberTools 18, was later converted in a GROMACS format using the acpype algorithm [67]. MD simulations were performed with GROMACS 5.0.4 [68]. An integration time step of 2 fs was used and all covalent bonds involving hydrogen atoms were constrained with the LINCS algorithm. The Particle Mesh Ewald algorithm [69] was used in order to account for electrostatic interactions. Simulations were done in the isothermal-isobaric NPT ensemble, at a temperature of 300 K, using a velocity-rescaling thermostat [70]. Preliminary energy minimization was done with the steepest descend algorithm.
An initial equilibration of the membrane was performed for 100 ns with the protein atoms harmonically restrained with a force constant of 1000 kJ mol−1 nm−2, reaching a constant value (92 × 92 × 151 Å3) of the simulation box size. Constraints were then slowly released, and the system was thermalized to the target temperature of 300 K in about 10 ns. Then, the aromatase in complex with glyphosate was relaxed by performing a 100 ns MD simulation rescaling the motion of the center of mass of aromatase and the ligand, followed by an unbiased 100 ns MD simulation.

2.4. MELN cell Culture

MELN cells were provided by Dr. P. Balaguer (INSERM, Montpellier, France). They are estrogen-sensitive human breast cancer cells (MCF-7) transfected with the ERE-βGlob-Luc-SVNeo plasmid (ERE-βGlob-Luc-SVNeo) [71,72]. The integrated plasmid contains a luciferase reporter gene, the estrogen-responsive elements (ERE) and an antibiotic resistance selection gene (SVNeo). MELN cells were cultured at 37 °C and 5% CO2 in Dulbecco’s Modified Eagle’s Medium Nutrient Mixture F12-Ham (DMEM-F12), supplemented with phenol red, fetal bovine serum (FBS) (5% v/v), l-glutamine (4 mM), penicillin-streptomycin (100 U/mL–100 µg/mL), and G418 (1 mg/mL).

2.5. MELN Gene Reporter Assay

The assay was carried out as described by Balaguer et al., 1999 [73] with slight modifications [74]. For three days the cells were adapted to a test medium: DMEM-F12 without phenol red and supplemented with dextran-coated charcoal-treated FBS (5% v/v), l-glutamine (4 mM), and penicillin-streptomycin (100 U/mL-100 µg/mL). Then, the cells were seeded at a density of 40,000 cells/well, in 96-well plates (100 μL/well). After 24 h, the test medium of each well was replaced with a test medium containing pesticides (100 μL/wells), and the cells were incubated for 16 h. After the incubation, the luciferase activity was assessed adding 100 μL/well of the One Glo Reagent (One-Glo Luciferase Assay System, Promega, Madison, USA), mixing (5 min) and measuring the luminescence of each well by a luminometer (Infinite Reader M200 Pro, Tecan, Männedorf, Switzerland).
The stock solutions of thiacloprid and imidacloprid were prepared in DMSO, while the stock solution of glyphosate was prepared in a test medium. The stock solutions were stored at −20 °C and, shortly before exposure, different concentrations of pesticides were prepared in a test medium (glyphosate and imidacloprid: From 10−8 to 10−3 M; thiacloprid from 10−8 to 5 × 10−4 M, due to lower solubility). The final DMSO concentration was less than 0.1%. Cells exposed to the test medium were used as a negative control and five concentrations of E2 (from 10−12 to 10−8 M) were tested to obtain a standard positive curve of the reference compound (E2).
The estrogenic activity was calculated as the ratio of the activity induced by the treatment over the activity induced by the positive control with 17-βestradiol (E2). It was expressed in percentage considering the relative luciferase activity of E2 (10−8 M) as 100%. Since all experiments were performed in quadruplicate (four wells for each experimental condition), the estrogenic activity was expressed as the mean and standard deviation of four values. The estrogenic activity of pesticides was also evaluated by the determination of the relative potency of each pesticide in comparison with the reference compound (E2) and it was expressed as the E2 equivalency factor (EEF) [75]. The EEF was calculated using the concentrations of E2 and pesticides at which 50% of biological effect is achieved (EC50) through the formula: EEF = E2 EC50/pesticide EC50.
Three concentrations of pesticides (10−5, 2.5 × 10−4, and 5 × 10−4 M) were also tested: In combination with an ER-antagonist (tamoxifen 10−6 M), in order to confirm whether the observed effects were due to the ER activation, and in combination with E2 (10−10 M), in order to investigate the interaction between pesticides and E2 in MELN cells. The estrogenic activity of these treatments was expressed as relative luciferase activity and it was calculated as percentage of activity induced by the treatment with respect to the activity induced by the E2 10−10 M (relative luciferase activity of E2 10−10 M = 100%). The stock solutions of E2 and tamoxifen were prepared in ethanol and stored at −20 °C.

2.6. Data Analysis

Statistical analysis was performed using IBM SPSS Statistics 25.0 (IBM, Armonk, USA). The EC50 of E2 and pesticides was calculated by dose–response curves, which were estimated through a probit regression between the relative luciferase activity and Log transformed-concentrations of E2 or pesticides.
Data collected with the MELN gene reporter assay were not normally distributed, so the non-parametric Kruskal-Wallis test followed by the post-hoc Dunnett test was used to assess significant differences among the different experimental conditions. The differences were considered significant with p-value < 0.05.

3. Results

3.1. Effect of Pesticides on Aromatase Activity

In order to study the effect of pesticide compounds on aromatase activity, a direct competitive estrone ELISA was performed using the purified cytochrome P450 reductase (CPR), as an electron donor from NADPH, and aromatase. The aromatase activity was evaluated by measuring the estrone production in the absence and presence of three different concentrations of pesticides (0.5, 1, and 5 μM). As a control, anastrozole, a known aromatase inhibitor, was applied at a concentration of 1000 nM and the residual aromatase activity detected was 0.7%. As can be seen in Table 1, glyphosate partially reduced the aromatase activity at the concentrations tested. The enzyme activity decreased with the increase of glyphosate concentration. When adding 5 μM of glyphosate, the residual aromatase activity was 36%. Unlike glyphosate, imidacloprid and thiacloprid did not inhibit the enzyme activity (Table 1).

3.2. Effect of Glyphosate Concentration on Aromatase Activity

The effect of glyphosate on aromatase activity was further studied by exploring the concentration range of glyphosate applied from 50 to 1500 nM. Such concentrations of glyphosate were selected since they resemble the ones detected in human urine samples [76]. The experiment was carried out by the ELISA assay at the concentration of 50 and 400 nM androstenedione, respectively. The two different concentrations were chosen on the basis of the kinetic parameters of aromatase: The first one (50 nM) is close to the enzyme KM and the second one (400 nM) is saturating the enzyme (see next paragraph).
As shown in Figure 1, the activity of Aro is inhibited by 30% when the glyphosate concentration is ≥1000 nM. However, when the substrate concentration is 50 nM (black squares in Figure 1), the maximal inhibitory effect is already achieved when the glyphosate concentration is 100 nM. Therefore, the inhibitory effect of glyphosate strongly depends on the substrate concentration and it is only partial, indicating that this compound can be considered as a weak inhibitor.

3.3. Effect of Glyphosate on the Catalytic Parameters of Aro

In order to investigate the mechanism of aromatase inhibition by glyphosate, the kinetic parameters of the enzyme were evaluated using the estrone ELISA assay in the absence and presence of two different concentrations of the pesticide (1000 and 5000 nM). Different substrate concentrations were applied and the product formation rate was plotted as a function of the substrate concentration (Figure 2A). The plot showed hyperbolic trends and the catalytic parameters, shown in Table 2, were obtained by fitting the experimental data to the Michaelis-Menten equation.
When glyphosate is not present, the resulting KM and Vmax are 41.3 ± 8.2 nM and 0.018 ± 0.001 min−1, respectively. When 1000 nM of glyphosate was added, the Vmax value was significantly decreased whereas the KM value was not significantly affected. Interestingly, when the glyphosate concentration was increased to 5000 nM, both KM and Vmax were affected. Compared to the reaction without glyphosate, KM was increased by 2.2 folds, while Vmax was decreased to 0.011 ± 0.002 min−1 (Table 2).
The kinetic parameters show that the type of inhibition of aromatase by glyphosate depends on the herbicide concentration applied. Indeed, when using 1000 nM of glyphosate, the Vmax was decreased while KM did not change, indicating a non-competitive inhibition mechanism, meaning that, at this concentration, glyphosate does not compete with the substrate and binds to a site different from that where the substrate binds. When the concentration of glyphosate was increased to 5000 nM, both KM and Vmax were affected, and the Lineweaver-Burk plot shows a trend typical of a mixed inhibition mechanism (Figure 2B). Mixed inhibition is considered a more general case of non-competitive inhibition, in which the inhibitor exhibits unequal affinity for the free enzyme and for the enzyme-substrate complex.

3.4. Molecular Dynamics Simulations on Aromatase

Grounding on recent evidence demonstrating the existence of allosteric binding sites [50] and their possible exploitation for a non-competitive/mixed inhibition mechanism [54] we docked glyphosate into the two allosteric cavities. Namely, we docked it to Site 1, which lies along the most relevant access channel to the enzyme active site [77] and to Site 2, which instead lies at the interface with the cytochrome P450 reductase (CPR), supplying the electrons necessary for catalysis (Figure 3) [78]. In the docking pose in Site 1 and during MD simulations glyphosate engages a salt bridge interaction with its phosphate group and Arg192, as well as the formation of a hydrogen (H)-bond between Gln218 and the carboxylic group of the pesticide. Most importantly, the phosphate group of glyphosate makes up to two simultaneous H-bonds with Asp309 (Figure 3), which normally is engaged in stabilizing the binding of aromatase substrates.
Due to these interactions the molecule remains stably bound in the pocket for the whole MD simulation, in line with its inhibitory activity in the µM range. Remarkably, it was recently suggested that the binding of a small molecule in Site 1 triggered the displacement of the water molecule needed for the catalytic activity, which are normally H-bonding with the Asp309 and Arg192 residues, both being critical residues for the catalytic activity [41,57], thus inhibiting estrogen biosynthesis. Conversely, the docking pose obtained in Site 2, did not establish any relevant H-bond/salt bridge. As a result, the glyphosate dissociated from the pocket within the first few ns of MD simulations.

3.5. Detection of Estrogenic Activity with the MELN Gene Reporter Assay

The MELN gene reporter assay was carried out to evaluate the estrogenic activity of glyphosate, imidacloprid, and thiacloprid on MELN cells. In this study, different concentrations of each pesticide were tested. Concentrations similar to the pesticide levels measured in human urine [76] were selected as the lowest concentrations (10−8 and 10−7 M), while concentrations up to 10−3 M were selected as the highest concentrations, in order to assess the effect induced by each pesticide in a wide range of concentrations.
Our results showed that glyphosate did not increase the relative luciferase activity with respect to the negative control; therefore no estrogenic activity was detected testing this pesticide on MELN cells (Figure 4). On the contrary, the highest concentrations of glyphosate induced a small decrease of the relative luciferase activity, which may be due to a toxic effect of the pesticide on cells.
The null estrogenic activity induced by glyphosate was confirmed also by the exposure of cells to glyphosate in combination with tamoxifen (ER-antagonist) or with E2. Indeed, the wells treated with glyphosate and tamoxifen showed a relative luciferase activity equal to the wells treated with tamoxifen alone, and the wells treated with glyphosate and E2 showed a relative luciferase activity equal to the wells treated with E2 alone (data not shown). These results suggest that glyphosate does not interfere with the binding between ER and E2.
Regarding the neonicotinoid pesticides, in the present study, imidacloprid and thiacloprid significantly increased the relative luciferase activity with respect to the negative control, starting from 6.3 × 10−5 M (Log[nM] = 4.79) and 10−6 M (Log[nM] = 3), respectively (Kruskal-Wallis test followed by the post-hoc Dunnett test, p < 0.05). Since both the neonicotinoid pesticides induced a dose-dependent increase of the relative luciferase activity, in particular from 10−4 M (Log[nM] = 5) to the highest tested dose (Figure 4), a significant estrogenic activity of these two pesticides was detected on MELN cells.
The estrogenic activity of imidacloprid and thiacloprid was also quantitatively evaluated by the estimate of the concentrations of E2 and pesticides at which 50% of biological effect is achieved (EC50) and the E2 equivalency factor (EEF). The EC50 of E2 and pesticides was calculated by dose-response curves whereas the EEF was calculated through the formula: EEF = E2 EC50/pesticide EC50. The EC50 of imidacloprid and thiacloprid was 1.0 × 10−2 M (IC 95% 1.7 × 10−3–2.2 × 10−1 M) and 1.5 × 10−3 M (IC 95% 2.5 × 10−4–3.6 × 10−2 M), respectively, while the EEF was 5.4 × 10−10 (IC 95% 3.3 × 10−9–2.5 × 10−11) and 3.7 × 10−9 (IC 95% 2.2 × 10−8–1.5 × 10−10), respectively.
The exposure of cells to imidacloprid and thiacloprid in combination with tamoxifen confirmed that the estrogenic activity of the two pesticides was induced by the activation of ER. Indeed, the wells treated with the neonicotinoid pesticides and tamoxifen showed a relative luciferase activity that was lower compared to the wells treated with these pesticides alone. Furthermore, the relative luciferase activity of wells treated with the neonicotinoid pesticides and tamoxifen was similar to the relative luciferase activity measured in the negative control (Figure 5A,C).
Finally, the exposure of cells to imidacloprid and thiacloprid in combination with E2 induced an increase of the relative luciferase activity with respect to the E2 alone (Figure 5B,D). The increase was slight for imidacloprid while it was stronger for thiacloprid, suggesting a possible additive effect exerted by these pesticides in combination with E2.

3.6. Docking Calculation on ERα

In order to provide a rationale for the estrogenic activity exerted by the two neonicotinoids, thiacloprid and imidacloprid, we have performed docking calculations on ERα. First, the two molecules were docked into the ERαs estrogen binding site, using the crystal structure of 17-β-estradiol bound to the ERα dimer (PDB id: 1qku) [52]. Both neonicotinoids fit inside the estrogen binding pocket (Figure 6). In particular, imidacloprid forms a H-bond with the backbone of Gly521, while in thiacloprid, the Cl atom makes halogen bonds with the guanidinium group of Arg394 and the aromatic rings of Phe404 and Trp393. Halogen bonds are attractive interactions between the electrophilic region associated with the Cl halogen atom and the nucleophilic regions of the surrounding protein residues [79].
Next, in order to disclose if and how the two neonicotinoids exert an additive effect to estrogen binding, by occupying an allosteric cavity, we looked for the presence of druggable allosteric pockets in the protein. Interestingly, a high-ranking binding pocket was found in the proximity of the estrogen binding site of one monomer of ERα (Figure 6). The docking calculation performed on this pocket, strikingly revealed that imidacloprid can H-bond with Lys449 and Glu323. As well thiacloprid H-bonds with Lys449 and Trp393.

4. Discussion

In this study, the effect of three pesticides on two key targets of the endocrine system was evaluated using a combination of experimental and in silico methods that allowed an investigation at the molecular level.
The first target considered was aromatase, the key enzyme for estrogen production that resulted in being partially inhibited by glyphosate. The increase obtained in the aromatase activity in the presence of imidacloprid and thiacloprid, even if not statistically significant, can be interesting to be further investigated. Indeed, these two compounds could directly act on aromatase as allosteric activators or they could exert their action indirectly on CPR that has an essential role in catalysis and acts as an effector on human aromatase conformation [80,81].
Glyphosate is the most widely used active compound among herbicides and was already reported to affect aromatase expression and activity in cells [33,35]. The impact for human health has already been demonstrated: Indeed, an alteration in androgens/estrogens balance due to a lower aromatase expression as a consequence of the glyphosate presence changed the sperm nuclear quality impacting mammalian reproduction [82].
However, the direct interaction between glyphosate and the enzyme together with the inhibition mechanism was not reported yet. In this study, the ELISA assay revealed that the inhibition of aromatase by glyphosate is partial and weak and strongly depends on the substrate concentration. Moreover, at the lower concentration of glyphosate tested (1 μM), the inhibition was found to be non-competitive, while at the higher concentration used (5 μM), the inhibition turned into a mixed inhibition mode. These data suggest that glyphosate binds to an allosteric site both when the enzyme is free and in complex with the substrate. Classical MD simulations supply a structural model for glyphosate binding to aromatase, explaining at the atomic-level that this pollutant may exert its inhibitory activity at a low concentration by binding to the allosteric Site 1 previously identified by Magistrato et al. [50] and demonstrated to bind small-molecules inhibitors [54]. When the concentration increases, a mixed inhibition is observed again compatible with the presence of an allosteric site. The binding of glyphosate probably becomes stronger or another allosteric site is occupied by the inhibitor. However, MD simulations showed that glyphosate did not stably bind to the second allosteric pocket (Site 2) identified in previous studies.
The second important target studied was the estrogen receptor that is responsible for estrogen binding and signal transduction in cells. MCF-7 cells, stably transfected with an estrogen-regulated luciferase gene (MELN cells), were used to assess the estrogenic activity of the pesticides. Glyphosate induced no estrogenic activity on MELN cells, moreover glyphosate did not change the effect induced by tamoxifen or E2. These results are in accordance with the results of the studies of Kojima et al., 2004 [46] and Gasnier et al., 2009 [35], who did not find any agonistic or antagonistic effect of glyphosate on ER. On the contrary, our results are different than the results of Thongprakaisang et al., 2013 [44] and Mesnage et al., 2017 [45] who found a significant increase of ER-induced estrogenic activity. A possible explanation of this discrepancy could be that in our study the estrogenic activity was evaluated on MELN cells while the other two studies applied T47D-KBluc cells, which might be more sensitive to the glyphosate activity than MELN cells. Moreover, our results are consistent with the results of the Endocrine Disruptor Screening Program (EDSP) conducted by the United States Environmental Protection Agency (US EPA, 2015) which concluded that “there is no convincing evidence of a potential interaction with the estrogen pathway for glyphosate”.
In the present study, both imidacloprid and thiacloprid induced a dose-response estrogenic activity mediated by the ER activation, in particular starting from 10−4 M. A consistent result was found in the study of Kojima et al., 2004 [46] in which no estrogenic activity was detected exposing transfected-cells to imidacloprid at concentrations lower than 10−5 M; a contradicting result was found by Westlund and Yargeau, 2017 [48], who did not find any estrogenic activity testing thiacloprid at concentrations comparable to ours. The conflicting results of thiacloprid should be interpreted considering the different cell models applied in our study with respect to the study of Westlund and Yargeau, 2017 (i.e., the present study was performed on mammalian cells while the previous one on yeast cells) [48]. Indeed, yeast cells are characterized by a different membrane permeability, transport proteins, and signal transduction pathways with respect to mammalian cells which may have influenced the results [42], as reported before also for other nuclear receptors [83].
The estrogenic activity of imidacloprid and thiacloprid was observed at concentrations higher than the pesticide levels measured in human biological samples [76,84]. Although such high levels of pesticides are not found in human biological samples, low doses of these pesticides should not be considered harmless. Indeed, in biological fluids these pesticides may be present in combination with other EDCs with a similar action mode and, thus, these pesticides together with other EDCs might cause an overall estrogenic effect. Moreover, these high concentrations of pesticides are similar to levels in some environmental matrices, where they could cause adverse health effects on wildlife animals. For example, in guttation water or morning dew of plants, imidacloprid was found to be present at a concentration up to 346 mg/L [85].
Interestingly, the results of the present study showed that, when the ER-antagonist tamoxifen was added to these pesticides, the estrogenic activity was still higher than the negative control. This result can be explained by the competition between each pesticide and tamoxifen for ER binding. Moreover, in the presence of high concentrations of imidacloprid and thiacloprid, an additive effect with E2 was also observed. The docking calculation suggested that these two neonicotinoids may bind to both the orthosteric and allosteric pockets of ERα, suggesting a putative mechanism to rationalize their observed estrogenic activity.

5. Conclusions

In conclusion, this study provides further evidence about the action of some pesticides as endocrine disrupting chemicals (EDCs) targeting important proteins of the endocrine system. In particular, it shows that the inhibitory effects of the three compounds tested on aromatase are partial and their estrogenic effects occur at relatively high concentrations. However, possible additive estrogenic effects with the physiological hormone 17-βestradiol are present. Furthermore, it has to be taken into account that pesticides are usually introduced in the environment with their co-formulants that can be also biologically active as EDCs. Previous studies already showed that aromatase is inhibited by co-formulants of glyphosate-based herbicides [33,86]. Moreover, more studies are needed to investigate a possible additive effect of different pesticides that can be contemporarily present in the environment.
Our study also provides an integrated approach based on different assays and computational methods that allowed gaining new information about the possible interaction of pesticides and key targets at a molecular level. Such information can be exploited to predict the possible impact of other compounds on estrogen production and signaling in order to develop safer compounds for human health and environment.

Author Contributions

Conceptualization, T.S., A.S., A.M., G.G., and G.D.N.; data curation, C.Z., T.S., M.G., A.M., and G.D.N.; formal analysis, C.Z., T.S., M.G., A.S., A.M., and G.D.N.; funding acquisition, A.M., G.G., and G.D.N.; investigation, C.Z., T.S., M.G., S.B., A.S., A.M., G.G., and G.D.N.; methodology, C.Z., M.G., S.B., A.S., and A.M.; project administration, G.D.N.; resources, T.S., A.M., and G.G.; software, C.Z., M.G., S.B., A.S., and A.M.; supervision, T.S., A.M., G.G., and G.D.N.; validation, C.Z., M.G., and S.B.; visualization, C.Z., M.G., and A.S.; writing—original draft, C.Z. and G.D.N.; writing—review and editing, C.Z., T.S., M.G., A.S., A.M., G.G., and G.D.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the CRT Foundation grant (project “Exposome” RF = 2016.2780) to G.D.N., by the Italian Association for Cancer Research (MFAG 17134), “Mario and Valeria Rindi” fellowship for Italy, and by the “Against bRain cancEr: Finding personalized therapies with in silico and in vitro strategies” (ARES) CUP: D93D19000020007 POR FESR 2014 2020—1.3.b—Friuli Venezia Giulia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mie, A.; Andersen, H.R.; Gunnarsson, S.; Kahl, J.; Kesse-Guyot, E.; Rembiałkowska, E.; Quaglio, G.; Grandjean, P. Human health implications of organic food and organic agriculture: A comprehensive review. Environ. Health 2017, 16, 111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Prüss-Ustün, A.; Vickers, C.; Haefliger, P.; Bertollini, R. Knowns and unknowns on burden of disease due to chemicals: A systematic review. Environ. Health 2011, 10, 9. [Google Scholar] [CrossRef] [Green Version]
  3. Kortenkamp, A.; Faust, M.; Scholze, M.; Backhaus, T. Low-level exposure to multiple chemicals: Reason for human health concerns? Environ. Health Perspect. 2007, 115, 106–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Yang, M.; Park, M.S.; Lee, H.S. Endocrine disrupting chemicals: Human exposure and health Risks. J. Environ. Sci. Health C Environ. Carcinog. Ecotoxicol. Rev. 2006, 24, 183–224. [Google Scholar] [CrossRef] [PubMed]
  5. Schug, T.T.; Janesick, A.; Blumberg, B.; Heindel, J.J. Endocrine disrupting chemicals and disease susceptibility. J. Steroid Biochem. Mol. Biol. 2011, 127, 204–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Vandenberg, L.N.; Colborn, T.; Hayes, T.B.; Heindel, J.J.; Jacobs, D.R.; Lee, D.H.; Shioda, T.; Soto, A.M.; Vom, S.F.; Welshons, W.V.; et al. Hormones and endocrine-disrupting chemicals: Low-dose effects and nonmonotonic dose responses. Endocr. Rev. 2012, 33, 378–455. [Google Scholar] [CrossRef]
  7. Takayanagi, S.; Tokunaga, T.; Liu, X.; Okada, H.; Matsushima, A.; Shimohigashi, Y. Endocrine disruptor bisphenol A strongly binds to human estrogen-related receptor γ (ERRγ) with high constitutive activity. Toxicol. Lett. 2006, 167, 95–105. [Google Scholar] [CrossRef]
  8. Mnif, W.; Hassine, A.I.H.; Bouaziz, A.; Bartegi, A.; Thomas, O.; Roig, B. Effect of endocrine disruptor Pesticides: A Review. Int. J. Environ. Res. Public Health 2011, 8, 2265–2303. [Google Scholar] [CrossRef] [Green Version]
  9. Bernhardt, E.S.; Rosi, E.J.; Gessner, M.O. Synthetic chemicals as agents of global change. Front. Ecol. Environ. 2017, 15, 84–90. [Google Scholar] [CrossRef]
  10. Storck, V.; Karpouzas, D.G.; Martin-Laurent, F. Towards a better pesticide policy for the European Union. Sci. Total Environ. 2017, 575, 1027–1033. [Google Scholar] [CrossRef]
  11. Amrhein, N.; Deus, B.; Gehrke, P.; Steinrücken, H.C. The Site of the Inhibition of the Shikimate Pathway by glyphosate: ii. interference of glyphosate with chorismate formation in vivo and in vitro. Plant Physiol. 1980, 66, 830–834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Herrmann, K.M.; Weaver, L.M. The Shikimate Pathway. Annu. Rev. Plant Phys. 1999, 50, 473–503. [Google Scholar] [CrossRef] [PubMed]
  13. Simon-Delso, N.; Amaral-Rogers, V.; Belzunces, L.P.; Bonmatin, J.M.; Chagnon, M.; Downs, C.; Furlan, L.; Gibbons, D.W.; Giorio, C.; Girolami, V.; et al. Systemic insecticides (neonicotinoids and fipronil): Trends, uses, mode of action and metabolites. Environ. Sci. Pollut. Res. Int. 2015, 22, 5–34. [Google Scholar] [CrossRef] [PubMed]
  14. Matsuda, K.; Ihara, M.; Sattelle, D.B. Neonicotinoid Insecticides: Molecular Targets, Resistance, and Toxicity. Annu. Rev. Pharmacol. Toxicol. 2020, 60, 241–255. [Google Scholar] [CrossRef] [PubMed]
  15. Sánchez-Bayo, F.; Tennekes, H.A. Time-Cumulative Toxicity of Neonicotinoids: Experimental Evidence and Implications for Environmental Risk Assessments. Int. J. Environ. Res. Public Health 2020, 17, 1629. [Google Scholar] [CrossRef] [Green Version]
  16. Mensah, P.K.; Palmer, C.G.; Odume, O.N. Ecotoxicology of glyphosate and glyphosate-based herbicides — toxicity to wildlife and humans. In Toxicity and Hazard of Agrochemicals; Larramendy, M.L., Soloneski, S., Eds.; Intech: Rijeka, Croatia, 2015; pp. 281–304. [Google Scholar]
  17. Gibbons, D.; Morrissey, C.; Mineau, P. A review of the direct and indirect effects of neonicotinoids and fipronil on vertebrate wildlife. Environ. Sci. Pollut. Res. 2015, 22, 103–118. [Google Scholar] [CrossRef] [Green Version]
  18. Thompson, D.A.; Lehmler, H.-J.; Kolpin, D.W.; Hladik, M.L.; Vargo, J.D.; Schilling, K.E.; LeFevre, G.H.; Peeples, T.L.; Poch, M.C.; LaDuca, L.E.; et al. A critical review on the potential impacts of neonicotinoid insecticide use: Current knowledge of environmental fate, toxicity, and implications for human health. Environ. Sci. Processes Impacts 2020, 22, 1315–1346. [Google Scholar] [CrossRef]
  19. Bass, C.; Denholm, I.; Williamson, M.S.; Nauen, R. The global status of insect resistance to neonicotinoid insecticides. Pestic. Biochem. Physiol. 2015, 121, 78–87. [Google Scholar] [CrossRef] [Green Version]
  20. Jeschke, P.; Nauen, R.; Schindler, M.; Elbert, A. Overview of the status and global strategy for neonicotinoids. J. Agric. Food Chem. 2011, 59, 2897–2908. [Google Scholar] [CrossRef]
  21. The European Commission. Legislation 132. Off. J. Eur. Union 2018, 61, 31–40. [Google Scholar]
  22. Zhang, Q.; Li, Z.; Chang, C.H.; Lou, J.L.; Zhao, M.R.; Lu, C. Potential human exposures to neonicotinoid insecticides: A review. Environ. Pollut. 2018, 236, 71–81. [Google Scholar] [CrossRef] [PubMed]
  23. Annett, R.; Habibi, H.R.; Hontela, A. Impact of glyphosate and glyphosate-based herbicides on the freshwater environment: Impact of glyphosate-based herbicides. J. Appl. Toxicol. 2014, 34, 458–479. [Google Scholar] [CrossRef] [PubMed]
  24. Bai, S.H.; Ogbourne, S.M. Glyphosate: Environmental contamination, toxicity and potential risks to human health via food contamination. Environ. Sci. Pollut. Res. 2016, 23, 18988–19001. [Google Scholar] [CrossRef] [PubMed]
  25. Van Bruggen, A.H.C.; He, M.M.; Shin, K.; Mai, V.; Jeong, K.C.; Finckh, M.R.; Morris, J.G. Environmental and health effects of the herbicide glyphosate. Sci. Total Environ. 2018, 616–617, 255–268. [Google Scholar] [CrossRef]
  26. Morrissey, C.A.; Mineau, P.; Devries, J.H.; Sanchez-Bayo, F.; Liess, M.; Cavallaro, M.C.; Liber, K. Neonicotinoid contamination of global surface waters and associated risk to aquatic invertebrates: A review. Environ. Int. 2015, 74, 291–303. [Google Scholar] [CrossRef]
  27. Lu, C.; Chang, C.-H.; Palmer, C.; Zhao, M.; Zhang, Q. Neonicotinoid residues in fruits and vegetables: An integrated dietary exposure assessment approach. Environ. Sci. Technol. 2018, 52, 3175–3184. [Google Scholar] [CrossRef]
  28. Hladik, M.L.; Main, A.R.; Goulson, D. Environmental risks and challenges associated with neonicotinoid insecticides. Environ. Sci. Technol. 2018, 52, 3329–3335. [Google Scholar] [CrossRef] [Green Version]
  29. Vencill, W.K. Herbicide Handbook, 8th ed.; Vencill, W.K., Ed.; Weed Science Society of America: Lawrence, KS, USA, 2002; pp. 231–234. [Google Scholar]
  30. Agostini, L.P.; Dettogni, R.S.; dos Reis, R.S.; Stur, E.; dos Santos, E.V.W.; Ventorim, D.P.; Garcia, F.M.; Cardoso, R.C.; Graceli, J.B.; Louro, I.D. Effects of glyphosate exposure on human health: Insights from epidemiological and in vitro studies. Sci. Total Environ. 2020, 705, 135808. [Google Scholar] [CrossRef]
  31. Di Nardo, G.; Gilardi, G. Human aromatase: Perspectives in biochemistry and biotechnology: Human Aromatase. Biotechnol. Appl. Biochem. 2013, 60, 92–101. [Google Scholar] [CrossRef]
  32. Benachour, N.; Sipahutar, H.; Moslemi, S.; Gasnier, C.; Travert, C.; Séralini, G.E. Time- and dose-dependent effects of roundup on human embryonic and placental Cells. Arch. Environ. Contam. Toxicol. 2007, 53, 126–133. [Google Scholar] [CrossRef]
  33. Richard, S.; Moslemi, S.; Sipahutar, H.; Benachour, N.; Seralini, G.-E. Differential effects of glyphosate and roundup on human placental cells and aromatase. Environ. Health Perspect. 2005, 113, 716–720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Benachour, N.; Séralini, G.E. Glyphosate formulations induce apoptosis and necrosis in human umbilical, embryonic, and placental cells. Chem. Res. Toxicol. 2009, 22, 97–105. [Google Scholar] [CrossRef] [PubMed]
  35. Gasnier, C.; Dumont, C.; Benachour, N.; Clair, E.; Chagnon, M.C.; Séralini, G.E. Glyphosate-based herbicides are toxic and endocrine disruptors in human cell lines. Toxicology 2009, 262, 184–191. [Google Scholar] [CrossRef]
  36. Caron-Beaudoin, E.; Viau, R.; Hudon-Thibeault, A.A.; Vaillancourt, C.; Sanderson, J.T. The use of a unique co-culture model of fetoplacental steroidogenesis as a screening tool for endocrine disruptors: The effects of neonicotinoids on aromatase activity and hormone production. Toxicol. Appl. Pharmacol. 2017, 332, 15–24. [Google Scholar] [CrossRef] [PubMed]
  37. Caron-Beaudoin, É.; Denison, M.S.; Sanderson, J.T. Effects of neonicotinoids on promoter-specific expression and activity of aromatase (CYP19) in human adrenocortical carcinoma (H295R) and primary umbilical vein endothelial (HUVEC) cells. Toxicol. Sci. 2016, 149, 134–144. [Google Scholar] [CrossRef] [Green Version]
  38. Caron-Beaudoin, É.; Viau, R.; Sanderson, J.T. Effects of neonicotinoid pesticides on promoter-specific aromatase (CYP19) expression in Hs578t breast cancer cells and the role of the VEGF pathway. Environ. Health Perspect. 2018, 126, 047014. [Google Scholar] [CrossRef]
  39. Di Nardo, G.; Breitner, M.; Sadeghi, S.J.; Castrignanò, S.; Mei, G.; Di Venere, A.; Nicolai, E.; Allegra, P.; Gilardi, G. Dynamics and flexibility of human aromatase probed by FTIR and time resolved fluorescence spectroscopy. PLoS ONE 2013, 8, e82118. [Google Scholar] [CrossRef] [Green Version]
  40. Baravalle, R.; Ciaramella, A.; Baj, F.; Di Nardo, G.; Gilardi, G. Identification of endocrine disrupting chemicals acting on human aromatase. Biochim. Biophys. Acta - Proteins Proteomics 2018, 1866, 88–96. [Google Scholar] [CrossRef]
  41. Lo, J.; Di Nardo, G.; Griswold, J.; Egbuta, C.; Jiang, W.; Gilardi, G.; Ghosh, D. Structural basis for the functional roles of critical residues in human cytochrome P450 aromatase. Biochemistry 2013, 52, 5821–5829. [Google Scholar]
  42. Kiyama, R.; Wada-Kiyama, Y. Estrogenic endocrine disruptors: Molecular mechanisms of action. Environ. Int. 2015, 83, 11–40. [Google Scholar] [CrossRef]
  43. Pavlin, M.; Spinello, A.; Pennati, M.; Zaffaroni, N.; Gobbi, S.; Bisi, A.; Colombo, G.; Magistrato, A. A computational assay of estrogen receptor α antagonists reveals the key common structural traits of drugs effectively fighting refractory breast Cancers. Sci. Rep. 2018, 8, 649. [Google Scholar] [CrossRef] [PubMed]
  44. Thongprakaisang, S.; Thiantanawat, A.; Rangkadilok, N.; Suriyo, T.; Satayavivad, J. Glyphosate induces human breast cancer cells growth via estrogen receptors. Food Chem. Toxicol. 2013, 59, 129–136. [Google Scholar] [CrossRef] [PubMed]
  45. Mesnage, R.; Phedonos, A.; Biserni, M.; Arno, M.; Balu, S.; Corton, J.C.; Ugarte, R.; Antoniou, M.N. Evaluation of estrogen receptor alpha activation by glyphosate-based herbicide constituents. Food Chem. Toxicol. 2017, 108, 30–42. [Google Scholar] [CrossRef] [Green Version]
  46. Kojima, H.; Katsura, E.; Takeuchi, S.; Niiyama, K.; Kobayashi, K. Screening for estrogen and androgen receptor activities in 200 pesticides by in vitro reporter gene assays using Chinese hamster ovary cells. Environ. Health Perspect. 2004, 112, 524–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Mesnage, R.; Biserni, M.; Genkova, D.; Wesolowski, L.; Antoniou, M.N. Evaluation of neonicotinoid insecticides for oestrogenic, thyroidogenic and adipogenic activity reveals imidacloprid causes lipid accumulation. J. Appl. Toxicol. 2018, 38, 1483–1491. [Google Scholar] [CrossRef]
  48. Westlund, P.; Yargeau, V. Investigation of the presence and endocrine activities of pesticides found in wastewater effluent using yeast-based bioassays. Sci. Total Environ. 2017, 607–608, 744–751. [Google Scholar] [CrossRef]
  49. Huang, N.; Pandey, A.V.; Agrawal, V.; Reardon, W.; Lapunzina, P.D.; Mowat, D.; Jabs, E.W.; Vliet, G.V.; Sack, J.; Flück, C.E.; et al. Diversity and Function of Mutations in P450 Oxidoreductase in Patients with Antley-Bixler Syndrome and Disordered Steroidogenesis. Am. J. Hum. Genet. 2005, 76, 729–749. [Google Scholar] [CrossRef] [Green Version]
  50. Sgrignani, J.; Bon, M.; Colombo, G.; Magistrato, A. Computational approaches elucidate the allosteric mechanism of human aromatase inhibition: A novel possible route to small-molecule regulation of CYP450s activities? J. Chem. Inf. Model. 2014, 54, 2856–2868. [Google Scholar] [CrossRef]
  51. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision Glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar]
  52. Gangloff, M.; Ruff, M.; Eiler, S.; Duclaud, S.; Wurtz, J.M.; Moras, D. Crystal structure of a mutant hERα ligand-binding domain reveals key structural features for the mechanism of partial agonism. J. Biol. Chem. 2001, 276, 15059–15065. [Google Scholar] [CrossRef] [Green Version]
  53. Halgren, T. New method for fast and accurate binding-site identification and analysis. Chem. Biol. Drug Des. 2007, 69, 146–148. [Google Scholar] [CrossRef] [PubMed]
  54. Spinello, A.; Martini, S.; Berti, F.; Pennati, M.; Pavlin, M.; Sgrignani, J.; Grazioso, G.; Colombo, G.; Zaffaroni, N.; Magistrato, A. Rational design of allosteric modulators of the aromatase enzyme: An unprecedented therapeutic strategy to fight breast cancer. Eur. J. Med. Chem. 2019, 168, 253–262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef] [PubMed]
  56. Anandakrishnan, R.; Aguilar, B.; Onufriev, A.V. H++ 3.0: Automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations. Nucl. Aci. Res. 2012, 40, W537–W541. [Google Scholar] [CrossRef] [Green Version]
  57. Spinello, A.; Pavlin, M.; Casalino, L.; Magistrato, A. A dehydrogenase dual hydrogen abstraction mechanism promotes estrogen biosynthesis: Can we expand the functional annotation of the aromatase enzyme? Chem. Eur. J. 2018, 24, 10840–10849. [Google Scholar] [CrossRef]
  58. Liu, B.; Dong, L.; Yu, Q.; Li, X.; Wu, F.; Tan, Z.; Luo, S. Thermodynamic study on the protonation reactions of glyphosate in aqueous solution: Potentiometry, calorimetry and NMR spectroscopy. J. Phys. Chem. B 2016, 120, 2132–2137. [Google Scholar] [CrossRef]
  59. Wickstrom, L.; Okur, A.; Simmerling, C. Evaluating the performance of the ff99SB force field based on NMR scalar coupling data. Biophys. J. 2009, 97, 853–856. [Google Scholar] [CrossRef] [Green Version]
  60. Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J.L.; Dror, R.O.; Shaw, D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 2010, 78, 1950–1958. [Google Scholar] [CrossRef] [Green Version]
  61. Dickson, C.J.; Madej, B.D.; Skjevik, Å.A.; Betz, R.M.; Teigen, K.; Gould, I.R.; Walker, R.C. Lipid14: The Amber Lipid Force Field. J. Chem. Theory Comput. 2014, 10, 865–879. [Google Scholar] [CrossRef]
  62. Shahrokh, K.; Orendt, A.; Yost, G.S.; Cheatham, T.E. Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle. J. Comput. Chem. 2012, 33, 119–133. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  64. Bayly, C.I.; Cieplak, P.; Cornell, W.; Kollman, P.A. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: The RESP model. J. Phys. Chem. 1993, 97, 10269–10280. [Google Scholar] [CrossRef]
  65. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision A.02; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  66. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic atom type and bond type perception in molecular mechanical calculations. J. Mol. Graphics Model. 2006, 25, 247–260. [Google Scholar] [CrossRef] [PubMed]
  67. Sousa da Silva, A.W.; Vranken, W.F. ACPYPE - Antechamber Python Parser Interface. BMC Res. Notes 2012, 5, 367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  69. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  70. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [Green Version]
  71. Balaguer, P.; Boussioux, A.-M.; Demirpence, E.; Nicolas, J.C. Reporter cell lines are useful tools for monitoring biological activity of nuclear receptor ligands. Luminescence 2001, 16, 153–158. [Google Scholar] [CrossRef]
  72. Berckmans, P.; Leppens, H.; Vangenechten, C.; Witters, H. Screening of endocrine disrupting chemicals with MELN cells, an ER-transactivation assay combined with cytotoxicity assessment. Toxicol. In Vitro 2007, 21, 1262–1267. [Google Scholar] [CrossRef]
  73. Balaguer, P.; François, F.; Comunale, F.; Fenet, H.; Boussioux, A.-M.; Pons, M.; Nicolas, J.C.; Casellas, C. Reporter cell lines to study the estrogenic effects of xenoestrogens. Sci. Total. Environ. 1999, 233, 47–56. [Google Scholar] [CrossRef]
  74. Schilirò, T.; Porfido, A.; Longo, A.; Coluccia, S.; Gilli, G. The E-screen test and the MELN gene-reporter assay used for determination of estrogenic activity in fruits and vegetables in relation to pesticide residues. Food Chem. Toxicol. 2013, 62, 82–90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Körner, W.; Hanf, V.; Schuller, W.; Kempter, C.; Metzger, J.; Hagenmaier, H. Development of a sensitive E-screen assay for quantitative analysis of estrogenic activity in municipal sewage plant effluents. Sci. Total Environ. 1999, 225, 33–48. [Google Scholar] [CrossRef]
  76. Gillezeau, C.; Van Gerwen, M.; Shaffer, R.M.; Rana, I.; Zhang, L.; Sheppard, L.; Taioli, E. The evidence of human exposure to glyphosate: A review. Environ. Health 2019, 18, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Magistrato, A.; Sgrignani, J.; Krause, R.; Cavalli, A. Single or multiple access channels to the CYP450s active Site? An answer from free energy simulations of the human aromatase enzyme. J. Phys. Chem. Lett. 2017, 8, 2036–2042. [Google Scholar] [CrossRef]
  78. Ritacco, I.; Saltalamacchia, A.; Spinello, A.; Ippoliti, E.; Magistrato, A. All-atom simulations disclose how cytochrome reductase reshapes the substrate access/egress routes of its partner CYP450s. J. Phys. Chem. Lett. 2020, 11, 1189–1193. [Google Scholar] [CrossRef]
  79. Scholfield, M.R.; Vander Zanden, C.M.; Carter, M.; Ho, P.S. Halogen bonding (X-bonding): A biological perspective: Halogen Bonding (X-Bonding). Protein Sci. 2013, 22, 139–152. [Google Scholar] [CrossRef]
  80. Parween, S.; Di Nardo, G.; Baj, F.; Zhang, C.; Gilardi, G.; Pandey, A.V. Differential effects of variations in human P450 oxidoreductase on the aromatase activity of CYP19A1 polymorphisms R264C and R264H. J. Steroid Biochem. Mol. Biol. 2020, 196, 105507. [Google Scholar] [CrossRef]
  81. Zhang, C.; Catucci, G.; Di Nardo, G.; Gilardi, G. Effector role of cytochrome P450 reductase for androstenedione binding to human aromatase. Int. J. Biol. Macromol. 2020, 164, 510–517. [Google Scholar] [CrossRef]
  82. Cassault-Meyer, E.; Gress, S.; Séralini, G.-É.; Galeraud-Denis, I. An acute exposure to glyphosate-based herbicide alters aromatase levels in testis and sperm nuclear quality. Environ. Toxicol. Pharmacol. 2014, 38, 131–140. [Google Scholar] [CrossRef]
  83. Bureik, M.; Brück, N.; Hübel, K.; Bernhardt, R. The human mineralocorticoid receptor only partially differentiates between different ligands after expression in fission yeast. FEMS Yeast Res. 2005, 5, 627–633. [Google Scholar] [CrossRef] [Green Version]
  84. Wang, L.; Liu, T.; Liu, F.; Zhang, J.; Wu, Y.; Sun, H. Occurrence and profile characteristics of the pesticide imidacloprid, preservative parabens, and their metabolites in human urine from rural and urban China. Environ. Sci. Technol. 2015, 49, 14633–14640. [Google Scholar] [CrossRef] [PubMed]
  85. Tapparo, A.; Giorio, C.; Marzaro, M.; Marton, D.; Soldà, L.; Girolami, V. Rapid analysis of neonicotinoid insecticides in guttation drops of corn seedlings obtained from coated seeds. J. Environ. Monit. 2011, 13, 1564. [Google Scholar] [CrossRef] [PubMed]
  86. Defarge, N.; Takács, E.; Lozano, V.; Mesnage, R.; Spiroux de Vendômois, J.; Séralini, G.E.; Székács, A. Co-formulants in glyphosate-based herbicides disrupt aromatase activity in human cells below toxic levels. Int. J. Environ. Res. Public Health 2016, 13, 264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijerph 17 05664 g001 550
Figure 1. Aromatase activity in the presence of different concentrations of glyphosate and 50 nM (black squares) or 400 nM (red circles) of the substrate androstenedione. Statistical significance *: p-value < 0.05 versus C+.
Figure 1. Aromatase activity in the presence of different concentrations of glyphosate and 50 nM (black squares) or 400 nM (red circles) of the substrate androstenedione. Statistical significance *: p-value < 0.05 versus C+.
Ijerph 17 05664 g001
Ijerph 17 05664 g002 550
Figure 2. (A) Michaelis-Menten plots and (B) Lineweaver-Burk linearization for aromatase activity in the absence (black circles) and presence of 1000 nM (red squares) and 5000 nM (blue triangles) of glyphosate. In Panel (A), the data were fitted to the Michaelis-Menten equation using the Sigma Plot software to obtain the kinetic parameters.
Figure 2. (A) Michaelis-Menten plots and (B) Lineweaver-Burk linearization for aromatase activity in the absence (black circles) and presence of 1000 nM (red squares) and 5000 nM (blue triangles) of glyphosate. In Panel (A), the data were fitted to the Michaelis-Menten equation using the Sigma Plot software to obtain the kinetic parameters.
Ijerph 17 05664 g002
Ijerph 17 05664 g003 550
Figure 3. Representative structure of aromatase embedded in a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane, with phosphorous and oxygen atoms shown as tan and red van der Waals (vdw) spheres, and cholesterol (yellow surface) membrane. Sites 1 and 2 are shown as orange and purple transparent surfaces, respectively. The heme and androstenedione (ASD) are displayed in a vdw representation. The protein is shown as green new cartoons. The inset reports a close view of structure of aromatase in complex with glyphosate, as obtained from the most representative cluster of the molecular dynamics simulation trajectory. The heme moiety, ASD, and glyphosate are depicted in balls and sticks. The key catalytic residue Asp309, lining the binding cavity, is shown in licorice and colored by the atom name.
Figure 3. Representative structure of aromatase embedded in a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane, with phosphorous and oxygen atoms shown as tan and red van der Waals (vdw) spheres, and cholesterol (yellow surface) membrane. Sites 1 and 2 are shown as orange and purple transparent surfaces, respectively. The heme and androstenedione (ASD) are displayed in a vdw representation. The protein is shown as green new cartoons. The inset reports a close view of structure of aromatase in complex with glyphosate, as obtained from the most representative cluster of the molecular dynamics simulation trajectory. The heme moiety, ASD, and glyphosate are depicted in balls and sticks. The key catalytic residue Asp309, lining the binding cavity, is shown in licorice and colored by the atom name.
Ijerph 17 05664 g003
Ijerph 17 05664 g004 550
Figure 4. Estrogenic activity of pesticides measured with the MELN gene reporter assay. Data are expressed as means and standard deviations of the relative luciferase activity (% of C+, E2 10−8 M). The relative luciferase activity of the C+ is 100.0 ± 5.9%, while the relative luciferase activity of the negative control (test medium) is 30.5 ± 1.1%.
Figure 4. Estrogenic activity of pesticides measured with the MELN gene reporter assay. Data are expressed as means and standard deviations of the relative luciferase activity (% of C+, E2 10−8 M). The relative luciferase activity of the C+ is 100.0 ± 5.9%, while the relative luciferase activity of the negative control (test medium) is 30.5 ± 1.1%.
Ijerph 17 05664 g004
Ijerph 17 05664 g005 550
Figure 5. Estrogenic activity, measured with the MELN gene reporter assay, of imidacloprid (A,B) and thiacloprid (C,D) in combination with E2 (10−10 M) (B,D) or in combination with tamoxifen (10−6 M) (A,C). Data are expressed as the relative luciferase activity (% of E2 10−10 M). C−: Negative control; E2: E2 10−10 M; T: Tamoxifen 10−6 M. * p < 0.05 vs. C−; # p < 0.05 vs. E2; Kruskal-Wallis test followed by the post-hoc Dunnett test.
Figure 5. Estrogenic activity, measured with the MELN gene reporter assay, of imidacloprid (A,B) and thiacloprid (C,D) in combination with E2 (10−10 M) (B,D) or in combination with tamoxifen (10−6 M) (A,C). Data are expressed as the relative luciferase activity (% of E2 10−10 M). C−: Negative control; E2: E2 10−10 M; T: Tamoxifen 10−6 M. * p < 0.05 vs. C−; # p < 0.05 vs. E2; Kruskal-Wallis test followed by the post-hoc Dunnett test.
Ijerph 17 05664 g005
Ijerph 17 05664 g006 550
Figure 6. Model of estrogen receptor α dimer (PDB ID 1qku [52]) in complex with the neonicotinoids, imidacloprid and thiacloprid. The allosteric pocket is shown as a green transparent surface. Imidacloprid, thiacloprid, and 17-β-estradiol are displayed in a van der Waals representation and colored by the atom name. The protein is shown as violet new cartoons. The insets report a close view of docking poses of imidacloprid and thiacloprid inside the estrogen binding site (red circles) and onto the newly identified allosteric pocket (dark blue circles). The ligand and the residues establishing the most important interactions are depicted in balls and sticks and licorice representations, respectively, and colored by the atom name.
Figure 6. Model of estrogen receptor α dimer (PDB ID 1qku [52]) in complex with the neonicotinoids, imidacloprid and thiacloprid. The allosteric pocket is shown as a green transparent surface. Imidacloprid, thiacloprid, and 17-β-estradiol are displayed in a van der Waals representation and colored by the atom name. The protein is shown as violet new cartoons. The insets report a close view of docking poses of imidacloprid and thiacloprid inside the estrogen binding site (red circles) and onto the newly identified allosteric pocket (dark blue circles). The ligand and the residues establishing the most important interactions are depicted in balls and sticks and licorice representations, respectively, and colored by the atom name.
Ijerph 17 05664 g006
Table
Table 1. Effect of pesticides on the aromatase activity.
Table 1. Effect of pesticides on the aromatase activity.
PesticideRelative Activity (%)
500 nM1000 nM5000 nM
1Glyphosate76.6 ± 11.3 *74.5 ± 7.6 *36.0 ± 19.5 *
2Imidacloprid100.1 ± 8.892.9 ± 28.0120.5 ± 17.6
3Thiacloprid100.9 ± 11.9153.6 ± 56.9146.9 ± 42.1
Statistical significance *: p-value < 0.05 versus positive control (C+).
Table
Table 2. Kinetic parameters obtained from the fitting of the data in Figure 2A to a Michaelis-Menten curve. The kinetic parameters are calculated for aromatase activity in the absence and presence of 1000 and 5000 nM of glyphosate.
Table 2. Kinetic parameters obtained from the fitting of the data in Figure 2A to a Michaelis-Menten curve. The kinetic parameters are calculated for aromatase activity in the absence and presence of 1000 and 5000 nM of glyphosate.
Glyphosate (nM)KM (nM)Vmax (min−1)
041.3 ± 8.20.018 ± 0.001
100035.3 ± 3.50.013 ± 0.001 *
500092.3 ± 20.7 *0.011 ± 0.002 *
Statistical significance*: p-value < 0.05 versus the values obtained in the absence of glyphosate.

Share and Cite

MDPI and ACS Style

Zhang, C.; Schilirò, T.; Gea, M.; Bianchi, S.; Spinello, A.; Magistrato, A.; Gilardi, G.; Di Nardo, G. Molecular Basis for Endocrine Disruption by Pesticides Targeting Aromatase and Estrogen Receptor. Int. J. Environ. Res. Public Health 2020, 17, 5664. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17165664

AMA Style

Zhang C, Schilirò T, Gea M, Bianchi S, Spinello A, Magistrato A, Gilardi G, Di Nardo G. Molecular Basis for Endocrine Disruption by Pesticides Targeting Aromatase and Estrogen Receptor. International Journal of Environmental Research and Public Health. 2020; 17(16):5664. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17165664

Chicago/Turabian Style

Zhang, Chao, Tiziana Schilirò, Marta Gea, Silvia Bianchi, Angelo Spinello, Alessandra Magistrato, Gianfranco Gilardi, and Giovanna Di Nardo. 2020. "Molecular Basis for Endocrine Disruption by Pesticides Targeting Aromatase and Estrogen Receptor" International Journal of Environmental Research and Public Health 17, no. 16: 5664. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17165664

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop