Next Article in Journal
Bioactive Compounds and Sensory Analysis of Freeze-Dried Prickly Pear Fruits from an Inter-Andean Valley in Peru
Next Article in Special Issue
Comparison between Cashew-Based and Petrochemical Hydroxyoximes: Insights from Molecular Simulations
Previous Article in Journal
Enhanced Solubility and Stability of Aripiprazole in Binary and Ternary Inclusion Complexes Using Hydroxy Propyl Beta Cyclodextrin (HPβCD) and L-Arginine
Previous Article in Special Issue
Comparison of Phytochemical Profiles of Wild and Cultivated American Ginseng Using Metabolomics by Ultra-High Performance Liquid Chromatography-High-Resolution Mass Spectrometry
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 9
10.3390/molecules28093861
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elagolix Sodium Salt and Its Synthetic Intermediates: A Spectroscopic, Crystallographic, and Conformational Study

by
Samuele Ciceri
1,2,
Diego Colombo
2,
Enrico M. A. Fassi
1,
Patrizia Ferraboschi
2,
Giovanni Grazioso
1,
Paride Grisenti
3,
Marco Iannone
4,
Carlo Castellano
5 and
Fiorella Meneghetti
1,*
1
Department of Pharmaceutical Sciences, University of Milan, Via L. Mangiagalli 25, 20133 Milano, Italy
2
Department of Medical Biotechnology and Translational Medicine, University of Milan, Via C. Saldini 50, 20133 Milano, Italy
3
Chemical-Pharmaceutical Consulting and IP Management, Viale G. da Cermenate 58, 20141 Milano, Italy
4
Tecnomed Foundation, University of Milano-Bicocca, Via Pergolesi 33, 20900 Monza, Italy
5
Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(9), 3861; https://doi.org/10.3390/molecules28093861
Submission received: 3 April 2023 / Revised: 28 April 2023 / Accepted: 28 April 2023 / Published: 3 May 2023
(This article belongs to the Special Issue Development of Analytical Methodologies for Selective Extraction, Separation, Identification, and Purification of Valuable Compounds)
Browse Figures
Versions Notes

Abstract

:
Elagolix sodium salt is the first marketed orally active non-peptide gonadotropin-releasing hormone receptor antagonist (GnRHR-ant) for the management of hormone dependent diseases, such as endometriosis and uterine fibroids. Despite its presence on the market since 2018, a thorough NMR analysis of this drug, together with its synthetic intermediates, is still lacking. Hence, with the aim of filling this literature gap, we here performed a detailed NMR investigation, which allowed the complete assignment of the 1H, 13C, and 15N NMR signals. These data allowed, with the support of the conformational analysis, the determination of the stereochemical profile of the two atropisomers, detectable in solution. Moreover, these latter were also detected by means of cellulose-based chiral HPLC, starting from a sample prepared through an implemented synthetic procedure with respect to the reported ones. Overall, these results contribute to further understanding of the topic of atropisomerism in drug discovery and could be applied in the design of safe and stable analogs, endowed with improved target selectivity.

Graphical Abstract

1. Introduction

Elagolix sodium salt (CAS number 832720-36-2, 1, Figure 1) is the first non-peptide orally active gonadotropin-releasing hormone receptor antagonist (GnRHR-ant), approved since 2018 by FDA for the management of moderate-to-severe endometriosis-associated pain [1] and, more recently, against the symptoms related to uterine fibroids [2]. The oral bioavailability of 1 represents a breakthrough with the past, ensuring a rapid, dose-dependent, and reversible suppression of the hypothalamic–pituitary–gonadal axis, not achievable with the peptide analogs marketed before elagolix sodium salt 1 [3,4].
Chemically, compound 1 is an uracil-based derivative substituted at positions 1,3,5, and 6, bearing a stereocenter in (R)-configuration. For its close analogue 2, the literature data [5] reported a decreased rotation of the C-C bond evidenced in red in Figure 1, due to the interaction of the o-fluorine of the 5-aryl group with the methyl at the 6-position of the uracil moiety and the electronegative oxygen atom of the carbonyl at the 4-position. Therefore, two atropisomers were detectable under physiological conditions and, with a half-life (t1/2) of ~45 min, 2 belongs to the atropisomeric class 2, which includes atropisomers with a rate of interconversion ranging between minutes and a few months, as classified by LaPlante et al. [6]. Atropisomerism arises from the hindered rotation about single bonds; several drugs present this type of isomerism [7,8,9] and can be resolved or not, depending on the extent of the hindered rotation. For example, the atropisomers of the oral uric acid transporter 1 (URAT1) inhibitor lesinurad [10], belonging to class 3, were resolved through a semipreparative enantioselective supercritical fluid chromatography (SFC). On the contrary, the atropoisomers of drugs belonging to class 1, such as sildenafil (Viagra®) [11], cannot be resolved. For class 2 [6,12,13], the moderate interconversion rate between atropisomers causes issues in the drug manufacturing and quality control processes, such as a low batch-to-batch reproducibility, and inconsistencies in the safety/efficacy evaluation of the drug. A valuable strategy in avoiding the obtainment of class 2 atropisomers can be the elimination of the source of atropisomerism or the rigidification of the molecular structure. These strategies have been already applied to some elagolix analogues [12,14]. Moreover, the two atropisomers could have different pharmacological behaviors [14]. For 1, although it shares with 2 the same biaryl system responsible for atropisomerism, data related to atropisomers are absent in the literature.
The crystal structure of elagolix bound to human GnRHR was recently reported [15], pointing out the shallow non-peptide antagonist binding site and the key regions for ligand recognition: the pyrimidine ring, the benzyl group at position 1 of the uracil moiety, and the phenyl groups at position 5, without mentioning atropisomerism.
Although the synthesis of 1 was previously described in a 2005 patent [16] and in a 2008 paper [17], the literature requires a comprehensive 1H, 13C, and 15N NMR study of this drug and its synthetic intermediates. Firstly, only the list of 1H NMR resonances of 1 and some intermediates was reported [16,17]. Then, in an international patent application [18], the previously non-assigned 1H NMR signals of the most advanced intermediate of 1 (the corresponding ethyl ester) were added. In 2018, two international applications, the first regarding novel solid forms of 1 and its ethyl ester [19], and the second concerning a modified synthetic pathway [20], listed the 1H NMR data of 1 and its primary intermediates, without their assignments. In 2021, the lists of 1H or 13C resonances of 1 and some intermediates, without the relative assignments, were published [21,22,23]. More recently [24], the degradation products of 1 were studied; although the 1H and 13C NMR resonances of 1 and of two oxidation products were completely assigned, the atropisomerism was not mentioned.
Here, in a continuing effort to fully characterize molecules exerting important therapeutic activities [25,26,27], a complete NMR (1H, 13C, 15N) characterization of 1, and its main synthetic intermediates was carried out, taking in account not only that 1 has a proved pharmacological activity as GNRHR-ant, but also one of its advanced synthetic intermediates [17]. In addition, we obtained single-crystals of the first intermediate endowed with the biaryl system responsible for atropisomerism, suitable for the investigation of its solid-state structure by X-ray diffraction (SC-XRD). Furthermore, a detailed conformational analysis was performed on 1, to shed some light onto the geometry of low energy conformers. Taken together, these outcomes allowed the determination of the stereochemical profile of the two atropisomers of 1 detectable in solution.
The samples of 1 and its intermediates, required for the analytical characterization, were prepared in order to implement some modifications concerning literature procedures [16,17,28].
Overall, the results here obtained could be considered, both from the synthetic and analytical points of view, as a valuable starting point for the design of safe and stable analogs, endowed with improved target selectivity.

2. Results and Discussion

2.1. Chemistry

The synthesis of 1, performed following the literature procedure reported in Ref. [17], led to the obtainment of a poor-quality product (86% HPLC purity) and in a low overall yield (8%). These outcomes prompted us to introduce some modifications, in order to improve the quality and yield of 1, using cheaper and safer reagents.
Starting from the commercially available compound 3 (Scheme 1), the suitable precursor of the uracil moiety 4 was prepared, in 98% yield.
To build the 6-methyl uracil ring, instead of hazardous and unstable diketene [16,17], we chose the safer tert-butyl acetoacetate, according to a 2009 patent [28], obtaining 5 in 70% yield (Scheme 1). Indeed, unlike tert-butyl acetoacetate, diketene readily polymerizes on standing [29], with a high risk of evolution of toxic gases and explosion [30]. Our attempt to use in the synthetic route the cheaper ethyl acetoacetate afforded a mixture of 5 and its isomer (1:0.6 by NMR). In contrast, by using the more stable diketene-acetone adduct, if compared with the diketene, the undesired isomer of 5 was obtained almost quantitatively, contrary to what is reported in Ref. [20]. (Scheme 2, i and ii, respectively).
The 1H two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiment allowed understanding of the regiochemical outcome of the synthesis of 5, starting from 4. The presence of a cross peak between the methyl group at the 6-position of the uracil moiety (2.15 ppm) and the benzyl protons (5.36 ppm) confirmed their proximity (Figure 2A), whereas, for the undesired isomer of 5, this correlation is not present (Figure 2B).
The synthetic steps from 5 to 1 were carried out as described in the 2009 patent [28], but employing intermediate 6 instead of its iodo-analog, to avoid the use of the highly toxic and expensive iodine monochloride. Therefore, we performed a careful optimization of the bromination of 5 to obtain 6. Indeed, under the experimental conditions reported in the 2008 paper [17] (slow addition of 2.0 eq of bromine to a 1.1 M solution of 5 in acetic acid), we obtained 6 from 5 in a poor yield (43%). On the contrary, the slow addition of a 1.1 M solution of bromine (1.1 eq) to a 0.25 M solution of 5, both in acetic acid, allowed us to obtain 6 in a higher yield (96%). From intermediate 6, through a Suzuki cross coupling with 2-fluoro-3-methoxyphenylboronic acid, catalyzed by bis(tri-tert-butyl-phosphine)palladium(0), derivative 7 was directly recovered in 88% yield by precipitation from the reaction mixture; pure 8 was obtained in 81% yield through nucleophilic substitution by 7 on (R)-2-((tert-butoxy-carbonyl)amino)-2-phenylethyl methane-sulfonate 10, followed by N-Boc deprotection with methane-sulfonic acid. Otherwise, reproducing the 2008 synthetic procedure [17], 8 was obtained through a Mitsunobu-like reaction in 6 followed by Suzuki cross coupling catalyzed by tetrakis(triphenylphosphine)palladium(0), but each of these two steps required a careful chromatographic purification to remove the by-products triphenylphosphine oxide, tert-butyl hydrazine dicarboxylate, and the palladium catalyst. At the end, the ethyl ester 9, obtained in 50% yield by N-alkylation of 8 with ethyl 4-bromobutyrate, was hydrolyzed with sodium hydroxide to afford 1, isolated in 88% yield by extraction with methyl isobutyl ketone, followed by precipitation from n-heptane. The overall yield of this process was 21% and the HPLC purity of 1 was improved to 99.9%. In comparison with the salt formation described in the original synthesis [17], in which the isolated elagolix acid form was converted into 1 employing an ion exchange resin, the direct extraction of the sodium salt from the alkaline hydrolytic reaction mixture [28] was chosen to avoid the formation of the lactam impurity 11 (Scheme 3). We isolated and characterized 11 when we followed the literature procedure [11] (Supplementary Materials, Figures S43 and S44).

