AERAdvances in Enzyme Research2328-4846Scientific Research Publishing10.4236/aer.2022.104006AER-122263ArticlesBiomedical&Life Sciences Engineering Medicine&Healthcare Molecular Docking Studies of Estrone-Coumarin Derivatives as Aromatase and 17<i>β</i>-HSD1 Inhibitors Related to Hormone Receptor Positive (HR+) Breast Cancer SilviaAlejandra Meza-Ireta1BlancaColin-Lozano1PenélopeMerino-Montiel1*JoséLuis Vega-Báez1SaraMontiel-Smith1Facultad de Ciencias Químicas, Ciudad Universitaria, Benemérita Universidad Autónoma de Puebla, Puebla, México0912202210048310011, October 202227, December 2022 30, December 2022© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Hormone Receptor positive (HR+) breast cancer is the most common malignancy in women. New strategies in the treatments have targeted the estrogen biosynthesis pathways including the inhibition of the aromatase and 17 β-HSD1 enzymes. The present work, describes the study of a new family of 9 hybrid compounds derived from estrone attached to a coumarin fragment, linked through different lengths of hydrocarbon chains. The activity of these compounds was evaluated by molecular docking with two relevant enzymes in breast cancer (HR+). It has been proposed nine compounds as 17 β-HSD1 inhibitors and six as aromatase inhibitors. We found important interactions with key amino acids at the orthosteric site of each enzyme and their score values compared to the crystallographic ligand. The in silico analysis showed good score values in the proposed compounds, where the steroidal portion presented important interactions with Met374 and Tyr155 in aromatase and in 17 β-HSD1 respectively. Highlighting Compounds 2, 5 and 8 with an aromatic ring at the C4 position of the coumarin moiety, which favored arene-H type interactions essential for protein-ligand recognition. In addition, the results related to the 17 β-HSD1 enzyme demonstrated how the length of the linker influences the interaction; the best score was found for derivative 8 with a chain of 8 methylenes.

Estrone Coumarins Aromatase Inhibitors 17<i>β</i>-HSD1 Inhibitors Molecular Docking1. Introduction

Hormone Receptor positive (HR+) breast cancer is the most aggressive and frequent malignant neoplasm in women, representing the main cause of oncological mortality [1]. In this regard, aromatase inhibitors are compounds that suppress estrogen biosynthesis through reversible blockade of this enzyme and are mainly used as aids in the treatment of this condition [2]. Aromatase belongs to the group of the cytochrome P450 family, and in breast tissue, it catalyzes the demethylation of carbon 19 of androstenedione, producing an 18-carbon phenolic estrogen, estrone E1 [3]. Furthermore, estrone sulfate is the most abundant circulating estrogen in postmenopausal women and it can be converted to estrone via sulfatase [4]. In the same way, the 17β-Hydroxysteroid Dehydrogenase Type 1 (17β-HSD1) enzyme is involved in the catalytic reduction of estrone to estradiol in both pathways and it is overexpressed in breast cancer cells [5]. Moreover, the inhibition of this enzyme prevents the formation of estradiol (E2), a hormone that stimulates the growth of breast tumors and endometriosis through the activation of the Estrogen Receptor α (ERα) [6], therefore, it was proposed as an interesting pharmacological target.

It has been described steroidal and non-steroidal inhibitors for these molecular targets, steroid inhibitors usually bind reversibly to the catalytic site of the enzyme, mimicking either the substrate or the product. In this context, aromatase and 17β-HSD1 enzymes have in common that estrone bind to both [7].

However, these antiestrogen compounds have limitations: their mixed agonist, antagonist profile (ERα), and the development of tumor resistance, while the non-steroidal inhibitors showed low selectivity and inhibits other CYP450 isoforms. Based on the above, it is evident the necessity of new molecules with dual action, the selective inhibition of both targets and simultaneously to reduce the risk of adverse effects, drug interactions and resistance (Figure 1) [8] - [14].

