Next Article in Journal
3D-Printing of Drug-Eluting Implants: An Overview of the Current Developments Described in the Literature
Previous Article in Journal
Chiroptical Sensing of Amino Acid Derivatives by Host–Guest Complexation with Cyclo[6]aramide
Previous Article in Special Issue
Insight into the Reaction of Alexidine with Sodium Hypochlorite: A Potential Error in Endodontic Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel 1,2,3-Triazole Derivatives as Mimics of Steroidal System—Synthesis, Crystal Structures Determination, Hirshfeld Surfaces Analysis and Molecular Docking

1
Department of Inorganic Chemistry, Faculty of Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
2
Department of Organic Chemistry, Faculty of Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(13), 4059; https://doi.org/10.3390/molecules26134059
Submission received: 11 June 2021 / Revised: 30 June 2021 / Accepted: 1 July 2021 / Published: 2 July 2021

Abstract

:
Herein, we present the synthesis and crystal structures determination of five 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives containing halogen atoms, 6ae, which may be used as an excellent mimic of steroids in the drug development process. Good quality crystals obtained for all of the synthesized compounds allowed the analysis of their molecular structures. Subsequently, the determined crystal structures were used to calculate the Hirshfeld surfaces for each of the synthesized compounds. Furthermore, results of our docking studies indicated that synthesized derivatives are able to bind effectively to the active sites of selected enzymes and receptors involved in the hormone biosynthesis and signaling pathways, analogously to the native steroids.

Graphical Abstract

1. Introduction

Triazoles are a class of compounds showing very interesting properties, e.g., hydrogen bond (HB) formation, π–π stacking interaction, large dipole moments, bioisosteric effects, and therefore, they have been successfully used as scaffolds in the synthesis of antimicrobial, antiviral, and antitumor agents [1]. Importantly, they do not undergo hydrolysis under acidic or basic conditions, and they withstand metabolic degradation, which is desired in the design of new pharmaceuticals. The above advantages indicate that the derivatization of the triazole ring may lead to compounds possessing interesting biological properties. For example, the introduction of two additional phenyl rings with diverse substituents to the 1,2,3-triazole ring should allow the production of 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives demonstrating similarities to some natural compounds, e.g., hormones.
The hormone signaling pathway is a well-established target for the development of hormone-dependent cancer drugs (e.g., breast cancer) [2]. For example, one of the clinically used drugs—Tamoxifen 1 (Figure 1) acts as a selective estrogen receptor modulator (SERM). On the other hand, chemotherapeutics, which may influence the hormone formation process, are of high therapeutic importance. The biosynthesis of active steroids (e.g., estradiol (E2) and androstenediol (Adiol)) in cancer tissues mainly depends on three enzymatic pathways: aromatase (AROM), 17β-hydroxysteroid dehydrogenase (17β-HSD) and steroid sulfatase (STS) [3]. For example, currently used Letrozole 2 and Anastrozole 3 (Figure 1) block the conversion of androgens to estrogens via the inhibition of the AROM complex. In light of recent research indicating that the disorders in sulfation/desulfation processes may be responsible for numerous pathologies [4], the other enzyme implicated in the steroidogenesis process—STS—is becoming a new, interesting molecular target in the development of novel and effective hormone-dependent cancer treatment methods. Recently, we have developed a series of tricyclic compounds, featuring the 1,2,3-triazole unit, as potent STS inhibitors [5,6]. In the course of our investigation, we have found that the most active analog, MD77 4 (Figure 1), inhibited the STS enzyme with an IC50 value of 36 nM when evaluated in an enzymatic assay. Our studies indicated that the highest inhibitory activities were exhibited by derivatives containing a fluorine atom at the meta position of the terminal aromatic ring. According to the molecular docking calculations, it was noticed that the fluorine atoms presented in MD77 may interact with the Arg98 residue located in the STS active site. This additional interaction may stabilize the inhibitor–enzyme complex, resulting in improved inhibitory activity [5].
In the present paper, we have reported the crystal structures and Hirshfeld surfaces determination of five 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives. The crystal structures provide insights into the possible conformations and the character of the intermolecular interactions of different derivatives. Additionally, we have described preliminary docking studies indicating effective binding of the synthesized compounds to the active sites of AROM, STS, 17β-HSD, and estrogen receptors (ER). The collected data indicate that the reported triazole derivatives could be utilized in the development of novel inhibitors of proteins involved in the hormone biosynthesis and signaling pathways and thus become a perfect starting point in the development of novel anticancer agents.

2. Results and Discussion

2.1. Chemistry

Three compounds 6ac were resynthesized according to the two-step synthetic protocol that we previously described (Scheme 1) [5]. Two novel compounds 6d and 6e were synthesized according to the same synthetic pathway using chlorinated aniline derivatives. In the first step, 4-[(trimethylsilyl)ethynyl]phenol 5 was obtained by the Sonogashira reaction between the 4-iodophenol and ethynyltrimethylsilane in the presence of palladium(II) chloride (PdCl2), triphenyl phosphine (Ph3P), copper(I) iodide (CuI) and triethylamine (NEt3). The Sonogashira coupling is a widely used cross-coupling reaction applied in organic practice to generate carbon–carbon bonds [7]. Next, the appropriate aniline derivatives were transformed into the corresponding azides with tert-butyl nitrite (t-BuONO) and azidotrimethylsilane (TMSN3). Then, to the obtained solution, 4-[(trimethylsilyl)ethynyl]phenol 5, a 1 M solution of tetrabutylammonium fluoride (TBAF) in tetrahydrofuran (THF), copper(II) sulfate pentahydrate (CuSO4·5H2O) and a 1 M aqueous solution of sodium ascorbate were added. After a brief work-up of the reaction mixture, the corresponding 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives 6ae were isolated. Recrystallization in acetonitrile (ACN) (6ce), methanol (6b), or acetone (6a) allowed to obtain crystals of compounds suitable for X-ray diffraction measurement.

