QSAR models using a large diverse set of estrogens

Endocrine disruptors (EDs) have a variety of adverse effects in humans and animals. About 58,000 chemicals, most having little safety data, must be tested in a group of tiered assays. As assays will take years, it is important to develop rapid methods to help in priority setting. For application to large data sets, we have developed an integrated system that contains sequential four phases to predict the ability of chemicals to bind to the estrogen receptor (ER), a prevalent mechanism for estrogenic EDs. Here we report the results of evaluating two types of QSAR models for inclusion in phase III to quantitatively predict chemical binding to the ER. Our data set for the relative binding affinities (RBAs) to the ER consists of 130 chemicals covering a wide range of structural diversity and a 6 orders of magnitude spread of RBAs. CoMFA and HQSAR models were constructed and compared for performance. The CoMFA model had a r2 = 0.91 and a q2LOO = 0.66. HQSAR showed reduced performance compared to CoMFA with r2 = 0.76 and q2LOO = 0.59. A number of parameters were examined to improve the CoMFA model. Of these, a phenol indicator increased the q2LOO to 0.71. When up to 50% of the chemicals were left out in the leave-N-out cross-validation, the q2 remained significant. Finally, the models were tested by using two test sets; the q2pred for these were 0.71 and 0.62, a significant result which demonstrates the utility of the CoMFA model for predicting the RBAs of chemicals not included in the training set. If used in conjunction with phases I and II, which reduced the size of the data set dramatically by eliminating most inactive chemicals, the current CoMFA model (phase III) can be used to predict the RBA of chemicals with sufficient accuracy and to provide quantitative information for priority setting.     

Comparative molecular field analysis (CoMFA) model using a large diverse set of natural, synthetic and environmental chemicals for binding to the androgen receptor

A large number of natural, synthetic and environmental chemicals are capable of disrupting the endocrine systems of experimental animals, wildlife and humans. These so-called endocrine disrupting chemicals (EDCs), some mimic the functions of the endogenous androgens, have become a concern to the public health. Androgens play an important role in many physiological processes, including the development and maintenance of male sexual characteristics. A common mechanism for androgen to produce both normal and adverse effects is binding to the androgen receptor (AR). In this study, we used Comparative Molecular Field Analysis (CoMFA), a three-dimensional quantitative structure-activity relationship (3D-QSAR) technique, to examine AR-ligand binding affinities. A CoMFA model with r2 = 0.902 and q2 = 0.571 was developed using a large training data set containing 146 structurally diverse natural, synthetic, and environmental chemicals with a 10(6)-fold range of relative binding affinity (RBA). By comparing the binding characteristics derived from the CoMFA contour map with these observed in a human AR crystal structure, we found that the steric and electrostatic properties encoded in this training data set are necessary and sufficient to describe the RBA of AR ligands. Finally, the CoMFA model was challenged with an external test data set; the predicted results were close to the actual values with average difference of 0.637 logRBA. This study demonstrates the utility of this CoMFA model for real-world use in predicting the AR binding affinities of structurally diverse chemicals over a wide RBA range.     

A reactivity pattern for discrimination of ER agonism and antagonism based on 3-D molecular attributes

Various models have been developed to predict the relative binding affinity (RBA) of chemicals to estrogen receptors (ER). These models can be used to prioritize chemicals for further tiered biological testing to assess the potential for endocrine disruption. One shortcoming of models predicting RBA has been the inability to distinguish potential receptor antagonism from agonism, and hence in vivo response. It has been suggested that steroid receptor antagonists are less compact than agonists; thus, ER binding of antagonists may prohibit proper alignment of receptor protein helices preventing subsequent transactivation. The current study tests the theory of chemical bulk as a defining parameter of antagonism by employing a 3-D structural approach for development of reactivity patterns for ER antagonists and agonists. Using a dataset of 23 potent ER ligands (16 agonists, 7 antagonists), molecular parameters previously found to be associated with ER binding affinity, namely global ( E HOMO ) and local (donor delocalizabilities and charges) electron donating ability of electronegative sites and steric distances between those sites, were found insufficient to discriminate ER antagonists from agonists. However, parameters related to molecular bulk, including solvent accessible surface and negatively charged Van der Waal's surface, provided reactivity patterns that were 100% successful in discriminating antagonists from agonists in the limited data set tested. The model also shows potential to discriminate pure antagonists from partial agonist/antagonist structures. Using this exploratory model it is possible to predict additional chemicals for their ability to bind but inactivate the ER, providing a further tool for hypothesis testing to elucidate chemical structural characteristics associated with estrogenicity and anti-estrogenicity.     