Coumarins are molecules formed by the union of a benzene and an α-pyrone [15]. These molecules have shown a wide variety of pharmacological activities, mainly antiproliferative effects. In addition, it has been observed that coumarins show in vitro inhibitory activities of some enzymes involved in cancer, which could secondarily attenuate the clinical progress and the pathophysiological manifestations observed in this disease [16] [17] [18] [19]. Moreover, estrone derivatives or other steroidal compounds have greater specificity than non-steroidal ones for binding to enzymes involved in estrogen biosynthesis [20] [21]. Considering the described above, we chose a model of study derivatives of estrone bearing a coumarin fragment.

Leonetti et al. have reported that the chain length of carbon atoms that function as linkers, directly influences the potency of the molecules. The absence or the presence of very short or very large linkers results in the decrease of aromatase inhibition and selectivity (molecular promiscuity). Therefore, the length of the linker can affect the inhibitory potency by changing the position of the H-bond acceptors in the ligand binding site [22] [23] [24].

The drug development process requires a multidisciplinary effort comprising main stages: 1) target selection and validation; 2) compound screening and lead optimization; 3) preclinical studies; and 4) clinical trials [25]. Therefore, this whole process takes an average of 800 million dollars and 15 years to complete [26]; for this reason, tools such as Computer-Aided Drug Design (CADD) have gained popularity in recent years, playing an important role in drug development. Automated Molecular Docking is the most widely used instrument in CADD, its purpose is to propose a binding mode between two molecules, with overall minimum energy, considering the principle of complementarity of the structures and interaction forces involved in molecular recognition [27] [28]. In this sense, there are different types of interactions that can be essential in the affinity of a ligand with a protein target and in the stability of a ligand-protein complex. The nature of the interactions is diverse, for example: hydrogen-bridging type interactions (caused by steric effects due to entropy and solvent-related forces), electrostatic forces (charges), van der Waals interactions (electrodynamics forces) and hydrophobic interactions. The type of interaction observed is dependent on the size, shape, and physicochemical properties of the amino acid residues of the binding pocket and the class and flexibility of the ligands [29] [30].

Thus, Molecular Docking is very useful to predict the binding mode of a ligand to a protein and its possible affinity to it, optimizing the development and discovery of new bioactive molecules and/or drugs. It also improves the designation of the leading compounds through the structure-activity relationship, and facing the synthesis stage, it helps to filter and prioritize certain compounds from a large collection of molecules, saving time and costs.

Accordingly, it was decided to carry out molecular docking only of Compounds (1 - 6) with linkers of 4 and 6 carbon atoms on the aromatase enzyme.

The proposed molecules 1 - 6 (Figure 1) were designed based on the following pharmacophoric fragments: three coumarins as a polar group (non-steroidal fragment) attached through hydrocarbon chains of different lengths to estrone as an enzyme recognition group. In this work, we have designed and screened in silico a series of new estrone derivatives with different groups of coumarins against aromatase and 17β-HSD1 enzymes to find new chemical entities that could achieve benefits in the treatment of breast cancer (HR+).

2. Materials and Methods2.1. Molecular Docking

