Next Article in Journal
Structural Characterization, Cytotoxicity, and the Antifungal Mechanism of a Novel Peptide Extracted from Garlic (Allium sativa L.)
Previous Article in Journal
Impact of Ultrasonication on the Self-Assembly Behavior and Gel Properties of Bovine Bone Collagen I
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Insights and Intermolecular Energy for Some Medium and Long-Chain Testosterone Esters

1
National Institute For R&D of Isotopic and Molecular Technologies, 67-103 Donat, 400293 Cluj-Napoca, Romania
2
Department of Physics and Chemistry, Technical University of Cluj-Napoca, 400114 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(7), 3097; https://doi.org/10.3390/molecules28073097
Submission received: 21 February 2023 / Revised: 17 March 2023 / Accepted: 28 March 2023 / Published: 30 March 2023

Abstract

:
Testosterone (17β-Hydroxyandrost-4-en-3-one) is the primary male anabolic-androgenic steroid. The crystal structures of two medium and two long esterified forms of testosterone, including enanthate, cypionate, decanoate and undecanoate, were determined by X-ray single crystal diffraction. The samples were also characterized by powder X-ray diffraction, FT-IR spectroscopy and thermal analysis (DTA, TG). Crystal packings and supramolecular features were described. The analysis of structural features was accomplished by computational methods in terms of the type of intermolecular interactions, crystal energies and Hirshfeld surfaces analysis. From a pharmaceutical point of view, the solubility of compounds was investigated.

Graphical Abstract

1. Introduction

Testosterone (17β-Hydroxyandrost-4-en-3-one) is chemically a cholesterol derivative and the primary naturally occurring anabolic steroid. It can be regarded as an androstane derivative and the primary male sex hormone. It fulfills an important role in the development and maintenance of the male reproductive system, promoting the appearance of secondary male characteristics, including increased lean muscle mass, bone mass, body hair growth [1], and prevents osteoporosis [2,3]. The hormone is also present in the body of women in much smaller quantities than in men, being associated with certain physiological roles such as sexual function, musculoskeletal, cardiovascular and vulvovaginal health and cognitive function [4]. As an androgen receptor agonist, it exerts the anabolic-androgenic properties which are common to all of its derivatives/analogues which belong to this family [5]. Medically, testosterone is used to relieve symptoms related to low testosterone in men, breast cancer in women and in hormone therapy for transgender men [6]. Targeting androgen receptors has been shown to have a positive overall impact on health and the state of well-being, and it triggers an increased protein synthesis in certain tissues which contain androgen receptors, among which are bones and muscle tissues [7]. Optimal testosterone levels in elderly men induce an overall positive impact on general health; they help to decrease body and visceral fat mass, help increase lean muscle mass, improve cholesterol levels, and enhance carbohydrate metabolism [8]. Anabolic agents such as testosterone and its analogues/derivatives are often used by athletes, weightlifters and bodybuilders in order to boost performance and aid muscle recovery [9]. It can also be used to relieve or treat protein breakdown in catabolic conditions [10]. It should be emphasized that testosterone is used illicitly by athletes and its use is prohibited by the World Anti-Doping Agency.
Testosterone (Figure 1a) is characterized by a short half-life (a few hours), being subject to esterification, which leads to a considerable increase in the half-life depending on the length of the ester (of the carbon chain). Esterified forms offer the advantage of increasing the half-life when administered by intramuscular or subcutaneous injection, thus avoiding the need for daily administration [11]. Various esterified forms of testosterone show half-lives from less than 24 h for testosterone acetate [12] up to 36 days for testosterone undecanoate [13].
The current paper approaches and characterizes two medium and two long esterified testosterone prodrugs:
(i)
Testosterone enanthate: Androst-4-en-17β-ol-3-one 17β-heptanoate, (TEna, Figure 1b);
(ii)
Testosterone cypionate: Androst-4-en-17β-ol-3-one 17β-cyclopentylpropionate, (TCyp, Figure 1c);
(iii)
Testosterone decanoate: Androst-4-en-17β-ol-3-one 17β-decanoate, (TDec, Figure 1d);
(iv)
Testosterone undecanoate: Androst-4-en-17β-ol-3-one 17β-undecanoate, (TUnd, Figure 1e).
The numbering scheme for atoms and steroid backbone rings is presented in accordance with established notations for derivatives and analogues within this family. (Figure 1a) [14].
The literature reports the crystal structures of other shorter testosterone esters, along with their corresponding CCDC Refcodes, as follows: propionate (ZZZRCG01) [15], buciclate (ETEVII) [16], non-esterified form (TESTON10) [17], isocaproate, phenylpropionate and one propionate polymorph (with deposition numbers 2192707, 2192708 and 2192706) [18] and also acetate (TELYEO) [19]. Apart from testosterone esters, the crystal structures of other testosterone derivatives were also investigated [20,21,22]. This manuscript aims to analyze the structural features of various medium and long esters. Structural characterization was achieved by experimental methods including X-ray single crystal diffraction, powder X-ray diffraction, FT-IR spectroscopy and DTA/TG thermal analysis, and a quantitative evaluation of intermolecular interactions was completed by computational methods exploring lattice energies and intermolecular interactions in the solid state.
From a pharmaceutical perspective, it is known that certain pharmaceutical active ingredients are characterized by poor solubility in water, but at the same time they act as lipophilic compounds. Therefore, these compounds are better suited to dissolution in various lipid preparations [23]. Various agents such as hormones (progesterone, estradiol, testosterone esters and their derivatives), deoxycorticosterone and certain vitamins (K and E) are included in formulations which use oils as carriers [24,25]. The solubilities in various lipid-based mixtures were measured in light of this aspect.
Considering the very large number of scientific papers about this hormone, we considered that it would be useful to supplement it with studies related to the structure, intermolecular interactions and solubility.

