Applications of GnRH in the control and management of fertility in female animals

Gonadotropin releasing hormone (GnRH) has long been recognized as a potential target for the control and management of fertility in female animals. Attempts to apply GnRH-based technology to manage fertility have focussed on the development of GnRH agonists, antagonists and vaccines. All of these methods have potential, but the widespread application of these technologies has been limited to date. The greatest advance in the use of GnRH-based technology for long-term fertility control in recent years has been the development and commercialization of depot formulations that release GnRH agonists for periods of up to 1 year. These products have a broad range of potential applications in production and domestic animal management. The further development and commercialization of GnRH vaccines has been hampered by the variability of response between individual animals. The need to use adjuvant and multiple boosters also make this a less attractive option than the current GnRH agonist technology. However, GnRH vaccines have the advantage that they do not induce the initial stimulatory response that follows GnRH agonist administration. GnRH antagonists and GnRH-toxin conjugates show promise but are in an earlier phase of development. To date, no depot or long-acting formulations of antagonists have been developed. GnRH-toxin conjugates have yet to achieve permanent sterilization, but further dose-response trials may advance this approach.     

Regulation of primate testicular function by GnRH analogues

The potential of GnRH analogues for regulating testicular function is reviewed. Our experiments showed that constant infusion of GnRH agonists effectively suppressed testicular function in monkeys. In men, however, spermatogenesis could not be suppressed to achieve azoospermia uniformly. GnRH antagonists, although at much higher dosages than agonists, caused a more rapid and uniform inhibition of testis function. Spermatogenesis was reversibly disrupted at the spermatogonial level. Concomitant testosterone supplementation, used to maintain libido and potency, attenuated the antitesticular effects of GnRH analogues. In monkeys testosterone appears to stimulate spermatogenesis directly on the testicular level, while evidence has been obtained that in rats testosterone can also stimulate the release and synthesis of FSH under antagonist mediated blockage of pituitary GnRH receptors. When extrapolating to human studies special care has to be exerted in the selection of testosterone substitution regimens. Although the agonistic and antagonistic analogues of GnRH ultimately exert their antireproductive effects via inhibition of gonadotropin secretion the antagonists may have the greater potential for male fertility regulation due to quicker pituitary and testicular suppression.     

Gonadotropin-releasing hormone (GnRH) and its natural analogues: a review

The pivotal role of gonadotropin-releasing hormone (GnRH) during the hormonal regulation of reproductive processes is indisputable. Likewise, many factors are known to affect reproductive function by influencing either GnRH release from hypothalamus or pituitary gland responsiveness to GnRH. In veterinary medicine, GnRH and its agonists (GnRHa) are widely used to overcome reduced fertility by ovarian dysfunction, to induce ovulation, and to improve conception rate. GnRHa are, moreover, integrative part of other pro-fertility treatments, e.g. for synchronization of the estrous cycle or stimulation for embryo transfer. Additionally, continuous GnRH which shows desensitizing effects of the pituitary-ovarian axis has been recommended for implementation in anti-fertility treatments like inhibition of ovulation or reversible blockade of the estrous cycle. Just as much, another group of GnRH analogues, antagonists, are now in principle disposable for use. For a few decades, GnRH was thought to be a unique structure with a primary role in regulation gonadotropins. However, it became apparent that other homologous ligands of the GnRH receptor (GnRHR) exist. In the meantime, more than 20 natural variants of the mammalian GnRH have been identified in different species which may compete for binding and/or have their own receptors. These GnRH forms (GnRHs) have apparently common and divergent functions. More studies on GnRHs should contribute to a better understanding of reproductive processes in mammals and interactions between reproduction and other physiological functions. Increased information on GnRHs might raise expectations in the application of these peptides in veterinary practice. It is the aim of this review to discuss latest results from evolutionarily based studies as well as first experimental tests and to answer the question how realistic might be the efforts to develop effective and animal friendly practical applications for endogenous GnRHs and synthetic analogues.     

New GnRH analogs in canine reproduction

