Next Article in Journal
The HMOX1 Pathway as a Promising Target for the Treatment and Prevention of SARS-CoV-2 of 2019 (COVID-19)
Next Article in Special Issue
G Protein-Coupled Estrogen Receptor Mediates Cell Proliferation through the cAMP/PKA/CREB Pathway in Murine Bone Marrow Mesenchymal Stem Cells
Previous Article in Journal
GANT61 and Lithium Chloride Inhibit the Growth of Head and Neck Cancer Cell Lines Through the Regulation of GLI3 Processing by GSK3β
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens

by
Anna Helena Mazurek
1,
Łukasz Szeleszczuk
1,*,
Thomas Simonson
2 and
Dariusz Maciej Pisklak
1
1
Chair and Department of Physical Pharmacy and Bioanalysis, Department of Physical Chemistry, Medical Faculty of Pharmacy, University of Warsaw, Banacha 1 str., 02-093 Warsaw Poland
2
Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique, 91-120 Palaiseau, France
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(17), 6411; https://doi.org/10.3390/ijms21176411
Submission received: 21 July 2020 / Revised: 30 August 2020 / Accepted: 1 September 2020 / Published: 3 September 2020
(This article belongs to the Special Issue The Role of Estrogen Receptors in Health and Disease)

Abstract

:
In this review, applications of various molecular modelling methods in the study of estrogens and xenoestrogens are summarized. Selected biomolecules that are the most commonly chosen as molecular modelling objects in this field are presented. In most of the reviewed works, ligand docking using solely force field methods was performed, employing various molecular targets involved in metabolism and action of estrogens. Other molecular modelling methods such as molecular dynamics and combined quantum mechanics with molecular mechanics have also been successfully used to predict the properties of estrogens and xenoestrogens. Among published works, a great number also focused on the application of different types of quantitative structure–activity relationship (QSAR) analyses to examine estrogen’s structures and activities. Although the interactions between estrogens and xenoestrogens with various proteins are the most commonly studied, other aspects such as penetration of estrogens through lipid bilayers or their ability to adsorb on different materials are also explored using theoretical calculations. Apart from molecular mechanics and statistical methods, quantum mechanics calculations are also employed in the studies of estrogens and xenoestrogens. Their applications include computation of spectroscopic properties, both vibrational and Nuclear Magnetic Resonance (NMR), and also in quantum molecular dynamics simulations and crystal structure prediction. The main aim of this review is to present the great potential and versatility of various molecular modelling methods in the studies on estrogens and xenoestrogens.

1. Introduction

Successful applications of molecular modelling methods can be found in almost every branch of modern physics, chemistry, and biology. This versatility and popularity results from the constantly increasing computing power of both personal computers and specialized servers as well as the availability of molecular modelling software. The number of properties that can be accurately predicted and phenomena that can be explained as well as problems that can be solved using such calculations are enormous. Therefore, in this review, recent advances in molecular modelling applications to study the chemistry and biochemistry of estrogens [1] and xenoestrogens [2] will be presented. The aim of this article is not only to present the large volume of information concerning title compounds that have been obtained in recent years using in silico methods but also to present to researchers who are not specialized in molecular modelling methods the possible applications of these compounds and, in that way, encourage them to use such calculations in their own studies.
This review is organized as follows: first the title compounds, estrogens, and xenoestrogens are briefly summarized, with particular focus on those that have already been objects of computational studies. Then, the most important biomolecules that are involved in the metabolism and action of estrogens and xenoestrogens are described. In the main part of this review, the computational methods that have been employed in the studies on estrogens and xenoestrogens are presented. Each of the methods is briefly described, without too many details, as there are many very good reviews referenced in this article focusing on the basics of particular methods. Wherever possible, the computational results were compared to the corresponding experimental ones; however, in many cases, such comparison was impossible either due to the lack of experimental results or the purely theoretical character of the published work. In the next section, a critical analysis of the reviewed methods is presented, supported by the presentation of some technical aspects of the reviewed studies. This was done to facilitate the choice of a certain method or its properties. From the authors’ perspective, the number of studies on estrogens and xenoestrogens is constantly increasing; however, only in some of them are the experiments supported by suitable theoretical studies. Therefore, our aim was to present to researchers working with estrogens and xenoestrogens selected computational tools that can be used to facilitate and improve their studies.

1.1. Estrogens and Xenoestrogens: Types, Main Representatives, and Their Toxicity

Estrogens are a group of natural steroid sex hormones. There are four of them: estrone (E1), estradiol (E2), estriol (E3), and estretrol (E4) [1] (see the Supplementary Materials: Figure S1). The last one is produced only during pregnancy by the fetus liver [3]. Among these four compounds, E2 plays the most important role in the human organism and, therefore, is of high importance in breast or ovarian cancer progression. With regard to their relative binding affinity (RBA) to Estrogen Receptor (ER), with the exclusion of estretrol, estrogens are ranked in the following order: estradiol > estrone > estriol. In comparison to estradiol, the activity and potency of estrone and estriol are, respectively, 10 and 100 times lower than that of E2 [4]. E1 mainly functions as estradiol’s metabolite and, at the same time, serves as its precursor (the estrone-to-estradiol transformation is reversible) [5]. On the other hand, in non-pregnant women, estriol levels in the blood are hardly detectable, whereas during pregnancy its amount distinctly grows because it is produced by the placenta as well [6]. All estrogens are used as medication in menopausal hormone therapy, although estradiol is the most applied [7,8]. In such external applications, they should be considered as xenoestrogens.
Apart from estrogens, other non-endogenous substances, called endocrine-disrupting chemicals (EDCs), can bind to the ER. According to the European Commission Regulation from 2018 [9], EDCs are substances that have adverse effects in non-target organisms, have an endocrine mode of action, and exhibit these adverse effects as a consequence of an endocrine mode of action. In the Guidance on the identification and studies regarding ECDs, published by the European Chemical Agency (ECA) and the European Food Safety Authority (EFSA) [10], the importance of in silico studies in the process of EDC research is clearly pointed out.
More detailed information on in silico examination regarding EDCs is provided in OECD (Organization for Economic Co-operation and Development) Guidance (Update v3, 2017) [11], where Quantitative Structure–Activity Relationship (QSAR) methodologies, prediction of the metabolic transformation (ADME), and CYP450 metabolism investigations are listed as Level 1 methods in the identification and study process. It is emphasized that applications of such methods can be used to identify the groups of chemicals and structural characteristics that are responsible for the observed in vivo effects and can serve as a tool to explain possible differences between in vitro and in vivo results, clarifying EDCs’ mechanism of action. In the Conceptual Framework for Testing and Assessment of EDCs [12], at investigation Level 1, not only QSAR and ADME but also ‘other in silico model predictions’ have been listed. Additively, special QSAR Guidance on the topic has been published.
A large group of EDCs that are well-known and have been investigated for decades are xenoestrogens, which are xenobiotics that cause either an estrogenic or an antiestrogenic effect [13]. Nowadays, they are widespread and originate from different sources (see the Supplementary Materials: Figure S2). In the in silico studies, the main emphasis has been put, so far, on the investigation of pharmaceuticals, phytoestrogens, bisphenol A, and phthalates.

1.2. Estrogen-Related Biomolecules as the Molecular Modelling Study Objects

Numerous studies involving various in silico methods, but mainly ligand docking, have been performed to look for new inhibitors of the enzymes metabolizing estradiol. These proteins are a possible target for new drugs and are depicted in Table 1. Most of the research deals with aromatase [14,15,16,17,18,19,20,21,22], as its inhibitors, like letrozole, are already in medicinal use [23]. However, those inhibitors exhibit some detrimental side effects, which is the reason why research continues on this topic. 17β-Hydroxysteroid dehydrogenase (17β-HSD) [24,25], sulfatase (STS) [26], as well as sulfotransferase (SULT) [27] have also been taken under close consideration using molecular modelling in order to contribute to anti-cancer drug development. The first of the two enzymes transforms estrone into estradiol. The second inactivates estradiol by transforming it into a sulfated form.
Even if conceptually estradiol is central to the above-mentioned studies, none of them used it as a target ligand. Nevertheless, in order to model estradiol activity metabolism in molecular studies, the most common objects are estradiol itself, its natural receptor—the estrogen receptor (ER)—and its main hydroxylating enzymes: CYP1A1 [28], CYP1A2 [29], CYP1B1 [30], and SHBG (sex hormone-binding globulin). The last one is a protein that can bond estradiol (E2); thus, it has direct impact on the amount of free E2 in plasma and, consequently, on the hormone’s bioavailability [31]. Taking into account the interest of in silico research for the enzymes listed in this paragraph, only topics associated with these systems are covered in the main part of this review.
Two subtypes of the estrogen receptor are known, ERα and ERβ, with tissue-specific expression [32]. ERα is found mainly in the mammary gland, uterus, ovary (thecal cells), male reproductive organs, prostate, liver, and adipose tissue. ERβ is present in the prostate, bladder, ovary (granulosa cells), colon, adipose tissue, and immune system [33]. Despite being encoded by different genes, both estrogen receptors show high homology, and in both of them the E domain contains the ligand-binding domain (LBD) and C domain, which is the DNA binding domain (DBD) [34]. The homology between ERα and ERβ in DBD is more than 95% [35]. In LBD the homology is about 55% [35]. In ER structures, two transactivation functions are present, called AF-1 (located in N-terminal domain) and AF-2 (located in LBD) [36]. They contain the nuclear location signals and, after proper exposure of their surface, are responsible for incorporation of the co-activators, which is a necessary step to induce activation of the intercellular signaling pathways.
Binding of 17β-estradiol and any agonist to the ER requires creating a hydrogen bond with His524 (in ERα) [37] or His475 (in ERβ) [38]. This leads to a unique agonist-bound conformation of the receptor’s LBD, characterized by a specific repositioning of the H12 helix, which is the most C-terminal helix of the LBD (molecular switch) [39]. On the contrary, selective ER modulators (SERMs) such as raloxifene or tamoxifen induce relocation of H12 into the co-activator binding cleft, which blocks AF-2 activity. Finally, pure antagonists completely destabilize H12 [40].
All the above actions concerning estrogens’ binding to ER are genomic actions. This means that they require translocation of the estrogen–ER complex to the nucleus and interaction with chromatin at specific sequences, known as estrogen response elements [41]. There are other estrogen signaling paths that are non-genomic and involve indirect regulation of gene expression [42]. They include activation of various protein-kinase cascades after binding of an estrogen molecule to a membrane receptor, usually GPER1 (G protein-coupled ER1) [43]. It has been proven that estrogen binding to GPER1 shows similarity to estrogen–ER binding; however, estrogen’s affinity for GPER1 is significantly lower than for ER [43]. Nevertheless, all other steroid hormones are characterized by even lower affinities towards GPER1.
Experimental analysis of the GPER1 structure has been limited, as the protein is rather refractory to both X-ray crystallography and NMR due to its relatively high lipophilicity [44]. This explains the importance of computational homology modelling [45]. Homology models are later used to simulate binding with estrogens. There are only a few studies dealing directly with this subject [44,46]. Application of homology modelling, molecular docking, and molecular dynamics provides insight into the process of estrogen binding and helps to explain the induced non-genomic effects. As non-genomic estrogen actions through GPER1 and kinase cascades alter the cell membrane shape, further research on cell membranes and estrogens is also performed (see Section 2.1.3).

2. Application of Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens

Molecular modelling could be of a great help to experimentalists. Ligand docking directs research to the most probable hit molecules; QSAR, an officially accepted OECD method, predicts toxicity and often can facilitate research; and MD allows one to observe the time-dependent mechanisms (e.g., in membranes) and, therefore, helps to explain the experimentally observed phenomena. What is more, two types of theories on which calculations are based deliver two levels of accuracy. These are QM and MM. The former is more accurate, but as a consequence, calculations are more time-consuming. The latter, on the contrary, allows one to obtain a general view on the topic in a shorter time. The diversity of methods highly contributes to more successful experimental studies as it saves time and points out these research approaches that have a high probability of success.

2.1. Application to Estrogens

2.1.1. Ligand Docking Using Force Field Methods

  • Principles of docking and re-docking
Protein–ligand docking is a technique used to predict the orientation and active conformation of a molecule in an active center. Based on this prediction, the binding energy between protein and ligand is calculated [47]. Ligand docking can identify the chemical bonds crucial for activity and the specific atoms/residues that are responsible for them. Target proteins are usually acquired from the RCSB Protein Data Bank (PDB) [48] (Table 2). Presence of a ligand in the crystallographic receptor structure simplifies and speeds up researchers’ work, as a proper region for docking is already plainly indicated.
Accessibility of the X-ray protein structures with docked ligands gives also possibility to prove the quality of the chosen simulating docking method and parameters [49]. This is performed via extracting the ligand from the available in the PDB structure and re-docking it to this crystallographic measurements-based protein. The results are obtained in a form of the root mean square deviation (RMSD) of atom positions [49]. It is calculated as differences between the original experimental atom positions in the crystallized structure and the theoretically obtained molecular docking results [50]. It is acknowledged that the used modelling tools are able to identify the correct ligand pose, are repeatable and reliable when RMSD is smaller than 2Å [51].
Molecular modelling performed on estrogen-related receptors is not an exception. As the studies are performed mostly to search for new agonists or antagonists, the ligands re-docked into the ERs or SHBG are mainly E2 [52,53,54,55,56,57] and 4-hydroxytamoxifen [52,58,59]. RMSD value in most of the studies varies from 0.26 to 1.4Å, which proves the correctness of the docking methods in finding the proper orientation of the ligand in the active site.
  • Enzymes and receptors used as targets in estrogen-related docking studies
A great many proteins have been used in molecular docking studies of estradiol. They include 17β-hydroxysteroid dehydrogenase [60], progesterone receptor [61], protein disulfide isomerase [62], SHBG [63], CYP1B1 [64], steroid sulfatase [65], and even the voltage- and Ca2+-activated K+ channel β1 subunit [66]. In all studies, an emphasis is put on the importance of OH-hydrogen binding [67,68]. The relevance of this emphasis is confirmed by experimental results. One example is a study where the estrogen analog with the highest affinity, measured in a fluorescence polarization displacement assay, appears to have the second highest predicted affinity [67].
Docking serves either to investigate the binding of estradiol to the ER or to dock other molecules, potentially new drugs, such as potent and highly selective estriol analogues [69]. In the second case, information from previous studies on estradiol behavior in an active site serves as reference data [70,71,72,73]. This is a starting point for the commonly applied research sequence: investigation of E2 binding mode, comparison with the calculation results for potential drugs, and confirmation of the hypothesis by analytical techniques. Such a three-step process has been performed, for example, for estrogen-dependent MCF-7 cancer cells [70]. Eighteen compounds with antiproliferative activity have been docked into the cavity where E2 normally binds. It was predicted that, compared to E2, an additional aromatic ring is involved in the binding mechanism. The prediction was confirmed by site-directed mutagenesis.
In some cases, allosteric modulators have been studied. This means that estradiol must be present in the binding site during the docking process [74,75]. To test the prediction of allosteric activity, in vivo experiments are often performed. For example, in one study [74] it was shown that one particular compound can properly fit into the region of the binding pocket, along with E2. Afterwards, this ligand was investigated in vivo and was demonstrated, indeed, to be a new ER-beta-selective, negative allosteric modulator of E2 binding.
Including an estradiol molecule in the docking studies helps to properly score the resulting data [76] and rank the tested molecules according to their binding affinity [77]. For large ligand sets, high-throughput screening with a pure agonist (estradiol), an antagonist (tamoxifen), and decoys (known non-binders) can be performed. This can aid in setting a laboratory’s experimental testing priority, reducing the cost and time of its research, and boosting its effectivity [77]. This explains why in vivo and in vitro investigations are often combined with in silico ones [78,79]. Altogether, it enables the discovery of new possible drug molecules that could influence estradiol’s signaling pathways [80].
  • Docking studies of plant-derived potential xenoestrogens