2. Results

2.1. Crystal Structures Analysis

In Figure S1 (Supplementary Materials), the match between the experimental powder X-ray diffraction patterns and simulated ones obtained based on the CIF files is presented. The data show a good overall agreement, which is indicative of purity, structural homogeneity and that the analyzed single crystals are representative of the entire samples. However, certain reflections are missing or barely visible. This feature can be attributed to a certain degree to the preferred orientation of crystallites in some samples (TDec and TUnd are pronounced in this regard).
Table 1 provides detailed single-crystal data and refinement information for the studied compounds.

2.1.1. TEna (Testosterone Enanthate)

The enanthate ester of testosterone was found to crystallize in the monoclinic P21 space group, with two individual steroid molecules in the asymmetric unit (denoted with A and B suffixes, Figure 2a). They were bridged via C-H⋯O bonding between the O3B carbonyl acceptor of enanthate tail and CH2 donor (C21A-H21B⋯O3B interaction with an interaction magnitude of 45.1 kJ/mol). The values of interaction energies are further explored in detail in Section 2.3: Intermolecular energies evaluation. The ketone O1B oxygen was further building bifurcated C-H⋯O interactions, one to the CH3 group (C19A-H19B⋯O1B with an interaction energy of −29.3 kJ/mol) and one weakly linking the carbonyl O1B oxygen with the six membered A ring of molecule B (C4B-H4B⋯O1B, having a magnitude of −14.0 kJ/mol). On the other hand, O1A ketone oxygen was involved in the C4A-H4A⋯O1A hydrogen bond, which bound two A type molecules with an interaction magnitude of −14.6 kJ/mol. Overall packing shows that steroid molecules were assembled roughly diagonal with respect to the aoc plane (Figure 2b).

2.1.2. TCyp (Testosterone Cypionate)

The cypionate ester of testosterone crystallized in the rare non-centrosymmetric P21212 orthorhombic space group with two distinct molecules in the asymmetric unit (Figure 3a). O1A ketone oxygen was bounded in bifurcated hydrogen bonding with two neighboring A rings (C2B-H2BA⋯O1A with an energy of −20.3 kJ/mol and −14.8 kJ/mol for C4B-H4B⋯O1A). Further, both O3A and O3B carbonyl oxygens were involved in intermolecular interactions with other B steroid rings (C7A-H7AB⋯O3B was characterized by a high intercontact energy of −50.4 kJ/mol and weakly bounded energy of −18.8 kJ/mol for C7B-H8BB⋯O3A). The geometries of intermolecular interactions are detailed in Table S1 (Supplementary Materials), and in Figure 3b the crystal packing perspective seen along the c-axis is presented.

2.1.3. TDec (Testosterone Decanoate)

The decanoate ester of testosterone is one of the longest available and was found to crystallize in the orthorhombic P212121 non-centrosymmetric space group, and consisted of two independent molecules in the asymmetric unit (Figure 4a). Crystal structure stability was held by a combination of C-H⋯O intermolecular bonding with O1B ketone oxygen, with O3A and O3B carbonyl oxygen playing the important roles of acceptors (a high value of interaction energy of −53.1kJ/mol for C21B-H21B⋯O3A; Etot = −43.6 kJ/mol for C26B-H26B⋯O1B and C12B-H12A⋯O3B). The supramolecular arrangement seen along the a-axis shows that testosterone decanoate molecules adopt a W-shaped arrangement (Figure 4b).

2.1.4. TUnd (Testosterone Undecanoate)

This is the longest testosterone ester approached in the actual study. From a crystallographic point of view, it was found to crystallize in the same P21 space group as the enanthate ester. In a manner similar to the previous esters, its asymmetric unit contains two individual steroid molecules (Figure 5a). By the use of both O3A and O3B carbonyl oxygens in the asymmetric unit, the steroid molecules adopt chain-like arrangements in the a-axis direction (the asymmetric unit molecules are tightly bounded with an energy of −60.8 kJ/mol for selected contact C21B-H21A⋯O3A, and for the second interaction in the a-axis direction C21A-H21D⋯O3B, Etot = −51.1 kJ/mol). Once again, O1B oxygen from ketone serves as a hydrogen bonding acceptor, forming bifurcated interactions (C4B-H4B⋯O1B weakly connects a neighboring A ring, Etot = −13.7 kJ/mol, while C19A-H19D⋯O1B with Etot = −28.7 kJ/mol links the methyl CH3 of a neighboring A molecule). In a similar way, with testosterone enanthate the molecules are packed roughly diagonally with respect to the aoc plane (Figure 5b).
After analyzing the four crystal structures, the following conclusions can be drawn:
(i)
All four esters crystallize in various non-centrosymmetric monoclinic and orthorhombic space groups;
(ii)
Asymmetric units in all four crystals are characterized by two individual molecules, as opposed to the other short testosterone esters reported in the literature, such as acetate [19], propionate, isocaproate, phenylpropionate [18] and buciclate [16], which contain only one molecule;
(iii)
Although C-H⋯O interactions participate in the formation of supramolecular 3D assemblies, their weight is quite small compared to dispersion effects (see the crystal energies analysis section);
(iv)
C-H⋯O bonds are characterized by hydrogen⋯carbonyl donor-acceptor distances (Table S3, Supplementary Materials), which fall into the same range as other structures of the steroid class [26,27,28,29];
(v)
The six membered A rings depict an intermediate sofa-half-chair conformation, B and C exhibit chair geometry, and the five membered D backbone rings display an intermediate envelope-half-chair geometry. The C17 methylated form [30] and other short esters have been shown to possess similar geometries [16,18].

2.2. Crystal Lattice Energies Evaluation