2.2. Crystal Structures and Hirshfeld Surfaces Analysis

Experimental X-ray crystal structures of 6ae may be helpful in the identification of the intermolecular interactions that can arise between the studied molecules and their environment. Good quality crystals obtained for all of the compounds allowed the analysis of their molecular structures (for the parameters of crystallographic data, see Table 3 in the Materials and Methods Section). The asymmetric units contained either one or two molecules of 6ae; typical examples and numbering schemes are illustrated in Figure 2a–c for compounds 6a, 6c, and 6d, respectively (the molecular structures of the remaining two compounds can be found as Figures S1 and S2 in ESI). Bond lengths between the atoms, which are collected in Table S1 in ESI, were very similar in all molecules and, in our opinion, did not depend much on the substitution pattern; the mean values of C–C bond lengths in the phenol rings were very close when compared between the derivatives ranging from 1.390 to 1.393 Å. The mean values of C–C bonds in the halogenated rings were invariably shorter and varied from 1.380 Å to 1.387 Å (Table S1, ESI).
The substituted phenyl rings in 6ae can adopt different mutual orientations. The torsion angles between the six-membered rings varied from 0.90° in 6e to 50.41° in 6a. It can be observed that the same molecule can exhibit different torsion angles within the same crystal—as in 6a and 6e (Table S1, ESI). Moreover, the mutual rotation of the rings may introduce a helical type of chirality, and we do observe two helical enantiomers of 6a in crystals. In the examples presented in Figure 2 for 6a, 6c, and 6d compounds, atoms of halogens are asymmetrically attached to the phenyl ring (with regard to triazole substitution), which leads to two possibilities of their orientation with regard to the triazole ring. The mono-substituted rings in 6a and 6d realize common conformation with the torsions N2-N3-C11-F1/Cl1 17.46°/20.32°, respectively, whereas in 6c, the additional intramolecular contact between H8 and F1 (2.415 Å) may impose the opposite arrangement with the torsion N2-N3-C11-F2 equal to 154.62° (Figure 2).
Geometrical parameters of intermolecular interactions characteristic for compounds 6ae are collected in Table S2 (ESI) and the crystal packings are illustrated in Figure 3a,b and in Figures S3–S5 (ESI). The most typical interactions between the molecules included classical hydrogen bonds, which usually formed between the phenol group of one molecule and N1/N2 atoms of the triazole ring of the adjacent molecule (Figure 3a). The only exception to this pattern was compound 6d—in the crystal packing of 6d the OH⋯O supported chains are observed, and such change probably allowed the formation of additional halogen bonding interaction, as indicated in Figure 3b and Table S2 (ESI). The intermolecular interactions of derivatives containing fluorine atoms always include CH⋯F contacts; the shortest one, which is observed in 6c, exhibits the donor-acceptor distance below 3 Å, implicating quite strong interaction. As shown in Figure 3a, the intermolecular forces are usually strongest in the approximate plane of each molecule, which leads to the formation of hydrogen-bonded 2-D arrangements.
Based on the determined crystal structures of compounds 6ae, the Hirshfeld surfaces (HSs) were mapped with a dnorm function, which are illustrated in Figure 4 with fingerprint decomposition. The calculated HSs of molecules 6ae indicated the presence of strong O–H⋯N contacts between molecules, which are represented by the red spots shown in Figure 4. Moreover, the blue regions correspond to weak interactions, such as C–H⋯H contacts. According to the decomposed fingerprint plots of compounds 6ae, it was noticed that the most important interactions between molecules are van der Waals forces. The diagram of 6a shows that the C⋯H (33.4%) bonds are major factors in the crystal packing with hydrophobic H⋯H (24.7%) interactions, which are the next highest contribution. The decomposed fingerprint plots confirmed that with the increase in the number of fluorine atoms in the molecule of 6b and 6c, the F⋯X interactions become very important (28.4% in 6b) or even prevailing (38.7% in 6c) intermolecular forces. Based on the fingerprint plots of 6d, it was detected that the major factors in the crystal packing are hydrophobic H⋯H (22.2%) and C⋯H (34.7%) bonds. As in 6a6c series, with the increase in the number of halogen atoms, their intermolecular interactions grow in strength and become an important structure-building factor. Accordingly, in 6e, Cl⋯X (30.1%) becomes comparable with the hydrophobic H⋯H plus C⋯H (35.1%).

2.3. Computational Studies

2.3.1. The Lipinski’s Rule of Five Calculations

The free access web tool SwissADME server (Swiss Institute of Bioinformatics, Lausanne, Switzerland) was used to predict drug-like physicochemical (PC) properties of compounds 6ae based on Lipinski’s Rule of Five (molecular weight less than 500, log P or coefficient partition between −5 and 5, less than five HB donors, and less than ten HB acceptors) [8]. Our calculations indicated that compounds 6ae demonstrate desired drug-like PC features (the collected data are summarized in Table 1).

2.3.2. Molecular Docking

Our previous research indicated that some of the sulfamoylated derivatives of compounds 6ae (e.g., MD77 4) demonstrated very high STS inhibitory properties, and therefore, they might be recognized as drug candidates in the treatment of hormone-dependent types of cancers. In the present studies, we performed molecular docking calculations to the active sites of several molecular targets (AROM, STS, 17β-HSD1, ERα, and ERβ) for compounds 6ae using AutoDock Vina 1.1.2 software (Molecular Graphics Laboratory, The Scripps Research Institute, LaJolla, CA, USA). The obtained data indicated that the 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol core may be useful in the development of novel inhibitors of the abovementioned proteins. The summarized results of docking calculations are presented in Table 2. Subsequently, the visualizations of the examples of docked compounds using VMD 1.9 software (University of Illinois at Urbana—Champaign, Urbana, IL, USA) and detailed identification of plausible interactions using BIOVIA software (Dassault Systémes, Discovery Studio Visualiser, San Diego, CA, USA) were performed.