Docking plant-derived substances into ER and comparing their binding energies and interactions (above all hydrogen bonding) with those of estradiol is a common practice. Such studies deliver information on conformational flexibility [81], the ability of the investigated molecules to selectively modulate ERα/β ability [82,83], their possible reductive influence on breast cancer risk [84], or their applicability to reduce menopausal symptoms [85]. Most importantly, data derived in this way very often show good consistency with experiments [86]. For the obtained data, a correlation with estradiol, but sometimes also with tamoxifen [75] (ER antagonist) or genistein [87,88], is found. This last substance serves as an important reference, as genistein has a higher affinity for ERβ than 17β-estradiol.
Even if docking into ER is most widely used, modulation of CYP450 activity by plant-derived substances with regard to estrogenic effects has also been examined [89]. Indeed, the same signaling pathways are regulated by CYP450–estradiol interactions, and there are plenty of data available on this topic. Moreover, other calculation methods are also applied, including ADMET and molecular dynamics (MD) [90,91] (for MD description and examples, see Section 2.1.3). Nowadays, in silico methods are a standard tool in plant xenoestrogens studies. Thus, molecular modelling is often performed in parallel to either in vitro [92] or in vivo toxicity studies [93], and in most cases an agreement between data obtained from both sources is found. This was the case in the docking-based binding affinities of compounds derived from C. elegans and their measured reproductive toxicity [93]. Wide examination of many substances of a natural origin is possible thanks to the thorough structural knowledge on the estradiol molecule bound to the ER (Table 2).

2.1.2. Quantitative Structure–Activity Relationship (QSAR)

QSAR methods use mathematical models to correlate structural characteristics with the biological activity of a set of compounds that have closely related structures [94]. Empirical and theoretical molecular descriptors are used. Conceptually, three main types of QSARs are known: ones based on fragment analysis of a system (here, a pharmacophore [95] concept is used), ones that consider the given system as a whole (descriptors are computed from scalar quantities), and 3D-QSAR. In this last type, descriptors are obtained by application of a force field (3D approach). To achieve a 3D target structure, either software-based alignment or manual superimposition on the crystallographic data must be performed [96]. This is a necessary step, as different ligand-binding modes and bioactive conformations are possible. Examples of 3D-QSAR are Comparative Molecular Field Analysis (CoMFA) [97] and Comparative Molecular Similarity Indices Analysis (CoMSIA) [98]. From the created QSAR models, predictions on the bio-activity and toxicity of other molecules are made.
A smaller group of structure–activity studies are ones that are non-quantitative, namely SAR ones (structure–activity relationship). They include virtual screening with a pharmacophore and a large set of molecules. Afterwards, in vitro evaluation of the data is performed. This methodology has been used to determine potential 17β-HSD inhibitors [99,100] that could be applied in the treatment of osteoporosis provoked by estradiol deficiency.
Among all molecular modelling approaches, QSAR is one of the most commonly used to examine estradiol’s structure and activity. Often, data on estradiol serve only as a reference for information gathered on new possible drugs [14,101] such as raloxifene derivatives [102] or estradiol metabolites. Through comparison with experiments, it has been shown that QSAR models have a high sensitivity and specificity in providing relative binding affinities (RBAs), where estradiol’s RBA equals 100% [103]. Among the most important descriptors are ones calculated with quantum mechanics at the DFT level [104,105].
While QSAR is mostly applied to ER-binding, other proteins like CYP1A2 [106] and CYP1B1 [107] have also been targeted. One study used the CoMFA approach to model estriol’s influence on CYP1A1 [106]. It focused on the inhibition of estradiol to mutagenic 4-OH estradiol transformation. Out of 90 steroid candidates, thioestrone was selected and shown to have the desired inhibitory ability. Its mechanism of action was revealed by molecular modelling. It is desired because thioestrone’s -SH group is closer to the iron atom in the CYP1B1 heme than the -OH group in natural ligands of the enzyme, namely estradiol and estrone.
OECD Guidance identifies QSAR as an important element in toxicity evaluations. As a result, QSAR plays an important role in substance risk assessment [108]. It has been applied in studies dealing with mutations in enzymes that metabolize estradiol [109] and to study estradiol oxidation and emerging contaminants [110]. In the latter study, the most accurate DFT descriptors were used. This enabled a better understanding of the degradation mechanisms.

2.1.3. Advanced Docking Using Combined Quantum Mechanics/Molecular Mechanics (QM/MM) or Molecular Dynamics (MD) Methods

Molecular modelling with QM is much more accurate than MM methods. However, the direct use QM approaches in drug design is limited due to the cost and size of protein structures. In recent years, QM/MM has been gaining attention, as it allows one to consider a whole protein–ligand complex and not only a binding site [111,112]. Due to computational limitations, MM calculations on the outer part of a receptor deliver only approximate data. Nevertheless, QM/MM provides more knowledge on the protein’s influence than is obtained when only the LBD region is considered.
More commonly applied methodology is MD. It simulates time-dependent processes and provides data that are otherwise unavailable [52,54]. For example, in case of ERs or SHBG ligands MD can point out whether the examined substance is receptor’s agonist or antagonist. This assumption is based on the RMSF (root mean square fluctuation) value which is extracted from the MD trajectories. RMSF represents the flexibility of the amino acid residues [50]. Both RMSF and RMSD depend on the interactions between the protein and the ligand and are the result of the ligand movements in the active site trying to achieve the appropriate position [63]. If the same ligand undergoes the MD process in both ERα and ERβ and its RMSF and RMSD values significantly differ for these two receptors, it indicates that the investigated molecule probably occupies more favorably one of the investigated receptors. Such data leads to the hypothesis of the ligand’s selectivity.
Moreover, in terms of the ERs, RMSF and RMSD values suggest whether the analyzed molecule is the receptor’s agonist or antagonist [50,52,53,54,59]. As already mentioned in Section 1.2, the positioning of the H12 helix is differently influenced by agonists and antagonists. If the estrogenicity of the compound is known, comparison of RMSF values with the known molecule’s estrogenic effect, can serve as an evaluation of the applied docking parameters.
Nevertheless, MD requires a lot of computation time, which increases with the system’s size and simulation length. In return, it describes the dynamics of the system as well as entropic effects associated with the protein–ligand interaction.
MD can also provide quantitative estimates of relative binding affinities through techniques known as “free energy calculations”. The most rigorous is the “free energy perturbation,” or FEP family of methods [113,114,115,116,117]. To compare two ligands, say A and B, one introduces a model where both ligands are present, each with a partial occupancy. This is closely analogous to a crystallographic refinement where a particular group (ligand or side chain, say) has two possible conformations. Each ligand is assigned a weight between 0 and 1, say wA and wB = 1 − wA. These multiply terms in the energy function involving either ligand. By varying the weights gradually, one can effectively remove one ligand and introduce the other. Thus, when wB = 0, ligand A is fully weighted while interactions of ligand B with its surroundings have a zero weight: ligand B is “invisible” to its environment. Usually, interactions within the ligand are not weighted. As wB changes from zero to one, B is introduced and A is removed. Intermediate weight values correspond to an “alchemical” mixture, where both ligands are partially present. MD simulations (or Monte Carlo simulations) are performed for a few wA values, typically around 10. This series of simulations mimics a gradual, reversible replacement of A by B. Energy statistics are collected from all simulations. From these, a free energy difference between A and B can be obtained. The same process is carried out for the unbound ligands, solvated by a box of water. Subtracting the bound and unbound free energy changes yields the binding free energy difference. The method requires force field parameters for each ligand but has no other adjustable parameters. The tradeoff is that not one, but several MD simulations are required, and these should be sufficiently long. Indeed, FEP accuracy is limited by the MM force field, but also the amount of conformational sampling that is carried out. The method can also be used to compute the binding free energy changes due to point mutations of the protein. In this case, a particular residue is modeled with two side chains, each having a partial occupancy. Nowadays, FEP can be applied to one ligand or mutation per day on a medium-sized computer cluster.
A simplified version of FEP is to use only two simulations per ligand: one bound to the protein and one in solution. From these, binding free energies can also be obtained, if one is willing to extrapolate from a wA value of unity to a value of zero. To counter the use of such a large extrapolation, one introduces empirical weighting factors that multiply the interactions between the ligand and its surroundings (protein or solution). Usually one is applied to the electrostatic interactions and one to the Lennard–Jones interactions. The extrapolation and use of interaction energies have led to the name Linear Interaction Energy, or LIE method [118,119,120]. For a thorough review of its theoretical basis, see [121]. The tradeoff for its speed is that experimental data are needed to adjust the values of the empirical weights, which are not very transferable between different proteins and classes of ligands. Once the weights are optimized for the molecules of interest, predictions can be made.
In many studies, estradiol has been simulated in a complex with ER. The QM/MM approach seems to be the most accurate. It uses a QM description of the ligand and its binding pocket, but thanks to the MM description of more distant protein regions, it preserves information on the whole enzyme’s impact [122]. Very recently, QM/MM elucidated the important role of estradiol’s D-ring in the active site of ERα [123]. What is more, it helped to understand the influence of each enzyme segment on the ERα-agonist (estradiol, diethylsilbestrol) binding [124]. Most importantly, the binding energies of E2 and DES correlated well with experimental agonist binding affinities for the ER.
To analyze the changes in the investigated complexes upon DNA binding, not only MM, as with estradiol and DNA [125], but also MD has been applied [126]. The latter study identified specific bases within the aptamer (short-stranded DNA/RNA, binds only specific molecules [127]) and demonstrated the importance of water-mediated hydrogen bonds in the aptamer–estradiol complex.
MD is also often used to describe the effects of enzyme mutations on ligand binding. This methodology has been applied in estradiol studies, to compare wildtype and mutated CYP1B1 [128], which is mostly responsible for the 2-OH-hydroxylation of estradiol. MD serves also as a tool to explain the results of molecular docking. For example, it has been stated that the interactions between human α-fetoprotein and agonists (estradiol, estrone, diethylsilbestrol) were caused by van der Waals forces, whereas binding of antagonists (tamoxifen and its analogues) was equally based on hydrophobic and electrostatic interactions [129]. Moreover, information from MD simulations can be used to construct a pharmacophore in order to screen protein databases for a desired type of ligand. This methodology was used to search for substrates and inhibitors of the estrone-SULT [130]. Nine selected molecules were consistent with the ones indicated by the experiment.
Another nuclear receptor, the farnesoid X receptor (FXR), was used as a test case for free energy calculations in 2018 [116,117,120]. FXR is involved in regulating bile acid, lipid, and glucose homeostasis [131], and it has been linked to hepatocarcinogenesis [132]. Its hormone binding site is hydrophobic with few conserved interaction motifs and strong induced fit effects. With FEP, mean errors were about 1.5 kcal/mol for relative binding free energies of around 30 ligands, and the largest errors were about 2.5 kcal/mol. In one of the studies [117], ligands were first docked to the receptor, then compared using FEP. LIE gave similar errors for the same protein and 47 ligands [120], at a lower computational cost, but required optimization of the two adjustable LIE parameters using a subset of the ligands.
An earlier LIE study [119] considered the ER binding of estradiol and a series of xenoestrogens. A training set of 19 ligands was used to optimize the LIE parameters. A mean unsigned error of 0.6 kcal/mol was then obtained for a test set of 13 ligands. Several binding poses (3–4) were considered for each ligand; this was not too expensive because only short MD simulations were run.

2.1.4. Other MD-Based Studies of Estrogens

  • Membranes
As ERs are located in either the cytoplasm or the lipid bilayer (mER, membrane ER) [27,133], and estradiol itself is a steroid hormone, closer insight into the ligand’s interaction with the cell membrane is an obvious research target. Since transfer through this cellular barrier is not a stable state, but a dynamic process, MD calculations could be seen as a method of choice. Recent studies revealed [134] that, regarding estradiol’s long axis and the lipid acyl chains, E2 adopts a perpendicular position in the membrane. By having four rings located near the membrane interface, participation of the hormone’s 3-OH and 17β-OH groups in hydrogen bonds and electrostatic interactions with the lipids are possible.
Combining MD and QM enables further research into E2 membrane crossing. It provides information not only about the E2 orientation in the membrane (MD) but also about the strength of the electrostatic potential mapped on the electron density surface (QM) [135]. These data, derived from a HDL disc model, enable deeper insight into the mechanism of E2 incorporation into lipid membranes and is an important step forward to develop tissue-specific discs encircled by a membrane, which would serve as transporters for E2 or its derivatives.
Similar methodology (QM, MD) was used to study the removal of hormonal pollutants from water. It has been shown that by using high levels of salinity, which increases the strength of hydrogen bonding and hydrophobic interactions, one can perform a membrane-based sorption of 17α-ethinyl estradiol on the polyethersulfone membrane [136].
  • Nanotubes
Nowadays, estradiol is a relatively common water pollutant, and numerous studies have been performed to find a reliable tool for its removal. One possibility is to use nanotubes. For that purpose, within the molecular modelling approach, mostly single-walled carbon nanotubes (SWNTs) [137] are used. Free energies of adsorption have been calculated with a QM approach. In some cases, MD was also implemented [138]. The target ligands for these studies were 17β-estradiol and its medically useful derivative 17α-ethinyl-estradiol [139,140,141]. These simulations revealed the adsorption energy, a preferential sorption among different nanotubes and estradiol derivatives, and provided a molecular explanation for the observed results [142].

2.1.5. Density Functional Theory (DFT) Calculations in the Study of Estrogens

DFT is a QM approach that determines the ground-state properties of a many-body system by applying the electron density concept. The underlying concept is the Hohenberg–Kohn theorem [143], later developed into the Kohn–Sham theory (KS-DFT) [144]. Firstly, the energy of the system for a non-degenerate stationary state is uniquely determined by its electron density, which depends on three spatial coordinates. For this reason, the energy is expressed as a functional of the (scalar) electron density function. Secondly, according to the H-K theorems, the minimum energy occurs for a unique, precise electron density in the ground state. KS-DFT includes the Coulomb interactions between electrons and considers the energy of the exchange and correlation interactions. For a long time, the dispersion energy (the energy of the long-ranged electron correlation) [145] represented a difficult problem, as it is a time-dependent phenomenon, and it was not included in KS-DFT. However, nowadays, dispersion corrections are available and can be included in DFT functions [146]. Therefore, application of DFT leads to the highest obtainable accuracy in calculations. The only existing drawbacks are the risk of underestimating the energy [147] and the time needed to acquire the results. DFT-based calculations are especially widely used in solid-state studies.
Although application of DFT is mostly concentrated on the investigation of single molecules (see the paragraphs below), it is also used to determine total binding energies between a ligand and protein. This is the case for systems composed of ER, SERMs, and two widely used ER antagonists: 4-hydroxy tamoxifen (4OH-T) and raloxifene (RAL) [148]. The results show that the 4OH-T-SERMs set binds more strongly to the ER than the RAL-SERMs set. This is fully in agreement with the experimental data and, once more, as many other studies, confirms the high accuracy of the DFT calculations.
  • Crystal structure prediction
DFT is the theoretical basis for periodic calculations performed on solids, often pharmaceuticals, in order to find and depict new polymorphic forms of drugs or potentially bioactive molecules [149]. It can also be a part of the Crystal Structure Prediction (CSP) approach. Such methodology has been recently used to study crystal structures of 17β-estradiol. As a result, an estradiol hemihydrate has been computationally determined [150]. DFT calculations are also often necessary to refine a crystal structure obtained from powder X-ray diffraction (PXRD) experiments. These calculations are mostly consistent with experimental data, as in the case of estradiol ethinyl cocrystals [151].
DFT-based methodology has also been applied to examine the dissolution process in a study of estradiol cocrystals [152] and to calculate the free energy of solvation in the estradiol–ER complex [153].
  • NMR and vibrational properties calculations