Total crystal energies and the breakdown in four energy components were evaluated using the atom–atom CLP model (Table 2). Further, the energies of the four medium and long esterified forms were compared with their non-esterified base form (TBas) and other short esters reported in the literature: acetate (TAce), propionate (TPro), isocaproate (Tiso) and phenylpropionate (TPhp) [18].
All nine crystals were characterized by dominant values of dispersion energies, the values increasing with the increase in ester length: TAce, being the shortest ester, shows a value of −126.1 kJ/mol, whereas TUnd, which possesses an eleven-carbon-longer ester chain, has a value of −175.4 kJ/mol.
Regarding the Coulombian energy, because the structures lack strong hydrogen bondings, it contributes the least weight in the crystal packings, and their magnitudes are similar in all eight esters with an energy range of −15.5 kJ/mol for acetate (TAce) to −22.5 kJ/mol for enanthate (TEna). On the other hand, the Coulombic energy (−33.3 kJ/mol) in the non-esterified testosterone base form (TBas), which possesses classical hydroxyl-carbonyl O-H⋯O hydrogen bonding, shows an increased contribution to the total.
Polarization terms display similar values in six out of the eight esters (TAce, TPro, Tiso, TPhp, Tena, TCyp), but a slight increase in the longest two esters: −70.4 kJ/mol in decanoate (TDec) and −70.0 kJ/mol in undecanoate, showing that these two are slightly polarized.
Repulsion components do not show a trend related to ester length, and significant values are found in TBas and TEna.
With regard to total lattice energies, it can be noticed that crystal became more tightly bounded with the increase in ester length, from TBas −150.8 kJ/mol to Tund −220.7 kJ/mol, showing the highest stability of all nine.
Previous reports in the literature show comparable results of other agents from the steroid family, where dispersion terms play the dominant role and the total energy increases with the increase in ester length [31,32,33,34].

2.3. Intermolecular Energies Evaluation

In order to have a complete perspective of the packing of the crystals, an evaluation of the intermolecular energies of the molecules for which the atoms are located at a smaller distance from each other than the sum of the Van der Waals radii was carried out. The attraction energies (electrostatic, polarization, dispersion), respectively, the repulsion energy and the total energy are presented in Table S4 (Supplementary Materials).
The following conclusions result from Table S4:
Due to the long carbon tails, the molecules have an elongated shape. The interactions between parallel placed molecules are dominated by dispersion energy, due to the greater number of contacts and the small distances between them. For molecules that are placed end to end, the electrostatic interaction energy becomes important.
(i)
Since the crystalline structures do not possess classic strong hydrogen bonds but C-H⋯O interactions, the electrostatic terms have small values in all cases;
(ii)
The values of the polarization energies are very low, which indicates that the molecules are not polarized;
(iii)
The values of the total energies are relatively low and vary in a wide range (−10.4 to −60.8 kJ/mol) due to the different orientations of the molecules relative to each other; between the neighboring molecules, there are rather weak interactions due to the lack of strong hydrogen bonds;
(iv)
Table S4 (Supplementary Materials) shows that the dispersion energy is dominant.

2.4. Hirshfeld Surfaces and Fingerprint Plots Analysis

Molecular Hirshfeld surfaces of investigated crystals were generated via dnorm and compared to each other (Figure S2, Supplementary Materials). Based on the fact that the asymmetric units of all esters are comprised by two individual molecules, they required treating separately.
Intermolecular contacts, which are shorter than the sum of van der Waals radii, were exemplified by arrows and listed in Table S1 (Supporting Information). Based on color coding (red, white, blue), the surfaces can be explained as follows: red dots denote strong intermolecular interactions with distances shorter than the sum of vdW radii, white patches represent the intermolecular contacts close to the sum of vdW radii, while weak interactions with longer distances are illustrated in blue. The fingerprint plots for each crystal, which are related to the 2D representation of the Hirshfeld surfaces, were generated and analyzed as well (Figure S3, Supplementary Materials).
Analyzing the Hirshfeld diagrams and fingerprint plots, certain structural features can be summarized:
(i)
Fingerprint plots of esterified forms (Figure S3, Supporting Information) are asymmetric, which is characteristic of crystals with two or more molecules in the asymmetric unit and appears as a result of a distinct molecular environment in crystal;
(ii)
The plots are illustrating H⋯O/O⋯H spikes, which denote the existence of C-H⋯O hydrogen bonds, but for TDec the H⋯O/O⋯H spikes are less protruding as a consequence of longer distances in the C-H⋯O bonds, which are closer to the sum of vdW radii;
(iii)
A quantitative breakdown of fingerprint diagrams (Table S2, Supplementary Materials) shows the similarities in individual contributions in all four esters, with the highest percentage in H⋯H contacts, followed by a medium percentage of O⋯H/H⋯O intercontacts and a smaller one for C⋯H/H⋯C, respectively;
(iv)
Based on the large percentages of H⋯H contacts for all crystals (fingerprint breakdown in Table S2, Supplementary Materials) corroborated crystal and intermolecular energies (Table 2 and Table S4 from Supplementary Materials) validate that dispersion components play the major role in overall packing.

2.5. Solubility Evaluation

The solubility measurements are represented graphically in Figure 6 and listed in Table S3 (Supporting information). Abbreviations were used in Figure 6: MCT for medium-chain triglyceride, and GSO for grape seed oil. Unfortunately, the evaluation of the solubility of the enanthate ester was not possible because the solutions became transparent, gelatinous mixtures. Analyzing the values, it can be noted that there is no correlation between solubility and the length of the ester. The shortest ester (TCyp) showed an average solubility of roughly 255 mg/mL, the longest one (TUnd) was characterized by the lowest solubility of roughly 197 mg/mL, and a considerably greater solubility of about 2.5–3 fold greater (534 mg/mL) was recorded in the case of TDec.
What is worth noting is the fact that all six mixtures have similar solubilities.
The solubilities of some short esters have been previously reported, averaging 110 mg/mL for the acetate ester, 175.5 mg/mL for the propionate, 139 mg/mL for the phenylpropionate and 447 mg/mL for the isocaproate [18]. These show lower solubilities compared to long esters, but are comparable in orders of magnitude.
The solubility in castor oil of hydroxyprogesterone caproate polymorphs was previously reported in the literature as being 278 mg/mL and 301 mg/mL, respectively, at a temperature of 20° [35], which are comparable to the prodrugs investigated.