Gonadotropin-releasing hormone (GnRH) analogs, which include agonists and antagonists, have been produced by amino acid substitutions of the native GnRH molecule to create greater potency and longer duration of effectiveness. The aim of this article was to review the pharmacological effects and the existing clinical literature of new GnRH analogs, namely agonists released from long-term formulations and third generation antagonists, in domestic dogs. Long-term administration of agonists functions through desensitization and down-regulation of GnRH pituitary receptors inhibiting gonadotropin production and release after an initial stimulation. Conversely, GnRH antagonists bind to gonadotrope GnRH receptors and compete successfully with endogenous GnRH for occupancy, thereby inhibiting the pituitary-gonadal axis immediately. There is a promising place for both agonists and antagonists in future canine reproduction. They can be used in the control of estrous cycle, hormone-dependent diseases as well as in contraception. Some information on the effectiveness and safety of these new analogs in canine reproduction is already available, yet further work is needed before they could be widely recommended. The increase in gonadotropins and gonadal steroids following administration of agonists might have adverse effects when used on hormone-dependent diseases. This initial "flare" should also be carefully managed in anestrous and prepubertal bitches. At present, the main application of antagonists seems limited to situations where an acute endocrine, inhibitory effect is required, e.g. proestrus or pregnancy termination. Future commercial availability of long-acting, single-dose antagonists could go far towards controlling pet population.     

Role, mechanism of action and application of gonadoliberins in reproductive processes

Gonadoliberin (gonadotropin releasing hormone, GnRH) plays a central role in the regulation of reproductive functions as it regulates the release of both luteinizing hormone (LH) and follicle stimulating hormone (FSH). The isolation and structure determination of GnRH opened the possibility of its use for influencing reproductive processes. This possibility initiated a rapid development in the design of potent and long-acting GnRH agonists and antagonists. The most important structural modifications of GnRH leading to superagonists are the D-amino acid substitutions in position 6 combined with Pro9-ethylamide or azaGly10 at the C-terminus. We have synthesized several superagonists of GnRH according to these substitution principles. Furthermore, our L-isoaspartyl modification in position 6, as a new approach to GnRH agonist design, also resulted in superactive analogs. The recently discovered sequences of non-mammalian GnRH-s opened new routes for us to synthesize species specific GnRH agonists. All three groups of the above mentioned GnRH analogs have been successfully used for the treatment of sexual disorders of different animals (cattle, pigs, rabbits, etc.). Ovulation synchronization and a 30% increase in the fertility rate could be achieved by using GnRH agonists in cattle breeding. Analogs derived from species specific sequences could be applied for the induced artificial propagation of fish even out of the spawning season. It is known that superactive GnRH analogs can suppress the growth of certain hormone-dependent tumours. In vitro and in vivo tests of our analogs showed promising antitumour activity in breast cancer which might be explained by the mechanism of desensitization. Almost a hundred antagonist analogs of GnRH have been developed in our laboratory. The most effective ones contain 4 or 5 D-amino acids, and one of them is even orally active. The inhibition of ovulation can also be achieved by the administration of GnRH superagonists. This phenomenon might also be explained by the desensitization of LH-release. Radioactive analogs specifically labeled with tritium in different amino acid residues have been synthesized and used for studying tissue distribution and biodegradation of gonadoliberins. Analogs containing a photoreactive group have been prepared and applied for the trials of GnRH receptor isolation.     

GnRHs and GnRH receptors

GnRH is the pivotal hypothalamic hormone regulating reproduction. Over 20 forms of the decapeptide have been identified in which the NH2- and COOH-terminal sequences, which are essential for receptor binding and activation, are conserved. In mammals, there are two forms, GnRH I which regulates gonadotropin and GnRH II which appears to be a neuromodulator and stimulates sexual behaviour. GnRHs also occur in reproductive tissues and tumours in which a paracrine/autocrine role is postulated. GnRH agonists and antagonists are now extensively used to treat hormone-dependent diseases, in assisted conception and have promise as novel contraceptives. Non-peptide orally-active GnRH antagonists have been recently developed and may increase the flexibility and range of utility. As with GnRH, GnRH receptors have undergone co-ordinated gene duplications such that cognate receptor subtypes for respective ligands exist in most vertebrates. Interestingly, in man and some other mammals (e.g. chimp, sheep and bovine) the Type II GnRH receptor has been silenced. However, GnRH I and GnRH II still appear to have distinct roles in signalling differentially through the Type I receptor (ligand-selective-signalling) to have different downstream effects. The ligand-receptor interactions and receptor conformational changes involved in receptor activation have been partly delineated. Together, these findings are setting the scene for generating novel selective GnRH analogues with potential for wider and more specific application.     

Two mutations in extracellular loop 2 of the human GnRH receptor convert an antagonist to an agonist