One of the most common DFT applications is to predict NMR and vibrational data. For NMR properties, GIAO [154] and GIPAW [155] methods are implemented. Estradiol being relatively complex (many atoms, presence of both aromatic and non-aromatic rings) could cause computational problems. For this reason, it has been part of a many-ligand study to prove the applicability of the GIAO method as well as its ability to calculate JHH, JHC, and JCC NMR coupling constants [156,157]. Generally, the results obtained were satisfying. However, nowadays it is more common to study the solid state with the second method. GIPAW NMR calculations were also done recently for a new polymorph of 17β-estradiol [158]. The calculations helped to improve the assignment accuracy of chemical shifts obtained from the experiment and, therefore, to elucidate the structure of a new anhydrous estradiol form Figure 1.
Another study showed the applicability of DFT-NMR to explain the nature of a more complex system: the transversal distribution of 17β-estradiol in lipid membranes [159]. NOESY 2D NMR spectra contained cross-peaks between the hormone and lipids. Here, too, DFT calculations helped to properly interpret the experimental data and, as a consequence, greatly aided in describing the position of the estradiol aromatic ring in the membrane. An implication of such study is to increase our understanding of estradiol’s transfer through a lipid bilayer [159,160]. A second area where DFT methodology is widely applied is in calculating vibrational properties. The direct usefulness of the method is manifested through its contribution to accurate assignment of the vibrational modes in IR or Raman studies. This has been implemented, for example, in investigations of estrogens [160] and estradiol-17 valerate [161]. In most reported cases, computationally generated spectra were in very good agreement with the experimental data.
From a wider perspective, DFT calculations enable one to properly describe the examined subject. Additionally, in case of estrogens and estrogen derivatives DFT-derived spectra (both vibrational [162,163] and NMR [158]) stay in a good agreement with the experimental data and are often the only way to properly assign bands (Figure 2). In order to obtain theoretical spectra which ideally meet the experimental ones, scaling factors must be implemented [164]. Thanks to such combination of the theoretical and experimental approach, the first full interpretation of estriol IR spectrum could have been published [164]. The analogical situation has been reported for E1, E2, E3, and ethynylestradiol Raman spectra [163]. Simulation was necessary to identify unique marker peaks in the finger-print region what was useful to differentiate between very similar estrogen structures.
Another example shows that the vibrational frequency calculations of estradiol alone and in monohydrated form [165] has given insight into hydrogen bond formation by estradiol’s D-ring. This has been used to discuss the relationship between the stability of hydrated clusters and the estradiol conformation. In another study [166], DFT-based IR and Raman frequencies were used to investigate estradiol and tamoxifen structures. This helped to understand the intermolecular interactions made by these two molecules and to interpret opposite estrogenic effects. However, it should be mentioned that in the studies of estrogens, Raman spectroscopy is used less frequently than IR. This is due to the inherently weaker signals and common presence of fluorescence interference from the contaminants [163]. Apart from IR and Raman, low-frequency vibrations could also be calculated with the help of DFT. A good example is the assignment of vibrational modes in terahertz spectra for testosterone, estradiol, and estrone [167].
  • Removal of estrogenic pollutants
With regard to estrogen removal from the environment, QM calculations are performed not only on nanotubes but also on other sorbents (e.g., lignocellulosic material) [168]. The adsorption energy shows that this adsorbent could be used to remove all three main estrogens, E1, E2, and E3, from a solution by means of liquid phase extraction. Another extended study revealed the applicability of reduced graphene oxide modified with silver nanoparticles in electrochemical detection of estriol [169]. Firstly, MD simulations at 1000 K were performed in order to obtain 100 conformers of estriol. Later, these structures were used in a semi-empirical Hartree–Fock geometry that was pre-optimized with solvent simulated via the conductor-like Screening Model COSMO [170]. Afterwards, the most stable conformers were fully optimized by the DFT-based software. In this step, the solvent was included via application of the Polarizable Continuum Model [171]. Then, on the structure with the lowest energy, MD at 300 K was performed. Electric properties of the newly obtained conformers were determined with DFT calculations. This study is a good example of the wide range of in silico methods that can be applied to molecularly model the investigated subject (i.e., to develop a method to detect estriol in tap water and urine samples).

2.2. Application of Various Molecular Modelling Methods in the Study of Xenoestrogens

2.2.1. Various Molecular Modelling Methods Applied in Xenoestrogen Studies

The same methods as for estradiol, above, are also used to investigate xenoestrogens. They include QM/MM [172] and QSAR approaches [173,174,175] and confirm the importance of the hydrogen bond between xenoestrogen molecules and the His524 residue in the active site. A variety of xenoestrogens, such as bisphenol A and C, butylparaben, 4-octylphenol, DDE, phthalate, zaeranalol, estradiol, 4-OH-tamoxifen [176], and several proteins (ER [177], SHBG [178]), have been studied. These studies either looked for structural similarities between different ligands or focused on one specific molecule, such as zearalenone [69]. The latter study suggested a lack of agonistic activity against the ER due to the lack of any stable, functionally active conformation of the tested molecule in the LBD. On the other hand, according to QSAR analyses performed on zearalenone analogues [179] and metabolites [180], these zearalenone-related compounds show some estrogenicity due to the presence of a keto/hydroxyl group, a trans double bond in the macrolide ring, and two hydroxyl groups in the aromatic ring, which participate in binding to ER.
As endocrine-disrupting chemicals are present in the environment, research into their removal is constantly being performed. For example, to study dioxin adsorption on graphene, DFT [181,182] and MD [183] calculations have been undertaken.

2.2.2. Bisphenol A

  • Bisphenol A (BPA)–ER complex studies
Bisphenol analogues are among the most studied xenoestrogens. Molecular modelling helps show how these closely related structures adopt agonist/antagonist orientations in the estradiol binding pocket [184] and delineate the binding modes of each bisphenol molecule [185]. A deeper study, concentrating not only on the pure ligand–receptor binding but also including the allosteric effects and application of MD calculations, showed that BPA causes changes in a full-size receptor, and its effect is not limited to the separate domains [186].
A separate set of studies investigated BPAs with halogen substituents on the phenolic rings. All studies showed that a hydrogen bond with His524 for an agonist and with Thr347 for an antagonist was created, exactly as in the case of estradiol and 4-OH-tamoxifen, respectively [187,188]. What is more, just as for estradiol, the stability of helix H12 is crucial in halogenated BPAs–ER complexes [188,189]. The in silico results have been confirmed by experimentally measured affinities.
Although BPA mimics estradiol’s action in the LBD, the QM/MM study revealed that, in comparison to other tested EDCs, it exhibited lower estrogenic activity, probably due to the lack of interaction with His524 [190]. In turn, application of MD helped to elucidate mechanisms driving BPA–ER binding. According to that study, direct hydrogen bonds and hydrophobic interactions are responsible for the binding [173]. Like the previously mentioned experiment, this one confirmed that the ER binding affinity is slightly lower for BPA than for estradiol. MD not only helps to elucidate bound conformations and binding energies between LBD and BPA in the ER [191], but it also gives insight into the influence of that binding on the whole protein, including the DBD. One of the studies [185] showed that the allosteric effects in the LBD due to BPA binding could cause relaxation of the DBD and, therefore, alter ER’s function. Other researchers reported the influence of bisphenol compounds on the protein’s allosteric modulation, altering the Helix12 stability and reducing the recruitment potency of co-activators [187]. This knowledge can be useful in the process of estimating the toxicity of compounds.
  • Risk assessment and removal attempts
In 2012, a protocol for in silico risk assessment of BPA on the ER was proposed [189]. Later, many studies dealing with BPA removal from water were performed [192,193]. Applied molecular modelling techniques include DFT and MD. Since BPA results from the depolymerization of, for example, polycarbonates, both BPA and the initial polymer have been studied under periodic boundary conditions (pbc) [194]. The same DFT-pbc approach has been used to model the photocatalytic degradation of BPA caused by cobalt-doped BiOCl nanosheets [195] and by the effect of humidity combined with UV irradiation [196]. Another example is the evaluation of BPA’s binding to microextraction coatings [197]. In this case, dispersion-corrected DFT was applied.

2.2.3. Phthalates

The second group of xenoestrogens most widely examined by molecular modelling are the phthalates [198,199]. In silico investigations showed that estradiol has a lower ER binding affinity than phthalates. The highest RBA is exhibited by monophthalates [200]. In the case of SHBG binding, score values suggest that short-chain phthalates are more potent than long-chain ones [201]. This agrees with known experimental data.
To look for an effective method of phthalate removal from water, as with estradiol, DFT and MD calculations on SWNT–pollutant complexes have been performed [202]. The adsorption energy has been calculated, and the adsorbent’s chemical groups responsible for the binding have been determined. MD has also been used to examine polymer–solvent interactions while looking for a new, more eco-friendly substitute for plasticizers [203,204].

2.2.4. Technical Aspects of Calculations Performed on (xeno)Estrogens

The computational method most commonly applied in the analysis of (xeno)estrogens is molecular docking. Available publications show that, for this purpose, the most common software packages are Maestro Schrödinger and AutoDockTools. The former is also widely applied in Virtual Screening Workflow (Table 3, N° 9–11), which includes high-throughput virtual screening (HTVS) [205] and molecular docking with either standard or extra precision (SP, XP). In both cases, the OPLS 2005 force field is used. The mentioned (N° 9) consensus score is an effective score that enables ranking of the investigated ligands. It is a combination of different scores, like DockingScore (GlideScore + state penalties for protonation) [206], MM/GBSA Score (binding free energy calculations based on the MD trajectories) [207], and QSAR Score.
Virtual Screening Workflow makes it possible to screen large sets of ligands. It helps to differentiate between ‘actives’ (compounds active against the target protein), ‘decoys’ (compounds of known non-activity against the target protein), and ‘inhibitors/activators’ (potential bio-active substances). As a consequence, the Virtual Screening Workflow approach guides future ligand synthesis and helps in setting a priority for in vitro testing.
Often, a computational step following molecular docking is MD. Here, the most applied codes are GROMACS and AMBER using CHARMM and AMBER force fields, respectively. In most cases, the TIP3 (transferable intermolecular potential with three points) solvent model is applied [242].
Computational methods applied in (xeno)estrogens studies that deal not with a solvent environment but with a solid state are Crystal Structure Prediction (CSP) (N° 26, 27) [235,236,237] and DFT-based calculation of spectroscopic (IR, Raman, NMR) properties (N° 28–34). For the latter, the most commonly applied codes are CASTEP and Gaussian with GGA PBE [239] and B3LYP [241] functionals, respectively. These two functionals seem to be the most reasonable for the investigated subjects. GGA PBE establishes the non-homogeneity in the electron density and leads to more precise results, which in turn is of high importance in NMR spectra calculations. B3LYP is a hybrid functional based on combining DFT and Hartree–Fock theories and finds its application in spectroscopic spectra generation.
On the contrary, CSP [235,236,237] is a multi-step and much more complicated methodology, as it implements both MM and QM and sometimes even MD. Firstly, MM calculations are performed to generate and rank possible compound conformations. Afterwards, the selected conformers are subjected either to ab initio calculations on a molecule or to DFT-D (dispersion corrected) [146] calculations performed on the whole crystal structure. This enables one to observe conformational polymorphisms in the first case and packing polymorphisms in the second. The lattice energies obtained could be adjusted if kinetic factors (like temperature) are included. For that purpose, time- and computational power-consuming MD must be applied.
In contrast to molecular docking or MD simulations, one of the most important calculation methodologies for (xeno)estrogens, QSAR, is independent of the protein and based solely on the ligand structure. This explains why different codes must be applied for QSAR. Their application to (xeno)estrogens has already been described in detail, and the available codes have been compared in [14].
The above-mentioned codes and parameters cover the most common calculations. However, it is impossible to point out the best ones due to the insufficient number of studies. We can only observe that, out of the gathered data, some standard methodologies emerge.
However, a comparison of the applied methodologies in terms of their usage as well as their advantages and drawbacks is possible. The most important aspects have been gathered in Table 4.
To predict the binding affinities or the interactions between the (xeno)estrogens and biomacromolecules, either simple molecular docking or more sophisticated methods such as QM/MM, MD/MM, or FEP can be used. Notably, the more sophisticated methods require not only more specialized software but more computational time and power. Since, to the best of our knowledge, no study has been reported comparing the accuracy of various ligand docking methods applied to the particular group of (xeno)estrogens, no specific indications can be provided.
When focusing on the structural and physicochemical properties of estrogens and xenoestrogens, DFT-based methods have proven their high accuracy and reasonable calculation time. Therefore, such computations can be performed to obtain structural, spectroscopic (IR, Raman, NMR), and thermodynamic data of estrogens, xenoestrogens, their complexes, and solid-state forms such as solvates, salts, and co-crystals. Standard DFT functions (B3LYP for isolated compounds and PBE for periodic structures) have been found to be accurate in multiple studies.

3. Conclusions

In this review, it was clearly shown that molecular modelling methods are valuable tools in studies on estrogens and xenoestrogens. Their relatively low cost, requiring only certain specialized software licenses and computing servers, their increased personal and environmental safety, and their reasonable accuracy make molecular modelling methods unique and modern tools for these studies. In this article, the most common biomolecules studied using molecular modelling were presented, together with appropriate references to the published results. This group of molecules is composed mostly of enzymes participating in the metabolism of estrogens, along with estrogen receptors and even specific nucleic acid domains. While most studies focused on predicting the affinities of small molecules (ligands) to the chosen receptors, the computational research is not limited to this aspect. Another important role for modelling is to explain the conformational changes resulting from binding. Such in silico studies would not be possible without the very large number of already deposited, high-quality crystal structures of estrogen-related proteins that can be easily accessed and used in molecular modelling studies. An overview of those structures was presented in this review. The oldest studies in which molecular modelling was used to study the biochemistry of estrogens focused on molecular docking with molecular mechanics. More recent studies have used more sophisticated methods such as molecular dynamics or combinations of quantum mechanics and molecular mechanics. Further, in this review, it was shown that computational studies concern not only interactions between estrogens and biomacromolecules, but they also can be used to describe phenomena such as migration of estrogens through lipid bilayers or their adsorption on various materials. This can help predict the most efficient way to remove them from the environment when treated as pollutants. Further, it has been shown by multiple examples that quantum molecular modelling methods, such as those based on density functional theory, can be successfully used in structural studies on new solid forms of estrogens such as salts, co-crystals, hydrates, and polymorphs as well as on the complexes of estrogens with other molecules (i.e., cyclodextrins). In addition, the possibility to accurately calculate vibrational and NMR properties can be very helpful to explain spectroscopic results. Finally, we presented similar molecular modelling studies on xenoestrogens such as Bisphenol A and phthalates. Therefore, taking into consideration the versatility and confirmed accuracy of molecular modelling methods, it is not surprising that they have been listed in the specific guidance for studies on EDC published by international organizations such as ECA, EFSA, and OECD.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/17/6411/s1.

Author Contributions

Conceptualization, A.H.M. and Ł.S.; investigation, A.H.M. and Ł.S.; writing—original draft preparation, A.H.M. and T.S.; writing—review and editing, A.H.M., Ł.S., T.S., and D.M.P.; visualization, A.H.M., Ł.S., and D.M.P.; supervision, A.H.M. and Ł.S.; project administration, Ł.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

17β-HSD17β-hydroxysteroid dehydrogenase
ADMETAdsorption distribution metabolism elimination toxicity
BPABisphenol A
CoMFAComparative molecular field analysis
DBDDNA binding domain
DFTDensity functional theory
E1Estrone
E2Estradiol
E3Estriol
E4Estretrol
EREstrogen receptor
EDCsEndocrine-disrupting chemicals
FEPFree energy perturbation
H-KHohenberg–Kohn theorems
K-SKohn–Sham theorems
LBDLigand-binding domain
LIELinear interaction energy
MDMolecular dynamics
MMMolecular mechanics
PDBProtein Data Bank
QMQuantum mechanics
QSARQuantitative structure–activity relationship
RBARelative binding affinity
SERMsSelective estrogen receptor modulators
SHBGSex hormone-binding globulin
STSSulfatase
SULTSulfotransferase