AROM

The calculated binding free energies of compounds 6ad to AROM were comparable and in the range of −7.5 to −8.3 kcal mol−1 (Table 2). Only the binding free energy for compound 6e was significantly higher (−5.5 kcal mol−1). In general, the binding free energies calculated for fluorinated derivatives were slightly better than their chlorinated analogs. However, all of the calculated values were less favorable than the binding free energy of androstenedione (−12.4 kcal mol−1), indicating a slightly worse match to the enzyme’s active site than the natural substrate. Analysis of docking modes of compounds 6ae indicated that triazole rings are close to the heme Fe2+ (e.g., 4.14 Å for compound 6b, Figure 5) analogously as a methyl group of androstenedione (3.94 Å). Furthermore, the halogenated rings of compounds 6ae occupied the same region of the enzyme’s binding site like a five-membered ring of androstenedione. For example, the distance between one of the fluorine atoms of compound 6b and Asp309 was 2.36 Å (in the case of androstenedione, the distance between the carbonyl group and Asp309 was 3.02 Å). A more detailed list of the plausible interactions between compounds 6ae and the AROM enzyme was obtained using BIOVIA and presented in Table S3 (ESI).

STS

Initially, the docking calculations to the STS protein were performed for compounds 6ae and E1 (used as a reference). The calculated binding free energies of compounds 6ae were in the range of −8.1 to −8.3 kcal mol−1 (Table 2) and were slightly less favorable than the binding free energy value of E1 (−8.9 kcal mol−1). There were no significant differences between binding free energy values for fluorinated and chlorinated analogs. Analysis of docking modes of compounds 6ae and E1 indicated their similar binding manner to the STS active site. As it is presented in Figure 6A for representative 6e derivative, we found that the –OH groups of compound 6e and E1 are in a short distance to the catalytic amino acid residue fGly75 (3.13 and 2.72 Å, respectively) coordinated to the Ca2+ ion. Furthermore, we detected that the halogen atoms of compounds 6ae occupied the same region of the STS active site as the carbonyl group of E1, indicating the presence of possible interactions with Arg98.
However, taking into consideration the mechanism of the enzymatic reaction catalyzed by the STS, we prepared and docked a set of sulfated derivatives of compounds 6ae. Their calculated binding free energies were in the range of −6.6 to −7.8 kcal mol−1 (Table 2) and were more favorable than the binding free energy value of E1S (−6.3 kcal mol−1), used as a reference. Figure 6B shows the structure of the sulfated derivative of 6c and E1S docked into the active site of STS. We found that the sulfate groups of the sulfated derivative of 6c and E1S were located very close to the fGly75 residue (3.32 and 2.99 Å, respectively). Analogously, as it was detected in the case of phenolic derivatives 6ae, the halogen atoms of their sulfated derivatives were also in close distances to the Arg98 residue indicating the presence of interactions, which may stabilize the enzyme–ligand complexes.
Importantly, more favorable binding free energy values were calculated for sulfates of compounds 6ae in comparison with E1S, indicating better binding of ligands based on the 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol core before catalysis than the natural substrate. On the other hand, less favorable binding free energies were calculated for phenolic derivatives 6ae in comparison with E1, indicating easier dissociation of the compounds after catalysis. Both of the mentioned observations suggest that the enzymatic reaction for compounds containing 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol core may be more effective in comparison with the hydrolysis of E1S to E1. A more detailed list of the plausible interactions between compounds 6ae in their phenolic and sulfate forms and STS enzyme was obtained using BIOVIA and presented in Table S4 (ESI).

17β-HSD1

The calculated binding free energies of compounds 6ae to 17β-HSD1 were comparable (in the range of −7.9 to −8.5 kcal mol−1 Table 2) and slightly higher than the binding free energy value of E1 (−8.9 kcal mol−1). Fluorinated and chlorinated derivatives both demonstrated similar binding free energy values. Compounds 6ae docked to the 17β-HSD1’s active site analogously as it was determined for E1 (the halogenated rings of compound 6ae occupied the same region as a five-membered ring of E1). As it is presented for compound 6c (Figure 7), we found that the fluorine atoms present in the meta and para positions of the terminal aromatic ring were in close distances to the Ser142 and Tyr155 amino acid residues (the distances of meta-F–HO-Ser142, meta-F–HO-Tyr155 and para-F–HO-Tyr155 were 3.13, 2.95 and 3.31 Å, respectively). In comparison, the distances between the carbonyl group of the E1 and Ser142 and Tyr155 were 5.29 and 4.45 Å, respectively. The –OH groups of compound 6c and E1 occupied strictly the same region of the enzyme’s active site and were in short distances to the His221 and Gly282 amino acid residues (2.97 and 3.06 Å for 6c). All plausible interactions between compounds 6ae and the 17β-HSD1 enzyme were obtained using BIOVIA and collected in Table S5 (ESI).

ERα and ERβ

The calculated binding free energies of compounds 6ae to ERα and ERβ were similar (in the range of −7.9 to −8.8 kcal mol−1 and −7.6 to −8.9 kcal mol−1, respectively, Table 2), suggesting their effective association to the ERα and ERβ binding sites (the binding free energy values of E2 were −10.7 and −11.1 kcal mol−1, respectively). In both cases, the binding free energy values calculated for fluorinated derivatives were slightly better than their chlorinated analogs; however, in our opinion, the differences are too low to indicate them as significant. Compounds 6ae docked to the ERα’s and ERβ’s binding sites in a similar manner as reported for E2. As it is presented in Figure 8A, we found that the –OH groups of compound 6c and E2 were in short distances to the Glu353 and Arg394 amino acid residues of ERα (2.70 and 2.80 Å for 6c, respectively; 2.92 and 3.21 Å for E2, respectively). On the opposite side of the binding region of ERα, we detected that the halogens atoms of compounds 6ae occupied the same region as the -OH group of E2, indicating the presence of interactions with His524. For example, the distance between the para-substituted fluorine atom of 6c and the nitrogen atom of His524 was 2.72 Å (analogously, the distance between the –OH group of E2 and His524 was 2.93 Å). In case of docking to ERβ’s binding site (Figure 8B), we detected that –OH groups of compounds 6ae were in short distances to the Glu305 and Arg346 amino acid residues analogously to E2 (for example, 2.90 and 3.63 Å for 6c, respectively; 2.93 and 2.91 Å for E2, respectively). Furthermore, halogen atoms of compounds 6ae occupied the same region of ERβ like –OH group of E2. The para-substituted fluorine atom of compound 6c was at a distance of 3.05 Å to the nitrogen atom of His475 (for comparison, the distance between the –OH group of E2 and His475 was 3.17 Å). Tables S6 and S7 (ESI) summarizes all of the plausible interactions between compounds 6ae and both ERs detected using BIOVIA.