2.6. FT-IR Spectroscopy Analysis

The IR spectra of esters are shown and compared in Figure S4 (Supplementary Materials).
As can be seen from this figure, the absorption bands that appear in the 2800–3050 cm−1 wavenumbers region are similar and can be attributed to the symmetric and asymmetric C-H stretching within the CH, CH2, CH3 functional groups.
The two strong absorption bands that appear between wavenumbers 1670–1736 cm−1 are due to the carbonyl C=O groups present in the ketone and ester groups. At lower wavenumbers included in the 1611–1617 cm−1 range, sharp but less intense peaks are observed, which are attributed to C=C stretching of the carbons within the steroid A rings. In the fingerprint area (below 1500 cm−1), the spectra for TEna, TDec, TUnd are similar because only the number of carbon atoms in the tails differentiates them, while TCyp shows differences due to the presence of the cyclopentyl ring.

2.7. DTA/TG Analysis

DTA/TG curves of analyzed compounds are shown in Figure S5 (Supporting Information). Small endothermic peaks appear at (42 °C for Tena), (97 °C for TCyp), (49 °C for TDec) and (65 °C for TUnd), and are assigned to the melting points of the compounds. It should be mentioned that TEna has the lowest melting point, and in the solubility evolution process it became gelatinous even though it was assessed at 14 °C.
With the increase in temperature, no thermal events occurred until the onset temperature of 320 °C for long esters (TDec, TUnd) and up to 350 °C for medium esters (TEna, TCyp).
Further, two broad exothermic peaks were recorded. The first peak had its maximum at the following temperatures: (381 °C, TEna), (393 °C, TCyp), (378 °C, TDec) and (387 °C, TUnd). This peak can be attributed to the process of degradation and oxidation of organic matter, accompanied by a significant loss of mass. The second exothermic peak, with maxima (544 °C, TEna), (537 °C, TCyp), (540 °C, TDec), (540 °C, TUnd), can be associated with the decomposition and oxidation of the remaining material simultaneously with a less significant loss of mass. Similar thermal behavior has been reported for other active pharmaceutical ingredients such as riboflavin and norfloxacin [36].

3. Materials and Methods

3.1. Materials and Recrystallization Experiments

Crystalline powders of esters for scientific research purposes were received from Wuhan Shu Mai Technology Co. (Wuhan, China), while the solvents were received from Merck. Various oils of USP grade were received from Sigma-Aldrich (St. Louis, MO, USA), Tex Lab supply and Med Lab supply (USA).
Good quality single crystals for X-ray data collection were successfully grown in alcohols: isopropilic alcohol (TEna), or mixtures such as methanol–ethanol mixture (TDec, TUnd). The exception is TCyp, which was recrystallized in hexane solution.

3.2. Powder X-ray Diffraction

The diffraction patterns were recorded with a Bruker D8 Advance diffractometer with the X-ray tube operating at 40 kV and 40 mA. It was equipped with a germanium (1 1 1) monochromator used to obtain only the CuKα1 radiation, and a LINXEYE detector. X-ray diffraction patterns for Tena and TUnd were recorded in the 2–40° range, and for TCyp and TDec they were recorded in the 3–40° (2θ) range using the DIFFRAC plus XRD Commander Program and a scanning speed of 0.02°/s.

3.3. Single Crystal X-ray Diffraction and Structures Refinement

Single crystal X-ray diffraction data were collected using a SuperNova diffractometer equipped with dual micro-sources (Cu and Mo) with the X-ray tube operating at 50 kV and 0.8 mA. The collection, reduction and correction of Lorentz, polarization and absorption effects were completed with CrysAlis PRO software [37]. The crystal structures of TCyp and TDec were solved with direct methods by SHELXS [38], TEna with SHELXD [38] and TUnd with SHELXT [39]. All were further refined with Least Squares minimization SHELXL [40], all mentioned programs being included in the Olex2 software [41].
A riding model was considered for hydrogen atoms treatment: the isotropic displacement parameter Uiso(H) = 1.2Ueq(C) for ternary CH groups [C-H = 0.93 Å], secondary CH2 groups [C-H = 0.97 Å] and 1.5Ueq(C) for all methyl CH3 groups [C-H = 0.96 Å].

3.4. Computational Methods

Computational results were achieved based on the fractional coordinates of atoms in the unit cell, determined by single crystal X-ray diffraction analysis. C-H bond distances for CH, CH2 and CH3 groups were normalized to C-H = 1.083 Å.
The crystal energies were calculated using the atom–atom potential Coulomb–London–Pauli (CLP) included in the CLP-Pixel package [42]. The method implies the evaluation of total lattice energies and the breakdown in four terms: Coulombic, polarization, and dispersion, which are attractive, and the fourth repulsive term.
CrystalExplorer software [43] was used to calculate the magnitudes and nature of intermolecular interaction energies. The interactions with distances shorter than the sum of van der Waals radii were computed as pairwise and imply the summation of four energy components: electrostatic (Eele), polarization (Epol), dispersion (Edis) and exchange repulsion term (Erep) [44]. The B3LYP/6-31G(d,p) exchange correlation functional and basis set was used with the scale factors kele = 1.019, kpol = 0.651, kdisp = 0.901, and krep = 0.811.
The dnorm function was used to generate the Hirshfeld surfaces and fingerprint plots in CrystalExplorer software [43].