The crystal structures of the aromatase and 17β-HSD1 were obtained from the Protein Data Bank (https://www.rcsb.org/) with the codes PDB: 5 JKW and 1 FDT, respectively. The molecular docking was carried out with the Molecular Operating Environment software (https://www.chemcomp.com/MOE%20v2020.0901).

The unnecessary molecules were removed (H2O and SO4), the hydrogen atoms, charges and energy minimization were adjusted with the AMBER10: EHT force field and R-field solvation from the MOE suite. The ligands (1 - 6) were built with ChemDraw Professional 16.0 software. The minimization and ionization of the compounds were performed with AMBER10: EHT as the force field. As a placement function Triangle matcher was selected, the scores were calculated with the London ΔG function and the selection of the best poses was made using the GBVI/WSA ΔG (Generalized-Born Volume Integral/Weighted Surface area) scoring function. These parameters were established for docking validation. The binding site was defined using the crystallographic ligand (orthosteric site). For each ligand it was generated 100 repetitions and carried out 100 poses. The best binding poses were visually inspected, the pose that most overlapped with the crystallographic ligand was chosen and its scores (ΔG = kcal/mol) was reported. The graphical representations of ligand interactions were created in MOE.

2.2. Docking Validation

The Docking protocols were validated, by a redocking into the aromatase and 17β-HSD1 ligand-binding sites with its co-crystal-ligands Testosterone (TST) and estradiol (E2) respectively. The Root Mean Square Deviation (RMSD) between the co-crystal ligand and the redocked molecule of aromatase was 0.2912 ± 0.018 Å and to the17β-HSD1 was 0.1381 ± 0.0086 Å. Validation was performed in triplicate.

3. Results and Discussion3.1. Aromatase (CYP19A1)

To validate our model preparation, the Root Mean Square Deviation (RMSD) was included as a measure of re-docked success. In which a value of 0.2912 ± 0.018 Å was obtained. Figure 2 shows the pose with the lowest RMSD and a score of −8.3711 kcal/mol. It can be observed that the hydrophobic side chains Trp224, Phe221, Phe134, Ile133, Val370 and Met374 are placed around the testosterone molecule (TST, crystallographic ligand), shaping the ligand binding pocket (Figure 2). In addition, TST bonds with its 17β-hydroxyl oxygen and receives a hydrogen bond from the amide backbone -NH of Met374 [31] [32] [33].

To explain the presumable binding mode with human aromatase, Compounds 1 - 6 were subjected to an in silico molecular docking. The results revealed that Compounds 1 - 6 fit into the binding pocket and the carbonyl group (C17) of the steroid moiety form a hydrogen bond interaction with Met374. This portion achieved good overlap with the TST at the binding site and exhibited hydrophobic interactions with the amino acids that make up this region (Figure 2). It is important to mention that Compounds 2 and 5 obtained the best scores of −12.0826 kcal/mol and −12.0233 kcal/mol, respectively (Figure 3). It is known that aromatic rings can act as hydrogen bond acceptors stablishing meaningful interactions with the amino acid residues [34], the benzene ring attached to the coumarin fragment in Compounds 2 and 5 forms a H-arene interaction type with the amino acids Val369 and Val370.

Unlike most cytochrome P450 enzymes, near the active site of the aromatase there is a cavity large enough to shelter a secondary interaction with the substrate or with an additional non-steroidal molecule, this interaction could influence the enzymatic activity of the aromatase [22]. In this context, the coumarin

moiety of the studied compounds showed interactions involving the aromatic rings, these interactions are critical for protein-ligand recognition [35]. Compounds (2 and 5) substituted with a benzene ring at C4 of the coumarin portion (Figure 4), showed arene-H type interactions, the electron-rich π cloud interacts with hydrogen atoms (Δ+) of the amino acids Val369 and Val370 in Compound 2 (-NH), and with Ala438 in Compound 5 (-CH3). These additional interactions could be related with the score values obtained by the different ligands (Table 1), and therefore, with their probably affinity to the aromatase.

Additional interactions with the aromatase enzyme
Ligand∆G BindingConventional Hidrogen BondHydrophobic InteractionsOthers
*TST−8.3711Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477
1−11.3929Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443 Phe430 Val369 Pro429 Met364 Ile132 Leu152 Met303 Ala438 Met446 Phe148 Ile442Gly439
2−12.0826Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443 Phe430 Val369 Pro429 Met364 Pro368Val369 Val370
3−11.3581Met374 Ala438Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443 Ile132 Leu152 Met303 Ala438 Met446 Phe148 Ile442
4−11.7687Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443 Phe430 Pro429 Met364 Leu152 Met303 Ala438 Met446 Met447
5−12.0233Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443 Phe430 Val369 Pro429 Met364Ile132 Leu152 Met303 Ala438 Met446 Phe148 Met447 Phe203Met303 Val370 Ala438
6−11.7957Met374Trp224 Phe221 Phe134 Ile133 Val370 Met374 Ile305 Ala306 Val373 Leu372 Leu477 Ala307 Met311 Ala443Ile132 Leu152 Met303 Ala438 Met446 Phe148 Ile442Met303 Cys437

Moreover, Compounds 1, 3 and 6 displayed hydrogen bond interactions with Met374, a key residue for the ligand recognition and aromatase inhibition; furthermore, they showed additional interactions at the coumarin fragment. Compound 1 showed π-H type interactions with the Gly439 residue and the α-pyrone ring of the coumarin skeleton, Compound 3 formed a hydrogen bond with Ala438, residue which acts as a backbone donor, while for Compound 6, the residues Met303 and Cys437 (amino acids with a sulfur atom), behaved as side chain acceptors (Figure 5).

According to the structure/score relationship analysis (binding energy) it was predicted for Compounds 1 - 6 better scores than those of the crystallographic ligand TST. Moreover, no relevant difference was observed in the binding energy between the compounds with a chain of 4 and 6 carbon atoms (Figure 6). However, Compounds 2 and 5 with the benzene substituent at position 4 of the coumarin fragment obtained the best scores of −12.0826 and −12.0233 kcal/mol, respectively. This observation indicates that the aromatic type substituents improve and even reinforce the interaction with the aromatase binding site.

3.2. 17β-Hydroxysteroid Dehydrogenase Type 1 (17β-HSD1)

In the validation of the 17β-HSD1 enzyme, an RMSD value of 0.2067 ± 0.0228 Å was obtained. Figure 7 shows the pose with the lowest RMSD and a score of

−7.2788 kcal/mol. In this image, estradiol (E2, crystallographic ligand) is located near the catalytic triad (Ser142, Tyr155 and Lys195), the 17β-hydroxyl oxygen accepts a hydrogen bond from the -OH group of Tyr155.

The lipophilic pocket of the 17β-HSD1 is formed by Gly92, Leu93, Gly94, Ala291, Phe192 and Val196, which is connected to a hydrophilic zone composed of Ser142, Tyr155, Glu194 and Lys195; this hydrophilic zone is crucial for the catalytic reaction of estradiol [36] [37]. It has been observed that some inhibitors of 17β-HSD1, with a steroidal moiety attached to a polar fragment through different methylene chains, establish interactions between the polar group and the cofactor binding site. Among the compounds tested, those with the best results are the derivatives with linkers of 6 or 8 methylenes. [11] [38] [39] [40].

To favor secondary interactions that might improve the scores observed for 17β-HSD1, it was added three compounds bearing eight carbon atoms as a linker. The nine Compounds 1 - 9 showed interactions of the arene-H type with Tyr155 (Figure 8 and Figure 9) and some of them with the hydrophobic residues Pro187 and Phe259 at the binding site.

Notably, the presence of the coumarin attached through variable chain lengths to the A ring of the steroidal moiety has cause a reverse orientation of the molecule in the active site, with the A ring facing the catalytic site (Tyr155), and D ring pointing to the recognition edge of the binding site (Figure 8 and Figure 9). This inverse orientation causes inhibition of the enzyme, forming a dead-end complex that cannot be catalyzed. This type of inhibition, also reported in other enzymes, has been related with high concentrations of substrate [41] [42] [43] [44] [45].

In the conversion of E1 (estrone) to E2 catalyzed by 17β-HSD1, NADPH behave as a co-factor; at the cofactor binding site, NADPH is surrounded by the amino acids Gly9, Cys10, Ser11, Ser12, Ile14, Gly15, Arg37, Asn90, Gly92, and Lys195, forming important hydrogen-bonding interactions (Figure 10) [32]. It can be observed that for the products 1 - 9, the aliphatic linkers together with the coumarin fragments are in the same NADPHs cofactor binding region. Furthermore, in Compounds 1, 4, 6, 8 and 9, the oxygens of the lactone portion act

as hydrogen bond acceptors of the residues Ser11, Ser12, Gly15, Asn90 and Lys195. Moreover, the electron-rich π cloud of the unsaturated lactone moiety of Compounds 3, 5 and 7 interacted with hydrogen atoms (Δ+) of the amino acids Ile14, Leu93 and Gly92 in an arene-H interaction type. Additionally, Compounds 2 and 8 bearing a Ph group at the C4 position of coumarin, show an arene-H type interaction with Gly92 and Arg37 respectively (Figure 11).

Although all nine compounds obtained good scores compared to E2, as the methylene chain increased from 4 to 8 atoms, better scores were obtained, highlighting Compounds 7 and 8 with scores of −10.2370 kcal/mol and −11.3970 kcal/mol respectively (Figure 11).

4. Conclusions

The in silico analysis showed good score values in the proposed compounds, the steroidal portion presented important interactions with Met374 and Tyr155 in the catalytic site of aromatase and 17β-HSD1 respectively. Notably, Compounds 2, 5 and 8 with an aromatic ring attached to C4 of the coumarin moiety, favored arene-H type interactions, essential for protein-ligand recognition. In addition, for the enzyme 17β-HSD1, it was observed how the length of the hydrocarbon chain influences in the number of binding points, being Compound 8 with a chain of eight methylenes, the one that presented the best score. This linker allows the steroidal portion to locate into the substrate binding site, while the coumarin portion interacts with the cofactor binding site with amino acids Arg37 and Ser11.

As a result of the analysis performed in this work, it can be concluded that the increment of the length of the linker up to 8 methylenes, optimizes the number of interactions at different sections of the enzyme. Moreover, it has been established an excellent starting point for the synthetic stage of these promising molecules as inhibitors of these two enzymes related to breast cancer (HR+). Molecular docking helped predict the binding mode of the designed Compounds 1 - 9 and their target proteins and predict the affinity and probable activity of the molecules. The synthesis and enzyme inhibition studies of these compounds will help to corroborate the most suitable structure against these molecular targets.

Acknowledgements

Authors thank CONACYT-México (Grant 320334 “Paradigmas y Controversias 2022”). S. A. M.-I. thanks CONACYT-México for the predoctoral grant.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Meza-Ireta, S.A., Colin-Lozano, B., Merino-Montiel, P., Vega-Báez, J.L. and Montiel-Smith, S. (2022) Molecular Docking Studies of Estrone-Coumarin Derivatives as Aromatase and 17β-HSD1 Inhibitors Related to Hormone Receptor Positive (HR+) Breast Cancer. Advances in Enzyme Research, 10, 83-100. https://doi.org/10.4236/aer.2022.104006

ReferencesDeSantis, C.E., Ma, J., Goding Sauer, A., Newman, L.A. and Jemal, A. (2017) Breast Cancer Statistics, 2017, Racial Disparity in Mortality by State. CA: A Cancer Journal for Clinicians, 67, 439-448. https://doi.org/10.3322/caac.21412Brueggemeier, R.W., Hackett, J.C. and Diaz-Cruz, E.S. (2005) Aromatase Inhibitors in the Treatment of Breast Cancer. Endocrine Reviews, 26, 331-345. https://doi.org/10.1210/er.2004-0015Ghosh, D., Griswold, J., Erman, M. and Pangborn, W. (2009) Structural Basis for Androgen Specificity and Oestrogen Synthesis in Human Aromatase. Nature, 457, 219-223. https://doi.org/10.1038/nature07614Roberts, K., Rochefort, J., Bleau, G. and Chapdelaine, A. (1980) Plasma Estrone Sulfate Levels in Postmenopausal Women. Steroids, 35, 179-187. https://doi.org/10.1016/0039-128X(80)90101-4Aka, J.A., Mazumdar, M. and Lin, S.-X. (2009) Reductive 17β-Hydroxysteroid Dehydrogenases in the Sulfatase Pathway: Critical in the Cell Proliferation of Breast Cancer. Molecular and Cellular Endocrinology, 301, 183-190. https://doi.org/10.1016/j.mce.2008.10.042Aka, J.A., Mazumdar, M., Chen, C.Q., Poirier, D. and Lin, S.X. (2010) 17β-Hydroxysteroid Dehydrogenase Type 1 Stimulates Breast Cancer by Dihydrotestosterone Inactivation in Addition to Estradiol Production. Molecular Endocrinology, 24, 832-845. https://doi.org/10.1210/me.2009-0468Thomas, M.P. and Potter, B.V. (2013) The Structural Biology of Oestrogen Metabolism. The Journal of Steroid Biochemistry and Molecular Biology, 137, 27-49. https://doi.org/10.1016/j.jsbmb.2012.12.014Ramsay, R.R., Popovic-Nikolic, M.R., Nikolic, K., Uliassi, E. and Bolognesi, M.L. (2018) A Perspective on Multi-Target Drug Discovery and Design for Complex Diseases. Clinical and Translational Medicine, 7, Article No. 3. https://doi.org/10.1186/s40169-017-0181-2Tremblay, M.R., Auger, S. and Poirier, D. (1995) Synthesis of 16-(Bromoalkyl)-Estradiols Having Inhibitory Effect on Human Placental Estradiol 17β-Hydroxysteroid Dehydrogenase (17-β-HSD Type 1). Bioorganic & Medicinal Chemistry, 3, 505-523. https://doi.org/10.1016/0968-0896(95)00041-EPelletier, J.D. and Poirier, D. (1996) Synthesis and Evaluation of Estradiol Derivatives with 16α-(Bromoalkylamide), 16α-(Bromoalkyl) or 16α-(Bromoalkynyl) Side Chain as Inhibitors of 17β-Hydroxysteroid Dehydrogenase Type 1 without Estrogenic Activity. Bioorganic & Medicinal Chemistry, 4, 1617-1628. https://doi.org/10.1016/0968-0896(96)00154-XPoirier, D., Boivin, R.P., Tremblay, M.R., Bérubé, M., Qiu, W., et al. (2005) Estradiol-Adenosine Hybrid Compounds Designed to Inhibit Type 1 17β-Hydroxysteroid Dehydrogenase. Journal of Medicinal Chemistry, 48, 8134-8147. https://doi.org/10.1021/jm058235eKarthikeyan, C., Solomon, V.R., Lee, H. and Trivedi, P. (2013) Design, Synthesis and Biological Evaluation of Some Isatin-Linked Chalcones as Novel Anti-Breast Cancer Agents: A Molecular Hybridization Approach. Biomedicine & Preventive Nutrition, 3, 325-330. https://doi.org/10.1016/j.bionut.2013.04.001Abu Almaaty, A.H., Elgrahy, N.A., Fayad, E., Abu Ali, O.A., Mahdy, A.R.E., et al. (2021) Design, Synthesis and Anticancer Evaluation of Substituted Cinnamic Acid Bearing 2-Quinolone Hybrid Derivatives. Molecules, 26, Article No. 4724. https://doi.org/10.3390/molecules26164724Ghosh, D., Lo, J., Morton, D., Valette, D., Xi, J., et al. (2012) Novel Aromatase Inhibitors by Structure-Guided Design. Journal of Medicinal Chemistry, 55, 8464-8476. https://doi.org/10.1021/jm300930nBalewski, U., Szulta, S., Jalińska, A. and Kornicka, A. (2021) A Mini-Review: Recent Advances in Coumarin-Metal Complexes with Biological Properties. Frontiers in Chemistry, 9, Article 781779. https://doi.org/10.3389/fchem.2021.781779Luna-Guevara, M.L., Luna-Guevara, J.J., Hernández-Carranza, P., Ruíz-Espinosa, H. and Ochoa-Velasco, C.E. (2018) Phenolic Compounds: A Good Choice against Chronic Degenerative Diseases. Studies in Natural Products Chemistry, 59, 79-108. https://doi.org/10.1016/B978-0-444-64179-3.00003-7He, Q., Liu, J., Lan, J.S., Ding, J., Sun, Y., et al. (2018) Coumarin-Dithiocarbamate Hybrids as Novel Multitarget AChE and MAO-B Inhibitors against Alzheimer’s Disease: Design, Synthesis and Biological Evaluation. Bioorganic Chemistry, 81, 512-528. https://doi.org/10.1016/j.bioorg.2018.09.010Sharma, A., Tiwari, M. and Supuran, C.T. (2013) Novel Coumarins and Benzocoumarins Acting as Isoform-Selective Inhibitors against the Tumor-Associated Carbonic Anhydrase IX. Journal of Enzyme Inhibition and Medicinal Chemistry, 29, 292-296. https://doi.org/10.3109/14756366.2013.777334Touisni, N., Maresca, A., McDonald, P.C., Lou, Y., Scozzafava, A., et al. (2011) Glycosyl Coumarin Carbonic Anhydrase IX and XII Inhibitors Strongly Attenuate the Growth of Primary Breast Tumors. Journal of Medicinal Chemistry, 54, 8271-8277. https://doi.org/10.1021/jm200983eAhmad, I. and Shagufta (2015) Recent Developments in Steroidal and Nonsteroidal Aromatase Inhibitors for the Chemoprevention of Estrogen-Dependent Breast Cancer. European Journal of Medicinal Chemistry, 102, 375-386. https://doi.org/10.1016/j.ejmech.2015.08.010Selcer, K.W., Jagannathan, S., Rhodes, M.E. and Li, P.K. (1996) Inhibition of Placental Estrone Sulfatase Activity and MCF-7 Breast Cancer Cell Proliferation by Estrone-3-Amino Derivatives. The Journal of Steroid Biochemistry and Molecular Biology, 59, 83-91. https://doi.org/10.1016/S0960-0760(96)00089-1Recanatini, M., Bisi, A., Cavalli, A., Belluti, F., Gobbi, S., et al. (2001) A New Class of Nonsteroidal Aromatase Inhibitors: Design and Synthesis of Chromone and Xanthone Derivatives and Inhibition of the P450 Enzymes Aromatase and 17α-Hydroxylase/C17,20-Lyase. Journal of Medicinal Chemistry, 44, 672-680. https://doi.org/10.1021/jm000955sLeonetti, F., Favia, A., Rao, A., Aliano, R., Paluszcak, A., et al. (2004) Design, Synthesis, and 3D QSAR of Novel Potent and Selective Aromatase Inhibitors. Journal of Medicinal Chemistry, 47, 6792-6803. https://doi.org/10.1021/jm049535jGhodsi, R. and Hemmateenejad, B. (2016) QSAR Study of Diarylalkylimidazole and Diarylalkyltriazole Aromatase Inhibitors. Medicinal Chemistry Research, 25, 834-842. https://doi.org/10.1007/s00044-016-1530-1Chan, H.S., Shan, H., Dahoun, T., Vogel, H. and Yuan, S. (2019) Advancing Drug Discovery via Artificial Intelligence. Trends in Pharmacological Sciences, 40, 592-604. https://doi.org/10.1016/j.tips.2019.06.004Schlander, M., Hernandez-Villafuerte, K., Cheng, C.Y., Mestre-Ferrandiz, J. and Baumann, M. (2021) How Much Does It Cost to Research and Develop a New Drug? A Systematic Review and Assessment. PharmacoEconomics, 39, 1243-1269. https://doi.org/10.1007/s40273-021-01065-yMedina-Franco, J.L., Martinez-Mayorga, K., Fernández-de Gortari, E., Kirchmair, J. and Bajorath, J. (2021) Rationality over Fashion and Hype in Drug Design. F1000Research, 10, Article No. 397. https://doi.org/10.12688/f1000research.52676.1Prieto-Martínez, F.D., Arciniega, M. and Medina-Franco, J.L. (2018) Acoplamiento Molecular: Avances Recientes y Retos. TIP Revista Especializada en Ciencias Químico-Biológicas, 21, 65-87. https://doi.org/10.22201/fesz.23958723e.2018.0.143Naqvi, A.A.T., Mohammad, T., Hasan, G.M. and Hassan, M.I. (2018) Advancements in Docking and Molecular Dynamics Simulations towards Ligand-Receptor Interactions and Structure-function Relationships. Current Topics in Medicinal Chemistry, 18, 1755-1768. https://doi.org/10.2174/1568026618666181025114157Pinzi, L. and Rastelli, G. (2019) Molecular Docking: Shifting Paradigms in Drug Discovery. International Journal of Molecular Sciences, 20, Article No. 4331. https://doi.org/10.3390/ijms20184331Ghosh, D., Egbuta, C. and Lo, J. (2018) Testosterone Complex and Non-Steroidal Ligands of Human Aromatase. The Journal of Steroid Biochemistry and Molecular Biology, 181, 11-19. https://doi.org/10.1016/j.jsbmb.2018.02.009Bhuvaneswari, K., Sivaguru, P. and Lalitha, A. (2020) Synthesis, Anticancer Evaluation, and Docking Studies of Some Novel Azo Chromene Derivatives. Journal of the Chinese Chemical Society, 67, 1877-1886. https://doi.org/10.1002/jccs.201900481Li, J., Tang, Y., Li, W. and Tu, Y. (2020) Mechanistic Insights into the Regio- and Stereoselectivities of Testosterone and Dihydrotestosterone Hydroxylation Catalyzed by CYP3A4 and CYP19A1. Chemistry—A European Journal, 26, 6214-6223. https://doi.org/10.1002/chem.201905272Egli, M. and Sarkhel, S. (2006) Lone Pair-Aromatic Interactions: To Stabilize or Not to Stabilize. Accounts of Chemical Research, 40, 197-205. https://doi.org/10.1021/ar068174uHayakawa, D., Sawada, N., Watanabe, Y. and Gouda, H. (2020) A Molecular Interaction Field Describing Nonconventional Intermolecular Interactions and Its Application to Protein-Ligand Interaction Prediction. Journal of Molecular Graphics and Modelling, 96, Article ID: 107515. https://doi.org/10.1016/j.jmgm.2019.107515Breton, R., Housset, D., Mazza, C. and Fontecilla-Camps, J.C. (1996) The Structure of a Complex of Human 17β-Hydroxysteroid Dehydrogenase with Estradiol and NADP+ Identifies Two Principal Targets for the Design of Inhibitors. Structure, 4, 905-915. https://doi.org/10.1016/S0969-2126(96)00098-6Hernández-López, H., Leyva-Ramos, S., Azael Gómez-Durán, C.F., Pedraza-Alvarez, A., Rodríguez-Gutiérrez, I.R., et al. (2020) Synthesis of 1,4-Biphenyl-Triazole Derivatives as Possible 17β-HSD1 Inhibitors: An in Silico Study. ACS Omega, 5, 14061-14068. https://doi.org/10.1021/acsomega.0c01519Tremblay, M.R., Lin, S.X. and Poirier, D. (2001) Chemical Synthesis of 16β-Propylaminoacyl Derivatives of Estradiol and Their Inhibitory Potency on Type 1 17β-Hydroxysteroid Dehydrogenase and Binding Affinity on Steroid Receptors. Steroids, 66, 821-831. https://doi.org/10.1016/S0039-128X(01)00116-7Qiu, W., Campbell, R.L., Gangloff, A., Dupuis, P., Boivin, R.P., et al. (2002) A Concerted, Rational Design of Type 1 17β-Hydroxysteroid Dehydrogenase Inhibitors: Estradiol-Adenosine Hybrids with High Affinity. The FASEB Journal, 16, 1695-1624. https://doi.org/10.1096/fj.02-0026fjeSalah, M., Abdelsamie, A.S. and Frotscher, M. (2019) Inhibitors of 17β-Hydroxysteroid Dehydrogenase Type 1, 2 and 14: Structures, Biological Activities and Future Challenges. Molecular and Cellular Endocrinology, 489, 66-81. https://doi.org/10.1016/j.mce.2018.10.001Li, T., Stephen, P., Zhu, D., Shi, R. and Lin, S. (2019) Crystal Structures of Human 17β-Hydroxysteroid Dehydrogenase Type 1 Complexed with Estrone and NADP+ Reveal the Mechanism of Substrate Inhibition. The FEBS Journal, 286, 2155-2166. https://doi.org/10.1111/febs.14784Chen, C., Joo, J.C., Brown, G., Stolnikova, E., Halavaty, A.S., et al. (2014) Structure-Based Mutational Studies of Substrate Inhibition of Betaine Aldehyde Dehydrogenase BetB from Staphylococcus aureus. Applied and Environmental Microbiology, 80, 3992-4002. https://doi.org/10.1128/AEM.00215-14Efimov, I., Basran, J., Sun, X., Chauhan, N., Chapman, S.K., et al. (2012) The Mechanism of Substrate Inhibition in Human Indoleamine 2,3-Dioxygenase. Journal of the American Chemical Society, 134, 3034-3041. https://doi.org/10.1021/ja208694gMazumdar, M., Fournier, D., Zhu, D.W., Cadot, C., Poirier, D., et al. (2009) Binary and Ternary Crystal Structure Analyses of a Novel Inhibitor with 17β-HSD Type 1: A Lead Compound for Breast Cancer Therapy. Biochemical Journal, 424, 357-366. https://doi.org/10.1042/BJ20091020Ghosh, D. and Vihko, P. (2001) Molecular Mechanisms of Estrogen Recognition and 17-Keto Reduction by Human 17β-Hydroxysteroid Dehydrogenase 1. Chemico-Biological Interactions, 130, 637-650. https://doi.org/10.1016/S0009-2797(00)00255-6