GnRH regulates the reproductive system through cognate G protein-coupled receptors in vertebrates. Certain GnRH analogs that are antagonists at mammalian receptors behave as agonists at Xenopus laevis and chicken receptors. This phenomenon provides the opportunity to elucidate interactions and the mechanism underlying receptor activation. A D-Lys(iPr) in position 6 of the mammalian GnRH receptor antagonist is required for this agonist activity (inositol phosphate production) in the chicken and X. laevis GnRH receptors. Chimeric receptors, in which extracellular loop domains of the human GnRH receptor were substituted with the equivalent domains of the X. laevis GnRH receptor, identified extracellular loop 2 as the determinant for agonist activity of one of the mammalian antagonists: antagonist 135-18. Site-directed mutagenesis of nine nonconserved residues in the C-terminal domain of extracellular loop 2 of the human GnRH receptor showed that a minimum of two mutations (Val(5.24(197))Ala and Trp(5.32(205))His) is needed in this region for agonist activity of antagonist 135-18. Agonist activity of antagonist 135-18 was markedly decreased by low pH (<7.0) compared with GnRH agonists. These findings indicate that D-Lys(iPr)(6) forms a charge-supported hydrogen bond with His(5.32(205)) to stabilize the receptor in the active conformation. This discovery highlights the importance of EL-2 in ligand binding and receptor activation in G protein-coupled receptors.     

Conserved amino acid residues that are important for ligand binding in the type I gonadotropin-releasing hormone (GnRH) receptor are required for high potency of GnRH II at the type II GnRH receptor

GnRH I regulates reproduction. A second form, designated GnRH II, selectively binds type II GnRH receptors. Amino acids of the type I GnRH receptor required for binding of GnRH I (Asp2.61(98), Asn2.65(102), and Lys3.32(121)) are conserved in the type II GnRH receptor, but their roles in receptor function are unknown. We have delineated their functions using mutagenesis, signaling and binding assays, immunoblotting, and computational modeling. Mutating Asp2.61(97) to Glu or Ala, Asn2.65(101) to Ala, or Lys3.32(120) to Gln decreased potency of GnRH II-stimulated inositol phosphate production. Consistent with proposed roles in ligand recognition, mutations eliminated measurable binding of GnRH II, whereas expression of mutant receptors was not decreased. In detailed analysis of how these residues affect ligand-dependent signaling, [Trp2]-GnRH I showed lesser decreases in potency than GnRH I at the Asp2.61(97)Glu mutant. In contrast, [Trp2]-GnRH II showed the same loss of potency as GnRH II at this mutant. This suggests that Asp2.61(97) contributes to recognition of His2 of GnRH I, but not of GnRH II. GnRH II showed a large decrease in potency at the Asn2.65(101)Ala mutant compared with analogs lacking the CO group of Gly10NH2. This suggests that Asn2.65(101) recognizes Gly10NH2 of GnRH II. GnRH agonists showed large decreases in potency at the Lys3.32(120)Gln mutant, but antagonist activity was unaffected. This suggests that Lys3.32(120) recognizes agonists, but not antagonists, as in the type I receptor. These data indicate that roles of conserved residues are similar, but not identical, in the type I and II GnRH receptors.     

A chicken gonadotropin-releasing hormone receptor that confers agonist activity to mammalian antagonists. Identification of D-Lys(6) in the ligand and extracellular loop two of the receptor as determinants

Mammalian receptors for gonadotropin-releasing hormone (GnRH) have over 85% sequence homology and similar ligand selectivity. Biological studies indicated that the chicken GnRH receptor has a distinct pharmacology, and certain antagonists of mammalian GnRH receptors function as agonists. To explore the structural determinants of this, we have cloned a chicken pituitary GnRH receptor and demonstrated that it has marked differences in primary amino acid sequence (59% homology) and in its interactions with GnRH analogs. The chicken GnRH receptor had high affinity for mammalian GnRH (K(i) 4.1 +/- 1.2 nM), similar to the human receptor (K(i) 4.8 +/- 1.2 nM). But, in contrast to the human receptor, it also had high affinity for chicken GnRH ([Gln(8)]GnRH) and GnRH II ([His(5),Trp(7),Tyr(8)]GnRH) (K(i) 5.3 +/- 0.5 and 0.6 +/- 0.01 nM). Three mammalian receptor antagonists were also pure antagonists in the chicken GnRH receptor. Another three, characterized by D-Lys(6) or D-isopropyl-Lys(6) moieties, functioned as pure antagonists in the human receptor but were full or partial agonists in the chicken receptor. This suggests that the Lys side chain interacts with functional groups of the chicken GnRH receptor to stabilize it in the active conformation and that these groups are not available in the activated human GnRH receptor. Substitution of the human receptor extracellular loop two with the chicken extracellular loop two identified this domain as capable of conferring agonist activity to mammalian antagonists. Although functioning of antagonists as agonists has been shown to be species-dependent for several GPCRs, the dependence of this on an extracellular domain has not been described.     