3.5. Solubility Evaluation

The esters’ solubility (mg/mL) was evaluated in mixtures which included different solvents and organic oils, such as medium-chain triglyceride (MCT), apricot, grape seed (GSO), castor, cottonseed, castor and sesame oil.
The solutions contained mixtures of benzyl alcohol, benzyl benzoate and oil, with a volumetric ratio of 2% benzyl alcohol, 20% benzyl benzoate and 78% oil. Numerous pharmaceuticals of various lipophilic agents (various steroids were included) use benzyl benzoate as a co-solvent/solubilizer; benzyl alcohol behaves as well as a solvent and inhibits the microbial growth, and meanwhile the oils act as carriers.
The evaluation of solubility was completed in successive steps at an ambient temperature of 14 °C by the addition of small quantities of raw compounds (roughly 2–3 mg per step), while the solutions were stirred for several hours until completely dissolved. In certain cases, an excess of raw samples remained undissolved (in suspension) and small quantities of mixtures were added, ensuring that the resulting solutions became homogenous, clear and transparent. Three measurements were taken to improve accuracy, and their average value was used.

3.6. FT-IR Spectroscopy

FT-IR spectra of the samples were obtained using a Jasco 6200 FT-IR spectrometer, with 256 scans and a resolution 4 cm−1, in the spectral range 4000–400 cm−1. The samples were pre-prepared in the form of KBr disks with the pellet technique.

3.7. Differential Thermal Analysis (DTA) and Thermogravimetric Analysis (TG)

DTA/TGA measurements were performed with a simultaneous thermogravimetric and differential thermal analyzer Shimadzu DTG-60H. The samples were heated with a heating rate of 10°C/min, using an alumina sample cell (diameter 5.8 × 2.5 mm2) under a dry nitrogen purge (70 mL/min).

4. Conclusions

The crystal structures of four medium and long testosterone-based esters, including enanthate, cypionate, decanoate and undecanoate, were determined and reported. They were shown to belong to non-centrosymmetric monoclinic and orthorhombic space groups. The skeleton rings’ geometries were found to be sofa-half-chair in A rings and chair conformation for B and C rings, while five membered D rings depicted envelope geometry. Supramolecular arrangements were driven by dominant dispersion components, and the hydrogen⋯carbonyl C-H⋯O interactions played a lesser role in crystal stability. The solubility was evaluated, and the undecanoate form had the lowest solubility, followed by cypionate, while the decanoate form had a considerably higher solubility of roughly two- to threefold.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28073097/s1, Figure S1: Powder X-ray diffraction patterns comparison Experimental and simulated XRPD: TEna (a), TCyp (b), TDec (c), TUnd (d); Figure S2: Hirshfeld surfaces mapped with dnorm ilustrating the contacts referred in Table S1. Surfaces were represented with the clour scale in the ranges as follows: TEna (a) −022 (red) to 174 (blue), −0.19 (red) to 1.93 (blue) for TCyp, −0.09 (red) to 1.68 (blue) for TDec, −0.20 (red) to 1.73 (blue) for TUnd; Figure S3: Fingerprint diagrams of analysed crystals; Table S1: Hydrogen bond geometry for studied crystals (Å, o); Table S2: Contributions to the Hirshfeld surfaces for various intercontacts; Table S3: Solubility of esters in various mixtures; Table S4: Nature and magnitudes of intermolecular interaction energies for selected contacts (kJ/mol); Figure S4: FT-IR spectra of analyzed esters; Figure S5: DTA/TG curves of analyzed compounds

Author Contributions

Conceptualization, G.B. and A.T.; methodology, G.B.; software, A.T.; validation, G.B., P.P., L.M. and V.P.; formal analysis, L.M.; investigation, A.T., G.B., P.P., L.M. and V.P.; resources, G.B.; data curation, G.B.; writing—original draft preparation, G.B., A.T., V.P., L.M.; writing—review and editing, A.T.; visualization, G.B.; supervision, V.P. and G.B.; project administration, V.P.; funding acquisition, G.B. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support from the MCID through the “Nucleu” Program within the National Plan for Research, Development and Innovation 2022–2027, project PN 23 24 01 01. The service for the Oxford SuperNova diffractometer was supported by the Ministry of Research, Innovation, and Digitization through Programme 1—Development of the National Research and Development System, Subprogramme 1.2—Institutional Performance—Funding Projects for Excellence in RDI, Contract No. 37PFE/30.12.2021.

Institutional Review Board Statement

Not applicable

Informed Consent Statement

Not applicable

Data Availability Statement