3. Materials and Methods

3.1. Synthesis

4-iodophenol, trimethylsilylacetylene, PdCl2, PPh3, CuI, NEt3, 1 M solution of TBAF in THF, sodium ascorbate, CuSO4·5H2O, t-BuONO, TMSN3, all of the used aniline derivatives, and solvents are commercially available from Merck (Merck KGaA, Darmstadt, Germany). Solvents were dried and distilled using standard procedures. Melting points (uncorrected) were determined with a Stuart Scientific SMP30 apparatus (Stuart, Stone, UK). NMR spectra were recorded on a Bruker Avance III HD 400 MHz spectrometer (Bruker, Billerica, MA, USA). Chemical shifts are reported in ppm relative to the residue solvent peak (DMSO-d6 2.49 ppm for 1H, 39.5 ppm for 13C). Coupling constants are given in Hertz. IR spectra were recorded on a Nicolet 8700 spectrometer (Thermo Fisher, Waltham, MA, USA). Mass spectra were recorded on an Agilent 6540 Accurate Mass Q-TOF LC/MS System (Agilent, Santa Clara, CA, USA). Column chromatography was performed using silica gel 60 (230-400 mesh, Merck KGaA, Darmstadt, Germany). Elemental analysis was performed using a CHNS-Carlo Erba EA-1108. Preparative thin-layer chromatography was performed with Polygram SIL G/UV254 silica gel (Macherey–Nagel GmbH and Co. KG, Düren, Germany). The detailed procedure for the synthesis of compounds 6ae was previously described (with the characterization of compounds 6a-c) [5]. Recrystallization in acetonitrile (ACN) (6ce), methanol (6b), or acetone (6a) allowed to obtain crystals of compounds suitable for X-ray diffraction measurement.
4-[1-(3-chlorophenyl)-1H-1,2,3-triazol-4-yl]-phenol 6d. Yield 70%; mp 208–209 °C; νmax (KBr)/cm−1 3458, 1616, 1591, 1466, 1222, 1175, 1040, 839, 681; 1H NMR δH (400 MHz, DMSO) 9.70 (1H, s, OH), 9.19 (1H, s, CH), 8.07 (1H, t, J 2.0 Hz, Ar-H), 7.98–7.93 (1H, m, Ar-H), 7.75 (2H, d, J 8.7 Hz, Ar-H), 7.66 (1H, t, J 8.1 Hz, Ar-H), 7.60–7.55 (1H, m, Ar-H), 6.89 (2H, d, J 8.7 Hz, Ar-H); 13C NMR δC (101 MHz, DMSO) 158.2, 148.3, 138.2, 134.7, 132.1, 128.8, 127.3, 121.4, 120.1, 118.9, 118.7, 116.2. Anal. calcd for: C14H10ClN3O: C, 61.89; H 3.71; N, 15.47. Found: C, 61.97; H, 3.60; N, 15.51%. HRMS (m/z) [M−H] calcd 270.0434, found 270.0547.
4-[1-(3,5-dichlorophenyl)-1H-1,2,3-triazol-4-yl]-phenol 6e. Yield 54%; mp 241–244 °C; νmax (KBr)/cm−1 3126, 1614, 1591, 1471, 1226, 1177, 1057, 841, 662; 1H NMR δH (400 MHz, DMSO) 9.72 (1H, s, OH), 9.24 (1H, s, CH), 8.08 (2H, d, J 1.8 Hz, Ar-H), 7.76 (1H, t, J 1.8 Hz, Ar-H), 7.73 (2H, d, J 8.7 Hz, Ar-H), 6.89 (2H, d, J 8.7 Hz, Ar-H); 13C NMR δC (101 MHz, DMSO) 158.3, 148.4, 138.8, 135.7, 128.3, 127.3, 121.2, 118.9, 116.3. Anal. calcd for: C14H9Cl2N3O: C, 54.92; H 2.96; N, 13.73. Found: C, 54.85; H, 2.91; N, 13.86%. HRMS (m/z) [M−H] calcd 304.0044, found 304.0156.

3.2. X-ray Diffraction Measurement

Single crystal X-ray diffraction data of the compounds 6ae were collected at 120(2) K on a Stoe IPDS-2T diffractometer with graphite-monochromated Mo-Kα radiation. Crystals were cooled using a Cryostream 800 open flow nitrogen cryostat (Oxford Cryosystems, Long Hanborough, Oxford, UK). Data collection and image processing were performed with X-Area 1.75 (STOE and Cie Gmbh, Darmstadt, Germany) [9]. Intensity data were scaled with LANA (part of X-Area) in order to minimize differences of intensities of symmetry-equivalent reflections (multi-scan method). Structures were solved by direct methods, and all non-hydrogen atoms were refined with anisotropic displacement parameters by a full-matrix least squares procedure based on F2 using the SHELX–2014 program package [10,11]. The Olex [12] and Wingx [13] suites were used to prepare the final version of the CIF files. Mercury [14] was used to prepare the figures and to calculate the planes of six-membered rings. Hydrogen atoms were refined using an isotropic model with Uiso(H) values fixed to be 1.2 times Ueq for –CH groups. The hydrogens of OH groups were refined freely. A summary of crystallographic data is shown in Table 3.
CCDC 2,063,930–2,063,934 contain the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre.
The Hirshfeld surfaces and 2D fingerprint plots were generated using Crystal Explorer 17.5. Crystal structures were imported from CIF files. Hirshfeld surfaces were received using a high surface resolution and mapped with the dnorm function.