A reactivity pattern for discrimination of ER agonism and antagonism based on 3-D molecular attributes

Various models have been developed to predict the relative binding affinity (RBA) of chemicals to estrogen receptors (ER). These models can be used to prioritize chemicals for further tiered biological testing to assess the potential for endocrine disruption. One shortcoming of models predicting RBA has been the inability to distinguish potential receptor antagonism from agonism, and hence in vivo response. It has been suggested that steroid receptor antagonists are less compact than agonists; thus, ER binding of antagonists may prohibit proper alignment of receptor protein helices preventing subsequent transactivation. The current study tests the theory of chemical bulk as a defining parameter of antagonism by employing a 3-D structural approach for development of reactivity patterns for ER antagonists and agonists. Using a dataset of 23 potent ER ligands (16 agonists, 7 antagonists), molecular parameters previously found to be associated with ER binding affinity, namely global (E(HOMO)) and local (donor delocalizabilities and charges) electron donating ability of electronegative sites and steric distances between those sites, were found insufficient to discriminate ER antagonists from agonists. However, parameters related to molecular bulk, including solvent accessible surface and negatively charged Van der Waal's surface, provided reactivity patterns that were 100% successful in discriminating antagonists from agonists in the limited data set tested. The model also shows potential to discriminate pure antagonists from partial agonist/antagonist structures. Using this exploratory model it is possible to predict additional chemicals for their ability to bind but inactivate the ER, providing a further tool for hypothesis testing to elucidate chemical structural characteristics associated with estrogenicity and anti-estrogenicity.     

An integrated "4-phase" approach for setting endocrine disruption screening priorities--phase I and II predictions of estrogen receptor binding affinity

Recent legislation mandates the US Environmental Protection Agency (EPA) to develop a screening and testing program for potential endocrine disrupting chemicals (EDCs), of which xenoestrogens figure prominently. Under the legislation, a large number of chemicals will undergo various in vitro and in vivo assays for their potential estrogenicity, as well as other hormonal activities. There is a crucial need for priority setting before this strategy can be effectively implemented. Here we report an integrated computational approach to priority setting using estrogen receptor (ER) binding as an example. This approach rationally integrates different predictive computational models into a "Four-Phase" scheme so that it can effectively identify potential estrogenic EDCs based on their predicted ER relative binding affinity (RBA). The system has been validated using an in-house ER binding assay dataset for 232 chemicals that was designed to have both broad structural diversity and a wide range of binding affinities. When applied to 58,000 chemicals identified by Walker et al. as candidates for endocrine disruption screening, some 9100 chemicals were predicted to bind to ER. Of these, only 3600 were expected to bind to ER at RBA values up to 100,000-fold less than that of 17beta-estradiol. The method ruled out 83% of the chemicals as non-binders with a very low rate of false negatives. We believe that the same integrated scheme will be equally applicable to endpoints of other endocrine disrupting mechanisms, e.g. androgen receptor binding.     

An integrated "4-phase" approach for setting endocrine disruption screening priorities--phase I and II predictions of estrogen receptor binding affinity