CIF files of studied crystals were deposited with the Cambridge Crystallographic Data Centre with the following deposition numbers: 2233487 (TEna); 2233486 (TCyp); 2233488 (TDec); 2233489 (TUnd). Copies can be obtained free of charge on written application to CCDC, 12 Union Road, Cambridge, CB2 1EZ, UK (fax: +44-1223-336033); on request via e-mail to [email protected]; or by access to http://www.ccdc.cam.ac.uk (accessed on 20 February 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sinha-Hikim, I.; Taylor, W.E.; Gonzalez-Cadavid, N.F.; Zheng, W.; Bhasin, S. Androgen Receptor in Human Skeletal Muscle and Clutured Muscle Satellite Cells: Up-Regulation by Androgen Treatment. JCEM 2004, 89, 5245–5255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mohamad, N.V.; Soelaiman, I.N.; Chin, K.Y. A concise review of testosterone and bone health. Clin. Interv. Aging 2016, 11, 1317–1324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Katznelson, L. Terapeutic role of androgens in the treatment of osteoporosis in men. Baillieres Clin. Endocrinol. Metab. 1998, 12, 453–470. [Google Scholar] [CrossRef] [PubMed]
  4. Davis, S.R.; Wahlin-Jacobsen, S. Testosterone in women-the clinical significance. Lancet Diabetes Endocrionl. 2015, 3, 980–992. [Google Scholar] [CrossRef] [PubMed]
  5. Roy, A.K.; Lavrovsky, Y.; Song, C.S.; Jung, M.H.; Velu, N.K.; Bi, B.Y.; Chatterjee, B. Regulation of androgen action. Vitam. Horm. 1999, 55, 309–352. [Google Scholar] [CrossRef] [PubMed]
  6. Irwig, M.S. Testosterone therapy for transgender men. Lancet Diabetes Endocrinol. Metab. 2017, 5, 301–311. [Google Scholar] [CrossRef]
  7. Zuloaga, D.G.; Puts, D.A.; Jordan, C.L.; Breedlove, S.M. The role of androgen receptors in the masculinization of brain and behavior: What we’ve learned from the testicular feminization mutation. Horm. Behav. 2008, 53, 613–626. [Google Scholar] [CrossRef] [Green Version]
  8. Kelly, D.M.; Jones, T.H. Testosterone: A metabolic hormone in health and disease. J. Endocrionol. 2013, 217, 25–45. [Google Scholar] [CrossRef] [Green Version]
  9. Vlahopoulos, S.; Zimmer, W.E.; Jenster, G.; Belaguli, N.S.; Balk, S.P.; Brinkmann, A.O.; Lanz, R.B.; Zoumpourlis, V.C.; Schwartz, R.J. Recruitment of the androgen receptor via serum response factor facilitates expression of a myogenic gene. JBC 2005, 280, 7786–7792. [Google Scholar] [CrossRef] [Green Version]
  10. Braunstein, G.D. The influence of anabolic steroids on muscular strength. Princ. Med. Biol. 1997, 8, 465–474. [Google Scholar] [CrossRef]
  11. Vermeulen, A. Longacting steroid preparations. Acta Clin. Belg. 1975, 30, 48–55. [Google Scholar] [CrossRef] [PubMed]
  12. Forsdahl, G.; Erced, D.; Geisendorfer, T.; Turkalj, M.; Plavec, D.; Thevis, M.; Tretzele, L.; Gmeinera, G. Detection of testosterone esters in blood. Drug Test. Anal. 2015, 7, 983–989. [Google Scholar] [CrossRef] [PubMed]
  13. Behre, H.M.; Abshagen, K.; Oettel, M.; Hubler, F.; Nieschlag, E. Intramuscular injection of testosterone undecanoate for the treatment of male hypogonadism: Phase I studies. Eur. J. Endocrinol. 1999, 140, 414–419. [Google Scholar] [CrossRef] [Green Version]
  14. Elks, J.; Ganellin, C.R. The Dictionary of Drugs: Chemical Data: Chemical Data, Structures and Bibliographies, 1st ed.; Springer: Easton, PA, USA, 1990; p. 652. [Google Scholar]
  15. Reisch, J.; Eki-Gucer, N.; Takacs, M.; Henkel, G. Photochemische Studien, 54. Photodimerisierung von Testosteronpropionat in kristallinem Zustand und Kristallstruktur von Testosteronpropionat. LACHDL 1989, 6, 595–597. [Google Scholar] [CrossRef]
  16. Alcock, N.W.; Sanders, K.J.; Rodger, A. Potential injectable contraceptive steroids: Testosterone buciclate. Acta Cryst. 2004, E60, 348–349. [Google Scholar] [CrossRef]
  17. Roberts, P.J.; Pettersen, R.C.; Sheldrick, G.M.; Isaacs, N.W.; Kennard, O. Crystal and molecular structure of 17β-hydroxyandrost-4-en-3-one (testosterone). J.Chem. Soc. Perkin Trans. 2 1973, 14, 1978–1984. [Google Scholar] [CrossRef]
  18. Turza, A.; Popescu, V.; Mare, L.; Borodi, G. Structural Aspects and Intermolecular Energy for Some Short Testosterone Esters. Materials 2022, 15, 7245. [Google Scholar] [CrossRef]
  19. Böcskei, Z.; Gérczei, T.; Bodor, A.; Schwartz, R.; Náray-Szabó, G. Three Testosterone Derivatives. Acta Cryst. 1996, C52, 2899–2903. [Google Scholar] [CrossRef]
  20. Silva, M.R.; Beja, A.M.; Moreira, V.M.; Santos, R.C.; Salvador, A.R. 3-Oxoandrost-4-en-17β-yl iH-imidazole-1-carboxylate. Acta Cryst. 2007, E63, o4824. [Google Scholar] [CrossRef]
  21. Hoyama, H.; Shiro, M.; Sato, T.; Tsukuda, Y. Crystal Structures of 8β-Methyltestosterone 17 β-Monobromoacetate and Testosterone 17 β-p-Bromobenzoate. J. Chem. Soc. B 1970, 2, 443–452. [Google Scholar] [CrossRef]
  22. Ohrt, J.M.; Haner, B.A.; Norton, D.A. Crystal data (I) for some halogenated steroids. Acta Cryst. 1965, 19, 280. [Google Scholar] [CrossRef]
  23. Land, L.M.; Li, P.; Baummer, P.M. The Influence of Water Content of Triglyceride Oils on the Solubility. Pharm. Res. 2005, 5, 784–788. [Google Scholar] [CrossRef] [PubMed]
  24. Remington, J.P. Remington: The Science and Practice of Pharmacy, 23rd ed.; Adejare, A., Ed.; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2020. [Google Scholar]
  25. Rowe, R.C.; Sheskey, P.J.; Quinn, M.E. (Eds.) Handbook of Pharmaceutical Excipients, 5th ed.; Libros Digitales-Pharmaceutical Press: London, Great Britain, 2006. [Google Scholar]
  26. Rendle, D.F.; Trotter, J. Crystal and molecular structure of 17β-hydroxy-17α-methyl-2oxa-5α-androstan-3-one. J. Chem. Soc. Perkin Trans. 1975, 2, 1361–1365. [Google Scholar] [CrossRef]
  27. Turza, A.; Borodi, G.; Pop, M.M.; Ulici, A. Polymorphism and β-cyclodextrin complexation of methyldrostanolone. J. Mol. Struct. 2022, 1250, 131852. [Google Scholar] [CrossRef]
  28. Turza, A.; Ulici, A.; Pop, M.M.; Borodi, G. Solid forms and β-cyclodextrin complexation of turinabol. Acta Cryst. 2022, C78, 305–313. [Google Scholar] [CrossRef] [PubMed]
  29. Courseille, C.; Precigoux, G.; Leroy, F.; Busetta, B. 5α-Androstan-17β-ol-3one, C19H30O2. Cryst. Struct. Commun. 1973, 2, 441. [Google Scholar]
  30. Gaedecki, Z.; Grochulski, P.; Wawrzak, Z. Structure of 17α-methyl-testosterone semihydrate, C20H30O2·1/2H2O. J. Crystallogr. Spectrosc. Res. 1989, 19, 577–587. [Google Scholar] [CrossRef]
  31. Turza, A.; Miclaus, M.O.; Pop, A.; Borodi, G. Crystal and molecular structures of boldenone and four boldenone steroid esters. Z. Kristallogr. Cryst. Mater. 2019, 234, 671–683. [Google Scholar] [CrossRef]
  32. Borodi, G.; Turza, A.; Camarasan, P.A.; Ulici, A. Structural studies of Trenbolone, Trenbolone Acetate, Hexahydrobenzylcarbonate and Enanthate esters. J. Mol. Struct. 2020, 1212, 128127. [Google Scholar] [CrossRef]
  33. Borodi, G.; Turza, A.; Bende, A. Exploring the Polymorphism of Drostanolone Propionate. Molecules 2020, 25, 1436. [Google Scholar] [CrossRef] [Green Version]
  34. Turza, A.; Borodi, G.; Pop, A.; David, M. Structural studies of some androstane based prodrugs. J. Mol. Struct. 2022, 1248, 131440. [Google Scholar] [CrossRef]
  35. Caplette, J.; Frigo, T.; Jozwiakowski, M.; Shea, H.; Mirmehrabi, M.; Müller, P. Characterization of new crystalline forms of hydroxyprogesterone caproate. Int. J. Pharm. 2017, 527, 42–51. [Google Scholar] [CrossRef] [PubMed]
  36. Ferreira, L.T.; Alarcon, R.T.; Perpetuo, G.L.; Bannach, G. Investigation and characterization by TG/DTG-DTA and DSC of the fusion of Riboflavin, and its interaction with the antibiotic norfloxacin in the screening of cocrystal. J. Therm. Anal. Calorim. 2019, 136, 581–588. [Google Scholar] [CrossRef] [Green Version]
  37. Rigaku, O.D. CrysAlis PRO; Rigaku Oxford Diffraction Ltd.: Yarnton, UK, 2015. [Google Scholar]
  38. Sheldrick, G.M. A short history of SHELX. Acta Cryst. 2008, A64, 112–122. [Google Scholar] [CrossRef] [Green Version]
  39. Sheldrick, G.M. SHELXT—Integrated space-group and crystal-structure determination. Acta Cryst. 2015, A71, 3–8. [Google Scholar] [CrossRef] [Green Version]
  40. Sheldrick, G.M. Crystal structure refinement with SHELXL. Acta Cryst. 2015, C71, 3–8. [Google Scholar] [CrossRef] [Green Version]
  41. Dolomanov, O.V.; Bourhis, L.J.; Gildea, R.J.; Howard, J.A.K.; Puschmann, H. OLEX2: A complete structure solution, refinement and analysis program. J. Appl. Cryst. 2009, 42, 339–341. [Google Scholar] [CrossRef]
  42. Gavezzotti, A. Efficient computer modeling of organic materials. The atom–atom, Coulomb–London–Pauli (AA-CLP) model for intermolecular electrostatic-polarization, dispersion and repulsion energies. New J. Chem. 2011, 35, 1360–1368. [Google Scholar] [CrossRef]
  43. Turner, M.J.; McKinnon, J.J.; Wolff, S.K.; Grimwood, D.J.; Spackman, P.R.; Jayatilaka, D.; Spackman, M.A. Spackman, CrystalExplorer17; University of Western Australia: The Nedlands, Australia, 2017. [Google Scholar]
  44. Mackenzie, C.F.; Spackman, P.R.; Jayatilaka, D.; Spackman, M.A. CrystalExplorer model energies and energy frameworks: Extension to metal coordination compounds, organic salts, solvates and open-shell systems. IUCrJ 2017, 4, 575–587. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Chemical perspective of testosterone (17β-Hydroxyandrost-4-en-3-one), showing the backbone labeling system (a) and the studied testosterone-based prodrugs: enanthate (b), cypionate (c), decanoate (d), undecanoate (e).