3.3. Computational Studies

3.3.1. The Lipinski’s Rule of Five Calculations

Calculation of the standard properties of compound 6ae based on Lipinski’s Rule of Five was performed using the SwissADME server (Swiss Institute of Bioinformatics, Lausanne, Switzerland).

3.3.2. Ligands Preparation for Molecular Docking

The 3D structure of ligands 6ae, their sulfated analogs, and reference compounds (androstenedione, E1S, E1, and E2) were prepared with the Portable HyperChem 8.0.7 Release (Hypercube, Inc., Gainesville, FL, USA). Prior to docking calculations, the structure of each ligand was optimized using an MM+ force field and the Polak–Ribière conjugate gradient algorithm (terminating at a gradient of 0.05 kcal mol−1 Å−1).

3.3.3. Protein Preparation for Molecular Docking

The X-ray structures of the AROM, STS, 17β-HSD1, ERα and ERβ used for molecular modeling studies were taken from the Protein Databank (Protein Data Bank accession codes: 3EQM, 1P49, 6MNC, 1A52, and 5TOA, respectively). After standard preparation procedures (including removal of water molecules and other ligands as well as addition hydrogen atoms and Gasteiger charges to each atom), docking analysis was carried out.

3.3.4. Molecular Docking

Docking studies were carried out using Autodock Vina 1.1.2 software (The Molecular Graphic Laboratory, The Scripps Research Institute, La Jolla, CA, USA) [15]. For the docking studies, the corresponding grid box parameters were used:
- AROM: a grid box size of 24 Å × 24 Å × 24 Å centered on Asp309 amino acid residue (x = 88.230, y = 49.522, z = 51.205);
- STS: a grid box size of 24 Å × 24 Å × 24 Å centered on the fGly75 amino acid residue (x = 72.135, y = −1.720, z = 28.464);
- 17β-HSD1: a grid box size of 24 Å × 24 Å × 24 Å centered on Ser142 amino acid residue (x = 20.760, y = −6.625, z = −22.837);
- ERα: a grid box size of 24 Å × 24 Å × 24 Å centered on Glu353 (x = 105.702, y = 19.417, z = 103.747);
- ERβ: a grid box size of 24 Å × 24 Å × 24 Å centered on Met336 (x = 16.635, y = 41.368, z = 18.491).
Graphic visualizations of the 3D model were generated using VMD 1.9 software (University of Illinois at Urbana—Champaign, Urbana, IL, USA). Identification of the ligand–protein interactions was performed using Discovery Studio Visualiser v20. 1. 0. 19,295 (BIOVIA, Dassault Systémes, San Diego, CA, USA).

4. Conclusions

In summary, we synthesized five 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives containing halogen atoms, 6ae. Single crystals of such compounds were obtained, and their crystal structures were successfully confirmed. The determined crystal structures were used to calculate the Hirshfeld surfaces for each of the synthesized compounds, which indicated that the structures of all compounds are mainly characterized by H⋯H, C⋯H, and X⋯H (where X is F or Cl) interactions. The decomposed fingerprint plots confirmed that with the increase in the number of halogen atoms in the molecules, the interactions involving fluorine or chlorine atoms become very important or even prevail the intermolecular structure-building factor. Furthermore, the performed docking studies proved potentially effective binding of compounds based on a 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol core to the active sites of several molecular targets (AROM, STS, 17β-HSD1, ERα, and ERβ). Importantly, we also noticed that the halogen atoms of compounds 6ae play an important role in the putative binding process, and their presence is crucial for the creation of stabilizing interactions in the active sites of examined molecular targets. This finding corresponds well with the crystallographic data and Hirshfeld surfaces analysis, which indicated that the interactions involving halogen atoms are a very important structure-building factor.
Our research indicated that the studied compounds might be a great starting point in the development of inhibitors of enzymes involved in the hormone biosynthesis and signaling pathways. Importantly, the calculated binding free energies for sulfates of compounds 6ae were better than the binding free energy of the STS native substrate (E1S). As we mentioned, we previously reported that sulfamoylated derivatives of halogenated 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenols are potent STS inhibitors, which proves our present findings. It is also worth noting that the binding free energies for compounds 6ae were comparable to the binding free energy of 17β-HSD1 native substrate (E1), indicating that such compounds may also potentially demonstrate 17β-HSD1 inhibitory properties. However, such hypothesis needs to be evaluated in further biological experiments. Additionally, our docking studies and calculations of standard PC properties indicated that compounds 6ae demonstrate favorable drug-like features and may be potentially used as drugs. However, due to the limited accuracy of the utilized AutoDock Vina, further computational studies using more accurate methods (e.g., free energy perturbation calculations) should be performed for finally designed compounds.

Supplementary Materials

The following are available online. Figure S1: Molecular structure of 6b. Thermal ellipsoids at 50%, Figure S2: Molecular structures of 6e. Thermal ellipsoids at 50%, Figure S3: Crystal packings and intermolecular interactions in 6b., Figure S4: Crystal packings and intermolecular interactions in 6c., Figure S5: Crystal packings and intermolecular interactions in 6e, Table S1: Bond lengths and angles for 6ae, Table S2: Structural parameters for intramolecular interactions in 6ae, Table S3: The ligand-protein (AROM) interactions (and distances [Å]) identified using BIOVIA, Dassault Systémes, Discovery Studio Visualiser, Table S4: The ligand-protein (STS) interactions (and distances [Å]) identified using BIOVIA, Dassault Systémes, Discovery Studio Visualiser, Table S5: The ligand-protein (17β-HSD1) interactions (and distances [Å]) identified using BIOVIA, Dassault Systémes, Discovery Studio Visualiser, Table S6: The ligand-protein (ERα) interactions (and distances [Å]) identified using BIOVIA, Dassault Systémes, Discovery Studio Visualiser, Table S7: The ligand-protein (ERβ) interactions (and distances [Å]) identified using BIOVIA, Dassault Systémes, Discovery Studio Visualiser.