2.2. NMR Spectroscopy

A detailed NMR study was carried out on 1 and its precursors 5–9 (Scheme 1). 1H NMR spectra of 1 were recorded in deuterated methanol (CD3OD), for a better comparison with the signals already reported in the literature [17], while the spectra of intermediates 5 and 7–9 were recorded in deuterated chloroform (CDCl3) at 298 K. A mixture of CDCl3/CD3OD 9:1 was used for 6, due to its poor solubility in CDCl3. Table 1, Table 2 and Table 3 summarizes the unambiguous assignments of all 1H, 13C and 15N signals, established by combining the information gathered from 1D NMR spectra and 2D homo-correlated (COSY and NOESY) and hetero-correlated (1H–13C HSQC, 1H–13C HMBC, and 1H-15N HMBC) NMR spectra. Atoms were numbered as reported in Figure 1. Most of the proton assignments were accomplished using general knowledge of chemical shift dispersion, with the aid of the proton–proton coupling pattern (1H NMR spectra), gs-COSY, and NOESY experiments. In ambiguous cases, gs-HSQC and gs-HMBC spectra were used as a definitive and unequivocal tool to make specific assignments, especially for quaternary carbons. For the attribution of 13C NMR resonances, in some cases, the analysis of the nJ13C-19F [31] turned out to be very useful, although this coupling increased the complexity of the spectrum.
In particular, for compounds 5–7, the resonances of all protons were assigned using general knowledge of chemical shift dispersion, with the aid of the proton–proton coupling pattern and the COSY experiment for the aromatic protons. In the case of H-10, an additional coupling constant is recognizable (J1H-19F = 12.3 Hz) due to the presence of fluorine at the 9-position. The correct structural features of 5 were verified by a NOESY experiment. Indeed, as shown in Scheme 2 two different isomers can be obtained under these reaction conditions. The presence of a cross peak between the three H-15 (2.15 ppm) and the two H-7 (5.36 ppm) in the NOESY spectrum of 5 (Figure 2A), absent in the spectrum of its isomer (Figure 2B), confirmed its correct chemical structure. Starting from compound 8, obtained from 7 by introduction of the (R)-phenylamino ethyl moiety at the 3-position of the uracil core, for all the following intermediates, two distinct resonances for the H-21 proton (6.85 and 6.78 ppm, 6.83 and 6.76 ppm, and 6.76 and 6.61 ppm for 8, 9, and 1, respectively), in a ratio 1:1.07, 1:1.04, and 1:1.13 (Figure 3A), were observed. The same phenomenon was detected, though less conspicuous (distance between the two signals ≤ 0.01 ppm), for H-15 and H-22 signals of 9 and 1 (Figure 3B,C). The presence of a (R)-stereocenter in compounds 8, 9, and 1, together with the atropisomerism, explained the two signals observed in their spectra [5] (Figure 3). Furthermore, the ratio of the two signal sets, being different from 1.00, indicated that one of the two atropisomers was favored.
The absence of any mention of atropisomerism in Ref. [24] could be explained by the choice of DMSO-d6 as NMR solvent. Indeed, the 1H-NMR spectrum registered in this solvent is less indicative of highlighting the simultaneous existence of two diastereoisomers of 1, due to the presence of both the (R)-stereocenter and atropisomerism. From the comparison between the 1H NMR spectrum registered in CD3OD and that in DMSO-d6, it is evident that the most significant signal shown in the CD3OD spectrum (H-15 protons) is only partially resolved in DMSO-d6 (see Table 3 of reference [24] and the enclosed spectra in the Supplementary Materials Figure S41 of the present article). As a consequence, the signals due to the H-21, well resolved in both spectra, were not correctly assigned and integrated [24].
NOESY experiments performed on 1 dissolved in CD3OD revealed the proximity, within 5.0 Å, of H-15 (2.085 and 2.076 ppm) to H-21 (6.76 and 6.61 ppm) and H-7 (5.43 ppm), of H-19 (7.10 ppm) to H-22 (3.883 and 3.880 ppm), and of H-30 (2.40 ppm) to H-23 (4.26–4.07 ppm), H-24 (4.11 ppm), and H-26/H-26’ (7.29 ppm), in both atropisomers (Supplementary Materials Figure S40).
The assignments of carbon atoms of CH, CH2 and CH3 groups were confirmed by the gs-HSQC experiment. The quaternary carbon atoms were unambiguously assigned using the information obtained from the 1H–13C gs-HMBC experiment (the observed couplings for 1 are listed in Table 4) and from nJ13C-19F [31]. Due to atropisomerism, the 13C NMR spectra of compounds 8, 9, and 1 showed, for some carbon atoms, two singlets very close to each other (Table 2). In this case also, the 13C signals indicated a relative ratio between atropisomers other than 1.00.
1H-15N HMBC experiments were performed to further support the proposed structures of the studied compounds. Resonances of N-1, and N-3 were assigned for all compounds, and also of N-29 for 8, 9 and 1. The value of the chemical shift of N-1 remained almost unchanged from intermediate 5 to 1. On the contrary, the value of the chemical shift of N-3 for 5, 6, and 7 changed from 154.3, 153.4, and 151.9 ppm to 160.1, 161.3, and for 8, 9, and 1 to 160.5 ppm as a result of the N-alkylation of this nitrogen [32]. The same behavior was observed for N-29, whose chemical shift value changed from 31.2 ppm for 8 (primary amine) to 43.6 and 46.3 for 9 and 1 (secondary amine), respectively. For all compounds, the 1H–15N gs-HMBC spectra showed cross peaks of H-15 and H-7 with N-1, while, in the case of N-3 and N-29 of 8, 9, and 1, the cross peaks H-24/N-3, H-23/N-29, H-24/N-29 and H-31/N-29 were observed.

2.3. HPLC Analyses

The chromatogram of the HPLC analysis performed on 1 using a non-chiral reverse stationary phase (RP-C18) showed a broad peak, indicating the presence of two very similar but not distinguishable molecules (Figure 4A, see Section 3.2.8 and Section 3.4 for HPLC conditions). On the contrary, the HPLC analysis on a chiral stationary phase of cellulose tris(3,5-dimethylphenylcarbamate) (Figure 4B, see Section 3.2.8 and Section 3.4 for HPLC conditions) clearly showed the presence of two distinct peaks, corresponding to the two atropisomers (diastereoisomers) of 1.
Since 7 is the first intermediate carrying the biaryl system involved in the atropisomerism, we decided to analyse 7 by chiral HPLC (see section Section 3.2.4 and Section 3.4 for HPLC conditions). As shown in Figure 5, two peaks are detectable, confirming the existence of two atropisomers (enantiomers) in this compound also.
As reported in some literature examples of gas chromatographic (GC) and HPLC analysis of stereo-labile chiral compounds [33,34,35,36], the observed chromatographic profile, with a plateau between the two chromatographic peaks that does not reach the baseline, indicates the existence of a dynamic interconversion of the two atropisomers within the time scale of the analysis, the so called “on-column interconversion”.

2.4. Structure Description