Figure 1. Chemical perspective of testosterone (17β-Hydroxyandrost-4-en-3-one), showing the backbone labeling system (a) and the studied testosterone-based prodrugs: enanthate (b), cypionate (c), decanoate (d), undecanoate (e).
Molecules 28 03097 g001
Figure 2. Asymmetric unit of Tena presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along b-axis (b).
Figure 2. Asymmetric unit of Tena presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along b-axis (b).
Molecules 28 03097 g002
Figure 3. Asymmetric unit of TCyp presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along c-axis (b).
Figure 3. Asymmetric unit of TCyp presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along c-axis (b).
Molecules 28 03097 g003
Figure 4. Asymmetric unit of TDec presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along a-axis (b).
Figure 4. Asymmetric unit of TDec presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along a-axis (b).
Molecules 28 03097 g004
Figure 5. Asymmetric unit of TUnd presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along the b-axis (b).
Figure 5. Asymmetric unit of TUnd presenting non-hydrogen atoms at 50% probability level (a); overall packing viewed along the b-axis (b).
Molecules 28 03097 g005
Figure 6. Graphical representations of esters solubility (standard deviation is represented by error bars).
Figure 6. Graphical representations of esters solubility (standard deviation is represented by error bars).
Molecules 28 03097 g006
Table 1. Crystal structures and refinement data of investigated compounds.
Table 1. Crystal structures and refinement data of investigated compounds.
Identification CodeTEnaTCypTDecTUnd
Empirical formulaC26H40O3C27H40O3C29H46O3C30H48O3
Formula weight400.58411.58442.66456.68
Temperature/K293(2)293(2)293(2)293(2)
Crystal systemmonoclinicorthorhombicorthorhombicmonoclinic
Space groupP21P21212P212121P21
a/Å10.5644(6)20.9920(6)8.3367(3)10.5820(2)
b/Å8.2875(6)32.0578(8)19.0481(6)8.23950(10)
c/Å27.6918(15)7.2292(2)33.7169(12)32.3890(4)
α/°90909090
β/°94.652(5)909099.3530(10)
γ/°90909090
Volume/Å32416.5(3)4864.9(2)5354.2(3)2786.47(7)
Z4884
ρcalc g/cm31.1011.1241.0981.089
μ/mm−10.5420.5530.5310.523
F(000)880.01800.01952.01008.0
RadiationCuKα (λ = 1.54184)CuKα (λ = 1.54184)CuKα (λ = 1.54184)CuKα (λ = 1.54184)
2Θ range/°6.404 to 141.126.938 to 145.17.002 to 148.1648.3 to 141.722
Index ranges−12 ≤ h ≤ 12, −10 ≤ k ≤ 9, −33 ≤ l ≤ 26−25 ≤ h ≤ 25, −38 ≤ k ≤ 36, −7 ≤ l ≤ 8−8 ≤ h ≤ 10, −23 ≤ k ≤ 23, −41 ≤ l ≤ 41−12 ≤ h ≤ 12, −10 ≤ k ≤ 9, −39 ≤ l ≤ 39
Reflections collected15,91332,74975,17740,701
Independent reflections8094 [Rint = 0.0422, Rsigma = 0.0471]9217 [Rint = 0.0620, Rsigma = 0.0556]10,249 [Rint = 0.1805, Rsigma = 0.0748]10,192 [Rint = 0.0333, Rsigma = 0.0213]
Data/restraints/parameters8094/1/5299217/0/54510,249/0/58310,192/1/601
Goodness-of-fit on F21.0011.0311.0851.044
Final R indexes [I>=2σ (I)]R1 = 0.0549, wR2 = 0.1372R1 = 0.0706, wR2 = 0.1833R1 = 0.0938, wR2 = 0.1762R1 = 0.0428, wR2 = 0.1141
Final R indexes [all data]R1 = 0.0896, wR2 = 0.1652R1 = 0.1199, wR2 = 0.2265R1 = 0.1406, wR2 = 0.2100R1 = 0.0462, wR2 = 0.1186
Largest diff. peak/hole/e Å−30.15/−0.180.25/−0.170.16/−0.220.17/−0.21
Flack parameter−0.21(17)0.04(18)−0.02(16)−0.09(8)
Table 2. Crystal lattice energies [18].
Table 2. Crystal lattice energies [18].
StructureMolar Mass g/molEcoul (kJ/mol)Epol (kJ/mol)Edisp (kJ/mol)Erep (kJ/mol)Elatt (kJ/mol)
TBas288.43−33/3−47.5−130.960.9−150.8
TAce330.46−15.5−55.0−126.137.1−159.5
TPro344.49−18.3−55.5−126.433.9−166.3
TIso386.57−21.3−57.7−141.952.7−168.2
TPhp420.59−19.3−52.3−149.036.1−184.5
TEna400.60−22.5−60.4−157.464.6−175.7
TCyp412.61−17.6−58.5−142.837.0−181.9
TDec442.68−17.4−70.4−167.340.1−215.1
TUnd456.71−19.6−73.0−175.447.3−220.7
Ecoul: Coulombic energy; Epol: polarization energy; Edisp: dispersion energy; Erep: repulsion energy; Elatt: total crystal lattice energy.
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

Turza, A.; Pascuta, P.; Mare, L.; Borodi, G.; Popescu, V. Structural Insights and Intermolecular Energy for Some Medium and Long-Chain Testosterone Esters. Molecules 2023, 28, 3097. https://doi.org/10.3390/molecules28073097

AMA Style

Turza A, Pascuta P, Mare L, Borodi G, Popescu V. Structural Insights and Intermolecular Energy for Some Medium and Long-Chain Testosterone Esters. Molecules. 2023; 28(7):3097. https://doi.org/10.3390/molecules28073097

Chicago/Turabian Style

Turza, Alexandru, Petru Pascuta, Liviu Mare, Gheorghe Borodi, and Violeta Popescu. 2023. "Structural Insights and Intermolecular Energy for Some Medium and Long-Chain Testosterone Esters" Molecules 28, no. 7: 3097. https://doi.org/10.3390/molecules28073097

APA Style

Turza, A., Pascuta, P., Mare, L., Borodi, G., & Popescu, V. (2023). Structural Insights and Intermolecular Energy for Some Medium and Long-Chain Testosterone Esters. Molecules, 28(7), 3097. https://doi.org/10.3390/molecules28073097

Article Metrics

Back to TopTop