Author Contributions

Conceptualization, methodology, writing—original draft preparation, funding acquisition, M.D.; investigation, M.D., A.D., M.S., K.B., O.C., S.D.; writing—review and editing, M.D., A.D., K.B., O.C., S.D., supervision, J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Centre for Research and Development (Poland), the TANGO program, grant agreement number TANGO-IV-A/0004/2019-00.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds 6ae are available from the authors.

Abbreviations

17β-HSD17β-hydroxysteroid dehydrogenase
ACNacetonitrile
Adiolandrostenediol
AROMaromatase
CuSO4·5H2Ocopper(II) sulfate pentahydrate
DHESdehydroepiandrosterone
DHEASdehydroepiandrosterone sulfate
E1estrone
E2estradiol
E1Sestrone sulfate
ERestrogen receptor
HBhydrogen bond
HSHirshfeld surface
PCphysicochemical
SERMselective estrogen receptor modulator
STSsteroid sulfatase
TBAFtetrabutylammonium fluoride
t-BuONOtert-butyl nitrite
THFtetrahydrofuran
TMSN3azidotrimethylsilane

References

  1. Bonandi, E.; Christodoulou, M.S.; Fumagalli, G.; Perdicchia, D.; Rastelli, G.; Passarella, D. The 1,2,3-triazole ring as a bioisostere in medicinal chemistry. Drug Discov. Today 2017, 22, 1572–1581. [Google Scholar] [CrossRef] [PubMed]
  2. Saha, T.; Makar, S.; Swetha, R.; Gutti, G.; Singh, S.K. Estrogen signaling: An emanating therapeutic target for breast cancer treatment. Eur. J. Med. Chem. 2019, 177, 116–143. [Google Scholar] [CrossRef] [PubMed]
  3. Shah, R.; Singh, J.; Singh, D.; Singh Jaggi, A.; Singh, N. Sulfatase inhibitors for recidivist breast cancer treatment: A chemical review. Eur. J. Med. Chem. 2016, 114, 170–190. [Google Scholar] [CrossRef] [PubMed]
  4. Muller, J.W.; Gilligan, L.C.; Idkowial, J.; Arlt, W.; Foster, P.A. The regulation of steroid action by sulfation and desulfation. Endocr. Rev. 2015, 36, 526–563. [Google Scholar] [CrossRef] [PubMed]
  5. Daśko, M.; Demkowicz, S.; Rachon, J.; Biernacki, K.; Aszyk, J.; Kozak, W.; Masłyk, M.; Kubiński, K. New potent STS inhibitors based on fluorinated 4-(1-phenyl-1H-[1,2,3]triazol-4-yl)-phenyl sulfamates. J. Asian Nat. Prod. Res. 2019, 22, 1037–1044. [Google Scholar] [CrossRef] [PubMed]
  6. Ciupak, O.; Daśko, M.; Biernacki, K.; Rachon, J.; Masłyk, M.; Kubiński, K.; Martyna, A.; Demkowicz, S. New potent steroid sulphatase inhibitors based on 6-(1-phenyl-1H-1,2,3-triazol-4-yl)naphthalen-2-yl sulphamate derivatives. Enzym. Inhib. Med. Chem. 2021, 36, 238–247. [Google Scholar] [CrossRef] [PubMed]
  7. Sonogashira, K. Development of Pd-Cu catalyzed cross-coupling of terminal acetylenes with sp2-carbon halides. J. Organomet. Chem. 2002, 653, 46–49. [Google Scholar] [CrossRef]
  8. Lipinski, C.A. Lead- and drug-like compounds. The rule-of-five revolution. Drug Discov. Today Technol. 2004, 1, 337–341. [Google Scholar] [CrossRef] [PubMed]
  9. X-Area 1.75, Software Package for Collecting Single-Crystal Data on STOE Area-Detector Diffractometers, for Image Processing, Scaling Reflection Intensities and for Outlier Rejection; STOE and Cie GmbH: Darmstadt, Germany, 2015.
  10. Sheldrick, G.M. Crystal structure refinement with SHELXL. Acta Crystallogr. Sect. C Struct. Chem. 2015, 71, 3–8. [Google Scholar] [CrossRef] [PubMed]
  11. Sheldrick, G.M. SHELXL-2014; University of Göttingen and Bruker AXS: Karlsruhe, Germany, 2014. [Google Scholar]
  12. Dolomanov, O.V.; Bourhis, L.J.; Gildea, R.J.; Howard, J.A.K.; Puschmann, H. OLEX2: A complete structure solution, refirement and analysis program. J. Appl. Cryst. 2009, 42, 339–341. [Google Scholar] [CrossRef]
  13. Farrugia, L.J. WinGX and ORTEP for Windows: An update. J. Appl. Cryst. 2012, 45, 849–854. [Google Scholar] [CrossRef]
  14. Macrae, C.F.; Edgington, P.R.; McCabe, P.; Pidcock, E.; Shields, G.P.; Taylor, R.; Towler, M.; van de Streek, J. Mercury: Visualization and analysis of crystal structures. J. Appl. Cryst. 2006, 39, 453–457. [Google Scholar] [CrossRef] [Green Version]
  15. 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. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Chemical structures of Tamoxifen 1, Letrozole 2, Anastrozole 3, and MD77 4.