For 7, the first synthetic intermediate bearing the biaryl system contributes to atropisomerism, so we also decided to analyze its solid-state structure through SC-XRD. Single crystals were obtained by the slow evaporation of a MeOH solution after 1 week. Crystallographic data and refinement details are given in Section 3.
Compound 7 crystallized in the orthorhombic achiral space group P bca; its structure is shown in Figure 6 as an ORTEP diagram [37], indicating the arbitrary atom-numbering scheme used in the following discussion.
The overall molecular structure of the racemate was characterized by a dihydro-pyrimidine nucleus bound to a 5-(2-fluoro-3-methoxyphenyl moiety and a (trifluoromethyl)benzyl group. The angles between the best mean plane calculated for the heterocyclic ring and these two aromatic portions were 87.4(1)° and 67.2(1)°, respectively. The dihydro-pyrimidine was nearly planar, with a maximum deviation of 0.021(1) Å. Moreover, the conformation of the molecule was also characterized by the dihedral angles C2-N1-C7-C8 of −105.3(1)°, N1-C7-C8-C13 of −144.7(3)° and C4-C5-C16-C21 of 110.4(3).
The crystal packing, shown in Figure 7, was consolidated by strong dimeric H-bonds involving N3-H∙∙∙O1I (I at 1 − x, −y, −z); the donor-acceptor (D∙∙∙A) distance is 2.87(1) Å, the hydrogen and the acceptor (D-H∙∙∙A) are at 2.063(2) Å, and the angle is 155.2(5)°. These H-bonds give rise to the formation of molecular chains along the a axis. Parallel π-π stacking interactions between (trifluoromethyl)benzyl moieties were present and further contributed to the crystal packing: the distance between the centroids was 3.20(6) Å, while the angle between the centroid-centroid vector and the plane normal was 3.62(4)°. Loose CH∙∙∙O and CH∙∙∙F contacts contributed to the stabilization of the crystal structure; these non-traditional H-bonds were established between C12-H∙∙∙O2II, D∙∙∙A = 3.26(8) Å, D-H∙∙∙A = 2.63(1) Å, angle = 124.3(9)° (II at 1 + x, y, z), C11-H∙∙∙F4III, D∙∙∙A = 3.42(8) Å, D-H∙∙∙A = 2.54(1) Å, angle = 163.5(9)° (III at ½ + x, ½ − y, z) and C19-H∙∙∙F5IV, D∙∙∙A = 3.37(8) Å, D-H∙∙∙A = 2.48(1) Å, angle = 160.3(9)° (IV at x, ½ + y, ½ − z). The Van der Waals interactions between the aromatic rings of adjacent molecules favored a two-dimensional-sheet structure, in which the molecules in the two layers are inclined as shown in Figure 7.
It can be hypothesized that the crystallization conditions (methanol + 5% water) facilitated the solubilization of the compound and promoted the formation of a dense network of interactions, which ultimately allowed the orderly organization of the molecules to form a crystal.

2.5. Conformational Analysis

With the aim of rationalizing the obtained results, the conformational behavior of 1 was investigated by means of computational tools. The structure of 1 was simplified, as reported in Figure 8A, shortening the lateral chain carrying the carboxylate moiety, to reduce the number of dihedral angles to be scanned. Then, the conformational analysis was performed through a mixed molecular mechanics-DFT/B3LYP-D3/6-31G(**) level approach in the gas phase, and later in water and methanol as solvents (see Section 3 for details).
The four torsional angles of the simplified structure of 1 were initially scanned by increments of 10 degrees using the MacroModel tool of Maestro (release 2021-2, Schrödinger, LLC, New York, NY, USA), applying the OPLS4 force field. For simplicity, these calculations were conducted in two independent steps to identify the combination of torsion angles leading to low energy conformers. Firstly, τ1 and τ2 were combined and scanned (Figure 8B), and successively τ3 and τ4 were considered (Figure 8C). The torsion angles/energy 3D plots shown in Figure 8B,C displayed the low energy conformers’ combinations. Of note, only four τ1/τ2 combinations were energetically favored, while the τ3/τ4 combinations led to six different low energy conformers. Then, the τ1–τ4 torsion angles’ values leading to the lowest energy conformers were systematically combined, obtaining a library of 24 different conformers, which were successively optimized in the gas phase (Table 5). At a later stage, the energy of each optimized structure was recalculated in water (as single point calculation), and finally in methanol, the solvent employed in the NMR experiments (see Supplementary Materials, Table S1). These calculations were performed applying the DFT/B3LYP-D3/6-31G(**) level of theory [38,39], using the Jaguar module of Maestro (release 2021-2, Schrödinger, LLC, New York, NY, USA). Notably, some of them converged into the same energy minima conformer; in particular, 1N converged to 1T, while 1P converged to 1R (Table 5).
Finally, the τ5 dihedral angle, the one capable of generating the atropisomerism phenomena, was scanned in the lowest energy conformer of the simplified structure of 1 (Figure 8A, conformer 1M, Table 6).
These calculations were performed scanning the τ5 torsion angle, by an increment of 10° (from −180° to +180°). We contextually optimized the conformers at DFT/B3LYP-D3/6-31G(**) level of theory, using the CPCM water solvent model (see Section 3 for details). The achieved results (Figure 9, top panel) suggested the presence of four different energy minima conformers. Their geometries were geometrically optimized at the same level of theory, in order to improve the accuracy of these results. The outcomes showed that 1M-B converged into the minimum energy 1M-A, while 1M-C and 1M-D represented another two relative minimum energy conformers (Table 7). Single point calculations using the CPCM methanol solvent model (which is the solvent employed in the NMR experiments), displayed results similar to those observed using the CPCM water solvent model.
To conclude, 1M-A, the atropisomer aR, was the most abundant (83.6%) at room temperature, while 1M-C and 1M-D, corresponding to the opposite atropisomer aS, respectively, composed the remaining 16.4%, in both solvent models (water and methanol). Since one atropisomer is more abundant than the other in methanol solution, and this finding was in agreement with the NMR data (see Section 2.2.), the diagnostic H-21 (6.61 ppm), H-22 (3.880 ppm), and H-15 (2.076 ppm) signals can be assigned to the most abundant atropisomer aR, while the H-21 (6.76 ppm), H-22 (3.883 ppm), and H-15 (2.085 ppm) signals to the least abundant atropisomer aS.
Additionally, since the methoxy group is one of the key elements of the uracil pharmacophore as GnRHR-ant [3,40], its position was investigated in the lowest energy conformers of 1M-A and 1M-D. In particular, the dihedral angle highlighted by the red arrow in 1M-D (Figure 9) was rotated, forming the two conformers I and II. Here, the H-22 atoms were oriented in the direction of C15 (t6 = −66°, distance CH22-H19 = 3.6 Å) in I, and in opposite direction (t6 = +66°, distance CH22-H19= 3.6 Å) in II. The structures of the latter were different from those of the energy minima 1M-A and 1M-D, characterized by t6 of 0° and having the CH22-H19 distance of 2.56 Å. Energy optimization calculations at DFT/B3LYP-D3/6-31G(**) level of theory [38,39], in the gas phase and in CPCM methanol solvent model [41], suggested that the conformers with the lowest distance between CH22 and H19 (t6 = 0°) were the preferred ones when the CPCM methanol solvent model was used [41]. This outcome agreed with the presence of a cross peak between the H-19 (7.10 ppm) and H-22 (3.883 and 3.880 ppm) protons of the two atropisomers in the NOESY spectrum of 1 (Supplementary Materials Figure S40) and with the conformation of the solid-state structure determined for intermediate 7. Conversely, in the gas phase, the population of 1M-D conformers with distance CH22-H19 = 3.6 Å was at 69.5%, while that for the same conformers of 1M-A was 55.5% (Table 7).

3. Materials and Methods

3.1. General

All reagents and solvents were purchased from Sigma-Aldrich (Merck Life Science S.r.L., Milano, Italy). TLC analyses were performed on silica gel 60 F254 plates, precoated with a fluorescent indicator (Merck Life Science S.r.L., Milano, Italy); spots were detected by UV lamp 254 nm, or by a 0.3% w/v ninhydrin solution in n-butanol/acetic acid (100:3) and heating at 110 °C.
Optical rotation values were registered on an Anton Paar instrument (Mod MCP 100; Anton Paar Strasse 10, 8054 Graz, Austria) at 589 nm and 25 °C.

3.2. Synthesis of 1 from 3

3.2.1. 1-(2-fluoro-6-(trifluoromethyl)benzyl) urea (4)

Compound 4 was synthesized, starting from 2-fluoro-6-(trifluoromethyl)benzylamine, as reported in Ref. [17] in 98% yield.
Rf 0.60 (CH2Cl2/MeOH 9:1)
1H NMR (500 MHz, DMSO) δ 7.61-7.56 (m, 3H), 6.15 (t, J = 5.3 Hz, 1H), 5.47 (s, 2H), 4.36 (d, J = 5.0 Hz, 2H).
The other physico-chemical properties were in agreement with the reported ones [17,42].

3.2.2. 1-(2-fluoro-6-(trifluoromethyl)benzyl)-6-methylpyrimidine-2,4(1H,3H)-dione (5)

Compound 5 was synthesized, starting from 4, as in Ref. [28], in 70% yield. The physico-chemical properties were in agreement with those reported [17,28].
Rf 0.27 (hexane/ethyl acetate 1:1)

3.2.3. 5-bromo-1-(2-fluoro-6-(trifluoromethyl)benzyl)-6-methylpyrimidine-2,4(1H,3H)-dione (6)

A solution of bromine (2 mL; 38.5 mmol) in glacial acetic acid (36 mL) was slowly added (in about 20 min) to a solution of 5 (10.58 g; 35 mmol) in glacial acetic acid (140 mL) under vigorous stirring. The reaction progress was monitored by TLC analysis (hexane/ethyl acetate 7:3) until starting material disappearance (2h). The excess of bromine was removed bubbling nitrogen through the reaction mixture until almost complete disappearance of its colour. The reaction mixture was concentrated at reduced pressure (60 °C, 30 mmHg), affording a yellow residue which was suspended in tert-butyl methyl ether (25 mL). The suspension was stirred at room temperature for 2 h, filtered by suction and washed with tert-butyl methyl ether (5 mL). The solid was dried (60 °C, 2 mmHg, 8 h) affording 6 (12.86 g, 33.7 mmol; 96%). The physico-chemical properties were in agreement with those reported [17].
Rf 0.69 (hexane/ethyl acetate 1:1)

3.2.4. 5-(2-fluoro-3-methoxyphenyl)-1-(2-fluoro-6-(trifluoromethyl)benzyl)-6-methylpyrimidine-2,4(1H,3H)-dione (7)

An aqueous solution (10 mL) of KOH (4.64 g; 82.70 mmol) was added to a mixture of 6 (7.7 g; 20.20 mmol) and 2-fluoro-3-methoxy-phenyl boronic acid (4.12 g; 24.24 mmol) in acetone (30 mL) and water (20 mL). The obtained solution was degassed (bubbling argon for about 30 min) and tri-tert-butyl-phosphonium tetrafluoroborate (62 mg; 0.21 mmol) was added. After 20 min at 45 °C, Pd(OAc)2 (23 mg; 0.10 mmol) was added. The reaction mixture was refluxed under argon atmosphere and stirred while the reaction progress was monitored by TLC analysis (hexane/ethyl acetate 1:1) until disappearance of the starting material (24 h). The reaction mixture was cooled to 55 °C and glacial acetic acid (3.5 mL) was added and further stirred for 30 min. A white precipitate formed. The reaction mixture was cooled to room temperature and stirred at this temperature for 2 h. The white precipitate was filtered by suction and washed with water (15 mL) and then with methanol (30 mL). The solid was dried (60 °C; 2 mmHg; 8 h) affording 7.38 g (86% yield). The physico-chemical properties were in agreement with those reported [28,43].
Rf 0.36 (hexane/ethyl acetate 1:1)
RP-HPLC analysis (elution in gradient mode: Mobile phase: A = water with 0.05% trifluoroacetic acid; B = acetonitrile with 0.05% trifluoroacetic acid. Gradient: 95%A/5%B to 5%A/95%B over 50 min, then 5%A/95%B to 1%A/99%B over 0.1 min, then hold 1%A/99% for 0.8 min and back to 95%A/5% over 0.2 min, hold such gradient for 4 min. Flow rate 2.0 mL/min): tR = 25.1 min.
HPLC on chiral stationary phase (elution in isocratic mode: hexane/IPA 1:1. Flow rate 0.5 mL/min): tR (I atropisomer) = 27.4 min, tR (II atropisomer) = 32.8 min.

3.2.5. (R)-2-[(tert-butoxy-carbonyl)amino]-2-phenylethyl methane-sulfonate (10)

Compound 10 was synthesized, starting from (R)-phenyl-glycinol, as in Ref. [44], in 73% yield. The physico-chemical properties were in agreement with those reported [44].
Rf 0.80 (CH2Cl2/MeOH 95:5)

3.2.6. (R)-3-(amino(phenyl)methyl)-5-(2-fluoro-3-methoxyphenyl)-1-(2-fluoro-6-(trifluoromethyl)benzyl)-6-methylpyrimidine 2,4(1H,3H)-dione (8)

Compound 8 was synthesized, starting from compounds 7 and 10, as in Ref. [28], in 80% yield. The physico-chemical properties were in agreement with those reported [28,45].
[α]D25 + 15.1 (c 1, CHCl3).
Rf 0.17 (hexane/ethyl acetate/ Et3N 1:1:0.1)

3.2.7. Ethyl (R)-4-(((5-(2-fluoro-3-methoxyphenyl)-3-(2-fluoro-6-(trifluoromethyl)benzyl)-4-methyl-2,6-dioxo-3,6-dihydropyrimidin-1(2H)-yl)(phenyl)methyl)amino)butanoate (9)

Compound 9 was synthesized, starting from 8, as in Ref. [28], in 50% yield. The physico-chemical properties were in agreement with those reported [28,45].
[α]D25 − 5.7 (c 1, CHCl3).
Rf 0.44 (hexane/ethyl acetate/ Et3N 1:1:0.1)

3.2.8. Elagolix Sodium Salt (1)

Compound 1 was synthesized, starting from 9, as reported in Ref. [28] in 88% yield.
Rf 0.74 (CHCl3/MeOH/H2O 12.5:4:0.5)
[α]D25 + 12.2 (c 1, MeOH).
RP-HPLC analysis (elution in gradient mode in agreement with Ref. [28]): tR = 9.0 min.
HPLC on chiral stationary phase (elution in isocratic mode: hexane/IPA 6:4 + 0.1% TFA. Flow rate 0.3 mL/min): tR (I atropisomer) = 38.0 min, tR (II atropisomer) = 46.6 min.
The other physico-chemical properties were in agreement with those reported [17,28].

3.3. NMR Spectroscopy

NMR spectra were recorded on a Bruker AVANCE 500 spectrometer (Bruker, Billerica, MA, USA) equipped with a 5 mm broadband inverse (BBI) detection probe with field z-gradient operating at 500.13, 125.76, and 50.69 MHz for 1H, 13C, and 15N, respectively. The spectra were recorded at 298 K for 5, and 7–9 in chloroform-d (CDCl3, isotopic enrichment 99.9 atom % D) and, in the case of 1 and 6, in methanol-d4 (CD3OD, isotopic enrichment 99.9 atom % D) and chloroform-d/ methanol-d4 9:1 mixture, respectively. Chemical shifts (d) of the 1H NMR and 13C NMR spectra are reported in ppm using the signal for residual solvent proton resonance as the internal standard (1H NMR: CDCl3 7.26, CD3OD 3.31 ppm; 13C NMR: CDCl3 77.0 (central line), CD3OD 49.00 (central line) ppm). For 15N nuclei, in the gs-1H-15N HMBC experiment, nitromethane was used as the external reference (15N at 381.7 ppm [32]) whereas, in the gs-1H-15N HSQC experiment, ammonia (15N at 0.0 ppm [32]) was used. The pulse widths were 8.00 μs (90°) for 1H, 13.00 μs (90°) for 13C, and 27.50 μs (90°) for 15N. Data were collected and processed by XWIN-NMR software (version 3.5, Bruker, Billerica, MA, USA) running on a PC with Microsoft Windows 7. The samples (20 mg) were dissolved in the appropriate solvent (0.75 mL) in a 5 mm NMR tube. The acquisition parameters for 1D were as follows: 1H spectral width of 5000 Hz and 32 K data points providing a digital resolution of ca. 0.153 Hz per point, relaxation delay 20 s, “zg” pulse sequence of the Bruker library was used; 13C spectral width of 29,499 Hz, and 32 K data points providing a digital resolution of ca. 0.900 Hz per point, relaxation delay 2 s, “zgpg” pulse sequence of the Bruker library was used. The experimental error in the measured 1H-1H coupling constants was ±0.5 Hz. The splitting pattern abbreviations were as follows: s, singlet; d, doublet; t, triplet; q, quartet; m, multiplet; and br, broad signal. Except for NOESY, standard Bruker microprograms using gradient selection (gs) were applied for two-dimensional experiments. Gs-COSY-45 (“cosygpqf” pulse sequence of the Bruker’s library) and phase-sensitive NOESY (“noesyph” pulse sequence of the Bruker’s library) experiments were acquired with 512 t1 increments, 2048 t2 points, and a spectral width of 10.0 ppm. The NOESY experiments were performed on samples degassed under a flush of argon in a screwcap sample tube. There were no significant differences in the results obtained at different mixing times (0.5–1.5 s). The acquisition data for gs-HSQC (“hsqcetgp” pulse sequence of the Bruker’s library) and gs-HMBC (“hmbcgplpndqf” pulse sequence of the Bruker’s library) experiments were acquired with 512 t1 increments, 2048 t2 points, and a spectral width of 10.0 ppm for 1H and 240 ppm for 13C. Delay values were optimized to 140 Hz for 1J13C,1H and 8.0 Hz for nJ13C,1H.
The gs-1H-15N HMBC (“hmbcgpndqf” pulse sequence of the Bruker’s library) and HSQC (“hsqcetgpsi” pulse sequence of the Bruker’s library) experiments, were performed with 256 t1 increments, 1024 t2 points, and a spectral width of 10.0 ppm for 1H and 600 ppm for 15N setting an acquisition time of 0.5 s, a relaxation delay of 2 s, a 1J15N,1H value of 90.0 Hz, and a nJ15N,1H value of 4.5 Hz. This last parameter was set after several attempts between 1 and 11 Hz. The total experimental time for 1H-15N gs-HMBC analyses was about 12 h.
The NMR spectra are available in the Supplementary Materials section.

3.4. HPLC Analyses

The employed HPLC units were:
−
Chiral HPLC analysis: a Merck-Hitachi (Hitachi Ltd., Tokyo, Japan), equipped with a UV detector model L-4250, pump system model L-6200 and a chromato-integrator model D-2500. The column employed in the analyses was a Phenomenex Lux-Cellulose 1 (Phenomenex, Torrance, CA, USA). The dimension of the column is 250 mm × 4.6 mm, 3 µm. The elution was in isocratic mode with the indicated eluant and flow. All the samples were measured at λ = 254 nm and 25 °C.
−
RP-HPLC analysis: Agilent 1100 system (Agilent Technologies, Waldbronn, Germany) equipped with a Zorbax SB-C18 column (150 mm × 3.0 mm, 3.5 µm) for 1 and with a Supelco Discovery C18 (250 mm × 4.6 mm, 5.0 µm) for 7.
All the samples were measured at λ = 254 nm and 25 °C.

3.5. X-ray Analysis

X-Ray analyses were performed on a Bruker SMART APEX II CCD Single Crystal X-ray Diffractometer (Bruker, Karlsruhe, Germany), equipped with graphite-monochromated Mo-Kα radiation (λ = 0.71073 Å) at 298(2) K.
X-ray data were acquired in the θ range 2–20° recording four sets of 360 bidimensional CCD frames with the following operative conditions: omega rotation axis, scan width 0.5°, acquisition time 50 s, sample-to-detector distance 50 mm, phi angle fixed at four different values (0°, 90°, 180° and 270°) for the four different sets. Omega-rotation frames were processed using the SAINT software [46] for data reduction (intensity integration, background, Lorentz, and polarization corrections) and for the determination of accurate unit-cell dimensions. Absorption effects were empirically evaluated by the SADABS software [47], applying an absorption correction to the data.
The crystal structure was solved by direct methods and refined on F2 by full-matrix least-squares using Bruker’s SHELXL-2018/1 [48].
Details are summarized in Table 8. Crystallographic data were deposited to the Cambridge Crystallographic Data Center under accession number CCDC 2248213. The refined structure was inspected using ORTEP-3 (v. 2020.1) [37] and analyzed by Mercury 4.0 (v. 2021.3.0) [49] and PARST [50], within the WinGX suite (v. 2021.3) [50]. Graphical representations were rendered with Mercury.

3.6. Conformational Analysis

These studies were performed using the tools available in Maestro (release 2021-2, Schrödinger, LLC, New York, NY, USA). The simplified structure of 1 (Figure 8A) was created by means of the “Build” tool; the potential energy of the conformers, resulting from the different combination of the τ1–τ4 values (as depicted in Figure 8A), was calculated with the “Coordinate scan” tool of Macromodel, using the OPLS4 force field and choosing the proper couple of torsional angles to be scanned in the calculations. The plots reported in Figure 8B,C were obtained by the “plot coordinate scan” tool of Macromodel, focusing the displayed energy values from 0 to 5 kcal/mol. Then, all the energy minima identified by Macromodel, or manually constructed rotating the dihedral angles indicated throughout the Results section, were optimized by using Jaguar, at the DFT/B3LYP-D3/6-31G(**) level of theory [38,39].

4. Conclusions

In this paper we described an efficient synthesis of elagolix sodium salt (1), bearing axial chirality, which was thoroughly explored through spectroscopic, analytical and theoretical techniques. The obtained results offered novel clues to recognize and characterize its atropisomers. To optimize the original synthetic route, we applied a new strategy capable of improving the reaction yields and the purity of the intermediates and of the final product 1, using less hazardous reagents. The optimized synthetic pathway developed for the obtainment of elagolix sodium salt could be also followed in the next future for the synthesis of the selected analogues. We planned to develop a protocol that could be applied for the high-yield synthesis of new non-peptide congeners of 1, and to fill a literature gap concerning the NMR data of 1 and its synthetic intermediates, which are essential for comprehensively studying the pharmacological activity and interconversion thermodynamics of atropisomers.
We reached our goal, using high field NMR spectroscopy and theoretical calculations to analyze the molecular mobility of the atropisomers of 1, obtaining results that allowed their stereochemical characterization. Actually, the agreement of the NMR outcomes with the conformations located through the modelling study led to the pursuing of the stereochemical profile of the two atropisomers aR and aS detectable in solution, useful for deepening comprehension of the drug–target interactions.
In addition, a new chiral-phase HPLC method for controlling the synthetic steps and determining the chemical and optical purity was developed. Moreover, the solid-state structure of intermediate 7, crystallized as a racemic compound and investigated by SC-XRD, provided insights into its overall conformation, crystal packing and molecular interactions. The structural analysis of key dihedral angles showed values in agreement with the most abundant modeled conformer 1M-A. These outcomes contributed to shedding light on the structural determinants involved in the control of the spatial arrangement of the substituents within this molecular framework, useful for future development of derivatives with higher activity. The interest in developing small molecules orally active as GnRHR modulators is still very high because the therapies employing peptide modulators are characterized by several drawbacks (i.e., subcutaneous injection). In conclusion, the synthetic, spectroscopic, and crystallographic results of this investigation provided a valuable starting point for future studies to obtain, with minor structural modifications, separable atropisomers. These results can be extended to other new potentially active candidates to establish the relationship existing between their stereochemical features and the pharmacological properties. As atropisomers allow a predictable relative arrangement of groups in space, with those parameters in hand, we could design new strategies to address the constraints of the controlled stereoselective synthesis of atropisomers, stable at ambient temperature, offering an increasing structural diversity with application in the medicinal chemistry research field. Overall, this work demonstrated that the complete understanding of the atropisomerism phenomena is a very promising source of information for the development of drugs against a variety of diseases. Although there are still many gaps in the knowledge of these systems, the fundamental understanding of atropisomerism enables new, well-reasoned approaches to develop better therapeutic strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28093861/s1, Figure S1: 1H NMR of compound 5, Figure S2: 13C NMR of compound 5, Figure S3: COSY of compound 5, Figure S4: 1H-13C HSQC of compound 5, Figure S5: 1H-13C HMBC of compound 5, Figure S6: 1H-15N HMBC of compound 5, Figure S7: 1H-15N HSQC of compound 5, Figure S8: 1H NMR of compound 6, Figure S9: 13C NMR of compound 6, Figure S10: COSY of compound 6, Figure S11: 1H-13C HSQC of compound 6, Figure S12: 1H-13C HMBC of compound 6, Figure S13: 1H-15N HMBC of compound 6, Figure S14: 1H-15N HSQC of compound 6, Figure S15: 1H NMR of compound 7, Figure S16: 13C NMR of compound 7, Figure S17: COSY of compound 7, Figure S18: 1H-13C HSQC of compound 7, Figure S19: 1H-13C HMBC of compound 7, Figure S20: 1H-15N HMBC of compound 7, Figure S21: 1H-15N HSQC of compound 7, Figure S22: 1H NMR of compound 8, Figure S23: 13C NMR of compound 8, Figure S24: COSY of compound 8, Figure S25: 1H-13C HSQC of compound 8, Figure S26: 1H-13C HMBC of compound 8, Figure S27: 1H-15N HMBC of compound 8, Figure S28: 1H NMR of compound 9, Figure S29: 13C NMR of compound 9, Figure S30: COSY of compound 9, Figure S31: 1H-13C HSQC of compound 9, Figure S32: 1H-13C HMBC of compound 9, Figure S33: 1H-15N HMBC of compound 9, Figure S34: 1H NMR of compound 1, Figure S35: 13C NMR of compound 1, Figure S36: COSY of compound 1, Figure S37: 1H-13C HSQC of compound 1, Figure S38: 1H-13C HMBC of compound 1, Figure S39: 1H-15N HMBC of compound 1, Figure S40: NOESY of compound 1, Figure S41: 1H NMR of compound 1 in DMSO-d6, Figure S42: 1H NMR of compound 1 in D2O, Figure S43: 1H NMR of compound 11, Figure S44: COSY of compound 11, Table S1: Relative energies and equilibrium percentages of the optimized conformations endowed with ΔE < 5 kcal/mol in the gas phase, in water and methanol solvent models. Table S2: Internuclear H-15/H-21 distance estimated by 2D NOESY experiments.

Author Contributions

Conceptualization, P.G., P.F., S.C. and F.M.; organic synthesis, S.C. and M.I.; NMR and HPLC analyses, S.C. and D.C.; crystallography: S.C., C.C. and F.M.; computational studies, E.M.A.F. and G.G.; writing—original draft preparation, P.G., P.F., S.C., E.M.A.F., G.G. and F.M.; writing—review and editing, P.G., P.F., S.C., E.M.A.F., G.G. and F.M.; supervision, P.G., P.F. and F.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

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 all compounds mentioned in this manuscript are available from the authors.

References

  1. Lamb, Y.N. Elagolix: First Global Approval. Drugs 2018, 78, 1501–1508. [Google Scholar] [CrossRef] [PubMed]
  2. Ali, M.; Sara, A.R.; Al Hendy, A. Elagolix in the treatment of heavy menstrual bleeding associated with uterine fibroids in premenopausal women. Expert Rev. Clin. Pharmacol. 2021, 14, 427–437. [Google Scholar] [CrossRef] [PubMed]
  3. Betz, S.F.; Zhu, Y.-F.; Chen, C.; Struthers, R.S. Non-peptide gonadotropin-releasing hormone receptor antagonists. J. Med. Chem. 2008, 51, 3331–3348. [Google Scholar] [CrossRef] [PubMed]
  4. Tukun, F.-L.; Olberg, D.E.; Riss, P.J.; Haraldsen, I.; Kaass, A.; Klaveness, J. Recent Development of Non-Peptide GnRH Antagonists. Molecules 2017, 22, 2188. [Google Scholar] [CrossRef]
  5. Tucci, F.C.; Hu, T.; Mesleh, M.F.; Bokser, A.; Allsopp, E.; Gross, T.D.; Guo, Z.; Zhu, Y.-F.; Struthers, R.S.; Ling, N.; et al. Atropisomeric property of 1-(2,6-difluorobenzyl)-3-(2R)-amino-2-phenethyl-5-(2-fluoro-3-methoxyphenyl)-6-methyluracil. Chirality 2005, 17, 559–564. [Google Scholar] [CrossRef]
  6. LaPlante, S.R.; Edwards, P.J.; Fader, L.D.; Jakalian, A.; Hucke, O. Revealing atropisomer axial chirality in drug discovery. ChemMedChem 2011, 6, 505–513. [Google Scholar] [CrossRef]
  7. Basilaia, M.; Chen, M.H.; Secka, J.; Gustafson, J.L. Atropisomerism in the Pharmaceutically Relevant Realm. Acc. Chem. Res. 2022, 55, 2904–2919. [Google Scholar] [CrossRef]
  8. Toenjes, S.T.; Gustafson, J.L. Atropisomerism in medicinal chemistry: Challenges and opportunities. Future Med. Chem. 2018, 10, 409–422. [Google Scholar] [CrossRef]
  9. Clayden, J.; Moran, W.J.; Edwards, P.J.; LaPlante, S.R. The challenge of atropisomerism in drug discovery. Angew. Chem. Int. Ed Engl. 2009, 48, 6398–6401. [Google Scholar] [CrossRef]
  10. Wang, J.; Zeng, W.; Li, S.; Shen, L.; Gu, Z.; Zhang, Y.; Li, J.; Chen, S.; Jia, X. Discovery and Assessment of Atropisomers of (±)-Lesinurad. ACS Med. Chem. Lett. 2017, 8, 299–303. [Google Scholar] [CrossRef]
  11. Wang, H.; Liu, Y.; Huai, Q.; Cai, J.; Zoraghi, R.; Francis, S.H.; Corbin, J.D.; Robinson, H.; Xin, Z.; Lin, G.; et al. Multiple conformations of phosphodiesterase-5: Implications for enzyme function and drug development. J. Biol. Chem. 2006, 281, 21469–21479. [Google Scholar] [CrossRef]
  12. Guo, Z.; Chen, Y.; Huang, C.Q.; Gross, T.D.; Pontillo, J.; Rowbottom, M.W.; Saunders, J.; Struthers, S.; Tucci, F.C.; Xie, Q.; et al. Uracils as potent antagonists of the human gonadotropin-releasing hormone receptor without atropisomers. Bioorg. Med. Chem. Lett. 2005, 15, 2519–2522. [Google Scholar] [CrossRef]
  13. Wang, Z.; Meng, L.; Liu, X.; Zhang, L.; Yu, Z.; Wu, G. Recent progress toward developing axial chirality bioactive compounds. Eur. J. Med. Chem. 2022, 243, 114700. [Google Scholar] [CrossRef]
  14. Zhao, L.; Guo, Z.; Chen, Y.; Hu, T.; Wu, D.; Zhu, Y.-F.; Rowbottom, M.; Gross, T.D.; Tucci, F.C.; Struthers, R.S.; et al. 5-Aryluracils as potent GnRH antagonists-Characterization of atropisomers. Bioorg. Med. Chem. Lett. 2008, 18, 3344–3349. [Google Scholar] [CrossRef]
  15. Yan, W.; Cheng, L.; Wang, W.; Wu, C.; Yang, X.; Du, X.; Ma, L.; Qi, S.; Wei, Y.; Lu, Z.; et al. Structure of the human gonadotropin-releasing hormone receptor GnRH1R reveals an unusual ligand binding mode. Nat. Commun. 2020, 11, 5287. [Google Scholar] [CrossRef]
  16. Guo, Z.; Chen, Y.; Wu, D.; Chen, C.; Wade, W.; Dwight, W.J.; Huang, C.Q.; Tucci, F.C. Preparation of Pyrimidine-2,4(1H,3H)-dione Derivatives as Gonadotropin-Releasing Hormone Receptor Antagonists. WO Patent WO2005007165, 27 January 2005. [Google Scholar]
  17. Chen, C.; Wu, D.; Guo, Z.; Xie, Q.; Reinhart, G.J.; Madan, A.; Wen, J.; Chen, T.; Huang, C.Q.; Chen, M.; et al. Discovery of sodium R-(+)-4-{2-5-(2-fluoro-3-methoxyphenyl)-3-(2-fluoro-6-trifluoromethylbenzyl)-4-methyl-2,6-dioxo-3,6-dihydro-2H-pyrimidin-1-yl-1-phenylethylamino}butyrate (elagolix), a potent and orally available nonpeptide antagonist of the human gonadotropin-releasing hormone receptor. J. Med. Chem. 2008, 51, 7478–7485. [Google Scholar] [CrossRef]
  18. Peddireddy, S.R.; Allam, S.K.; Kottur, M.K.; Oruganti, S.; Kandagatla, B. Process for the Preparation of Elagolix Sodium and Its Polymorph. WO Patent WO2017221144, 18 December 2017. [Google Scholar]
  19. Vlasakova, R.; Cerna, I.; Obadalova, I.; Krejcik, L.; Dammer, O.; Svobodova, J.; Sembera, F. Solid Forms of Elagolix. Wo Patent WO2018224063, 13 December 2018. [Google Scholar]
  20. Sulake, R.S.; Shinde, S.R.; Siyan, R.S.; Bhise, N.B.; Singh, G.P. Process for the Preparation of Elagolix and Pharmaceutically Acceptable Salts Thereof. WO Patent WO2018198086, 1 November 2018. [Google Scholar]
  21. Lenna, R.; Fasana, A.; Ortiz, J. Process for the Preparation of the Sodium Salt of Elagolix and Its Intermediates. Wo Patent WO2021044230, 11 March 2022. [Google Scholar]
  22. Siripragada, M.R.; Pendyam, K.; Kallepally, S.; Avula, S.R.; Gottapu, V.N.; Pappula, V.R. An Improved Process for the Preparation of Elagolix Sodium. WO Patent WO2021064561, 8 April 2021. [Google Scholar]
  23. Ballete, R.; Jimenez Alonso, O.; Garcia Garcia, E.; Dobarro Rodriguez, A. 3-((R)-2-(Amino-2-phenylethyl)-1-(2-fluoro-6-trifluoromethylbenzyl)-5-iodo-6-methyl-1H-pyrimidine-2,4-dione or a Salt Thereof, Process for its Preparation, and Its in the Synthesis of Elagolix. WO Patent WO2021083554, 6 May 2021. [Google Scholar]
  24. Zhong, X.; Lv, Q.; Yong, Q.; Hu, W.; Li, D.; Ji, S.; Zhan, L.; Chen, W.; Li, M.; Lin, J.; et al. Forced degradation studies of elagolix sodium with the implementation of high resolution LC-UV-PDA-MSn (n = 1,2,3…) and NMR structural elucidation. J. Pharm. Biomed. Anal. 2023, 224, 115198. [Google Scholar] [CrossRef]
  25. Ferraboschi, P.; Chiara Sala, M.; Stradi, R.; Ragonesi, L.; Gagliardi, C.; Lanzarotti, P.; Ragg, E.M.; Mori, M.; Meneghetti, F. Full spectroscopic characterization of two crystal pseudopolymorphic forms of the antiandrogen cortexolone 17α-propionate for topic application. Steroids 2017, 128, 95–104. [Google Scholar] [CrossRef]
  26. Meneghetti, F.; Ferraboschi, P.; Grisenti, P.; Reza Elahi, S.; Mori, M.; Ciceri, S. Crystallographic and NMR Investigation of Ergometrine and Methylergometrine, Two Alkaloids from Claviceps purpurea. Molecules 2020, 25, 331. [Google Scholar] [CrossRef]
  27. Ciceri, S.; Colombo, D.; Ferraboschi, P.; Grisenti, P.; Iannone, M.; Mori, M.; Meneghetti, F. Vecuronium bromide and its advanced intermediates: A crystallographic and spectroscopic study. Steroids 2021, 176, 108928. [Google Scholar] [CrossRef]
  28. Gallagher, D.; Treiber, L.; Hughes, R.; Campopiano, O.; Wang, P.; Zhao, Y.; Chou, S.; Ouellette, M.; Hettinger, D. Processes for the Preparation of Uracil Derivatives. WO Patent WO2009062087, 14 May 2009. [Google Scholar]
  29. Larrañaga, M.D.; Lewis, R.J.; Lewis, R.A.; Hawley, G.G. Hawley’s Condensed Chemical Dictionary, 16th ed.; John Wiley & Sons, Inc: Hoboken, NJ, USA, 2016; ISBN 1119312469. [Google Scholar]
  30. Lewis, R.J. Sax’s Dangerous Properties of Industrial Materials, 12th ed.; Wiley: Hoboken, NJ, USA, 2012; ISBN 9780471701347. [Google Scholar]
  31. Badertscher, M.; Bühlmann, P.; Pretsch, E. Tables of Spectral Data for Structure Determination of Organic Compounds; Springer: Berlin/Heidelberg, Germany, 1983; ISBN 978-3-662-22455-7. [Google Scholar]
  32. Marek, R.; Lycka, A. 15N NMR Spectroscopy in Structural Analysis. COC 2002, 6, 35–66. [Google Scholar] [CrossRef]
  33. D’Acquarica, I.; Gasparrini, F.; Pierini, M.; Villani, C.; Zappia, G. Dynamic HPLC on chiral stationary phases: A powerful tool for the investigation of stereomutation processes. J. Sep. Sci. 2006, 29, 1508–1516. [Google Scholar] [CrossRef] [PubMed]
  34. Trapp, O.; Schoetz, G.; Schurig, V. Determination of enantiomerization barriers by dynamic and stopped-flow chromatographic methods. Chirality 2001, 13, 403–414. [Google Scholar] [CrossRef] [PubMed]
  35. Wolf, C. Stereolabile chiral compounds: Analysis by dynamic chromatography and stopped-flow methods. Chem. Soc. Rev. 2005, 34, 595–608. [Google Scholar] [CrossRef] [PubMed]
  36. König, W.A.; Gehrcke, B.; Runge, T.; Wolf, C. Gas chromatographic separation of atropisomeric alkylated and polychlorinated biphenyls using modified cyclodextrins. J. High Resol. Chromatogr. 1993, 16, 376–378. [Google Scholar] [CrossRef]
  37. Farrugia, L.J. WinGX and ORTEP for Windows: An update. J. Appl. Crystallogr. 2012, 45, 849–854. [Google Scholar] [CrossRef]
  38. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A Gen. Phys. 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  39. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  40. Guo, Z.; Zhu, Y.-F.; Gross, T.D.; Tucci, F.C.; Gao, Y.; Moorjani, M.; Connors, P.J.; Rowbottom, M.W.; Chen, Y.; Struthers, R.S.; et al. Synthesis and structure-activity relationships of 1-arylmethyl-5-aryl-6-methyluracils as potent gonadotropin-releasing hormone receptor antagonists. J. Med. Chem. 2004, 47, 1259–1271. [Google Scholar] [CrossRef]
  41. 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]
  42. Lan, C.B.; Auclair, K. Ammonium Chloride-Promoted Rapid Synthesis of Monosubstituted Ureas under Microwave Irradiation. Eur. J. Org. Chem. 2021, 2021, 5135–5146. [Google Scholar] [CrossRef]
  43. Capasso, C.; Winum, J.-Y. Novel method of treating macular degeneration: A patent evaluation (WO2018/107005). Expert Opin. Ther. Pat. 2019, 29, 749–752. [Google Scholar] [CrossRef]
  44. O’Brien, P.M.; Sliskovic, D.R.; Blankley, C.J.; Roth, B.D.; Wilson, M.W.; Hamelehle, K.L.; Krause, B.R.; Stanfield, R.L. Inhibitors of acyl-CoA:cholesterol O-acyl transferase (ACAT) as hypocholesterolemic agents. 8. Incorporation of amide or amine functionalities into a series of disubstituted ureas and carbamates. Effects on ACAT inhibition in vitro and efficacy in vivo. J. Med. Chem. 1994, 37, 1810–1822. [Google Scholar] [CrossRef]
  45. Liu, Y.; Hao, Q.; Lin, K.; Zhou, W.; Pan, J.; Chen, L.; Zhou, T. Preparation of Elagolix Intermediate. China Patent CN110498770, 30 April 2021. [Google Scholar]
  46. SMART & SAINT Software Reference Manual, Version 6.45; Bruker Analytical X-Ray Systems, Inc.: Madison, WI, USA, 2003.
  47. Sheldrick, G.M. SADABS, Version 2008/1; Bruker AXS Inc.: Karlsruhe, Germany, 2008.
  48. Sheldrick, G.M. SHELXL-2018; Universität Göttingen: Göttingen, Germany, 2018. [Google Scholar]
  49. Macrae, C.F.; Sovago, I.; Cottrell, S.J.; Galek, P.T.A.; McCabe, P.; Pidcock, E.; Platings, M.; Shields, G.P.; Stevens, J.S.; Towler, M.; et al. Mercury 4.0: From visualization to analysis, design and prediction. J. Appl. Crystallogr. 2020, 53, 226–235. [Google Scholar] [CrossRef]
  50. Nardelli, M. Parst: A system of fortran routines for calculating molecular structure parameters from results of crystal structure analyses. Comput. Chem. 1983, 7, 95–98. [Google Scholar] [CrossRef]
Molecules 28 03861 g001 550
Figure 1. Chemical structures of 1 and its analogue 2; the atropisomerism in 2 arises from hindrance to rotation about the single bond indicated by the red arrow.
Figure 1. Chemical structures of 1 and its analogue 2; the atropisomerism in 2 arises from hindrance to rotation about the single bond indicated by the red arrow.
Molecules 28 03861 g001
Molecules 28 03861 sch001 550
Scheme 1. Reagents and conditions: (i) urea, HCl, H2O, reflux; (ii) (a) tert-butyl acetoacetate, toluene, reflux, Dean-Stark apparatus; (b) p-TsOH·H2O, reflux; (iii) Br2, acetic acid; (iv) 2-fluoro-3-methoxyphenylboronic acid, Pd(OAc)2, [(t-Bu)3PH]BF, KOH, H2O, acetone, under argon atmosphere, reflux; (v) (a) 10, K2CO3, DMF; (b) MsOH; (vi) (a) ethyl 4-bromobutyrate, DIPEA, DMF, 60 °C; (b) filtration on silica; (vii) (a) NaOH, H2O, ethanol; (b) extraction with methyl isobutyl ketone (MIBK); (c) precipitation from heptane.
Scheme 1. Reagents and conditions: (i) urea, HCl, H2O, reflux; (ii) (a) tert-butyl acetoacetate, toluene, reflux, Dean-Stark apparatus; (b) p-TsOH·H2O, reflux; (iii) Br2, acetic acid; (iv) 2-fluoro-3-methoxyphenylboronic acid, Pd(OAc)2, [(t-Bu)3PH]BF, KOH, H2O, acetone, under argon atmosphere, reflux; (v) (a) 10, K2CO3, DMF; (b) MsOH; (vi) (a) ethyl 4-bromobutyrate, DIPEA, DMF, 60 °C; (b) filtration on silica; (vii) (a) NaOH, H2O, ethanol; (b) extraction with methyl isobutyl ketone (MIBK); (c) precipitation from heptane.
Molecules 28 03861 sch001
Molecules 28 03861 sch002 550
Scheme 2. Reaction of intermediate 4 with (i) ethyl aceto-acetate and (ii) diketene-acetone adduct.
Scheme 2. Reaction of intermediate 4 with (i) ethyl aceto-acetate and (ii) diketene-acetone adduct.
Molecules 28 03861 sch002
Molecules 28 03861 g002 550
Figure 2. NOESY experiments of compound 5 (A) and its isomer (B).
Figure 2. NOESY experiments of compound 5 (A) and its isomer (B).
Molecules 28 03861 g002
Molecules 28 03861 sch003 550
Scheme 3. Formation of the main by-product 11, observed when the literature procedure [11] was applied.
Scheme 3. Formation of the main by-product 11, observed when the literature procedure [11] was applied.
Molecules 28 03861 sch003
Molecules 28 03861 g003 550
Figure 3. (A) H-21, (B) H-22, and (C) H-15 signals of 1. The integration values are reported in blue. The H-15 signal is overlapped with H-32 signals.
Figure 3. (A) H-21, (B) H-22, and (C) H-15 signals of 1. The integration values are reported in blue. The H-15 signal is overlapped with H-32 signals.
Molecules 28 03861 g003
Molecules 28 03861 g004 550
Figure 4. HPLC chromatogram of 1 performed on (A) RP-C18 and (B) chiral stationary phase.
Figure 4. HPLC chromatogram of 1 performed on (A) RP-C18 and (B) chiral stationary phase.
Molecules 28 03861 g004
Molecules 28 03861 g005 550
Figure 5. Chiral HPLC chromatogram of 7.
Figure 5. Chiral HPLC chromatogram of 7.
Molecules 28 03861 g005
Molecules 28 03861 g006 550
Figure 6. ORTEP [37] diagram of 7 with an arbitrary atom-numbering scheme. Thermal ellipsoids were drawn at the 40% probability level.
Figure 6. ORTEP [37] diagram of 7 with an arbitrary atom-numbering scheme. Thermal ellipsoids were drawn at the 40% probability level.
Molecules 28 03861 g006
Molecules 28 03861 g007 550
Figure 7. Graphical representation of the crystal packing along a axis, showing the H-bonds as dotted blue lines.
Figure 7. Graphical representation of the crystal packing along a axis, showing the H-bonds as dotted blue lines.
Molecules 28 03861 g007
Molecules 28 03861 g008 550
Figure 8. (A). Simplified structure of 1 used for conformational analysis studies. The yellow arrows highlight the dihedral angles investigated by means of the conformational analysis. (B,C) Scan plots of the τ1–τ4 torsional angles.
Figure 8. (A). Simplified structure of 1 used for conformational analysis studies. The yellow arrows highlight the dihedral angles investigated by means of the conformational analysis. (B,C) Scan plots of the τ1–τ4 torsional angles.
Molecules 28 03861 g008
Molecules 28 03861 g009 550
Figure 9. Coordinate scan plot related to the τ5 torsional angle variation and low energy conformers, resulting from the data reported in Table 7.
Figure 9. Coordinate scan plot related to the τ5 torsional angle variation and low energy conformers, resulting from the data reported in Table 7.
Molecules 28 03861 g009
Table
Table 1. 1H NMR chemical shifts (ppm) a and coupling constants (Hz) b of compounds 5–9 and 1.
Table 1. 1H NMR chemical shifts (ppm) a and coupling constants (Hz) b of compounds 5–9 and 1.
1H567891
1//////
2//////
38.97 (brs)Exchanged with CD3OD8.52 (brs)///
4//////
55.59 (s)/////
6//////
75.36 (s)5.38 (s)5.47 (s)5.50 (m)5.47 (m)5.43 (m)
8//////
9//////
107.25 (dd, J = 12.3 and 8.5 Hz)7.20 (dd, J= 11.7 and 8.3 Hz)7.28 (dd, J = 11.8 and 8.3 Hz)7.26 (m, overlapped with 28)7.29 (m, overlapped with 28)7.44 (m)
117.42 (m)7.37 (m)7.42 (m)7.41 (m, overlapped with 26 and 26’)7.41 (m)7.53 (td, J = 8.1, 5.0 Hz)
127.54 (d, J = 7.9 Hz)7.48 (d, J = 7.8 Hz)7.54 (d, J = 7.9 Hz)7.55 (d, J = 7.9 Hz)7.54 (d, J = 7.9 Hz)7.62 (d, J = 7.9 Hz)
13//////
14//////
152.15 (s)2.35 (s)2.05 (s)2.07 (s)2.07 and 2.06 *2.085 (s) and 2.076 (s) *
16//////
17//////
18//////
19//6.97 (m)6.97 (td, J = 8.2 and 1.3 Hz)6.97 (tt, J = 8.1 and 1.6 Hz)7.10 (n.d., overlapped with 20)
20//7.10 (dd, J = 9.0 and 7.9 Hz)7.11 (tdd, J = 8.0, 2.8, and 1.3 Hz)7.11 (tdd, J = 8.0, 4.6, and 1.3 Hz)7.11 and 7.14 (n.d., overlapped with 19)
21//6.81 (m)6.85 (td, J = 6.2 and 1.5 Hz) and 6.78 (td, J = 6.2 and 1.4 Hz) *6.83 (ddd, J = 7.7, 6.0, and 1.6 Hz) and 6.76 (ddd, J = 7.7, 6.0, and 1.6 Hz) *6.76 (m) and 6.61 (m) *
22//3.88 (s)3.89 (s)3.889 (s) and 3.886 (s) *3.883 (s) and 3.880 (s) *
23///4.29 (m o ddd, J = 16.2, 13.1, 9.8 Hz, Ha) and 4.11 (m, Hb)4.28 (m, Ha) and 4.04(m, Hb)4.26–4.07 (m)
24///4.42 (dd, J = 9.8 and 4.5 Hz)4.12 (m)4.11 (n.d.)
25//////
26 and 26′///7.41 (m, overlapped with 11)7.37 (m)7.29 (n.d., overlapped with 28)
27 and 27′///7.32 (m)7.30 (m, overlapped with 10)7.27 (n.d., overlapped with 26 and 26’)
28///7.27 (m, overlapped with 10)7.24 (m)7.22 (m)
29///2.131.63 (brs)/
30////2.46 (m, Ha) and 2.38 (m, Hb)2.40 (m)
31////1.68 (m)1.75 (m, Ha) and 1.68 (m, Hb)
32////2.29 (m)2.09 (m, overlapped with 15)
33//////
CH2CH3////4.084 (q, J = 7.1Hz,) and 4.082 (q, J = 7.1Hz) */
CH2CH3////1.21 (t, J = 7.1 Hz)/
nd = J(H,H) were not determined due to the overlapping. a Assignments from 1H–1H COSY, HSQC, and HMBC data in CDCl3 (5 and 7–9), CD3OD (1), and CDCl3/CD3OD 9:1 mixture (6) at 298K. b Coupling constants were obtained by direct inspection of the spectra. Experimental error in the measured 1H-1H coupling constants was ±0.5 Hz. * The two signals are due to the presence of two atropisomers.
Table
Table 2. 13C NMR chemical shifts (ppm) a and 13C-19F coupling constants (Hz) b data of 5–9 and 1.
Table 2. 13C NMR chemical shifts (ppm) a and 13C-19F coupling constants (Hz) b data of 5–9 and 1.
13C567891
1//////
2151.7150.8151.0152.3152.3 and 152.2 *153.3 and 153.2 *
3//////
4162.4159.6161.2161.8 and 161.6 *161.7 and 161.6 *163.5 and 163.3 *
5102.599.2108.8108.3 and 108.2 *108.3 and 108.2 *109.3 and 109.2 *
6153.6151.6151.9149.9 and 149.8 *149.7 and 149.6 * 152.6 and 152.5 *
741.1 (m)42.5 (m)41.8 (m)42.7 (brs)42.7 (brs)44.1 (brs)
8121.8 (d, J = 11.6 Hz)121.2 (d, J = 11.3 Hz)121.8 (d, J = 11.1 Hz)122.0 (d, J= 11.4 Hz)122.1 (d, J = 11.3 Hz)123.6 (overlapped with 12 and 16)
9161.4 (d, J = 250.0 Hz)161.1 (d, J = 249.6 Hz)161.4 (d, J = 250.0 Hz)161.3 (d, J = 249.5 Hz)161.3 (d, J = 250.0 Hz)162.8 (d, J =248.0 Hz)
10120.9 (d, J = 24.0 Hz)120.9 (d, J = 24.0 Hz)121.0 (d, J = 24.4 Hz)120.9 (d, J = 24.1 Hz)120.9 (d, J = 24.4 Hz)122.15 (d, J = 24.3 Hz) and 122.19 (d, J = 24.3 Hz) *
11129.5 (d, J = 10.0 Hz)129.6 (d, J = 9.7 Hz)129.5 (d, J = 9.8 Hz)129.3 (d, J = 9.8 Hz)129.3 (d, J = 9.6 Hz)130.9 and 130.8 *
12122.6 (m)122.5 (m)122.6 (m)122.6 (m)122.5 (m)123.7 (overlapped with 8 and 16)
13129.6 (dd, J = 30.0 and 3.6 Hz; partially hidden by 11)129.3 (dd, J = 30.9 and 3.7 Hz; 129.6 (dd, J = 30.9 and 4.0 Hz)129.5 (dd, J = 31.0 and 3.1, partially overlapped with 11)129.3 (dd, J = 30.2 and 3.0 Hz)130.6 (dd, J = 31.1 and 3.3 Hz)
14123.4 (dd, J = 273.6 and 4.0 Hz)123.3 (dd, J = 274.3 and 4.2 Hz)123.5 (dd, J = 273.8 and 3.9 Hz)123.5 (dd, J = 273.6 and 3.8 Hz)123.5 (dd, J = 271.1 and 3.7 Hz)125.1 (dd, J = 273.3 and 3.8 Hz)
1520.120.117.917.817.818.1
16//121.5 (d, J = 13.6 Hz)122.3 (d, J = 13.7 Hz) and 122.2 (d, J = 13.7 Hz) *122.37 (d, J = 13.8 Hz) and 122.35 (d, J = 13.8 Hz)123.8 (overlapped with 8 and 12)
17//149.9 (d, J = 246.4 Hz)150.2 (d, J = 246.3 Hz) and 150.1 (d, J = 246.3 Hz) *150.23 (d, J = 245.5) and 150.17 (d, J = 245.5)151.8 (d, J = 245.4 Hz)
18//148.0 (d, J = 11.0 Hz)148.0 (d, J = 11.0 Hz)148.0 (d, J = 11.0 Hz)149.5 (d, J = 11.2 Hz)
19//113.5113.3113.2114.8
20//123.8123.8 (overlapped with 21)123.80 (s) and 123.76 (s) *125.1 (s) and 125.0 (s) *
21//123.9 (d, J = 4.6 Hz)123.9 and 123.8 (overlapped with 20) *123.92 (s) and 123.85 (s) *125.0 and 124.9 *
22//56.3 56.256.256.8
23///49.10 and 49.14 *47.7 and 47.6 *48.1 and 47.9 *
24///54.1 and 54.2 *60.9 and 60.8 *62.2 and 62.0 *
25///143.7 and 143.6 *141.92 (s) and 141.89 (s) *142.1 and 141.9 *
26 and 26’///126.4 and 126.3 *127.2 and 127.1 (overlapped with 28) *128.7
27 and 27’///128.5 and 128.4 *128.4 and 128.3 *129.4
28///127.3 and 127.2 *127.2 (overlapped with 26 and 26′)128.5
29//////
30////46.51 (s) and 46.47 (s) *48.7
31///25.427.86 (s) and 27.83 (s) *
32////32.037.03 (s) and 37.00 (s) *
33////173.7
60.1		182.3
CH2CH3/////
CH2CH3////14.2/
a Assignments from HSQC and HMBC data in CDCl3 (5 and 7–9), CD3OD (1), and CDCl3/CD3OD 9:1 mixture (6) at 298K. b Coupling constants were obtained by direct inspection of the spectra. Experimental error in the measured 1H-1H coupling constants was ±0.5 Hz. * The two signals are due to the presence of two atropisomers.
Table
Table 3. 15N NMR chemical shifts (ppm) a of compounds 5–9 and 1.
Table 3. 15N NMR chemical shifts (ppm) a of compounds 5–9 and 1.
15N567891
1136.7137.9134.3134.3135.5137.1
3154.3 b153.4 b151.9 b160.1161.3160.5
29///31.243.646.3
a Assignments from 1H-15N HMBC data in CDCl3 (5 and 7–9), CD3OD (1), and CDCl3/CD3OD 9:1 mixture (6) at 298K using nitromethane as external reference (neat nitromethane = 381.7 ppm as reported in Ref. [32]). b Assignments from 1H-15N HSQC data in CDCl3 (5 and 7), and CDCl3 + 3% of DMSO-d6 (6) at 298K using ammonia as external reference (ammonia = 0.0 ppm as reported in Ref. [32]).
Table
Table 4. 1H-13C gs-HMBC long range couplings observed in the spectra of 1.
Table 4. 1H-13C gs-HMBC long range couplings observed in the spectra of 1.
13CHMBC (13C→1H)
C-2H-7 and H-23
C-4H-23 and H-15 (weak)
C-5H-21 and H-15
C-6H-15
C-8H-7, H-10, and H-12
C-9H-7, H-10, and H-11
C-13H-7 and H-11
C-14H-12
C-16H-20
C-17H-19 and H-21
C-18H-20 and H-22
C-25H-23, H-24, H-27, and H-27’
C-33H-31 and H-32
Table
Table 5. Geometrical features, relative energy values and equilibrium percentages of the predicted conformations in the gas phase, as resulting after geometry optimization.
Table 5. Geometrical features, relative energy values and equilibrium percentages of the predicted conformations in the gas phase, as resulting after geometry optimization.
Conformersτ1 (°)τ2 (°)τ3 (°)τ4 (°)Gas Phase ΔE (kcal/mol)Gas Phase (%)
1A−81.6−48.5−132.049.14.350.0
1B−70.0147.8−130.049.76.840.0
1C96.0−54.0−125.070.00.6720.5
1D111.8134.3−133.052.54.600.0
1E−80.5−47.0−61.3−52.44.670.0
1F−69.4150.3−55.0−55.57.090.0
1G113.0136.0−60.5−55.65.010.0
1H93.2−62.1−64.1−55.61.902.6
1I94.5−53.661.3−142.62.261.4
1J118.5147.459.4−140.84.860.0
1K−74.0−40.367.8−141.24.780.0
1L−69.1148.960.5−142.93.690.1
1M92.6−61.9−110.0−134.50.0063.9
1N113.6146.5−77.5147.43.980.1
1O−85.9−51.6−115.4−136.43.580.2
1P−70.9148.7−75.9146.35.330.0
1Q−81.0−46.4−75.4142.83.640.1
1R−71.1148.9−75.2147.05.400.0
1S92.7−55.0−76.1139.91.0910.1
1T113.9146.9−77.6147.53.930.1
1U−86.3−49.956.053.04.530.0
1V−70.1145.951.155.23.670.1
1W93.5−53.549.953.82.770.6
1X114.2142.751.353.65.320.0
Table
Table 6. Geometrical features, relative energies and equilibrium percentages after geometry optimization of the most populated conformations of a simplified structure of 1.
Table 6. Geometrical features, relative energies and equilibrium percentages after geometry optimization of the most populated conformations of a simplified structure of 1.
Codeτ1τ2τ3τ4τ5Water ΔE (kcal/mol)Water (%)
7
(crystal structure)		 −105.5−144.7110.4
1M−A87.7°−63.3°−104.7°−133.9°−111.8°0.0083.6
1M−C91.3°−60.2°−106.7°−132.6°75.2°2.072.6
1M−D92.9°−60.3°−110.0°−135.0°126.3°1.0613.8
Table
Table 7. Geometrical features, relative energies and equilibrium percentages after geometry optimization of the most populated conformations of 1M-A and 1M-D, created rotating the dihedral angle involving in the methoxy group (τ6).
Table 7. Geometrical features, relative energies and equilibrium percentages after geometry optimization of the most populated conformations of 1M-A and 1M-D, created rotating the dihedral angle involving in the methoxy group (τ6).
Conformert6Distance
CH22 – H19 (Å)		Gas Phase ΔE (kcal/mol)Gas Phase (%)Methanol ΔE (kcal/mol)Methanol (%)
1M-A02.560.0044.5086.0
1M-A (I)−663.600.2927.51.625.6
1M-A (II)+663.600.2728.01.388.4
1M-D02.560.2230.50.0084.5
1M-D (I)−663.600.3225.61.368.5
1M-D (II)+663.600.0043.91.487.0
Table
Table 8. Crystal data and structural parameters of 7.
Table 8. Crystal data and structural parameters of 7.
Crystal Data
Chemical formulaC20H15F5N2O3
Mr426.34
Crystal system, space groupOrthorhombic, P bca
a, b, c (Å)11.165 (2), 11.073 (2), 30.367(6)
V (Å3)3754.2(13))
Z8
F (000)1744
Density (g/cm3)1.509
Temperature (K)298
Radiation typeMo-Kα (λ = 0.71073 Å)
µ (mm−1)1.135
Crystal size (mm)0.06 × 0.03 × 0.02
Data collection
DiffractometerBruker-Axs Smart-Apex CCD
Tmin, Tmax0.893, 1.000
No. of measured, independent and observed [I > 2σ(I)] reflections29,138, 4175, 1715
Rint0.0298
Structure refinement
R, wR2, S0.0437 (I > 2σ(I)) and 0.1119 (all), 0.0973 (I > 2σ(I)) and 0.1126 (all), 0.760 (all)
No. of parameters271
No. of restraints0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ciceri, S.; Colombo, D.; Fassi, E.M.A.; Ferraboschi, P.; Grazioso, G.; Grisenti, P.; Iannone, M.; Castellano, C.; Meneghetti, F. Elagolix Sodium Salt and Its Synthetic Intermediates: A Spectroscopic, Crystallographic, and Conformational Study. Molecules 2023, 28, 3861. https://doi.org/10.3390/molecules28093861

AMA Style

Ciceri S, Colombo D, Fassi EMA, Ferraboschi P, Grazioso G, Grisenti P, Iannone M, Castellano C, Meneghetti F. Elagolix Sodium Salt and Its Synthetic Intermediates: A Spectroscopic, Crystallographic, and Conformational Study. Molecules. 2023; 28(9):3861. https://doi.org/10.3390/molecules28093861

Chicago/Turabian Style

Ciceri, Samuele, Diego Colombo, Enrico M. A. Fassi, Patrizia Ferraboschi, Giovanni Grazioso, Paride Grisenti, Marco Iannone, Carlo Castellano, and Fiorella Meneghetti. 2023. "Elagolix Sodium Salt and Its Synthetic Intermediates: A Spectroscopic, Crystallographic, and Conformational Study" Molecules 28, no. 9: 3861. https://doi.org/10.3390/molecules28093861

APA Style

Ciceri, S., Colombo, D., Fassi, E. M. A., Ferraboschi, P., Grazioso, G., Grisenti, P., Iannone, M., Castellano, C., & Meneghetti, F. (2023). Elagolix Sodium Salt and Its Synthetic Intermediates: A Spectroscopic, Crystallographic, and Conformational Study. Molecules, 28(9), 3861. https://doi.org/10.3390/molecules28093861

Article Metrics

Back to TopTop