Using gonadotropin-releasing hormone (GnRH) and GnRH analogs to modulate testis function and enhance the productivity of domestic animals

Gonadotropin releasing hormone (GnRH) controls the activity of the gonadotrope cells of the pituitary gland and, as a consequence, is a critical component of the endocrine cascade that determines the growth, development, and functional activity of testicular tissue. The use of GnRH and GnRH analogs is common in domestic animal production systems. Although GnRH and GnRH analogs are most commonly used to control the fertility and reproductive events in female animals, GnRH agonists and antagonists are increasingly used to modulate the fertility, behavior, and productivity of male animals as well. This review will focus on recent advances in this use of GnRH agonists and antagonists.     

Subtleties of structure-agonist versus antagonist relationships of opioid peptides and peptidomimetics

The development of novel delta opioid antagonists and delta opioid agonists structurally derived from the prototype delta antagonist TIPP (H-Tyr-Tic-Phe-Phe-OH), is reviewed. Both delta antagonists and delta agonists with extraordinary potency and unprecedented delta receptor selectivity were discovered. Some of them are already widely used as pharmacological tools and are also of interest as potential therapeutic agents for use in analgesia. The results of the performed structure-activity studies revealed that the delta antagonist versus delta agonist behavior of this class of compounds depended on very subtle structural differences in diverse locations of the molecule. These observations can be best explained with a receptor model involving a number of different inactive and active receptor conformations.     

2-phenyl-4-piperazinylbenzimidazoles: orally active inhibitors of the gonadotropin releasing hormone (GnRH) receptor

Antagonism of the gonadotropin releasing hormone (GnRH) receptor has shown positive clinical results in numerous reproductive tissue disorders such as endometriosis, prostate cancer and others. Traditional therapy has been limited to peptide agonists and antagonists. Recently, small molecule GnRH antagonists have emerged as potentially new treatments. This article describes the discovery of 2-phenyl-4-piperazinylbenzimidazoles as small molecule GnRH antagonists with nanomolar potency in in vitro binding and functional assays, excellent bioavailability (rat %F>70) and demonstrated oral activity in a rat model having shown significant serum leuteinizing hormone (LH) suppression.     

Conformationally constrained opioid peptide analogs with novel activity profiles

Summary Novel conformationally constrained opioid peptide analogs with δ antagonist, mixed μ agonist/δ antagonist or δ agonist properties were developed. TIP(P)-related δ antagonists showed unprecedented δ antagonist potency and δ receptor selectivity, and may have potential for use in analgesia in combination with μ agonists. A definitive model of their δ receptor-bound conformation was developed. Three prototype mixed μ agonist/δ antagonists were discovered. They represent the only known compounds with this pharmacological profile and, as expected, one of them was shown to be a potent analgesic and to produce no dependence and less tolerance than morphine. Novel dipeptide derivatives turned out to be potent and selective δ agonists. Because of their low molecular weight and lipophilic character, these compounds may cross the blood-brain barrier and, thus, may have potential as centrally acting analgesics.     

Effects of gonadotropin-releasing hormone agonists on meiotic maturation of follicle-enclosed oocytes in rabbits

The effects of GnRH agonists on in vitro maturation of rabbit follicle-enclosed oocytes were studied. Rabbit preovulatory follicles were cultured with or without hCG (10(2) ng/ml), buserelin (10(2)-10(5) ng/ml), or leuprolide (10(2)-10(5) ng/ml) for 14 hours in vitro. GnRH agonists induced the resumption of meiosis in the follicle-enclosed oocytes in a dose-dependent manner. The percentage of oocytes achieving GVBD following treatment with 10(5) ng/ml buserelin (87.9 +/- 6.3%) or 10(5) ng/ml leuprolide (86.0 +/- 4.1%) did not differ significantly from hCG-treated control (87.3 +/- 3.8%). Mature oocytes initially were detected within 2 hours of GnRH agonist exposure. Concomitant addition of a GnRH antagonist at 10(4) ng/ml significantly blocked the stimulatory effect of GnRH agonist on oocyte maturation. GnRH agonists significantly stimulated both prostaglandin (PG) E2 (PGE2) and PGF2 alpha production by preovulatory follicles (p less than 0.01), but secreted prostanoid levels did not differ significantly among different concentrations of GnRH agonists. Meiotic maturation of follicle-enclosed oocytes following GnRH agonist exposure began 2 hours earlier than production of PGs. PG production stimulated by GnRH agonists was reduced significantly by indomethacin. However, oocyte maturity in the presence of GnRH agonist plus indomethacin did not differ significantly from that of GnRH agonist alone. GnRH agonistic analogues induce the resumption of meiosis in follicle-enclosed oocytes in rabbits by a mechanism other than PG stimulation.     