Figure 1. Chemical structures of Tamoxifen 1, Letrozole 2, Anastrozole 3, and MD77 4.
Molecules 26 04059 g001
Scheme 1. Synthesis of 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives 6ae: (a) ethynyltrimethylsilane, PdCl2, Ph3P, CuI, NEt3, ACN; (b) t-BuONO, TMSN3, ACN; (c) TBAF, sodium ascorbate, CuSO4·5H2O.
Scheme 1. Synthesis of 4-(1-phenyl-1H-1,2,3-triazol-4-yl)phenol derivatives 6ae: (a) ethynyltrimethylsilane, PdCl2, Ph3P, CuI, NEt3, ACN; (b) t-BuONO, TMSN3, ACN; (c) TBAF, sodium ascorbate, CuSO4·5H2O.
Molecules 26 04059 sch001
Figure 2. Molecular structures of (a) 6a; (b) 6c; (c) 6d. Thermal ellipsoids at 50%.
Figure 2. Molecular structures of (a) 6a; (b) 6c; (c) 6d. Thermal ellipsoids at 50%.
Molecules 26 04059 g002
Figure 3. Crystal packings and intermolecular interactions in (a) 6a; (b) 6d.
Figure 3. Crystal packings and intermolecular interactions in (a) 6a; (b) 6d.
Molecules 26 04059 g003
Figure 4. Hirshfeld surfaces of (a) 6a, (b) 6b, (c) 6c, (d) 6d, (e) 6e. Red color: normalized contact distances for 6ae: dnorm shorter than the sum of van der Waals radii (for 6a: dnorm = −0.6166; 6b: dnorm = −0.5728; 6c: dnorm = −0.5901; 6d: dnorm = −0.3820; 6e: dnorm = −0.6177), white color: van der Waals contacts (for 6a: dnorm = 0.4140; 6b: dnorm = 0.3931; 6c: dnorm = 0.3738; 6d: dnorm = 0.3266; 6e: dnorm = 0.4040), blue color: normalized contact distances exceeding the sum of van der Waals radii (for 6a: dnorm = 1.3683; 6b: dnorm = 1.2271; 6c: dnorm = 1.1824; 6d: dnorm = 1.0046; 6e: dnorm = 1.1712); (fj), Hirshfeld surface fingerprint decompositions showing main types of interactions for compounds 6ae.
Figure 4. Hirshfeld surfaces of (a) 6a, (b) 6b, (c) 6c, (d) 6d, (e) 6e. Red color: normalized contact distances for 6ae: dnorm shorter than the sum of van der Waals radii (for 6a: dnorm = −0.6166; 6b: dnorm = −0.5728; 6c: dnorm = −0.5901; 6d: dnorm = −0.3820; 6e: dnorm = −0.6177), white color: van der Waals contacts (for 6a: dnorm = 0.4140; 6b: dnorm = 0.3931; 6c: dnorm = 0.3738; 6d: dnorm = 0.3266; 6e: dnorm = 0.4040), blue color: normalized contact distances exceeding the sum of van der Waals radii (for 6a: dnorm = 1.3683; 6b: dnorm = 1.2271; 6c: dnorm = 1.1824; 6d: dnorm = 1.0046; 6e: dnorm = 1.1712); (fj), Hirshfeld surface fingerprint decompositions showing main types of interactions for compounds 6ae.
Molecules 26 04059 g004
Figure 5. The binding mode of compound 6b (red) and androstenedione (blue) in the AROM active site.
Figure 5. The binding mode of compound 6b (red) and androstenedione (blue) in the AROM active site.
Molecules 26 04059 g005
Figure 6. The binding mode of compound: (A) 6e (red) and E1 (blue); (B) sulfate of 6c (green) and E1S (yellow) in the STS active site.
Figure 6. The binding mode of compound: (A) 6e (red) and E1 (blue); (B) sulfate of 6c (green) and E1S (yellow) in the STS active site.
Molecules 26 04059 g006
Figure 7. The binding mode of compound 6c (red) and E1 (blue) in the 17β-HSD1 active site.
Figure 7. The binding mode of compound 6c (red) and E1 (blue) in the 17β-HSD1 active site.
Molecules 26 04059 g007
Figure 8. The binding modes of compound 6c (red) and E2 (blue) in the ERα (A) and ERβ (B) binding sites.
Figure 8. The binding modes of compound 6c (red) and E2 (blue) in the ERα (A) and ERβ (B) binding sites.
Molecules 26 04059 g008
Table 1. Computed drug-like PC properties of compounds 6ae based on Lipinski’s Rule of Five.
Table 1. Computed drug-like PC properties of compounds 6ae based on Lipinski’s Rule of Five.
No.R1R2R3R4R5Molecular Weight (g/mol)Log PNumber of HB DonorsNumber of HB Acceptors
6aHFHHH255.252.6714
6bHFHFH273.242.9715
6cFFFHH291.233.2216
6dHClHHH271.702.9013
6eHClHClH306.153.4213
Table 2. Binding free energies of compounds 6ae and references (androstenedione, estrone sulfate (E1S), estrone (E1) and, estradiol (E2)) to the AROM, STS, 17β-HSD1, ERα, and ERβ.
Table 2. Binding free energies of compounds 6ae and references (androstenedione, estrone sulfate (E1S), estrone (E1) and, estradiol (E2)) to the AROM, STS, 17β-HSD1, ERα, and ERβ.
No.R1R2R3R4R5Binding Free Energy (kcal mol−1)
AROMSTS17β-HSD1ERαERβ
6aHFHHH−8.3−8.1/−7.1 *−7.9−8.8−8.5
6bHFHFH−7.9−8.3/−6.6 *−8.2−8.5−8.7
6cFFFHH−7.7−8.3/−7.9 *−8.5−8.5−8.9
6dHClHHH−7.5−8.1/−6.9 *−8.0−8.2−8.5
6eHClHClH−5.5−8.3/−6.9 *−8.4−7.9−7.6
Androstenedione-----−12.4----
E1S------−6.3---
E1------−8.9−8.9--
E2--------−10.7−11.1
* unsulfated/sulfated forms.
Table 3. Crystal data and structure refinement for 6ae.
Table 3. Crystal data and structure refinement for 6ae.
Identification Code6a6b6c6d6e
Empirical formulaC14H10FN3OC14H9F2N3OC14H8F3N3OC14H10ClN3OC14H9Cl2N3O
Formula weight (u)255.25273.24291.23271.70306.14
Temperature (K)120(2)120(2)120(2)120(2)120(2) K
Wavelength (Å)0.710730.710730.710730.710730.71073 Ĺ
Crystal systemTriclinicMonoclinicMonoclinicMonoclinicMonoclinic
Space groupP-1P21/nP21/cP21/cP21/c
Unit cell dimensions
a (Å)7.138(4)8.2384(16)9.162(3)27.8104(14)13.590(8)
b (Å)11.869(4)11.772(2)10.218(4)5.6525(2)14.898(6)
c (Å)14.148(4)12.358(3)12.658(4)7.3445(4)14.307(8)
α (°)89.54(2)90909090
β (°)75.50(3)102.211(17)93.66(3)94.728(4)116.70(4)
γ (°)88.42(3)90909090
Volume (Å3)1160.0(8)1171.4(4)1182.5(7)1150.61(9)2588(3)
Z44448
Z’21112
Density (calculated) (Mg/m3)1.4621.5491.6361.5681.571
Absorption coefficient (mm−1)0.1070.1230.1390.3260.499
F(000)5285605925601248
Crystal size (mm)0.154 × 0.116 × 0.0640.47 × 0.24 × 0.1050.646 × 0.383 × 0.2380.31 × 0.174 × 0.0920.269 × 0.209 × 0.151
Theta range for data collection (°)2.270 to 25.995.2.416 to 25.997.2.564 to 25.9862.204 to 25.9972.099 to 26.000
Index ranges−8 ≤ h ≤ 8, −14 ≤ k ≤ 14, −17 ≤ l ≤ 17−9 ≤ h ≤ 10, −14 ≤ k ≤ 14, −15 ≤ l ≤ 15−10 ≤ h ≤ 11, −12 ≤ k ≤ 12, −15 ≤ l ≤ 15−34 ≤ h ≤ 34, −6 ≤ k ≤ 5, −8 ≤ l ≤ 9−16 ≤ h ≤ 16, −18 ≤ k ≤ 18, −17 ≤ l ≤ 17
Reflections collected12,53814,02514,455587116,627
Independent reflections4556 [R(int) = 0.0543]2297 [R(int) = 0.0282]2325 [R(int) = 0.0268]2223 [R(int) = 0.0330]5074 [R(int) = 0.0311]
Completeness to theta = 25.242°99.9%100.0%99.9%98.6%99.7%
Refinement methodFull-matrix least-squares on F2Full-matrix least-squares on F2Full-matrix least-squares on F2Full-matrix least-squares on F2Full-matrix least-squares on F2
Data/restraints/parameters4556/0/3522297/0/1852325/0/1942223/0/1765074/0/369
Goodness-of-fit on F21.0501.0541.0631.1111.058
Final R indices (I > 2sigma(I))R1 = 0.0708, wR2 = 0.1924R1 = 0.0303, wR2 = 0.0830R1 = 0.0315, wR2 = 0.0825R1 = 0.0369, wR2 = 0.0869R1 = 0.0409, wR2 = 0.1051
R indices (all data)R1 = 0.0975, wR2 = 0.2167R1 = 0.0339, wR2 = 0.0848R1 = 0.0379, wR2 = 0.0858R1 = 0.0513, wR2 = 0.0948R1 = 0.0543, wR2 = 0.1121
Extinction coefficient0.064(10)n/an/an/an/a
Largest diff. peak and hole (e⋅Å−3)0.384 and −0.3180.266 and −0.2020.239 and −0.1950.299 and −0.3780.308 and −0.386
CCDC number2,063,9302,063,9312,063,9342,063,9322,063,933
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Daśko, M.; Dołęga, A.; Siedzielnik, M.; Biernacki, K.; Ciupak, O.; Rachon, J.; Demkowicz, S. Novel 1,2,3-Triazole Derivatives as Mimics of Steroidal System—Synthesis, Crystal Structures Determination, Hirshfeld Surfaces Analysis and Molecular Docking. Molecules 2021, 26, 4059. https://doi.org/10.3390/molecules26134059