References

  1. Bennink, H.C.; Verhoeven, C.; Zimmerman, Y.; Visser, M.; Foidart, J.-M.; Gemzell-Danielsson, K. Clinical effects of the fetal estrogen estetrol in a multiple-rising-dose study in postmenopausal women. Matur. Eur. Menopause J. 2016, 91, 93–100. [Google Scholar] [CrossRef] [Green Version]
  2. Singleton, D.W.; David, W.S. Xenoestrogen exposure and mechanisms of endocrine disruption. Front. Biosci. 2003, 8, 110–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bennink, F.C.; Holinka, C.F.; Visser, M.; Bennink, H.J.T.C. Maternal and fetal estetrol levels during pregnancy. Climacteric 2008, 11, 69–72. [Google Scholar] [CrossRef]
  4. Schreiner, W.E. The Ovary in Labhart, A. Clinical Endocrinology: Theory and Practice; Springer: Berlin/Heidelberg, Germany, 2012; p. 548. [Google Scholar]
  5. Kuhl, H. Pharmacology of estrogens and progestogens: Influence of different routes of administration. Climacteric 2005, 8 (Suppl. 1), 3–63. [Google Scholar] [CrossRef]
  6. Blackburn, S. Maternal, Fetal, & Neonatal Physiology; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar]
  7. Fait, T. Menopause hormone therapy: Latest developments and clinical practice. Drugs Context 2019, 8, 1–9. [Google Scholar] [CrossRef]
  8. Tofovic, S.P.; Jackson, E.K. Estradiol Metabolism: Crossroads in Pulmonary Arterial Hypertension. Int. J. Mol. Sci. 2019, 21, 116. [Google Scholar] [CrossRef] [Green Version]
  9. EC. Corrigendum to Commission Regulation (EU) 2018/605 of 19 April 2018 Amending Annex II to Regulation (EC) No 1107/2009 by Setting Out Scientific Criteria for the Determination of Endocrine Disrupting Properties. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32018R0605&rid=210 (accessed on 25 August 2020).
  10. Ciubotaru, R.M.; Oyedele, J.; Zancanaro, G. Annual assessment of Echinococcus multilocularis surveillance reports submitted in 2018 in the context of Commission Regulation (EU) No 1152/2011. EFSA J. 2018, 16, 33–36. [Google Scholar] [CrossRef] [Green Version]
  11. OECD. Guidance Document on Standardised Test Guidelines for Evaluating Chemicals for Endocrine Disruption, No. 150, Update v3, Series on Testing and Assessment, ENV/JM/MONO(2012)22; OECD: Paris, France, 2017; 988p. [Google Scholar]
  12. OECD. Conceptual Framework for Testing and Assessment of Endocrine Disrupters. Available online: https://www.oecd.org/env/ehs/testing/oecdworkrelatedtoendocrinedisrupters.htm (accessed on 25 August 2020).
  13. Safe, S.; Khan, S.; Wu, F.; Li, X. Chemical-induced estrogenicity in Xenoestrogens and phytoestrogens as SERMs and implications for risk assessment (Chapter 65). In Veterinary Toxicology Basic and Clinical Principles; Elsevier: Amsterdam, The Netherlands, 2007; pp. 811–822. [Google Scholar]
  14. Cotterill, J.; Palazzolo, L.; Ridgway, C.; Price, N.; Rorije, E.; Moretto, A.; Peijnenburg, A.; Eberini, I. Predicting estrogen receptor binding of chemicals using a suite of in silico methods—Complementary approaches of (Q)SAR, molecular docking and molecular dynamics. Toxicol. Appl. Pharmacol. 2019, 378, 114630. [Google Scholar] [CrossRef]
  15. Kerdivel, G.; Habauzit, D.; Pakdel, F. Assessment and Molecular Actions of Endocrine-Disrupting Chemicals That Interfere with Estrogen Receptor Pathways. Int. J. Endocrinol. 2013, 2013, 501851. [Google Scholar] [CrossRef] [Green Version]
  16. Park, J.; Czapla, L.; Amaro, R.E. Molecular Simulations of Aromatase Reveal New Insights into the Mechanism of Ligand Binding. J. Chem. Inf. Model. 2013, 53, 2047–2056. [Google Scholar] [CrossRef] [Green Version]
  17. Lokhande, K.; Mehere, A.S.; Yadav, A.K.; Swamy, K.V. Molecular Modeling and Docking Studies of Aromatase Inhibitors with Aromatase for ERP Breast Cancer. In Proceedings of the 86th Conference of Society of Biological Chemists: Emerging Discoveries in Health and Agricultural Sciences, New Delhi, India, 16–19 November 2017. [Google Scholar] [CrossRef]
  18. Kang, H.; Xiao, X.; Huang, C.; Yuan, Y.; Tang, D.; Dai, X.; Zeng, X. Potent aromatase inhibitors and molecular mechanism of inhibitory action. Eur. J. Med. Chem. 2018, 143, 426–437. [Google Scholar] [CrossRef] [PubMed]
  19. Awasthi, M.; Singh, S.; Pandey, V.P.; Dwivedi, U.N. Molecular docking and 3D-QSAR-based virtual screening of flavonoids as potential aromatase inhibitors against estrogen-dependent breast cancer. J. Biomol. Struct. Dyn. 2014, 33, 804–819. [Google Scholar] [CrossRef] [PubMed]
  20. Suvannang, N.; Nantasenamat, C.; Isarankura-Na-Ayudhya, C.; Prachayasittikul, V. Molecular Docking of Aromatase Inhibitors. Molecules 2011, 16, 3597–3617. [Google Scholar] [CrossRef] [Green Version]
  21. Rampogu, S.; Son, M.; Park, C.; Kim, H.-H.; Suh, J.-K.; Lee, K. Sulfonanilide Derivatives in Identifying Novel Aromatase Inhibitors by Applying Docking, Virtual Screening, and MD Simulations Studies. BioMed Res. Int. 2017, 2017, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Narayana, B.L.; Kishore, D.P.; Balakumar, C.; Rao, K.V.; Kaur, R.; Rao, A.R.; Murthy, J.N.; Ravikumar, M. Molecular Modeling Evaluation of Non-Steroidal Aromatase Inhibitors†. Chem. Boil. Drug Des. 2012, 79, 674–682. [Google Scholar] [CrossRef] [PubMed]
  23. Fischer, J.; Ganellin, C.R. Analogue-Based Drug Discovery; John Wiley & Sons: Hoboken, NJ, USA, 2006; p. 516. [Google Scholar]
  24. Maltais, R.; Ayan, D.; Trottier, A.; Barbeau, X.; Lagüe, P.; Bouchard, J.-E.; Poirier, D. Discovery of a Non-Estrogenic Irreversible Inhibitor of 17β-Hydroxysteroid Dehydrogenase Type 1 from 3-Substituted-16β-(m-carbamoylbenzyl)-estradiol Derivatives. J. Med. Chem. 2013, 57, 204–222. [Google Scholar] [CrossRef]
  25. Lespérance, M.; Barbeau, X.; Roy, J.; Maltais, R.; Lagüe, P.; Poirier, D. Chemical synthesis of C3-oxiranyl/oxiranylmethyl-estrane derivatives targeted by molecular modeling and tested as potential inhibitors of 17β-hydroxysteroid dehydrogenase type 1. Steroids 2018, 140, 104–113. [Google Scholar] [CrossRef]
  26. Maltais, R.; Djiemeny, A.N.; Roy, J.; Barbeau, X.; Lambert, J.-P.; Poirier, D. Design and synthesis of dansyl-labeled inhibitors of steroid sulfatase for optical imaging. Bioorg. Med. Chem. 2020, 28, 115368. [Google Scholar] [CrossRef]
  27. Rakers, C.; Schumacher, F.; Meinl, W.; Glatt, H.; Kleuser, B.; Wolber, G. In SilicoPrediction of Human Sulfotransferase 1E1 Activity Guided by Pharmacophores from Molecular Dynamics Simulations. J. Boil. Chem. 2015, 291, 58–71. [Google Scholar] [CrossRef] [Green Version]
  28. Kisselev, P.; Schunck, W.-H.; Roots, I.; Schwarz, D. Association of CYP1A1 Polymorphisms with Differential Metabolic Activation of 17 -Estradiol and Estrone. Cancer Res. 2005, 65, 2972–2978. [Google Scholar] [CrossRef] [Green Version]
  29. Hong, C.-C.; Tang, B.-K.; Hammond, G.L.; Tritchler, D.; Yaffe, M.; Boyd, N.F. Cytochrome P450 1A2 (CYP1A2) activity and risk factors for breast cancer: A cross-sectional study. Breast Cancer Res. 2004, 6, R352–R365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Nishida, C.R.; Everett, S.; De Montellano, P.R.O. Specificity Determinants of CYP1B1 Estradiol Hydroxylation. Mol. Pharmacol. 2013, 84, 451–458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Selby, C. Sex Hormone Binding Globulin: Origin, Function and Clinical Significance. Ann. Clin. Biochem. Int. J. Lab. Med. 1990, 27, 532–541. [Google Scholar] [CrossRef]
  32. Dechering, K.; Boersma, C.; Mosselman, S. Estrogen receptors alpha and beta: Two receptors of a kind? Curr. Med. Chem. 2000, 7, 561–576. [Google Scholar] [CrossRef] [PubMed]
  33. Paterni, I.; Granchi, C.; Katzenellenbogen, J.A.; Minutolo, F. Estrogen receptors alpha (ERα) and beta (ERβ): Subtype-selective ligands and clinical potential. Steroids 2014, 90, 13–29. [Google Scholar] [CrossRef] [Green Version]
  34. Ruff, M.; Gangloff, M.; Wurtz, J.; Moras, D. Estrogen receptor transcription and transactivation Structure-function relationship in DNA- and ligand-binding domains of estrogen receptors. Breast Cancer Res. 2000, 2, 353–359. [Google Scholar] [CrossRef] [Green Version]
  35. Kumar, R.; Zakharov, M.N.; Khan, S.H.; Miki, R.; Jang, H.; Toraldo, G.; Singh, R.; Bhasin, S.; Jasuja, R. The Dynamic Structure of the Estrogen Receptor. J. Amino Acids 2011, 2011, 1–7. [Google Scholar] [CrossRef] [Green Version]
  36. Billon-Galés, A.; Krust, A.; Fontaine, C.; Abot, A.; Flouriot, G.; Toutain, C.; Berges, H.; Gadeau, A.-P.; Lenfant, F.; Gourdy, P.; et al. Activation function 2 (AF2) of estrogen receptor-α is required for the atheroprotective action of estradiol but not to accelerate endothelial healing. Proc. Natl. Acad. Sci. USA 2011, 108, 13311–13316. [Google Scholar] [CrossRef] [Green Version]
  37. Gao, L.; Tu, Y.; Ågren, H.; Eriksson, L.A. Characterization of Agonist Binding to His524 in the Estrogen Receptor α Ligand Binding Domain. J. Phys. Chem. B 2012, 116, 4823–4830. [Google Scholar] [CrossRef]
  38. Hu, G.; Wang, J. Ligand selectivity of estrogen receptors by a molecular dynamics study. Eur. J. Med. Chem. 2014, 74, 726–735. [Google Scholar] [CrossRef]
  39. Gu, X. Helix 12 in the human estrogen receptor (hER) is essential for the hER function by overcoming nucleosome repression in yeast. J. Cell. Biochem. 2002, 86, 224–238. [Google Scholar] [CrossRef] [PubMed]
  40. Shiau, A.K.; Barstad, D.; Loria, P.M.; Cheng, L.; Kushner, P.J.; Agard, D.A.; Greene, G.L. The Structural Basis of Estrogen Receptor/Coactivator Recognition and the Antagonism of This Interaction by Tamoxifen. Cell 1998, 95, 927–937. [Google Scholar] [CrossRef] [Green Version]
  41. Fuentes, N.; Silveyra, P. Estrogen receptor signaling mechanisms. Adv. Protein Chem. Struct. Biol. 2019, 116, 135–170. [Google Scholar] [CrossRef] [PubMed]
  42. Stefkovich, M.L.; Arao, Y.; Hamilton, K.J.; Korach, K.S. Experimental models for evaluating non-genomic estrogen signaling. Steroids 2018, 133, 34–37. [Google Scholar] [CrossRef] [PubMed]
  43. Vrtacnik, P.; Ostanek, B.; Mencej-Bedrac, S.; Marc, J. The many faces of estrogen signaling. Biochem. Medica 2014, 24, 329–342. [Google Scholar] [CrossRef] [Green Version]
  44. Rosano, C.; Ponassi, M.; Santolla, M.F.; Pisano, A.; Felli, L.; Vivacqua, A.; Maggiolini, M.; Lappano, R. Macromolecular Modelling and Docking Simulations for the Discovery of Selective GPER Ligands. AAPS J. 2015, 18, 41–46. [Google Scholar] [CrossRef] [Green Version]
  45. Bruno, A.; Aiello, F.; Costantino, G.; Radi, M. Homology Modeling, Validation and Dynamics of the G Protein-coupled Estrogen Receptor 1 (GPER-1). Mol. Inform. 2016, 35, 333–339. [Google Scholar] [CrossRef]
  46. Méndez-Luna, D.; Martínez-Archundia, M.; Maroun, R.C.; Ceballos, G.; Fragoso-Vázquez, M.; González-Juárez, D.; Correa-Basurto, J. Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies, molecular dynamics, and docking simulations. J. Biomol. Struct. Dyn. 2015, 33, 1–12. [Google Scholar] [CrossRef]
  47. Taylor, R.; Jewsbury, P.; Essex, J.W. A review of protein-small molecule docking methods. J. Comput. Mol. Des. 2002, 16, 151–166. [Google Scholar] [CrossRef]
  48. RCSB PDB. Available online: https://www.rcsb.org/ (accessed on 25 August 2020).
  49. Landeros-Martinez, L.-L.; Glossman-Mitnik, D.; Orrantia-Borunda, E.; Flores-Holguin, N. A Combined Molecular Docking and Electronic Structure Study for a Breast Cancer Drug Design. In Molecular Docking; IntechOpen: London, UK, 2018. [Google Scholar] [CrossRef] [Green Version]
  50. Muchtaridi, M.; Megantara, S.; Dermawan, D.; Yusuf, M. Antagonistic mechanism of α-mangostin derivatives against human estrogen receptor α of breast cancer using molecular dynamics simulation. Rasayan J. Chem. 2019, 12, 1927–1934. [Google Scholar] [CrossRef]
  51. Cavasotto, C.; Aucar, M.G. High-Throughput Docking Using Quantum Mechanical Scoring. Front. Chem. 2020, 8, 246. [Google Scholar] [CrossRef] [PubMed]
  52. Spiriti, J.; Subramanian, S.R.; Palli, R.; Wu, M.; Zuckerman, D.M. Middle-way flexible docking: Pose prediction using mixed-resolution Monte Carlo in estrogen receptor α. PLoS ONE 2019, 14, e0215694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wolohan, P.; Reichert, D.E. CoMSIA and docking study of rhenium based estrogen receptor ligand analogs. Steroids 2007, 72, 247–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wierbowski, S.D.; Wingert, B.M.; Zheng, J.; Camacho, C.J. Cross-docking benchmark for automated pose and ranking prediction of ligand binding. Protein Sci. 2019, 29, 298–305. [Google Scholar] [CrossRef] [Green Version]
  55. Shtaiwi, A.; Adnan, R.; Khairuddean, M.; Khan, S.U. Computational investigations of the binding mechanism of novel benzophenone imine inhibitors for the treatment of breast cancer. RSC Adv. 2019, 9, 35401–35416. [Google Scholar] [CrossRef] [Green Version]
  56. Ling, J.; Oh, D.L.; Chia, A.Y.Y. Molecular Docking Studies of Glycyrrhizic Acid (GA), Glycyrrhetic Acid (GE) and Glabridin (GLA) with Estrogen Receptors (ERs). Biosci. Biotechnol. Res. Asia 2017, 14, 1–11. [Google Scholar] [CrossRef]
  57. Pang, X.; Fu, W.; Wang, J.; Kang, D.; Xu, L.; Zhao, Y.; Liu, A.; Du, G.-H. Identification of Estrogen Receptor α Antagonists from Natural Products via In Vitro and In Silico Approaches. Oxidative Med. Cell. Longev. 2018, 2018, 6040149. [Google Scholar] [CrossRef] [Green Version]
  58. Muchtaridi, M.; Dermawan, D.; Yusuf, M. Molecular Docking, 3D Structure-Based Pharmacophore Modeling, and ADME Prediction of Alpha Mangostin and Its Derivatives against Estrogen Receptor Alpha. J. Young Pharm. 2018, 10, 252–259. [Google Scholar] [CrossRef] [Green Version]
  59. Yu, E.; Xu, Y.; Shi, Y.; Yu, Q.; Liu, J.; Xu, L. Discovery of novel natural compound inhibitors targeting estrogen receptor α by an integrated virtual screening strategy. J. Mol. Model. 2019, 25, 278. [Google Scholar] [CrossRef]
  60. Abdelsamie, A.S.; Salah, M.; Siebenbürger, L.; Hamed, M.M.; Börger, C.; Van Koppen, C.J.; Frotscher, M.; Hartmann, R.W. Development of potential preclinical candidates with promising in vitro ADME profile for the inhibition of type 1 and type 2 17β-Hydroxysteroid dehydrogenases: Design, synthesis, and biological evaluation. Eur. J. Med. Chem. 2019, 178, 93–107. [Google Scholar] [CrossRef]
  61. Hasan, T.N.; B, L.G.; A Masoodi, T.; Shafi, G.; Alshatwi, A.A.; Sivashanmugham, P. Affinity of estrogens for human progesterone receptor A and B monomers and risk of breast cancer: A comparative molecular modeling study. Adv. Appl. Bioinform. Chem. 2011, 4, 29–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Fu, X.; Wang, P.; Zhu, B.T. Characterization of the Estradiol-Binding Site Structure of Human Protein Disulfide Isomerase (PDI). PLoS ONE 2011, 6, e27185. [Google Scholar] [CrossRef] [PubMed]
  63. Da Silva, A.J.; Dos Santos, E.S. Aqueous solution interactions with sex hormone-binding globulin and estradiol: A theoretical investigation. J. Boil. Phys. 2018, 44, 539–556. [Google Scholar] [CrossRef] [PubMed]
  64. Thomas, M.P.; Potter, B.V. The structural biology of oestrogen metabolism. J. Steroid Biochem. Mol. Boil. 2013, 137, 27–49. [Google Scholar] [CrossRef]
  65. Maltais, R.; Fournier, D.; Poirier, D.; Maltais, R. Quantitative Structure-Activity Relationship (QSAR) Study with a Series of 17α-Derivatives of Estradiol: Model for the Development of Reversible Steroid Sulfatase Inhibitors. QSAR Comb. Sci. 2009, 28, 1284–1299. [Google Scholar] [CrossRef]
  66. Granados, S.T.; Castillo, K.; Bravo-Moraga, F.; Sepúlveda, R.V.; Carrasquel-Ursulaez, W.; Rojas, M.; Carmona, E.; Lorenzo-Ceballos, Y.; González-Nilo, F.D.; González, C.; et al. The molecular nature of the 17β-Estradiol binding site in the voltage- and Ca2+-activated K+ (BK) channel β1 subunit. Sci. Rep. 2019, 9, 9965. [Google Scholar] [CrossRef]
  67. McCullough, C.; Neumann, T.S.; Gone, J.R.; He, Z.; Herrild, C.; (Nee Lukesh) Wondergem, J.; Pandey, R.K.; Donaldson, W.A.; Sem, D.S. Probing the human estrogen receptor-α binding requirements for phenolic mono- and di-hydroxyl compounds: A combined synthesis, binding and docking study. Bioorg. Med. Chem. 2013, 22, 303–310. [Google Scholar] [CrossRef] [Green Version]
  68. Kerdivel, G.; Le Guével, R.; Habauzit, D.; Brion, F.; Ait-Aissa, S.; Pakdel, F. Estrogenic Potency of Benzophenone UV Filters in Breast Cancer Cells: Proliferative and Transcriptional Activity Substantiated by Docking Analysis. PLoS ONE 2013, 8, e60567. [Google Scholar] [CrossRef] [Green Version]
  69. Sauvée, C.; Schäfer, A.; Sundén, H.; Ma, J.-N.; Gustavsson, A.-L.; Burstein, E.S.; Olsson, R. The A-CD analogue of 16β,17α-estriol is a potent and highly selective estrogen receptor β agonist. MedChemComm 2013, 4, 1439. [Google Scholar] [CrossRef]
  70. Kucinska, M.; Giron, M.D.; Piotrowska, H.; Lisiak, N.; Granig, W.H.; Lopez-Jaramillo, F.-J.; Salto, R.; Murias, M.; Erker, T. Novel Promising Estrogenic Receptor Modulators: Cytotoxic and Estrogenic Activity of Benzanilides and Dithiobenzanilides. PLoS ONE 2016, 11, e0145615. [Google Scholar] [CrossRef]
  71. Gonzalez, T.L.; Rae, J.M.; Colacino, J.A.; Richardson, R.J. Homology models of mouse and rat estrogen receptor-α ligand-binding domain created by in silico mutagenesis of a human template: Molecular docking with 17β-estradiol, diethylstilbestrol, and paraben analogs. Comput. Toxicol. 2019, 10, 1–16. [Google Scholar] [CrossRef] [PubMed]
  72. Yarger, J.G.; E Babine, R.; Bittner, M.; Shanle, E.; Xu, W.; Hershberger, P.; Nye, S.H. Structurally similar estradiol analogs uniquely alter the regulation of intracellular signaling pathways. J. Mol. Endocrinol. 2012, 50, 43–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Maltais, R.; Trottier, A.; Barbeau, X.; Lagüe, P.; Perreault, M.; Thériault, J.-F.; Lin, S.X.; Poirier, D. Impact of structural modifications at positions 13, 16 and 17 of 16β-(m-carbamoylbenzyl)-estradiol on 17β-hydroxysteroid dehydrogenase type 1 inhibition and estrogenic activity. J. Steroid Biochem. Mol. Boil. 2016, 161, 24–35. [Google Scholar] [CrossRef] [PubMed]
  74. E Starkey, N.J.; Li, Y.; Drenkhahn-Weinaug, S.K.; Liu, J.; Lubahn, D.B. 27-Hydroxycholesterol Is an Estrogen Receptor β-Selective Negative Allosteric Modifier of 17β-Estradiol Binding. Endocrinology 2018, 159, 1972–1981. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Singh, K.; Munuganti, R.S.N.; Leblanc, E.; Lin, Y.L.; Leung, E.; Lallous, N.; Butler, M.S.; Cherkasov, A.; Rennie, P.S. In silico discovery and validation of potent small-molecule inhibitors targeting the activation function 2 site of human oestrogen receptor α. Breast Cancer Res. 2015, 17, 27. [Google Scholar] [CrossRef] [Green Version]
  76. Durrant, J.D.; Carlson, K.E.; Martin, T.A.; Offutt, T.L.; Mayne, C.G.; Katzenellenbogen, J.A.; Amaro, R.E. Neural-Network Scoring Functions Identify Structurally Novel Estrogen-Receptor Ligands. J. Chem. Inf. Model. 2015, 55, 1953–1961. [Google Scholar] [CrossRef] [Green Version]
  77. Ng, H.W.; Zhang, W.; Shu, M.; Luo, H.; Ge, W.; Perkins, R.; Tong, W.; Hong, H. Competitive molecular docking approach for predicting estrogen receptor subtype α agonists and antagonists. BMC Bioinform. 2014, 15, S4. [Google Scholar] [CrossRef] [Green Version]
  78. Amr, A.E.-G.E.; Elsayed, E.; Al-Omar, M.; Eldin, H.O.B.; Nossier, E.S.; Abdalla, M.M. Design, Synthesis, Anticancer Evaluation and Molecular Modeling of Novel Estrogen Derivatives. Molecules 2019, 24, 416. [Google Scholar] [CrossRef] [Green Version]
  79. Stjernschantz, E.; Reinen, J.; Meinl, W.; George, B.J.; Glatt, H.; Vermeulen, N.P.; Oostenbrink, C. Comparison of murine and human estrogen sulfotransferase inhibition in vitro and in silico—Implications for differences in activity, subunit dimerization and substrate inhibition. Mol. Cell. Endocrinol. 2010, 317, 127–140. [Google Scholar] [CrossRef] [Green Version]
  80. Hanson, R.N.; McCaskill, E.; Tongcharoensirikul, P.; Dilis, R.; Labaree, D.; Hochberg, R.B. Synthesis and evaluation of 17α-(dimethylphenyl)vinyl estradiols as probes of the estrogen receptor-α ligand binding domain. Steroids 2012, 77, 471–476. [Google Scholar] [CrossRef] [Green Version]
  81. Grande, F.; Rizzuti, B.; Occhiuzzi, M.A.; Ioele, G.; Casacchia, T.; Gelmini, F.; Guzzi, R.; Garofalo, A.; Statti, G. Identification by Molecular Docking ofHomoisoflavones from Leopoldia comosa as Ligands of Estrogen Receptors. Molecules 2018, 23, 894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Powers, C.N.; Setzer, W.N. A molecular docking study of phytochemical estrogen mimics from dietary herbal supplements. Silico Pharmacol. 2015, 3, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Yuan, P.; Liang, K.; Ma, B.; Zheng, N.; Nussinov, R.; Huang, J. Multiple-Targeting and Conformational Selection in the Estrogen Receptor: Computation and Experiment. Chem. Boil. Drug Des. 2011, 78, 137–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Chen, Y.; Wang, J.; Hong, D.-Y.; Chen, L.; Zhang, Y.-Y.; Xu, Y.-N.; Pan, D.; Fu, L.-Y.; Tao, L.; Luo, H.; et al. Baicalein has protective effects on the 17β-estradiol-induced transformation of breast epithelial cells. Oncotarget 2017, 8, 10470–10484. [Google Scholar] [CrossRef] [PubMed]
  85. Fokialakis, N.; Alexi, X.; Aligiannis, N.; Boulaka, A.; Meligova, A.K.; Lambrinidis, G.; Kalpoutzakis, E.; Pratsinis, H.; Cheilari, A.; Mitsiou, D.J.; et al. Biological evaluation of isoflavonoids from Genista halacsyi using estrogen-target cells: Activities of glucosides compared to aglycones. PLoS ONE 2019, 14, e0210247. [Google Scholar] [CrossRef]
  86. Ayoub, N.M.; Siddique, A.B.; Ebrahim, H.Y.; Mohyeldin, M.M.; El Sayed, K.A. The olive oil phenolic (-)-oleocanthal modulates estrogen receptor expression in luminal breast cancer in vitro and in vivo and synergizes with tamoxifen treatment. Eur. J. Pharmacol. 2017, 810, 100–111. [Google Scholar] [CrossRef]
  87. Yuseran, H.; Hartoyo, E.; Nurseta, T.; Kalim, H. Molecular docking of genistein on estrogen receptors, promoter region of BCLX, caspase-3, Ki-67, cyclin D1, and telomere activity. J. Taibah Univ. Med. Sci. 2018, 14, 79–87. [Google Scholar] [CrossRef]
  88. Nanashima, N.; Horie, K.; Maeda, H. Phytoestrogenic Activity of Blackcurrant Anthocyanins Is Partially Mediated through Estrogen Receptor Beta. Molecules 2017, 23, 74. [Google Scholar] [CrossRef] [Green Version]
  89. Alam, S.; Khan, F. Virtual screening, Docking, ADMET and System Pharmacology studies on Garcinia caged Xanthone derivatives for Anticancer activity. Sci. Rep. 2018, 8, 5524. [Google Scholar] [CrossRef]
  90. Puranik, N.V.; Srivastava, P.; Bhatt, G.; Mary, D.J.S.J.; Limaye, A.M.; Sivaraman, J. Determination and analysis of agonist and antagonist potential of naturally occurring flavonoids for estrogen receptor (ERα) by various parameters and molecular modelling approach. Sci. Rep. 2019, 9, 7450. [Google Scholar] [CrossRef]
  91. Wang, T.; Wang, Y.; Zhuang, X.; Luan, F.; Zhao, C.; Cordeiro, M.N.D.S. Interaction of Coumarin Phytoestrogens with ERα and ERβ: A Molecular Dynamics Simulation Study. Molecules 2020, 25, 1165. [Google Scholar] [CrossRef] [Green Version]
  92. Wang, X.; Wang, G.-C.; Rong, J.; Wang, S.W.; Ng, T.B.; Zhang, Y.B.; Lee, K.F.; Zheng, L.; Wong, H.-K.; Yung, K.K.-L.; et al. Identification of Steroidogenic Components Derived From Gardenia jasminoides Ellis Potentially Useful for Treating Postmenopausal Syndrome. Front. Pharmacol. 2018, 9, 9. [Google Scholar] [CrossRef]
  93. Jeong, J.; Kim, H.; Choi, J. In Silico Molecular Docking and In Vivo Validation with Caenorhabditis elegans to Discover Molecular Initiating Events in Adverse Outcome Pathway Framework: Case Study on Endocrine-Disrupting Chemicals with Estrogen and Androgen Receptors. Int. J. Mol. Sci. 2019, 20, 1209. [Google Scholar] [CrossRef] [Green Version]
  94. Neves, B.J.; Braga, R.C.; Melo-Filho, C.C.; Moreira-Filho, J.T.; Muratov, E.N.; Andrade, C.H. QSAR-Based Virtual Screening: Advances and Applications in Drug Discovery. Front. Pharmacol. 2018, 9, 1275. [Google Scholar] [CrossRef] [Green Version]
  95. Dror, O.; Schneidman-Duhovny, D.; Inbar, Y.; Nussinov, R.; Wolfson, H.J. Novel Approach for Efficient Pharmacophore-Based Virtual Screening: Method and Applications. J. Chem. Inf. Model. 2009, 49, 2333–2343. [Google Scholar] [CrossRef] [Green Version]
  96. Fiser, A. Template-based protein structure modeling. Breast Cancer 2010, 673, 73–94. [Google Scholar] [CrossRef] [Green Version]
  97. Zhao, X.; Chen, M.; Huang, B.; Ji, H.; Yuan, M. Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) Studies on α1A-Adrenergic Receptor Antagonists Based on Pharmacophore Molecular Alignment. Int. J. Mol. Sci. 2011, 12, 7022–7037. [Google Scholar] [CrossRef]
  98. Sharma, R.; Dhingra, N.; Patil, S. CoMFA, CoMSIA, HQSAR and Molecular Docking Analysis of Ionone-based Chalcone Derivatives as Antiprostate Cancer Activity. Indian J. Pharm. Sci. 2016, 78, 54–64. [Google Scholar] [CrossRef] [Green Version]
  99. Vuorinen, A.; Engeli, R.T.; Leugger, S.; Kreutz, C.R.; Schuster, D.; Odermatt, A.; Matuszczak, B. Phenylbenzenesulfonates and -sulfonamides as 17β-hydroxysteroid dehydrogenase type 2 inhibitors: Synthesis and SAR-analysis. Bioorg. Med. Chem. Lett. 2017, 27, 2982–2985. [Google Scholar] [CrossRef]
  100. Vuorinen, A.; Engeli, R.; Meyer, A.; Bachmann, F.; Griesser, U.J.; Schuster, D.; Odermatt, A. Ligand-Based Pharmacophore Modeling and Virtual Screening for the Discovery of Novel 17β-Hydroxysteroid Dehydrogenase 2 Inhibitors. J. Med. Chem. 2014, 57, 5995–6007. [Google Scholar] [CrossRef]
  101. Chang, Y.-H.; Chen, J.-Y.; Hor, C.-Y.; Chuang, Y.-C.; Yang, C.-B.; Yang, C.-N. Computational Study of Estrogen Receptor-Alpha Antagonist with Three-Dimensional Quantitative Structure-Activity Relationship, Support Vector Regression, and Linear Regression Methods. Int. J. Med. Chem. 2013, 2013, 1–13. [Google Scholar] [CrossRef] [Green Version]
  102. Sodero, A.C.R.; Romeiro, N.C.; Da Cunha, E.F.F.; Magalhães, U.D.O.; De Alencastro, R.B.; Rodrigues, C.R.; Cabral, L.M.; Castro, H.C.; Albuquerque, M.G. Application of 4D-QSAR Studies to a Series of Raloxifene Analogs and Design of Potential Selective Estrogen Receptor Modulators. Molecules 2012, 17, 7415–7439. [Google Scholar] [CrossRef]
  103. Wang, P.; McInnes, C.; Zhu, B.T. Structural Characterization of the Binding Interactions of Various Endogenous Estrogen Metabolites with Human Estrogen Receptor α and β Subtypes: A Molecular Modeling Study. PLoS ONE 2013, 8, e74615. [Google Scholar] [CrossRef] [Green Version]
  104. Bhhatarai, B.; Wilson, D.M.; Price, P.S.; Marty, S.; Parks, A.K.; Carney, E. Evaluation of OASIS QSAR Models Using ToxCast™ in Vitro Estrogen and Androgen Receptor Binding Data and Application in an Integrated Endocrine Screening Approach. Environ. Health Perspect. 2016, 124, 1453–1461. [Google Scholar] [CrossRef] [Green Version]
  105. Bohari, M.; Srivastava, H.K.; Sastry, G.N. Analogue-based approaches in anti-cancer compound modelling: The relevance of QSAR models. Org. Med. Chem. Lett. 2011, 1, 3. [Google Scholar] [CrossRef] [Green Version]
  106. Dems, M.A.E.; Laib, S.; Latelli, N.; Ouddai, N. A DFT-based Quantitative structure activity relationship Study of organometallic estradiol derivatives. JCPS 2017, 10, 483–487. [Google Scholar]
  107. Zhang, T.; Wei, D.-Q.; Chou, K.-C. A pharmacophore model specific to active site of CYP1A2 with a novel molecular modeling explorer and CoMFA. Med. Chem. 2012, 8, 198–207. [Google Scholar]
  108. Poirier, D.; Roy, J.; Cortes-Benitez, F.; Dutour, R. Targeting cytochrome P450 (CYP) 1B1 with steroid derivatives. Bioorg. Med. Chem. Lett. 2016, 26, 5272–5276. [Google Scholar] [CrossRef]
  109. Kar, S.; Roy, K.; Leszczynski, J. Impact of Pharmaceuticals on the Environment: Risk Assessment Using QSAR Modeling Approach. Methods Mol. Biol. 2018, 1800, 395–443. [Google Scholar] [CrossRef]
  110. Colosi, L.M.; Huang, Q.; Weber, W.J. QSAR-assisted design of an environmental catalyst for enhanced estrogen remediation. Chemosphere 2010, 81, 897–903. [Google Scholar] [CrossRef] [Green Version]
  111. Rokhina, E.V.; Suri, R. Application of density functional theory (DFT) to study the properties and degradation of natural estrogen hormones with chemical oxidizers. Sci. Total Environ. 2012, 417, 280–290. [Google Scholar] [CrossRef] [PubMed]
  112. Warshel, A.; Levitt, M. Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J. Mol. Boil. 1976, 103, 227–249. [Google Scholar] [CrossRef]
  113. Jorgensen, W. The Many Roles of Computation in Drug Discovery. Science 2004, 303, 1813–1818. [Google Scholar] [CrossRef] [PubMed]
  114. Simonson, T.; Archontis, G.; Karplus, M. Free Energy Simulations Come of Age: Protein−Ligand Recognition. Acc. Chem. Res. 2002, 35, 430–437. [Google Scholar] [CrossRef] [Green Version]
  115. Fratev, F.; Steinbrecher, T.; Jónsdóttir, S.O. Prediction of Accurate Binding Modes Using Combination of Classical and Accelerated Molecular Dynamics and Free-Energy Perturbation Calculations: An Application to Toxicity Studies. ACS Omega 2018, 3, 4357–4371. [Google Scholar] [CrossRef] [Green Version]
  116. Schindler, C.; Rippmann, F.; Kuhn, D. Relative binding affinity prediction of farnesoid X receptor in the D3R Grand Challenge 2 using FEP+. J. Comput. Mol. Des. 2017, 32, 265–272. [Google Scholar] [CrossRef]
  117. Olsson, M.A.; García-Sosa, A.T.; Ryde, U. Binding affinities of the farnesoid X receptor in the D3R Grand Challenge 2 estimated by free-energy perturbation and docking. J. Comput. Mol. Des. 2017, 32, 211–224. [Google Scholar] [CrossRef]
  118. Åqvist, J.; Luzhkov, V.B.; Brandsdal, B.O. Ligand Binding Affinities from MD Simulations. Acc. Chem. Res. 2002, 35, 358–365. [Google Scholar] [CrossRef]
  119. Van Lipzig, M.M.H.; Ter Laak, A.M.; Jongejan, A.; Vermeulen, N.P.; Wamelink, M.; Geerke, D.P.; Meerman, J.H.N. Prediction of Ligand Binding Affinity and Orientation of Xenoestrogens to the Estrogen Receptor by Molecular Dynamics Simulations and the Linear Interaction Energy Method. J. Med. Chem. 2004, 47, 1018–1030. [Google Scholar] [CrossRef]
  120. Rifai, E.A.; Van Dijk, M.; Vermeulen, N.P.; Geerke, D.P. Binding free energy predictions of FXR agonists using LIE with reliability estimation: Application to the D3R Grand Challenge 2. J. Comput. Mol. Des. 2017, 32, 239–249. [Google Scholar] [CrossRef] [Green Version]
  121. Simonson, T. Protein: Ligand recognition: Simple models for electrostatic effects. Curr. Pharm. Des. 2013, 19, 4241–4256. [Google Scholar] [CrossRef] [PubMed]
  122. Gelpí, J.L.; Hospital, A.; Goñi, J.R.; Orozco, M. Molecular dynamics simulations: Advances and applications. Adv. Appl. Bioinform. Chem. 2015, 8, 37–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  123. Freindorf, M.; Furlani, T.R.; Kong, J.; Cody, V.; Davis, F.B.; Davis, P.J. Combined QM/MM Study of Thyroid and Steroid Hormone Analogue Interactions with αvβ3 Integrin. J. Biomed. Biotechnol. 2012, 2012, 1–12. [Google Scholar] [CrossRef] [PubMed]
  124. Kalaiarasi, C.; Manjula, S.; Kumaradhas, P. Combined quantum mechanics/molecular mechanics (QM/MM) methods to understand the charge density distribution of estrogens in the active site of estrogen receptors. RSC Adv. 2019, 9, 40758–40771. [Google Scholar] [CrossRef] [Green Version]
  125. Costa, A.H.L.; Clemente, W.S.; Bezerra, K.S.; Neto, J.X.L.; De Albuquerque, É.L.; Fulco, U.; Clemente, W.S., Jr. Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists. New J. Chem. 2018, 42, 19801–19810. [Google Scholar] [CrossRef]
  126. Hilder, T.A.; Hodgkiss, J.M. Molecular Mechanism of Binding between 17β-Estradiol and DNA. Comput. Struct. Biotechnol. J. 2016, 15, 91–97. [Google Scholar] [CrossRef] [PubMed]
  127. Eisold, A.; LaBudde, D. Detailed Analysis of 17β-Estradiol-Aptamer Interactions: A Molecular Dynamics Simulation Study. Molecules 2018, 23, 1690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Lakhin, A.V.; Tarantul, V.Z.; Gening, L.V. Aptamers: Problems, Solutions and Prospects. Acta Naturae 2013, 5, 34–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  129. Dutkiewicz, Z.; Mikstacka, R. Structure-Based Drug Design for Cytochrome P450 Family 1 Inhibitors. Bioinorg. Chem. Appl. 2018, 2018, 3924608. [Google Scholar] [CrossRef] [PubMed]
  130. Moldogazieva, N.T.; Ostroverkhova, D.S.; Kuzmich, N.N.; Kadochnikov, V.V.; Terentiev, A.A.; Porozov, Y. Elucidating Binding Sites and Affinities of ERα Agonists and Antagonists to Human Alpha-Fetoprotein by In Silico Modeling and Point Mutagenesis. Int. J. Mol. Sci. 2020, 21, 893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  131. Ma, K.; Saha, P.K.; Chan, L.; Moore, D.D. Farnesoid X receptor is essential for normal glucose homeostasis. J. Clin. Investig. 2006, 116, 1102–1109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  132. Yang, F.; Huang, X.; Yi, T.; Yen, Y.; Moore, D.D.; Huang, W. Spontaneous Development of Liver Tumors in the Absence of the Bile Acid Receptor Farnesoid X Receptor. Cancer Res. 2007, 67, 863–867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Levin, E.R. Cellular Functions of the Plasma Membrane Estrogen Receptor. Trends Endocrinol. Metab. 1999, 10, 374–377. [Google Scholar] [CrossRef]
  134. Soltysik, K.; Czekaj, P. Membrane estrogen receptors - is it an alternative way of estrogen action? J. Physiol. Pharmacol. Off. J. Pol. Physiol. Soc. 2013, 64, 129–142. [Google Scholar]
  135. Guo, J.J.; Yang, D.-P.; Tian, X.; Vemuri, V.K.; Yin, D.; Li, C.; Duclos, R.I.; Shen, L.; Ma, X.; Janero, D.R.; et al. 17β-estradiol (E2) in membranes: Orientation and dynamic properties. Biochim. Biophys. Acta (BBA) Biomembr. 2016, 1858, 344–353. [Google Scholar] [CrossRef]
  136. Vogel, A.; Scheidt, H.A.; Feller, S.E.; Metso, J.; Badeau, R.M.; Tikkanen, M.J.; Wähälä, K.; Jauhiainen, M.; Huster, D. The Orientation and Dynamics of Estradiol and Estradiol Oleate in Lipid Membranes and HDL Disc Models. Biophys. J. 2014, 107, 114–125. [Google Scholar] [CrossRef] [Green Version]
  137. Goh, J.Y.; Goh, K.S.; Yip, Y.M.; Ng, C.K. High salinity enhances adsorption of 17α-ethinyl estradiol by polyethersulfone membrane: Isotherm modelling and molecular simulation. engrXiv Preprints 2019. [Google Scholar] [CrossRef]
  138. Lakshmanan, S.; Kanwal, A.; Liu, S.; Patlolla, A.; Iqbal, Z.; Mitra, S.; Thomas, G.A.; Fagan, J.A.; Farrow, R.C. Improved Electrophoretic Deposition of Vertical Single Wall Carbon Nanotubes with Nanoscopic Electrostatic Lenses. Micromachines 2020, 11, 324. [Google Scholar] [CrossRef] [Green Version]
  139. Ulissi, Z.W.; Zhang, J.; Sresht, V.; Blankschtein, D.; Strano, M.S. 2D Equation-of-State Model for Corona Phase Molecular Recognition on Single-Walled Carbon Nanotube and Graphene Surfaces. Langmuir 2014, 31, 628–636. [Google Scholar] [CrossRef]
  140. Sun, W.; Li, M.; Zhang, W.; Wei, J.; Chen, B.; Wang, C. Sediments inhibit adsorption of 17β-estradiol and 17α-ethinylestradiol to carbon nanotubes and graphene oxide. Environ. Sci. Nano 2017, 4, 1900–1910. [Google Scholar] [CrossRef]
  141. Jiang, L.; Liu, Y.; Liu, S.; Zeng, G.; Hu, X.-J.; Hu, X.; Guo, Z.; Tan, X.; Wang, L.; Wu, Z. Adsorption of Estrogen Contaminants by Graphene Nanomaterials under Natural Organic Matter Preloading: Comparison to Carbon Nanotube, Biochar, and Activated Carbon. Environ. Sci. Technol. 2017, 51, 6352–6359. [Google Scholar] [CrossRef] [PubMed]
  142. Boateng, L.K.; Heo, J.; Flora, J.R.V.; Park, Y.-G.; Yoon, Y. Molecular level simulation of the adsorption of bisphenol A and 17α-ethinyl estradiol onto carbon nanomaterials. Sep. Purif. Technol. 2013, 116, 471–478. [Google Scholar] [CrossRef]
  143. Zaib, Q.; Khan, I.A.; Saleh, N.B.; Flora, J.R.V.; Park, Y.-G.; Yoon, Y. Removal of Bisphenol A and 17β-Estradiol by Single-Walled Carbon Nanotubes in Aqueous Solution: Adsorption and Molecular Modeling. Water Air Soil Pollut. 2012, 223, 3281–3293. [Google Scholar] [CrossRef]
  144. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864. [Google Scholar] [CrossRef] [Green Version]
  145. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef] [Green Version]
  146. Grimme, S. Density functional theory with London dispersion corrections. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 211–228. [Google Scholar] [CrossRef]
  147. Clark, S.J.; Segall, M.D.; Pickard, C.J.; Hasnip, P.J.; Probert, M.I.J.; Refson, K.; Payne, M.C. First principles methods using CASTEP. Z. Krist. Cryst. Mater. 2005, 220, 567–570. [Google Scholar] [CrossRef] [Green Version]
  148. Amusia, M.Y.; Msezane, A.Z.; Shaginyan, V.R. Density Functional Theory versus the Hartree–Fock Method: Comparative Assessment. Phys. Scr. 2003, 68, C133–C140. [Google Scholar] [CrossRef] [Green Version]
  149. Mota, K.; Neto, J.X.L.; Costa, A.L.; Oliveira, J.; Bezerra, K.; De Albuquerque, É.L.; Caetano, E.W.S.; Freire, V.; Fulco, U. A quantum biochemistry model of the interaction between the estrogen receptor and the two antagonists used in breast cancer treatment. Comput. Theor. Chem. 2016, 1089, 21–27. [Google Scholar] [CrossRef]
  150. Mazurek, A.H.; Szeleszczuk, Ł.; Pisklak, D.M. Periodic DFT Calculations—Review of Applications in the Pharmaceutical Sciences. Pharmaceutics 2020, 12, 415. [Google Scholar] [CrossRef]
  151. Stevenson, E.L.; Lancaster, R.W.; Buanz, A.B.M.; Price, L.S.; Tocher, D.A.; Price, S.L. The solid state forms of the sex hormone 17-β-estradiol. CrystEngComm 2019. [Google Scholar] [CrossRef] [Green Version]
  152. Du, R.; Xu, J.; Zhang, L.; Ning, L.; Li, S. Ethinyl estradiol cocrystals assembled by chain structures: Improvement in stability and solubility. New J. Chem. 2019, 43, 16889–16897. [Google Scholar] [CrossRef]
  153. Wang, J.-R.; Yang, Y.; Chen, X.; Mei, X. Solid-state characterization of 17β-estradiol co-crystals presenting improved dissolution and bioavailability. CrystEngComm 2016, 18, 3498–3505. [Google Scholar] [CrossRef]
  154. Watanabe, H.; Okiyama, Y.; Nakano, T.; Tanaka, S. Incorporation of solvation effects into the fragment molecular orbital calculations with the Poisson–Boltzmann equation. Chem. Phys. Lett. 2010, 500, 116–119. [Google Scholar] [CrossRef]
  155. Szeleszczuk, Ł.; Pisklak, D.M.; Zielińska-Pisklak, M.; Wawer, I. Effects of structural differences on the NMR chemical shifts in cinnamic acid derivatives: Comparison of GIAO and GIPAW calculations. Chem. Phys. Lett. 2016, 653, 35–41. [Google Scholar] [CrossRef]
  156. Charpentier, T. The PAW/GIPAW approach for computing NMR parameters: A new dimension added to NMR study of solids. Solid State Nucl. Magn. Reson. 2011, 40, 1–20. [Google Scholar] [CrossRef] [PubMed]
  157. Elyashberg, M.E. Identification and structure elucidation by NMR spectroscopy. TrAC Trends Anal. Chem. 2015, 69, 88–97. [Google Scholar] [CrossRef]
  158. Szeleszczuk, Ł.; Pisklak, D.M.; Zielińska-Pisklak, M.; Jurczak, E. A new polymorph of 17-β-estradiol and the application of different analytical techniques (ssNMR, PXRD, DSC, and FTIR) for its study. J. Mol. Struct. 2019, 1183, 274–280. [Google Scholar] [CrossRef]
  159. Singh, H.; Singh, S.; Srivastava, A.; Tandon, P.; Bharti, P.; Kumar, S.; Maurya, R. Conformational analysis and vibrational study of daidzein by using FT-IR and FT-Raman spectroscopies and DFT calculations. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2014, 120, 405–415. [Google Scholar] [CrossRef] [PubMed]
  160. Machado, N.; De Carvalho, L.A.E.B.; Otero, J.C.; Marques, M.P.M. A conformational study of hydroxyflavones by vibrational spectroscopy coupled to DFT calculations. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2013, 109, 116–124. [Google Scholar] [CrossRef] [Green Version]
  161. A Minaeva, V.; Minaev, B.F.; Hovorun, D.M. Vibrational spectra of the steroid hormones, estradiol and estriol, calculated by density functional theory. The role of low-frequency vibrations. Ukr. Biokhimichnyi Zhurnal 2009, 80, 82–95. [Google Scholar]
  162. Scheidt, H.A.; Badeau, R.M.; Huster, D. Investigating the membrane orientation and transversal distribution of 17β-estradiol in lipid membranes by solid-state NMR. Chem. Phys. Lipids 2010, 163, 356–361. [Google Scholar] [CrossRef] [PubMed]
  163. Vedad, J.; Mojica, E.-R.E.; Desamero, R.Z. Raman spectroscopic discrimination of estrogens. Vib. Spectrosc. 2018, 96, 93–100. [Google Scholar] [CrossRef] [PubMed]
  164. Oren, I.; Fleishman, S.J.; Kessel, A.; Tal, N.B. Free Diffusion of Steroid Hormones Across Biomembranes: A Simplex Search with Implicit Solvent Model Calculations. Biophys. J. 2004, 87, 768–779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  165. Ellena, J.; De Paula, K.; De Melo, C.C.; Da Silva, C.C.P.; Bezerra, B.P.; Venâncio, T.; Ayala, A.P. Temperature-Driven Isosymmetric Reversible Phase Transition of the Hormone Estradiol 17β Valerate. Cryst. Growth Des. 2014, 14, 5700–5709. [Google Scholar] [CrossRef]
  166. Morishima, F.; Inokuchi, Y.; Ebata, T. Laser Spectroscopic Study of β-Estradiol and Its Monohydrated Clusters in a Supersonic Jet. J. Phys. Chem. A 2012, 116, 8201–8208. [Google Scholar] [CrossRef] [Green Version]
  167. Borah, M.M.; Devi, T.G. The vibrational spectroscopic studies and molecular property analysis of Estradiol, Tamoxifen and their interaction by density functional theory. J. Mol. Struct. 2018, 1163, 205–220. [Google Scholar] [CrossRef]
  168. Cherkasova, O.; Nazarov, M.; Mankova, A.; Fedulova, E.; Volodin, V.; Minaeva, V.A.; Minaev, B.F.; Baryshnikov, G.V. Terahertz time-domain spectroscopy of testosterone, estradiol and estriol. In Proceedings of the 2010 International Kharkov Symposium on Physics and Engineering of Microwaves, Millimeter and Submillimeter Waves, Kharkiv, Ukraine, 21–26 June 2010; pp. 1–3. [Google Scholar]
  169. Hafizi, R.; Taheri, R.A.; Moghimi, H. Liquid phase extraction of nanosized biologically active estrogenic pollutants by using an efficient adsorbent. J. Mol. Liq. 2018, 266, 535–539. [Google Scholar] [CrossRef]
  170. Donini, C.A.; Silva, M.K.L.; Simões, R.; Cesarino, I. Reduced graphene oxide modified with silver nanoparticles for the electrochemical detection of estriol. J. Electroanal. Chem. 2018, 809, 67–73. [Google Scholar] [CrossRef] [Green Version]
  171. Klamt, A.; Schüürmann, G. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J. Chem. Soc. Perkin Trans. 2 1993, 2, 799–805. [Google Scholar] [CrossRef]
  172. Mennucci, B.; Tomasi, J.; Cammi, R.; Cheeseman, J.R.; Frisch, M.J.; Devlin, F.J.; Gabriel, S.; Stephens, P.J. Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules. J. Phys. Chem. A 2002, 106, 6102–6113. [Google Scholar] [CrossRef]
  173. Chinnasamy, K.; Kumaradhas, P. Intermolecular interactions and charge density distribution of endocrine-disrupting molecules (xenoestrogens) with ERα: QM/MM perspective. Struct. Chem. 2020, 31, 1013–1028. [Google Scholar] [CrossRef]
  174. Ruiz, P.; Ingale, K.; Wheeler, J.S.; Mumtaz, M. 3D QSAR studies of hydroxylated polychlorinated biphenyls as potential xenoestrogens. Chemosphere 2016, 144, 2238–2246. [Google Scholar] [CrossRef] [PubMed]
  175. Kim, M.; Li, L.Y.; Grace, J.R. Predictability of physicochemical properties of polychlorinated dibenzo-p-dioxins (PCDDs) based on single-molecular descriptor models. Environ. Pollut. 2016, 213, 99–111. [Google Scholar] [CrossRef]
  176. Eguchi, A.; Hanazato, M.; Suzuki, N.; Matsuno, Y.; Todaka, E.; Mori, C. Maternal–fetal transfer rates of PCBs, OCPs, PBDEs, and dioxin-like compounds predicted through quantitative structure–activity relationship modeling. Environ. Sci. Pollut. Res. 2015, 25, 7212–7222. [Google Scholar] [CrossRef]
  177. Delfosse, V.; Grimaldi, M.; Cavaillès, V.; Balaguer, P.; Bourguet, W. Structural and Functional Profiling of Environmental Ligands for Estrogen Receptors. Environ. Health Perspect. 2014, 122, 1306–1313. [Google Scholar] [CrossRef] [Green Version]
  178. Nwachukwu, J.C.; Srinivasan, S.; Bruno, N.E.; Nowak, J.; Wright, N.J.; Minutolo, F.; Rangarajan, E.S.; Izard, T.; Yao, X.-Q.; Grant, B.J.; et al. Systems Structural Biology Analysis of Ligand Effects on ERα Predicts Cellular Response to Environmental Estrogens and Anti-hormone Therapies. Cell Chem. Boil. 2017, 24, 35–45. [Google Scholar] [CrossRef]
  179. Cozzini, P.; Dellafiora, L. In silico approach to evaluate molecular interaction between mycotoxins and the estrogen receptors ligand binding domain: A case study on zearalenone and its metabolites. Toxicol. Lett. 2012, 214, 81–85. [Google Scholar] [CrossRef]
  180. Zhang, J.; Wu, W.; Song, Y.; Hou, L.; Li, T.; Guan, T.; Zhang, T.; Wang, Y. Homogeneous assay for zearalenone analogues and their docking studies with apo-/holo-estrogen receptors. Anal. Methods 2019, 11, 192–199. [Google Scholar] [CrossRef]
  181. Dellafiora, L.; Oswald, I.P.; Dorne, J.-L.; Galaverna, G.; Battilani, P.; Dall’Asta, C. An in silico structural approach to characterize human and rainbow trout estrogenicity of mycotoxins: Proof of concept study using zearalenone and alternariol. Food Chem. 2020, 312, 126088. [Google Scholar] [CrossRef]
  182. Yang, J.; Hu, C.T.; Zhu, X.; Zhu, Q.; Ward, M.D.; Kahr, B. DDT Polymorphism and the Lethality of Crystal Forms. Angew. Chem. 2017, 129, 10299–10303. [Google Scholar] [CrossRef]
  183. Zhang, H.; He, W.; Luo, X.; Lin, X.; Lu, X. Adsorption of 2,3,7,8-tetrochlorodibenzo-p-dioxins on intrinsic, defected, and Ti (N, Ag) doped graphene: A DFT study. J. Mol. Model. 2014, 20, 1–7. [Google Scholar] [CrossRef] [PubMed]
  184. Liu, C.; Li, H.; Johnston, C.T.; Boyd, S.A.; Teppen, B.J. Relating Clay Structural Factors to Dioxin Adsorption by Smectites: Molecular Dynamics Simulations. Soil Sci. Soc. Am. J. 2012, 76, 110–120. [Google Scholar] [CrossRef]
  185. Zhang, J.; Wu, W.; Wang, Y.; Xing, X.; Zhong, S.; Guan, T.; Zhang, T.; Hou, L.; Li, T. Estrogen receptor-based fluorescence polarization assay for bisphenol analogues and molecular modeling study of their complexation mechanism. Anal. Chim. Acta 2018, 1032, 107–113. [Google Scholar] [CrossRef] [PubMed]
  186. Sengupta, S.; Obiorah, I.; Maximov, P.; Curpan, R.; Jordan, V.C. Molecular mechanism of action of bisphenol and bisphenol A mediated by oestrogen receptor alpha in growth and apoptosis of breast cancer cells. Br. J. Pharmacol. 2013, 169, 167–178. [Google Scholar] [CrossRef] [PubMed]
  187. Liu, Y.; Qu, K.; Hai, Y.; Zhao, C. Bisphenol A (BPA) binding on full-length architectures of estrogen receptor. J. Cell. Biochem. 2018, 119, 6784–6794. [Google Scholar] [CrossRef] [PubMed]
  188. Zhang, J.; Li, T.; Wang, T.; Yuan, C.; Zhong, S.; Guan, T.; Li, Z.; Wang, Y.; Yu, H.; Luo, Q.; et al. Estrogenicity of halogenated bisphenol A: In vitro and in silico investigations. Arch. Toxicol. 2017, 92, 1215–1223. [Google Scholar] [CrossRef] [PubMed]
  189. Cao, H.; Wang, F.; Liang, Y.; Wang, H.; Zhang, A.; Song, M. Experimental and computational insights on the recognition mechanism between the estrogen receptor α with bisphenol compounds. Arch. Toxicol. 2017, 91, 3897–3912. [Google Scholar] [CrossRef]
  190. Zhuang, S.; Zhang, C.; Liu, W. Atomic Insights into Distinct Hormonal Activities of Bisphenol A Analogues toward PPARγ and ERα Receptors. Chem. Res. Toxicol. 2014, 27, 1769–1779. [Google Scholar] [CrossRef]
  191. Li, L.; Wang, Q.; Zhang, Y.; Niu, Y.; Yao, X.; Liu, H. The Molecular Mechanism of Bisphenol A (BPA) as an Endocrine Disruptor by Interacting with Nuclear Receptors: Insights from Molecular Dynamics (MD) Simulations. PLoS ONE 2015, 10, e0120330. [Google Scholar] [CrossRef]
  192. Delfosse, V.; Grimaldi, M.; Pons, J.-L.; Boulahtouf, A.; Le Maire, A.; Cavaillès, V.; Labesse, G.; Bourguet, W.; Balaguer, P. Structural and mechanistic insights into bisphenols action provide guidelines for risk assessment and discovery of bisphenol A substitutes. Proc. Natl. Acad. Sci. USA 2012, 109, 14930–14935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  193. Wei, D.; Li, J.; Chen, Z.; Liang, L.; Ma, J.; Wei, M.; Ai, Y.; Wang, X. Understanding bisphenol-A adsorption in magnetic modified covalent organic frameworks: Experiments coupled with DFT calculations. J. Mol. Liq. 2020, 301, 112431. [Google Scholar] [CrossRef]
  194. Bao, S.; Wu, S.; Huang, L.; Xu, X.; Xu, R.; Li, Y.; Liang, Y.; Yang, M.; Yoon, D.K.; Lee, M.; et al. Supramolecular Nanopumps with Chiral Recognition for Moving Organic Pollutants from Water. ACS Appl. Mater. Interfaces 2019, 11, 31220–31226. [Google Scholar] [CrossRef] [PubMed]
  195. Díaz, I.; Díez, E.; Camacho, J.; León, S.; Ovejero, G.; Cabanillas, S.L. Comparison between three predictive methods for the calculation of polymer solubility parameters. Fluid Phase Equilibria 2013, 337, 6–10. [Google Scholar] [CrossRef]
  196. Wang, C.-Y.; Zhang, Y.-J.; Wang, W.-K.; Pei, D.-N.; Huang, G.-X.; Chen, J.-J.; Zhang, X.; Yu, H.-Q. Enhanced photocatalytic degradation of bisphenol A by Co-doped BiOCl nanosheets under visible light irradiation. Appl. Catal. B Environ. 2018, 221, 320–328. [Google Scholar] [CrossRef]
  197. Motta, A.; La Mantia, F.P.; Ascione, L.; Mistretta, M. Theoretical study on the decomposition mechanism of bisphenol A polycarbonate induced by the combined effect of humidity and UV irradiation. J. Mol. Graph. Model. 2020, 99, 107622. [Google Scholar] [CrossRef]
  198. Liu, Y.; Liu, Y.; Liu, Z.; Du, F.; Qin, G.; Li, G.; Hu, X.; Xu, Z.; Cai, Z. Supramolecularly imprinted polymeric solid phase microextraction coatings for synergetic recognition nitrophenols and bisphenol A. J. Hazard. Mater. 2019, 368, 358–364. [Google Scholar] [CrossRef]
  199. Dvorakova, M.; Kejlová, K.; Rucki, M.; Jírová, D. Selected bisphenols and phthalates screened for estrogen and androgen disruption by in silico and in vitro methods. Neuro Endocrinol. Lett. 2018, 39, 409–416. [Google Scholar]
  200. Zhu, Q.; Liu, L.; Zhou, X.; Ma, M. In silico study of molecular mechanisms of action: Estrogenic disruptors among phthalate esters. Environ. Pollut. 2019, 255, 113193. [Google Scholar] [CrossRef]
  201. Josh, M.S.; Pradeep, S.; Adarsh, V.; Amma, K.V.; Devi, R.S.; Balachandran, S.; Sreejith, M.; Jaleel, U.A.; Benjamin, S. In silicoevidences for the binding of phthalates onto human estrogen receptor α, β subtypes and human estrogen-related receptor γ. Mol. Simul. 2013, 40, 408–417. [Google Scholar] [CrossRef]
  202. Sheikh, I.A.; Turki, R.F.; Abuzenadah, A.M.; Damanhouri, G.A.; A Beg, M. Endocrine Disruption: Computational Perspectives on Human Sex Hormone-Binding Globulin and Phthalate Plasticizers. PLoS ONE 2016, 11, e0151444. [Google Scholar] [CrossRef] [PubMed]
  203. Wang, S.; Wang, S.; Chen, M.; Xu, D.; Tang, L.; Wang, S. Elucidating Adsorption Mechanisms of Phthalate Esters upon Carbon Nanotubes/Graphene and Natural Organic Acid Competitive Effects in Water by DFT and MD Calculations. Bull. Korean Chem. Soc. 2015, 36, 1631–1636. [Google Scholar] [CrossRef]
  204. Liu, Y.; Zhang, R.; Wang, X.; Sun, P.; Chen, W.; Shen, J.; Xue, G. Hydrogenation induced deviation of temperature and concentration dependences of polymer-solvent interactions in poly(vinyl chloride) and a new eco-friendly plasticizer. Eur. Phys. J. Plus 2015, 130, 11. [Google Scholar] [CrossRef]
  205. Jacob, R.B.; Andersen, T.; McDougal, O.M. Accessible High-Throughput Virtual Screening Molecular Docking Software for Students and Educators. PLoS Comput. Boil. 2012, 8, e1002499. [Google Scholar] [CrossRef] [Green Version]
  206. Chaput, L.; Mouawad, L. Efficient conformational sampling and weak scoring in docking programs? Strategy of the wisdom of crowds. J. Chem. 2017, 9, 37. [Google Scholar] [CrossRef]
  207. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
  208. Software Website. GOLD—Protein Ligand Docking Software. Available online: www.ccdc.cam.ac.uk/solutions/csd-discovery/components/gold/ (accessed on 25 August 2020).
  209. Clark, M.; Cramer, R.D.; Van Opdenbosch, N. Validation of the general purpose tripos 5.2 force field. J. Comput. Chem. 1989, 10, 982–1012. [Google Scholar] [CrossRef]
  210. Software Website. Available online: https://www.charmm.org/ (accessed on 25 August 2020).
  211. Brooks, B.R.; Brooks, C.L.; Mackerell, A.D., Jr.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef]
  212. Software Website. Available online: http://www.swissdock.ch/ (accessed on 25 August 2020).
  213. Grosdidier, A.; Zoete, V.; Michielin, O. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011, 39, W270–W277. [Google Scholar] [CrossRef] [Green Version]
  214. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; Beer, T.A.P.D.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [Green Version]
  215. Dominguez, C.; Boelens, R.; Bonvin, A.M.J.J. HADDOCK: A Protein−Protein Docking Approach Based on Biochemical or Biophysical Information. J. Am. Chem. Soc. 2003, 125, 1731–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  216. Jian, Y.; He, Y.; Yang, J.; Han, W.; Zhai, X.; Zhao, Y.; Li, Y. Molecular Modeling Study for the Design of Novel Peroxisome Proliferator-Activated Receptor Gamma Agonists Using 3D-QSAR and Molecular Docking. Int. J. Mol. Sci. 2018, 19, 630. [Google Scholar] [CrossRef] [Green Version]
  217. Software Website. Available online: https://www.schrodinger.com/maestro (accessed on 25 August 2020).
  218. McAliley, J.H.; Bruce, D.A. Development of Force Field Parameters for Molecular Simulation of Polylactide. J. Chem. Theory Comput. 2011, 7, 3756–3767. [Google Scholar] [CrossRef] [PubMed]
  219. Halgren, T.A. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Compt. Chem. 1996, 17, 490–519. [Google Scholar] [CrossRef]
  220. Software Website. Available online: https://gaussian.com/ (accessed on 25 August 2020).
  221. Forli, S.; Huey, R.; Pique, M.E.; Sanner, M.F.; Goodsell, D.S.; Olson, A.J. Computational protein–ligand docking and virtual drug screening with the AutoDock suite. Nat. Protoc. 2016, 11, 905–919. [Google Scholar] [CrossRef] [Green Version]
  222. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2009, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  223. Software Website. Available online: http://autodock.scripps.edu/ (accessed on 25 August 2020).
  224. Lipparini, F.; Mennucci, B. Perspective: Polarizable continuum models for quantum-mechanical descriptions. J. Chem. Phys. 2016, 144, 160901. [Google Scholar] [CrossRef] [Green Version]
  225. Li, D.-D.; Meng, X.-F.; Wang, Q.; Yu, P.; Zhao, L.-G.; Zhang, Z.-P.; Wang, Z.; Xiao, W. Consensus scoring model for the molecular docking study of mTOR kinase inhibitor. J. Mol. Graph. Model. 2018, 79, 81–87. [Google Scholar] [CrossRef]
  226. Software Website. Available online: https://ambermd.org/ (accessed on 25 August 2020).
  227. Case, D.A.; Iii, T.E.C.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [Green Version]
  228. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C.L. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins Struct. Funct. Bioinform. 2006, 65, 712–725. [Google Scholar] [CrossRef] [Green Version]
  229. Available online: http://www.cosmo-model.org/ (accessed on 25 August 2020).
  230. Grubmüller, H.; Heller, H.; Windemuth, A.; Schulten, K. Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions. Mol. Simul. 1991, 6, 121–142. [Google Scholar] [CrossRef]
  231. Softwae Website. Available online: www.gromacs.org (accessed on 25 August 2020).
  232. Softwae Website. Available online: https://www.vasp.at/ (accessed on 25 August 2020).
  233. Softwae Website. Available online: https://www.ks.uiuc.edu/Research/namd/ (accessed on 25 August 2020).
  234. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  235. Karamertzanis, P.G.; Pantelides, C.C. Ab initio crystal structure prediction. I. Rigid molecules. J. Comput. Chem. 2005, 26, 304–324. [Google Scholar] [CrossRef] [PubMed]
  236. Karamertzanis, P.G.; Pantelides, C.C. Ab initio crystal structure prediction. II. Flexible molecules. Mol. Phys. 2007, 105, 273–291. [Google Scholar] [CrossRef]
  237. Vasileiadis, M.; Pantelides, C.C.; Adjiman, C.S. Prediction of the crystal structures of axitinib, a polymorphic pharmaceutical molecule. Chem. Eng. Sci. 2015, 121, 60–76. [Google Scholar] [CrossRef] [Green Version]
  238. Delley, B. DFT studies: From molecules and molecular environments to surfaces and solids. Comput. Mater. Sci. 2000, 17, 122–126. [Google Scholar] [CrossRef]
  239. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  240. Softwae Website. Available online: http://www.castep.org/ (accessed on 25 August 2020).
  241. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  242. Mark, P.; Nilsson, L. Structure and Dynamics of the TIP3P, SPC, and SPC/E Water Models at 298 K. J. Phys. Chem. A 2001, 105, 9954–9960. [Google Scholar] [CrossRef]