Sodium-induced alterations of opiate effects on the binding of 3H-dihydromorphine to mouse brain homogenates.

The inhibitory potency of opiate agonists on the stereo-specific binding of 3H-dihydromorphine in mouse brain homogenates was not affected by the presence of sodium ions. That of pure antagonists was greatly enhanced by NaCl whereas the inhibitory effects of mixed agonist-antagonists were reduced by NaCl, indicating that sodium ions might affect the agonist component more than the antagonist component of narcotic-antagonist analgesics. The inhibitory potency of the opiates tested in our system agrees with their potency in reducing the stereospecific binding of 3H-naloxone to rat membranes, the contractions of co-axially stimulated guinea pig ileum and their analgesic potency in animals and humans.     

Covalently cyclized agonist and antagonist analogues of bombesin and related peptides.

During a search for possible cyclization points in shortened, potent bombesin agonists and antagonists, it was found that the joining of amino acid residues in positions 6 and 14 by various means resulted in retention of significant binding affinity for rat pancreatic acini and murine Swiss 3T3 cells. In one series of analogues, Cys residues in these positions were used for bridging via a disulfide bond. (D)-C-Q-W-A-V-G-H-L-C-NH2 retained significant binding affinity for rat pancreatic acini cells and was a full amylase releasing agonist (EC50 187 nM). Potency was markedly increased by substituting D-Ala for Gly (EC50 67 nM compared to 10 nM for its linear counterpart) and was decreased by substituting L-Cys for D-Cys in this analogue (EC50 214 nM), thus strongly suggesting stabilization of peptide folding by the D residues. Elimination of the COOH-terminal amino acid produces competitive antagonists in the linear analogues; however, (D)-C-Q-W-A-V-G-H-C-NH2 was devoid of activity. Likewise, cyclization to position 13 with the 14 amino acids intact to give (D)-C-Q-W-A-V-G-H-C-L-NH2 resulted in an almost inactive peptide. On the other hand, as in the linear series, the reduced peptide bond analogue, (D)-C-Q-W-A-V-(D)-A-H-L-psi (CH2NH)-C-NH2, was a receptor antagonist (IC50 5.7 mM), albeit much weaker than the corresponding linear analogues, but with no residual agonist activity. Direct head-to-tail cyclization was also tried. Both cyclo[(D)-F-Q-W-A-V-G-H-L-L] (EC50 346 nM) and the shorter cyclo [Q-W-A-V-G-H-L-L] (EC50 1236 nM) were full agonists. Elimination of the COOH-terminal residue in cyclo[(D)-p-Cl-F-Q-W-A-V-(D)-A-H-L] produced an agonist (EC50 716 nM) rather than an antagonist. These results provide support for the proposal that both bombesin agonists and antagonists adopt a folded conformation at their receptor(s). Furthermore, the retention of appreciable potencies using several cyclization strategies and chain lengths suggests that further optimization of these structures both in terms of potency and ring size is possible. Since these peptides have increased conformational restriction, they should begin to serve as useful substrates for NMR and molecular modeling studies aimed at comparing the obviously subtle differences between agonist and antagonist structures.     

New protocol for commencing the GnRH antagonist in assisted conception treatment cycles: elimination of the premature LH surge with similar pregnancy rates