AMA Style

Daśko M, Dołęga A, Siedzielnik M, Biernacki K, Ciupak O, Rachon J, Demkowicz S. Novel 1,2,3-Triazole Derivatives as Mimics of Steroidal System—Synthesis, Crystal Structures Determination, Hirshfeld Surfaces Analysis and Molecular Docking. Molecules. 2021; 26(13):4059. https://doi.org/10.3390/molecules26134059

Chicago/Turabian Style

Daśko, Mateusz, Anna Dołęga, Magdalena Siedzielnik, Karol Biernacki, Olga Ciupak, Janusz Rachon, and Sebastian Demkowicz. 2021. "Novel 1,2,3-Triazole Derivatives as Mimics of Steroidal System—Synthesis, Crystal Structures Determination, Hirshfeld Surfaces Analysis and Molecular Docking" Molecules 26, no. 13: 4059. https://doi.org/10.3390/molecules26134059

APA Style

Daśko, M., Dołęga, A., Siedzielnik, M., Biernacki, K., Ciupak, O., Rachon, J., & Demkowicz, S. (2021). Novel 1,2,3-Triazole Derivatives as Mimics of Steroidal System—Synthesis, Crystal Structures Determination, Hirshfeld Surfaces Analysis and Molecular Docking. Molecules, 26(13), 4059. https://doi.org/10.3390/molecules26134059

Article Metrics

Back to TopTop