Figure 1. Experimental (top, green) Cross Polarization Magic Angle Spinning (CP MAS) and calculated (bottom, blue) GIPAW 13C solid-state NMR spectra of E2. Very good agreement between calculated and experimental values proves the usefulness of DFT calculations in solid-state analysis of estrogens. More details in [158]. Source: Author’s archive.
Figure 1. Experimental (top, green) Cross Polarization Magic Angle Spinning (CP MAS) and calculated (bottom, blue) GIPAW 13C solid-state NMR spectra of E2. Very good agreement between calculated and experimental values proves the usefulness of DFT calculations in solid-state analysis of estrogens. More details in [158]. Source: Author’s archive.
Ijms 21 06411 g001
Figure 2. Experimental (green) and calculated (blue) IR spectra of β-estradiol hemihydrate, selected range 700–1650 cm−1. Such calculations enable proper band assignments and thus facilitate the spectrum analysis. More information can be found in [158]. Source: Author’s archive.
Figure 2. Experimental (green) and calculated (blue) IR spectra of β-estradiol hemihydrate, selected range 700–1650 cm−1. Such calculations enable proper band assignments and thus facilitate the spectrum analysis. More information can be found in [158]. Source: Author’s archive.
Ijms 21 06411 g002
Table 1. Reagents of the main metabolic processes regarding estradiol that have been investigated with in silico methods.
Table 1. Reagents of the main metabolic processes regarding estradiol that have been investigated with in silico methods.
SubstrateEnzymeProductRef.
TestosteroneAromatase (CYP219A1)Estradiol[16]
Estradiol17β-OH-dehydrogenase (17β-HSD)Estrone[24,25]
EstradiolCYP1B14-OH-hydroxylated estradiol[30]
EstradiolCYP1A1, CYP1A22-OH-hydroxylated estradiol[28,29]
EstradiolSulfotransferase (SULT)Inactivated estradiol (sulfated)[27]
Inactivated (sulfated) estradiolSulfatase (STS)Activated estradiol[26]
Table 2. Selected proteins involved in the metabolism of estrogens present in the RCSB PDB (Protein Data Bank) [48]. All structures were obtained using X-ray diffraction.
Table 2. Selected proteins involved in the metabolism of estrogens present in the RCSB PDB (Protein Data Bank) [48]. All structures were obtained using X-ray diffraction.
ProteinRCSB PDB Reference CodeResolution (Å)Incorporated Ligands
βER5TOA2.5Estradiol
βER-LBD1QKM1.8Genistein
Phosphorylated βER-LBD3OLL1.5Estradiol, N-peptide linking
αER-LBD3UUC2.1Bisphenol C
αER-LBD mutant4Q503.074-hydroxytamoxifen
αER-LBD mutant2QXS1.7Raloxifene
αER-LBD2R6Y2.0SERM
17β-HSD1IOL2.3017β-estradiol
17β-HSD6MNC2.40Estrone
17β-HSD6MNE1.86Estrone, NADP+
17β-HSD3DHE2.30Dehydroepiandrosterone (DHEA)
17β-HSD4FJ02.23,7-dihydroxy flavone
17β-HSD4FJ12.3Genistein
SULT1E11AQU1.6Estradiol, PAP cofactor
SULT1E14JVM1.994Flame retardant, PAP cofactor
Placental E1/DHEA STS1P492.6l-octylglucoside, N-acetylo-d-glucosamine, Ca2+, PO43−
CYP1B16OyV3.101Estradiol
Table 3. Selected technical computation data in terms of ER and (xeno)estrogens regarding the publications cited in this article.
Table 3. Selected technical computation data in terms of ER and (xeno)estrogens regarding the publications cited in this article.
Code/Software UsedForce Field or
DFT Functional and Basis Set
Type of CalculationRef. MethodRef. in Article
1GOLD Molecular docking[208][57,69]
2-Ghemical 2.95
-Swiss Dock
-Tripos 5.2
-CHARMM
-Geometry optimization
-Molecular docking
[209,210,211,212][70]
3-Swiss model
-Hex 8.0, HADDOCK
-OPLS-Homology of receptors
-Molecular docking
[212,213,214,215][87]
4-Swiss model
-SybylX
-Tripos 5.2-Homology of receptors
-Molecular docking
[209,214,216][84]
5Maestro SchrödingerOPLS 2005, Glide SP, XPMolecular docking[217,218][83,85,86]
6Maestro SchrödingerMMFF94Geometry optimization, molecular docking[217,218,219]
7-Gaussian09W
-AutoDockTools
-B3LYP/6-31G(d)
-AutoDockZN
-Geometry optimization
-Molecular docking
[220,221,222,223][169,184]
8Gaussian03B3LYP/6-311++g**, PCMHydration enthalpy[220,224][183]
9Maestro SchrödingerZINC database
OPLS 2005, Glide SP
eHiTS docking module
consensus score
Energy minimization
HTVS
rank
[205,217,218,225][75]
10Maestro SchrödingerOPLS 2005, Glide SP, XPHTVS[205,217,218][76]
11Maestro SchrödingerOPLS 2005, Glide SP, XPSegregation: agonists/antagonists[217,218][77]
12Maestro Schrödinger-OPLS 2005, Grid (Glide)
-Desmond OPLS 2005
-Molecular docking, MD
-ADMET parameters
[217,218][90]
13-Maestro Schrödinger
-AMBER14
-AMBER14
-OPLS 2005, Glide
-FF03 (protein)
GAFF (ligand)
-MMPBSA, MMGBSA
-Docking
-MD
-Binding free energy, decomposition energy
[207,217,218,226,227,228][91]
14-MOPAC2016
-Gaussian09
-Gabedit package
-PM6 in HF, COSMO model
-B3LYP, PCM
-Verlet algorithm
-Pre-optimization, solvent model
-Optimization (DFT), solvent model
-MD
[156,220,224,229,230][169]
15-GOLD
-GROMACS
-Swiss Param Tool
-CHARMM27
-CHARMM27
-Molecular docking
-MD
-Ligand parametrization
[209,210,211,231][180]
16-SybylX
-AMBER11
-AutoDock 4.0
-Tripos 5.2
-AMBER
-AutoDockZN
-Geometry optimization
-MD
-Molecular docking
[216,221,222,223,226,227,228][186]
17-Gaussian09
-LeDock
-AMBER12
-AmberTools14
-B3LYP/-cc-pVTZ
-CHARMM
-AMBER
-MM/GBSA
-Geometry optimization
-Molecular docking
-MD
-Binding free energy
[207,220,221,222,223,226,227,228][188]
18-Gaussian09
-Molegro Virtual Docker
-AMBER Tools
-B3LYP/6-311++G(d,p)
-AMBER
-AMBER03
-Molecular electrostatic potential
-Molecular docking
-MD
[220,226,227,228][187]
19-Maestro Schrödinger
-Gaussian09
-AMBER10
-OPLS 2005
-HF, 6–31G*
-GAFF (ligand), ff03 (protein)
-Molecular docking
-Geometry optimization
-MD
[217,218,220,226,227,228][189]
20-VASP
-GROMACS
-PBE GGA (DFT-D3)
-GROMOS96
-Geometry optimization
-MD
[230,232][192]
21-GROMACS
-AutoDock Tools
-AutoDock Vina, Hex8.0.0
GROMACS
-AutoDockZN
-AutoDock Vina, GROMOS96
-Energy minimization
-Molecular docking
-MD
[221,222,223,231][199]
22-NAMD
-Spartan04
-Charm CMAP FF
-HF 3-21G
-MD
-QM
[233,234][125]
23-Gaussian 03
-AutoDock
-AMBBER
-B3LYP/6311**G
-AutoDockZN
-PM3/Amberff14SB FF
-Geometry optimization
-Molecular docking
-QM/MM
[220,221,222,223,226,227,228][114]
24-GROMACS
-Gaussian 09
-CHARMM (MM)
-GGA-D2 (QM)
-Geometry optimization
-DFT calculations
[210,211,220,221,222,223,229][115]
25-Maestro Schrödinger
-AMBER Tools
-OPLS 2005
Glide
-B3LYP/Amberff14SB
-Protein, ligand preparation (geometry optimization), molecular docking
-QM/MM
[217,218,226,227,228][124]
26-Crystal Predictor
-Crystal Optimizer (Gaussian)
-DMACRYS
-PBE0/631G(d,p)-Conformations
-Geometry optimization CSP
-Intermolecular lattice energies
[220,235,236,237][150]
27-GULP, DFTB+
-VASP
-optB88 level-Geometry pre-optimization CSP
-Geometry re-optimization
[232,235,236,237][181]
28DMol3DNP basis set, PBE GGAGeometry, energy optimization[238,239][182]
29CASTEPGGA PBEDFT, NMR[239,240][157]
30CASTEPGGA PBEDFT, structure parameters calculation[239,240][195]
31Gaussian09WB3LYP/631G(d)DFT, IR[238,241][164]
32Gaussian09M05-2X/6-311++G**DFT, IR[238,241][165]
33Gaussian09WB3LYP/6-31G (d,p)DFT, Raman[238,241][166]
34GaussianB3LYP/6-31G(d,p)DFT, IR[238,241][167]
The most widely used are commercial codes: Maestro Schrödinger [217], CASTEP [240], GOLD [208], Gaussian [220], AMBER [226], SybylX, VASP [220]; academic codes: AutoDock [209], CHARMM [210], GROMACS [232]. (AMBER and CHARMM are names of both the codes and the force fields.) AMBER (Assisted Model Building with Energy Refinement), CHARMM (Chemistry at HARvard Macromolecular Mechanics), CSP (Crystal Structure Prediction), GAFF (General AMBER Force Field), Glide (Grid-based Ligand Docking with Energetics), GOLD (Genetic Optimization for Ligand Docking), GROMACS (GROningen MAchine for Chemical Simulations), HADDOCK (High Ambiguity Driven protein-protein DOCKing), HTVS (high throughput virtual screening), MM/GBSA (Molecular Mechanics/Generalized Born Surface Area), PCM (polarizable continuum model), VASP (Vienna ab initio Simulation Package), Glide SP (Standard Precision), XP (Extra Precision). References in the last column refer to articles already cited in this review. These are examples of application of the listed methods in (xeno)estrogens research. References in the fourth column refer to articles that describe the theoretical basis of the listed software and calculation methods.
Table 4. Comparison of the calculation methods used in (xeno)estrogen investigations.
Table 4. Comparison of the calculation methods used in (xeno)estrogen investigations.
Calculation MethodPros and CapabilitiesCons and Limitations
Molecular docking-Explanation of a molecular basis for protein–ligand binding
-Relatively short calculation time
-Enables virtual screening for active compounds
-Lower accuracy when compared to QM methods
-Significant increase in time and complexity of calculations when combined with QM (QM/MM)
QSAR-Evaluation of estrogenicity
-No protein preparation needed
-No receptor–ligand binding data
-Large set of high-quality experimental data needed to obtain accurate results
QM (DFT-D)-High accuracy of calculations
-Simulation of IR, Raman, NMR spectra
-Thermodynamic calculations
-Long calculation time
-A lot of computational power needed
-Usually limited to small molecules and systems such as estrogen complexes, salts, co-crystals, etc.
QM/MM-High accuracy of calculations in the binding area (QM)
-Consideration of a whole complex (protein–ligand) with emphasis on the binding pocket
-Calculation time elongated due to QM
-Limitation of the QM-calculated area
MD-Simulation of dynamical processes
-Possibility to perform DFT-MD
-Significantly longer time required

Share and Cite

MDPI and ACS Style

Mazurek, A.H.; Szeleszczuk, Ł.; Simonson, T.; Pisklak, D.M. Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens. Int. J. Mol. Sci. 2020, 21, 6411. https://doi.org/10.3390/ijms21176411

AMA Style

Mazurek AH, Szeleszczuk Ł, Simonson T, Pisklak DM. Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens. International Journal of Molecular Sciences. 2020; 21(17):6411. https://doi.org/10.3390/ijms21176411

Chicago/Turabian Style

Mazurek, Anna Helena, Łukasz Szeleszczuk, Thomas Simonson, and Dariusz Maciej Pisklak. 2020. "Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens" International Journal of Molecular Sciences 21, no. 17: 6411. https://doi.org/10.3390/ijms21176411

APA Style

Mazurek, A. H., Szeleszczuk, Ł., Simonson, T., & Pisklak, D. M. (2020). Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens. International Journal of Molecular Sciences, 21(17), 6411. https://doi.org/10.3390/ijms21176411

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