Recent legislation mandates the US Environmental Protection Agency (EPA) to develop a screening and testing program for potential endocrine disrupting chemicals (EDCs), of which xenoestrogens figure prominently. Under the legislation, a large number of chemicals will undergo various in vitro and in vivo assays for their potential estrogenicity, as well as other hormonal activities. There is a crucial need for priority setting before this strategy can be effectively implemented. Here we report an integrated computational approach to priority setting using estrogen receptor (ER) binding as an example. This approach rationally integrates different predictive computational models into a "Four-Phase" scheme so that it can effectively identify potential estrogenic EDCs based on their predicted ER relative binding affinity (RBA). The system has been validated using an in-house ER binding assay dataset for 232 chemicals that was designed to have both broad structural diversity and a wide range of binding affinities. When applied to 58,000 chemicals identified by Walker et al. as candidates for endocrine disruption screening, some 9100 chemicals were predicted to bind to ER. Of these, only 3600 were expected to bind to ER at RBA values up to 100,000-fold less than that of 17 g -estradiol. The method ruled out 83% of the chemicals as non-binders with a very low rate of false negatives. We believe that the same integrated scheme will be equally applicable to endpoints of other endocrine disrupting mechanisms, e.g. androgen receptor binding.     

Experimental versus predicted affinities for ligand binding to estrogen receptor: iterative selection and rescoring of docked poses systematically improves the correlation

The computational determination of binding modes for a ligand into a protein receptor is much more successful than the prediction of relative binding affinities (RBAs) for a set of ligands. Here we consider the binding of a set of 26 synthetic A-CD ligands into the estrogen receptor ERα. We show that the MOE default scoring function (London dG) used to rank the docked poses leads to a negligible correlation with experimental RBAs. However, switching to an energy-based scoring function, using a multiple linear regression to fit experimental RBAs, selecting top-ranked poses and then iteratively repeating this process leads to exponential convergence in 4–7 iterations and a very strong correlation. The method is robust, as shown by various validation tests. This approach may be of general use in improving the quality of predicted binding affinities.     

Conformational studies on (17alpha,20Z)-21-(X-Phenyl)-19-norpregna-1, 3,5(10),20-tetraene-3,17beta-diols using 1D and 2D NMR spectroscopy and GIAO calculations of (13)C shieldings

Differences in agonist responses of the novel estrogen receptor ligands (17alpha,20Z)-(p-methoxyphenyl)vinyl estradiol (1), (17alpha, 20Z)-(o-alpha,alpha,alpha-trifluoromethylphenyl)vinyl estradiol (2), and (17alpha,20Z)-(o-hydroxymethylphenyl)vinyl estradiol (3) led us to investigate their solution conformation. In competitive binding assay studies, we observed that several phenyl-substituted (17alpha, 20E/Z)-(X-phenyl)vinyl estradiols exhibited significant estrogen receptor binding, but with variation (RBA (1) = 20; RBA (2) = 23; RBA (3) = 140 where estradiol RBA = 100) depending on the phenyl substitution pattern. Because the 17alpha-phenylvinyl substituent interacts with the key helix-12 of the ligand binding domain, we considered that differences in the preferred conformation of 1-3 could account for their varying binding affinity. 2D NMR experiments at 500 MHz allowed the complete assignment of the (13)C and (1)H spectra of 1-3. The conformations of these compounds in solution were established by 2D and 1D NOESY spectroscopy. A statistical approach of evaluating contributing conformers of 1-3 from predicted (13)C shifts correlated quite well with the NOE data. The 17alpha substituents of 1 and 2 exist in similar conformational equilibria with some differences in relative populations of conformers. In contrast, the 17alpha substituent of 3 exists in a different conformational equilibrium. The similarity in solution conformations of 1 and 2 suggests they occupy a similar receptor volume, consistent with similar RBA values of 20 and 23. Conversely, the different conformational equilibria of 3 may contribute to the significant binding affinity (RBA = 140) of this ligand.     

Reactivity profiles of ligands of mammalian retinoic acid receptors: A preliminary COREPA analysis