GnRH antagonists have been used with increasing frequency in assisted reproduction treatments over the past few years and have been associated with quicker and more profound LH suppression and shorter treatment cycles than conventional GnRH agonists. Usually, these are commenced on day 6 of FSH stimulation without allowing for patient variation in response to treatment. The study was aimed at individualising this protocol to the patients' ovarian response. The control group included 215 treatment cycles where the GnRH antagonist was commenced on day 6 of FSH stimulation. A new individualised protocol was formulated, applied to practice and 172 treatment cycles following that were analysed. The study group had no premature LH surges (LH > 10 iu x mL(-1)) compared to the control group who had a rate of 4.1%. There was also a higher fertilisation and clinical pregnancy rate in the study group (P = 0.06). It is concluded that the new individualised GnRH antagonist protocol eliminates premature LH surges in assisted conception treatment cycles and may improve clinical pregnancy rates compared to the conventional protocol of "day 6 commencement".     

Data-driven analysis in drug discovery

In the process of drug discovery for new chemical entities, application of appropriate pharmacological models often is not possible because the molecular mechanism of the compound is not yet elucidated. Therefore, a data-driven approach using generic tools designed to quantify characteristic patterns of concentration-response curves is required. This article outlines the options available for quantifying agonist and antagonist activity. Specifically, for agonists, the use of the Operational model for the determination of functional effects (equimolar potency ratios for full agonists, calculation of relative efficacy) is described. For antagonists, the measurement of pKB (-log of the equilibrium dissociation constant of the antagonist-receptor complex) for orthosteric antagonists that do not alter basal response (simple competitive antagonists), increase basal response (partial agonists), and decrease basal response (in constitutively active systems; inverse agonists) is discussed. In addition, this article considers methods to discern orthosteric receptor antagonism from allosteric antagonism whereby the agonist and antagonist bind to separate sites and interact through a conformational change in the receptor. Methods for the measurement of the pKB for allosteric modulators as well as co-operativity constants for these modulators is described.     

Structure-function studies of linear and cyclized peptide antagonists of the GnRH receptor.

Structurally new analogs of the peptidic GnRH receptor antagonist Cetrorelix as well as conformationally constrained cyclized deca- or pentapeptides were synthesized and selected peptides evaluated comprehensively. To understand how structural variations of the antagonistic peptide effect pharmacodynamic properties, binding affinities and antagonistic potencies toward the human and rat GnRH receptor were determined. Whereas large substituents in position 6 of linear peptides are compatible with high binding affinity (K(D) < 0.5 nM), all cyclized peptides except the cyclo[3-10] analog D-52391 depicted low binding affinity (K(D) > 10 nM). Binding affinity and antagonistic potency in vitro correlated for all peptides and surprisingly no discrimination between human and rat receptor proteins was observed. Since receptor residues W(101) and N(102) are involved in agonist and antagonist binding, equally potent but structurally different antagonists were tested for binding to the respective W(101)A and N(102)A mutants. In contrast to linear decapeptides, residues N(102) and W(101) are not involved in binding of D-23938 and W(101) is the critical residue for D-52391 binding. We conclude that although equally potent, peptidic GnRH receptor antagonists do have distinct interactions within the ligand binding pocket. Finally, selected antagonists were tested for testosterone suppression in male rats. The duration of testosterone suppression below castration levels differed largely from 1 day for Ganirelix to 27 days for D-23487. Systemic availability became evident as the most important parameter for in vivo efficacy.     

Identification of Thyroid Receptor Ant/Agonists in Water Sources Using Mass Balance Analysis and Monte Carlo Simulation

Some synthetic chemicals, which have been shown to disrupt thyroid hormone (TH) function, have been detected in surface waters and people have the potential to be exposed through water-drinking. Here, the presence of thyroid-active chemicals and their toxic potential in drinking water sources in Yangtze River Delta were investigated by use of instrumental analysis combined with cell-based reporter gene assay. A novel approach was developed to use Monte Carlo simulation, for evaluation of the potential risks of measured concentrations of TH agonists and antagonists and to determine the major contributors to observed thyroid receptor (TR) antagonist potency. None of the extracts exhibited TR agonist potency, while 12 of 14 water samples exhibited TR antagonistic potency. The most probable observed antagonist equivalents ranged from 1.4 to 5.6 µg di-n-butyl phthalate (DNBP)/L, which posed potential risk in water sources. Based on Monte Carlo simulation related mass balance analysis, DNBP accounted for 64.4% for the entire observed antagonist toxic unit in water sources, while diisobutyl phthalate (DIBP), di-n-octyl phthalate (DNOP) and di-2-ethylhexyl phthalate (DEHP) also contributed. The most probable observed equivalent and most probable relative potency (REP) derived from Monte Carlo simulation is useful for potency comparison and responsible chemicals screening.