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.     

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 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.     

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.     

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.     

Molecular modeling of class I and II alleles of the major histocompatibility complex in <Emphasis Type="Italic">Salmo salar</Emphasis>

Knowledge of the 3D structure of the binding groove of major histocompatibility (MHC) molecules, which play a central role in the immune response, is crucial to shed light into the details of peptide recognition and polymorphism. This work reports molecular modeling studies aimed at providing 3D models for two class I and two class II MHC alleles from Salmo salar (Sasa), as the lack of experimental structures of fish MHC molecules represents a serious limitation to understand the specific preferences for peptide binding. The reliability of the structural models built up using bioinformatic tools was explored by means of molecular dynamics simulations of their complexes with representative peptides, and the energetics of the MHC-peptide interaction was determined by combining molecular mechanics interaction energies and implicit continuum solvation calculations. The structural models revealed the occurrence of notable differences in the nature of residues at specific positions in the binding groove not only between human and Sasa MHC proteins, but also between different Sasa alleles. Those differences lead to distinct trends in the structural features that mediate the binding of peptides to both class I and II MHC molecules, which are qualitatively reflected in the relative binding affinities. Overall, the structural models presented here are a valuable starting point to explore the interactions between MHC receptors and pathogen-specific interactions and to design vaccines against viral pathogens.     

A rule-based expert system for chemical prioritization using effects-based chemical categories

A rule-based expert system (ES) was developed to predict chemical binding to the estrogen receptor (ER) patterned on the research approaches championed by Gilman Veith to whom this article and journal issue are dedicated. The ERES was built to be mechanistically transparent and meet the needs of a specific application, i.e. predict for all chemicals within two well-defined inventories (industrial chemicals used as pesticide inerts and antimicrobial pesticides). These chemicals all lack structural features associated with high affinity binders and thus any binding should be low affinity. Similar to the high-quality fathead minnow database upon which Veith QSARs were built, the ERES was derived from what has been termed gold standard data, systematically collected in assays optimized to detect even low affinity binding and maximizing confidence in the negatives determinations. The resultant logic-based decision tree ERES, determined to be a robust model, contains seven major nodes with multiple effects-based chemicals categories within each. Predicted results are presented in the context of empirical data within local chemical structural groups facilitating informed decision-making. Even using optimized detection assays, the ERES applied to two inventories of >600 chemicals resulted in only ~5% of the chemicals predicted to bind ER.     

Characterization of hydroxycinnamic acid derivatives binding to bovine serum albumin

Hydroxycinnamic acid derivatives (HCAs) are a group of naturally occurring polyphenolic compounds which possess various pharmacological activities. In this work, the interactions of bovine serum albumin (BSA) with six HCA derivatives, including chlorogenic acid (CHA), caffeic acid (CFA), m-coumaric acid (m-CA), p-coumaric acid (p-CA), ferulic acid (FA) and sinapic acid (SA) have been investigated by NMR spectroscopic techniques in combination with fluorescence and molecular modeling methods. Competitive STD NMR experiments using warfarin sodium and L-tryptophan as site-selective probes indicated that HCAs bind to site I in the subdomain IIA of BSA. From the analysis of the STD NMR-derived binding epitopes and molecular docking models, it was deduced that CHA, CFA, m-CA and p-CA show similar binding modes and orientation, in which the phenyl ring is in close contact with protein surface, whereas carboxyl group points out of the protein. However, FA and SA showed slightly different binding modes, due to the steric hindrance of methoxy-substituents on the phenyl ring. Relaxation experiments provided detailed information about the relationship between the affinity and structure of HCAs. The binding affinity was the strongest for CHA and ranked in the order CHA > CFA > m-CA ≥ p-CA > FA > SA, which agreed well with the results from fluorescence experiments. Based on our experimental results, we also conclude that HCAs bind to BSA mainly by hydrophobic interaction and hydrogen bonding. This study therefore provides valuable information for elucidating the mechanisms of BSA-HCAs interaction.     

Chemometrical classification of ephrin ligands and Eph kinases using GRID/CPCA approach

Eph receptor tyrosine kinases are divided on two subfamilies based on their affinity for ephrin ligands and play a crucial role in the intercellular processes such as angiogenesis, neurogenesis, and carcinogenesis. As such, Eph kinases represent potential targets for drug design, which requires the knowledge of structural features responsible for their specific interactions. To overcome the existing gap between available sequence and structure information we have built 3D models of eight ephrins and 13 Eph kinase ligand-binding domains using homology modeling techniques. The interaction energies for several molecular probes with binding sites of these models were calculated using GRID and subjected to chemometrical classification based on consensus principal component analysis (CPCA). Despite inherent limitations of the homology models, CPCA was able to successfully distinguish between ephrins and Eph kinases, between Eph kinase subfamilies, and between ephrin subfamilies. As a result we have identified several amino acids that may account for selectivity in ephrin-Eph kinase interactions. In general, although the difference in charge between ephrin and Eph kinase binding domains creates an attractive long-range electrostatic force, the hydrophobic and steric interactions are highly important for the short-range interactions between two proteins. The chemometrical analysis also provides the pharmacophore model, which could be used for virtual screening and de novo ligand design.     

Study of the interaction of aglycon of daunorubicin with human serum albumin by spectroscopy and modeling

The interaction between aglycon of daunorubicin (DNR-A) and human serum albumin (HSA) was investigated using fluorescence quenching and modeling. Results shown that fluorescence quenching of HSA by DNR-A resulted from the formation of DNR-A-HSA complex. The quenching constants were determined via measurement of the binding affinity between DNR-A and HSA using the Stern-Volmer equation. The thermodynamic parameters DeltaG, DeltaH, DeltaS and the binding distance r were calculated. Furthermore, SFS and UV spectra suggested that the complex changed the conformation of HSA and that hydrophobic interactions played a major role in DNR-A-HSA association, which was in good agreement with the results of the modeling study. Moreover, the SFS technique was successfully applied to determine the total proteins in biology samples with satisfactory results.     

Double affinity amplification of galectin-ligand interactions through arginine-arene interactions: synthetic, thermodynamic, and computational studies with aromatic diamido thiodigalactosides

A series of aromatic mono- or diamido-thiodigalactoside derivatives were synthesized and studied as ligands for galectin-1, -3, -7, -8N terminal domain, and -9N terminal domain. The affinity determination in vitro with competitive fluorescence-polarization experiments and thermodynamic analysis by isothermal microcalorimetry provided a coherent picture of structural requirements for arginine-arene interactions in galectin-ligand binding. Computational studies were employed to explain binding preferences for the different galectins. Galectin-3 formed two almost ideal arene-arginine stacking interactions according to computer modeling and also had the highest affinity for the diamido-thiodigalactosides (K(d) below 50 nM). Site-directed mutagenesis of galectin-3 arginines involved in binding corroborated the importance of their interaction with the aromatic diamido-thiodigalactosides. Furthermore, the arginine mutants revealed distinct differences between free, flexible, and solvent-exposed arginine side chains and tightly ion-paired arginine side chains in interactions with aromatic systems.