Retinoic acid and associated derivatives comprise a class of endogenous hormones that bind to and activate different families of retinoic acid receptors (RARs, RXRs), and control many aspects of vertebrate development. Identification of potential RAR and RXR ligands is of interest both from a pharmaceutical and toxicological perspective. The recently developed COREPA (COmmon REactivity PAttern) algorithm was used to establish reactivity profiles for a limited data set of retinoid receptor ligands in terms of activation of three RARs ( f , g , n ) and an RXR ( f ). Conformational analysis of a training set of retinoids and related analogues in terms of thermodynamic stability of conformers and rotational barriers showed that these chemicals tend to be quite flexible. This flexibility, and the observation that relatively small energy differences between conformers can result in significant variations in electronic structure, highlighted the necessity of considering all energetically reasonable conformers in defining common reactivity profiles. The derived reactivity patterns for three different subclasses of the RAR ( f , g , n ) were similar in terms of their global electrophilicity (nucleophilicity) and steric parameters. However, the profile of active chemicals with respect to interaction with the RXR- f differed qualitatively from that of the RARs. Variations in reactivity profiles for the RAR versus RXR families would be consistent with established differences in their affinity for endogenous retinoids, likely reflecting functional differences in the receptors.     

Reactivity profiles of ligands of mammalian retinoic acid receptors: a preliminary COREPA analysis

Retinoic acid and associated derivatives comprise a class of endogenous hormones that bind to and activate different families of retinoic acid receptors (RARs, RXRs), and control many aspects of vertebrate development. Identification of potential RAR and RXR ligands is of interest both from a pharmaceutical and toxicological perspective. The recently developed COREPA (COmmon REactivity PAttern) algorithm was used to establish reactivity profiles for a limited data set of retinoid receptor ligands in terms of activation of three RARs (alpha, beta, gamma) and an RXR (alpha). Conformational analysis of a training set of retinoids and related analogues in terms of thermodynamic stability of conformers and rotational barriers showed that these chemicals tend to be quite flexible. This flexibility, and the observation that relatively small energy differences between conformers can result in significant variations in electronic structure, highlighted the necessity of considering all energetically reasonable conformers in defining common reactivity profiles. The derived reactivity patterns for three different subclasses of the RAR (alpha, beta, gamma) were similar in terms of their global electrophilicity (nucleophilicity) and steric parameters. However, the profile of active chemicals with respect to interaction with the RXR-alpha differed qualitatively from that of the RARs. Variations in reactivity profiles for the RAR versus RXR families would be consistent with established differences in their affinity for endogenous retinoids, likely reflecting functional differences in the receptors.     

Androgen receptor binding affinity of pesticide "active" formulation ingredients. QSAR evaluation by COREPA method

The COREPA approach for identifying the COmmon REactivity PAttern of biologically similar chemicals was employed to upgrade the recently derived affinity pattern for high androgen receptor (AR) binding affinity. The training set consisted of 28 steroidal and nonsteroidal ligands whose AR binding affinity was determined in competitive binding assays (in terms of pKi). The interatomic distances between nucleophilic sites and their charges providing distinct and non-overlapping integral patterns for active and inactive chemicals were assumed that it was related with the endpoint, which was under study. These stereoelectronic characteristics were used to predict pKi values of pesticide "active" formulation ingredients in an attempt to identify chemicals with potential AR binding affinity.     

Androgen receptor binding affinity of pesticide "active" formulation ingredients. QSAR evaluation by COREPA method

The COREPA approach for identifying the COmmon REactivity PAttern of biologically similar chemicals was employed to upgrade the recently derived affinity pattern for high androgen receptor (AR) binding affinity. The training set consisted of 28 steroidal and nonsteroidal ligands whose AR binding affinity was determined in competitive binding assays (in terms of p K i ). The interatomic distances between nucleophilic sites and their charges providing distinct and non-overlapping integral patterns for active and inactive chemicals were assumed that it was related with the endpoint, which was under study. These stereoelectronic characteristics were used to predict p K i values of pesticide "active" formulation ingredients in an attempt to identify chemicals with potential AR binding affinity.