Docking of MBP and BPA to Human ERα and ERβ
Docking of MBP into human ERα and ERβ using AutoDock 4
[25],
[26] and AutoDock Vina
[27] gave two symmetric poses, which is not surprising because MBP has a phenolic ring at each end []. BPA also had two poses for one of the rings in ERα and ERβ. We analyzed both poses for MBP and BPA in human ERα and ERβ. In our analysis of the 3D models of MBP and BPA in both ERs, we use the term “first phenolic ring” to describe the ring that has contacts with ERα and ERβ that are similar to the A ring of E2.
Analysis of the crystal structures of ERα complexed with E2
[33],
[34] and other estrogens
[35] revealed that Glu-353 and Arg394 have important stabilizing contacts with the C3 hydroxyl on the A ring and His-524 with the 17β-hydroxyl on the D ring. Glu-305, Arg-346 and His-475 on ERβ have similar stabilizing contacts with estrogens. As reported below, the presence or absence of these contacts in the 3D models of ERα and ERβ with BPA and MBP is important analyzing the interaction between these chemicals and the ERs.
Analysis of MBP in Orientation 1 in Human ERα
In , we show the 3D model of MBP in human ERα in Orientation 1, in which C1 on MBP is closest to first phenolic ring. The distances between MBP and ERα are shown in and . For comparison, in and we show the distances between E2 and human ERα
[33],
[35],
[36],
[37],
[38].
| Table 1Distances between MBP and ERα. |
| Table 2Distances between E2 and ERα and ERβ. |
The first phenolic ring on MBP has contacts that are similar to that of the A ring on E2 with human ERα
[33],
[34],
[39]. The phenolic hydroxyl on MBP is 3.1 Å from Oε2 on Glu-353, 3.1 Å from Nη2 on Arg-394 and 2.8 Å from the backbone oxygen of Leu-387. MBP is 3.9 Å from Cε2 on Phe-404 [, ]. These contacts are similar to that for E2 with human ERα, except that Cδ2 on Leu-387 does not contact the phenolic hydroxyl on MBP, in contrast to the contact between Leu-387 and E2 in human ERα [, ].
The second phenolic hydroxyl in MBP is 3 Å, 2.8 Å and 4 Å from Nδ1 on His-524, the backbone oxygen on Gly-521 and Cβ on Leu-525, respectively, on ERα []. This phenolic hydroxyl also contacts Met-343 and Met-421 on ERα. These five residues stabilize the D ring on E2 in human ERα [].
There are, however, differences in some interactions between ERα and MBP compared to that with E2. While Gly-521 and Met-421 contact the second phenolic hydroxyl on MBP [], Gly-521 and Met-421 contact C16 on E2 in ERα []. While Leu-384 has two van der Waals contacts with MBP, Leu-384 does not contact E2 in ERα. While Thr-347 has a van der Waals contact with MBP, Thr-347 does not contact E2 in ERα. While Cβ on Leu-391 is 3.7 Å from the first phenolic hydroxyl on MBP, this contact is absent between ERα and E2. While Ala-350 contacts the linker between the two phenolic rings on MBP, Ala-350 contacts C1 on the A ring in E2 in ERα. Phe-404 and Met-421 have van der Waals contacts with C1 on MBP, which has no equivalent in E2 in ERα.
Analysis of MBP in Orientation 2 in Human ERα
As shown in and , analysis of ERα with MBP in Orientation 2 reveals that MBP has contacts with Glu-353, Arg-394, Phe-404, Met-343, Leu-384, Met-421, Gly-521, His-524 and Leu-525 that are similar to those found in Orientation 1 of MBP in ERα. Due to the reversed orientation of MBP in ERα, C1 on MBP has a van der Waals contact with Thr-347, and the other part of the linker contacts Ala-350.
Analysis of MBP in Orientation 1 in Human ERβ
shows MBP in Orientation 1 in human ERβ. For comparison, in , we show E2 in human ERβ
[28]. Many of the contacts between MBP and human ERβ shown in and are similar to that between MBP in Orientation 1 and human ERα [, ] and between E2 and ERβ []. Like the A ring in E2, the first phenolic ring on MBP has stabilizing contacts with Glu-305, Arg-346, Phe-356, Leu-339 and Leu-343 in ERβ. The second phenolic ring contacts His-475, Gly-472, Leu-476, Ile-373 and Ile-376 [, ].
| Table 3Distances between MBP and ERβ. |
Analysis of MBP in Orientation 2 in Human ERβ
shows MBP in Orientation 2 in human ERβ. Many of the contacts between MBP in Orientation 2 and human ERβ [] are similar to that between MBP in Orientation 1 and human ERβ [, ]] and between E2 and ERβ []. The backbone oxygen on Leu-339, Cβ on Ala-302 and the side chains on Leu-343 contact the first phenolic ring on MBP. These contacts are absent between MBP in Orientation 2 in ERα [, ].
Analysis of BPA in Orientation 1 in Human ERα
shows BPA in Orientation 1 in human ERα. The phenolic ring on BPA, corresponding to the A ring of E2, has stabilizing contacts with Oε2 on Glu-353, Nη2 on Arg-394, Cδ2 on Phe-404 and Cβ on Ala-350 [, ]. However, the backbone oxygen on Leu-387 is 6.4 Å from the phenolic hydroxyl and Leu-391 does not have a van der Waals contact with the phenolic ring.
| Table 4Distances between BPA and ERα. |
The second phenolic ring does not contact either Gly-521, His-524 or Leu-525 []. Instead, phenolic ring moves so that it contacts Cε2 and Cη on Phe-404, Cε on Met-421 and Cδ2 on Ile-424. Also, Ala-350 and Leu-384 and Thr-347 contact the linker on BPA.
Analysis of BPA in Orientation 2 in Human ERα
In , we show BPA in Orientation 2 in human ERα. The first phenolic ring on BPA contacts Oε2 on Glu-353, Nη2 on Arg-394, Cδ2 on Phe-404 and Cβ on Ala-350 [, ]. Leu-391 has a van der Waals contact with the phenolic ring. However, the backbone oxygen on Leu-387 does not contact BPA.
The second phenolic ring on BPA does not contact either Gly-521, His-524, Met-421 or Ile-424 []. Instead, the phenolic hydroxyl contacts Leu-525 and Thr-347. Leu-384 and Ala-350 also contact the second phenolic ring.
Analysis of BPA in Orientation 1 in Human ERβ
In , we show BPA in Orientation 1 in human ERβ. The first phenolic ring on BPA contacts Oε2 on Glu-305, Nη2 on Arg-346, Cε2 on Phe-404, Cβ on Leu339 and the backbone nitrogen on Ala-302 [, ]. The backbone oxygen on Leu-339 does not contact the phenolic hydroxyl. Cβ on Ala-302 and Cδ2 on Leu-343 do not contact the phenolic ring.
| Table 5Distances between BPA and ERβ. |
The second phenolic ring contacts Gly-472, His-475, Leu-476, Ile-373 and Met-336, but the second phenolic ring does not contact Met-295.
Analysis of BPA in Orientation 2 in Human ERβ
In , we show the minimized structure of BPA in Orientation 2 in human ERβ. The first phenolic ring on BPA contacts Oε2 on Glu-305, Nη2 on Arg-346, Cε2 on Phe-404, Cβ on Ala-350, Cδ2 on Leu-343 and the backbone oxygen on Leu-339 [, ].
The second phenolic ring does not contact either Gly-472, His-475, Leu-476 or Met-295 []. Instead, the phenolic hydroxyl contacts Phe-377 and Ile-373. Interestingly, Phe-356 contacts the second phenolic ring and Met-336 contacts the linker on BPA.
Docking Energy Analysis
We used X-Score
[31] and DSX
[31] to estimate the affinity of MBP and BPA in their different orientations in ERα and ERβ. and summarize these analyses for the X-Score and DSX. For both algorithms, MBP has an affinity for ERα and ERβ that is closer to that of E2 than is BPA for these receptors. This is consistent with previous assays of the activity of MBP and BPA
[22],
[23].
| Table 6Docking analysis of MBP and BPA in ERα and ERβ. |
| Table 7Docking analysis of MBP and BPA in ERα and ERβ. |