Next Article in Journal
Design of Downhole Robot Actuator System and Mechanical Behavior Analyses of the PRSM by Considering Elastic Errors and Radial Loads
Previous Article in Journal
Hybrid Memetic Algorithm to Solve Multiobjective Distributed Fuzzy Flexible Job Shop Scheduling Problem with Transfer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Degradation Process of Ethinylestradiol—Kinetic Study

1
Advanced Instrumental Screening Center, Faculty of Pharmacy, Victor Babeş University of Medicine and Pharmacy, 2 Eftimie Murgu Square, 300041 Timisoara, Romania
2
Research Centre for Pharmaco-Toxicological Evaluation, Faculty of Pharmacy, Victor Babeş University of Medicine and Pharmacy, 2 Eftimie Murgu Square, 300041 Timisoara, Romania
3
Faculty of Industrial Chemistry and Environmental Engineering, Politehnica University Timisoara, 2 Victoriei Square, 300006 Timisoara, Romania
4
Faculty of Medicine, Victor Babeş University of Medicine and Pharmacy, 2 Eftimie Murgu Square, 300041 Timisoara, Romania
*
Author to whom correspondence should be addressed.
Processes 2022, 10(8), 1518; https://doi.org/10.3390/pr10081518
Submission received: 21 June 2022 / Revised: 21 July 2022 / Accepted: 28 July 2022 / Published: 2 August 2022

Abstract

:
The present study reports the results obtained after the analysis of the thermal stability and decomposition kinetics of widely used synthetic derivative of estradiol, ethinylestradiol (EE), as a pure active pharmaceutical ingredient. As investigational tools, Fourier transformed infrared spectroscopy (FTIR), thermal analysis, and decomposition kinetics modeling of EE were employed. The kinetic study was realized using three kinetic methods, namely Kissinger, Friedman, and Flynn-Wall-Ozawa. The results of the kinetic study are in good agreement, suggesting that the main decomposition process of EE that takes place in the 175–375 °C temperature range is a single-step process, invariable during the modification of heating rate of the sample.

1. Introduction

Ethinylestradiol (17 alpha-ethinylestradiol, abbreviated EE) is a synthetic derivative of estradiol [1] and is an estrogen receptor agonist that binds to both forms of the receptor (ERα and ERβ) and to the G protein coupled estrogen receptor (GPER) [2]. If compared to estradiol, EE is characterized by a higher bioavailability when administered orally, increased resistance to metabolism, and stronger effects in certain parts of the body [1]. EE is commonly used in combination with progestin as a contraceptive, in the treatment of certain gynecological disorders, hormone-dependent cancers (such as prostate and breast cancer), and menopausal symptoms [1]. As a contraceptive, EE is usually available in low doses, in combination with progestin, to reduce the risk of side effects. The intensity of these side effects depends on the dose and route of administration and include weight gain, headaches, nausea, bloating, breast tenderness, and overall feminization (in males) [3]. Long-term administration may lead to an increased risk of blood clots, cardiovascular issues, liver damage, and high doses of EE have been associated with an increased risk of endometrial cancer [1]. More recently, a pilot clinical trial study was carried out by Cortés-Algara et al. regarding the use of EE along with norelgestromin as immunoregulators incorporated in transdermal patches, as an option for the treatment of COVID-19 disease [4]. Moreover, the role of several estrogens, including EE, in menstrual migraine was also studied [5]. The structural formula of EE is presented in Figure 1.
In the pharmaceutical market, worldwide, there are numerous dosage forms containing EE as the solitary active pharmaceutical ingredient (API), or alongside with other API (such as dienogest, levonorgestrel, chlormadinone, norethisterone, norelgestromin and desogestrel), designed for oral (tablets and coated tablets), vaginal (rings), or transdermal (patches with extended release) administration. Dependent of the administration route, the content of EE in each formulation varies commonly between 20 and 50 micrograms per tablet, up to 600 micrograms per patch and 2.7 mg per ring [6]. Since the amount of EE (micrograms) in each formulation is considerably lower than the amount of excipients (milligrams), the development of stable formulations is an imperious demand.
Kinetic analysis—as a direct implementation of thermal analysis—is an important investigational tool for characterization of drugs, including their behavior in pharmaceutical formulations, leading to crucial data regarding lifetime, shelf-life, and degradation mechanism of investigated compounds [7,8,9,10,11,12]. Taking into account the fact that the mechanism of thermolysis is not known for the majority of organic molecules lead to the main advantage of implementing isoconversional methods, since they allow the estimation of activation energy without knowing or assuming an explicit model for the differential or integral conversion functions [13]. Heterogeneous kinetics has developed over the years and during this period of time, a series of recommendations was elaborated by the International Confederation for Thermal Analysis and Calorimetry (ICTAC) Kinetics Committee, in order to improve the quality of these studies [14]. Numerous papers presents, also, the advantages and disadvantages of each type of kinetic modeling approaches, including model-fitting, model-free (isoconversional), and deconvolution analyses [15,16,17,18,19].
The kinetic method proposed by Kissinger and published initially in 1956 [20] and later in 1957 [21] is one of the most popular kinetic protocol that allows the estimation of the activation energy by using differential scanning calorimetry (DSC) data, differential thermal analysis (DTA) data, or derivative thermogravimetry (DTG) data. The simplicity of its use—however, it is tricky for non-expert users since the method can lead to false conclusions mainly for complex degradative processes—reveals a solitary activation energy value, based on the hypothesis of single-step kinetics. For this reason, methods such as Kissinger, classic Ozawa or ASTM E698 are always preliminary to isoconversional studies, which offer by far a more complete and objective perspective over the investigated processes [18].
Regarding the state-of-art for stability and degradation of EE, the literature data reveal several studies carried out between 1993 and 2022, but none of them refer to heterogeneous degradation of this API in solid state. Recent contributions were reported regarding the degradation of EE by different bacterial strains during wastewater treatment, including Enterobacter tabaci [22] and Acinetobacter [23], and as well the evaluation of kinetics of natural degradation and identification of the resulted products during photodegradation and oxidation were carried out [24]. Moreover, several studies were published regarding the risks and pollution characteristics of contaminants for aquatic ecosystems, including the endocrine disruptors such as EE [25,26,27,28]. Adsorption of EE from different aqueous environments was also studied [29,30,31,32], and as well of the effect over cellular and molecular variations of microalgae such as Chlorella pyrenoidosa [33].
Since the literature data reveal no information regarding the processes of heterogeneous decomposition of this API, we set our goal into carrying out an isoconversional kinetic study according to ICTAC 2000 protocol for this API, under thermal stress [34,35,36,37], using the derivative thermogravimetric (DTG) data collected for five different heating rates β = 2, 4, 6, 8, and 10 °C/min. The obtained results were processed according to Kissinger method, the differential method of Friedman, and the integral method of Flynn–Wall–Ozawa.

2. Materials and Methods

2.1. Samples and Preparation

Ethinylestradiol (EE), a commercial product of Sigma-Aldrich (St. Louise, MO, USA) was used without further purification. The purity of EE was according to European Pharmacopoeia (EP) Reference Standard. Up to the use, the sample was stored according to recommendations made by the supplier.

2.2. FTIR Investigations

Fourier-transform infrared spectroscopy studies (FT-IR), were performed using a Shimadzu Prestige-21 spectrometer (Duisburg, Germany), at 24 °C. The operating parameters set were: resolution of 4 cm−1 within the spectral range of 400–4000 cm−1, using KBr pellets. In the spectroscopic description of bands, the following abbreviations are used for different types of vibration: ν for stretching vibration, δ for internal deformation (bending).

2.3. Thermo-Analytical Investigations

The stability of the samples was performed by thermal behavior assess, using a Netzsch STA 449 C instrument (Netzsch-Gerätebau GmbH, Selb, Germany), in the range of 20–500 °C, air atmosphere, at 2, 4, 6, 8, and 10 °C/min heating rates. Each sample was exactly weighed in aluminum crucibles and the analysis was performed under artificial air at a flow rate of 20 mL/min. Air atmosphere was chosen since most of APIs are processed and stored under usual ambient atmosphere. The analysis was carried out in duplicate and the results are practically identical.

2.4. Kinetic Study

The kinetic processing of the data (Friedman and Flynn-Wall-Ozawa methods) was carried out on the main decomposition process of EE, using the AKTS—Thermokinetics Software (AKTS AG TechnoArk, Siders, Switzerland). The classical Kissinger method was employed using a template file created by our research team. All the aspects regarding the theoretical foundation and advantages of isoconversional kinetics are presented exhaustively in numerous papers [38,39,40].

3. Results and Discussion

3.1. FTIR Investigations

FTIR spectroscopy was used as an investigational technique for characterization of EE. The FTIR spectrum of EE on spectral range 4000–400 cm−1 is shown in Figure 2. Literature data present several characteristic bands for infrared spectral investigation of EE, especially for the compound adsorbed from aqueous medium [41]. However, a complete description of FTIR bands is not presented in the literature regarding the spectroscopic analysis of EE.
The stretching vibrations ν(O–H) for both OH moieties from the EE structure are observed in the spectral range 3650−3100 cm−1 as a broad band, suggesting the intense H-bonding between the molecules in solid state, overlapped with the sharp band of superficially adsorbed water: the bands are evidenced at 3606.89, 3500.80, and 3292.49 cm−1, respectively. The strong band at 3321.42 cm−1 can be assigned to ν(≡C–H) from alkyne moiety. The symmetric and asymmetric stretching vibrations for other C–H bonds, namely ν(C–H), including the ones from CH3 moiety and CH2 ones are represented by the bands at 2972.31, 2935.66, and 2866.22 cm−1. The sharp, weak bands at 2357.01 and 2322.29 cm−1 are characteristic bands, well individualized in the FTIR spectrum of compounds containing -C≡C- moieties, i.e., the ethinyl moiety of the API. The bands from 1614.42, 1585.49, 1500.62, 1471.69, 1446.61, and 1435.04 cm−1 are due to symmetric C-C stretching (νC-C and νC=C), as well for δas(CH3) and δ(CH2). The bands from 1373.32, 1357.89, 1298.09, and 1286.52 cm−1 are the consequence of symmetric methyl bending δs(CH3) and the latter two for the hydroxyls, δ(COH), respectively. The bands at 1255.66 and 1056.99 cm−1 are probably due to ν(C–O) vibration, while the other bands (recorded at 1182.36, 1134.14, 1111.00, 1020.34, 972.12, 929.69, 914.26, 879.54, 819.75, 788.89, 680.87, 644.22, 621.08, 569.00, 524.64, and 441.70 cm−1) from the fingerprint region are due to skeleton vibration and different combination bands, that cannot be correctly attributed to a certain bond without carrying theoretical simulations of vibrational spectra using density functional theory [42].

3.2. Thermoanalytical Investigations

Thermal stability of EE was carried out in dynamic oxidative atmosphere at a heating rate of 2 °C/min, as shown in Figure 3.
According to the thermoanalytical curves depicted in Figure 3, EE is thermally stable up to 71 °C, when a mass loss process begins, taking place between 71 and 101 °C (mass loss 3.18%, DTG process between 69 and 101 °C, DTGpeak at 87 °C, DSC process 61 and 101 °C, DSCpeak at 88 °C). This first process is attributed to the superficially adsorbed water removal from the sample. Anhydrous EE is stable up to 177 °C, when a decomposition process begins, overlapping the melting of the API (DSC endothermal process between 177–191 °C, DSCpeak at 185 °C), in good agreement with the melting interval mentioned on PubChem, suggesting the existence of polymorphic form II [43], instead of polymorphic form I, which melts at 146 °C [44].
The main decomposition process of EE takes place in the temperature range 187–324 °C (mass loss 59.93%, DTGpeak at 293 °C), accompanied by an exothermic effect on the DSC curve, with the peak at 287 °C. Further on, with the increase of temperature, the thermoanalytical profile becomes more complex, since the overlapping thermal degradation process takes place. At 500 °C, the residual mass is 14.18%, suggesting an advanced degradation of the structure under thermal stress.

3.3. Kinetic Study

The kinetic analysis was carried out over processed DTG, obtained in dynamic air atmosphere for the following five heating rates β: 2, 4, 6, 8, and 10 °C/min. The main decomposition process of EE that was subjected to kinetic analysis is the one that takes place after the formation of anhydrous form, which takes place in the following temperature ranges vs. heating rate: β = 2 °C/min (179–330 °C, Tmax at 293.8 °C), β = 4 °C/min (179–348 °C, Tmax at 311.7 °C), β = 6 °C/min (179–355 °C, Tmax at 319.8 °C), β = 8 °C/min (182–377 °C, Tmax at 326.9 °C), and β = 10 °C/min (187–378 °C, Tmax at 332.9 °C). A preliminary kinetic study was realized using the Kissinger method, which is based on the assumption that the degree of conversion is a constant and non-dependent of the heating rate at the DTG peak (Tmax) (Equation (1)):
ln   ( β T max 2 )   =   ln   ( A R E a 1 )   +   ln   [ n ( 1 α max ) n 1 ]     E a R 1 T max 1
The estimation of the activation energy of the decomposition (Ea) can be made by evaluating the slope of the linear plotting for experiments carried out at the five different heating rates (Figure 4), revealing that for the main stage of EE degradation, the value for Ea is 107.91 kJ/mol.
Even if current protocols regarding the kinetic studies highly recommend the use of isoconversional methods, as requested by international conventions established by the ICTAC protocols, the use of simple methods, such as Kissinger or ASTM E698 is valuable, since it permit a classification of kinetic mechanism for decomposition of organic molecules, and a separation of simple vs. complex processes. It is well-known that organic molecules possess numerous pathways of thermolysis since they contain numerous covalent bonds, and the intensity of thermal stress can determine the modification of the mechanism with increasing of heating rate, and as a consequence, the comparison of the results obtained by classical kinetic methods vs. isoconversional ones can lead to valuable information regarding the complexity of the degradation. Following this hypotheses, if the result obtained by Kissinger method is in good agreement with the ones suggested by differential and integral isoconversional methods, it can be said that the heating rate has no influence over modification of degradation mechanism.
Two isoconversional methods, namely the differential method of Friedman (Fr) and the integral method of Flynn–Wall–Ozawa (FWO), were used in order to evaluate the values of Ea of the decomposition vs. conversion degree α. The reaction progress vs. temperature is presented in Figure 5. As can be shown from Figure 5, the decomposition process is shifted to higher temperatures due to thermal inertia of the sample, as the heating rate increases.
The same tendency can be seen in Figure 6, where the plotting of reaction rate vs. temperature reveals the shifting of maximum heating rate at higher temperatures, as the hating rate increases. Moreover, it can be seen that the reaction rate is increased not only by temperature, but also by the heating rate of the sample.
The theoretical background of isoconversional methods is extensively reported in literature [45,46,47,48,49,50,51]; however, for a facile understanding of the results, we briefly present in this paper the final models implemented in kinetic analysis of EE.
The linearized mathematical equation of Friedman method (Fr) [45] is shown in Equation (2).
ln   ( β     d α / d T )   =   ln   [ A     f ( α ) ]     E a     R   1   T   1
For selected α at each heating rates, the plot of ln   ( β   d α / d T ) vs. 1 / T is linear. By evaluation of the slopes of these graphical representations (Figure 7), the values of the activation energy of the decomposition (Ea) are revealed (Table 1).
The Flynn–Wall–Ozawa (FWO) method [48,49] is represented, after Doyle linearization, by the mathematical model presented in Equation (3), where g(α) is the integral conversion function.
ln   β   =   ln   [ A     E a     R 1     g 1 ( α ) ]     5.331   1.052   E a     R 1     T 1
Similarly, the estimation of activation energy of the decomposition values (Ea) for all the conversion degree can be achieved by the plotting of ln β vs. T 1 , as can be seen in Figure 8 and Table 1, respectively.
As can be seen by comparing the results from the Kissinger kinetic method and the two isoconversional ones, the obtained values for the decomposition activation energies are in good agreement. Moreover, the fitting of Tmax points from the Kissinger method suggests a very good fitted regression model (R2 = 0.9968), leading to a possible conclusion that the degradative mechanism of EE under thermal stress is independent of the heating rate. However, this conclusion can be sustained only by the correlation of the result with the ones obtained by the use of the two isoconversional methods, the differential and the integral one, since there are situations where the Kissinger plot is almost perfectly linear, without detecting the complexity of the decomposition process(es) [18].
The necessity of estimating the value of activation energy for small variation of conversion degree (5%) is based on the fact that the individual Ea values falling outside the ± 10% interval around the medium value, clearly indicate a multi-step degradation process. Otherwise, if the variation of Ea vs. α values are inside the interval E a ¯ ± 0.1 E a ¯ , the mechanism of degradation consists in a single-step process, invariable with the modification of heating rate of the sample.
The isoconversional method of Friedman reveals a variation outside the ±10% limit, outside the medium value of Ea for EE solely at conversions superior to 85%, as seen in Table 1. This variation is not observed in the case of the integral method of Flynn–Wall–Ozawa, even if the individual Ea values have the tendency to decrease with the advance of the reaction.

4. Conclusions

In this paper, the results obtained after studying the process of heterogeneous degradation of the estrogen medication ethinylestradiol were reported, by means of thermal analysis and isoconversional kinetics. The analysis was also completed by FTIR spectroscopy, in order to prove the identity and purity of the compound.
Thermal analysis revealed that anhydrous EE has a good thermal stability (up to 177 °C), this fact being explained by the existence of the stable steroid moiety in the structure of this drug. The kinetic study was realized using three kinetic methods, namely Kissinger, Friedman and Flynn–Wall–Ozawa. The results of the kinetic study are in good agreement, suggesting that the main decomposition process of EE that take place in the 175–375 °C temperature range is a single-step process, invariable with the modification of heating rate of the sample.
However, due to differential process of the data according to Friedman method, at higher conversion degrees (α > 0.85), a clear indication for a modification of the degradative mechanism is revealed. This last observation can be explained by the fact that the main degradation process is immediately followed by another process, more complex, as can be also seen on the thermoanalytical profile of EE. Moreover, the complexity of the second degradative process can be explained by the fact that at higher temperatures, the thermolysis of steroid moiety is accentuated and also overlaps the degradation of the products resulted from the main process, that was kinetically investigated in this study.

Author Contributions

Conceptualization, S.S. and I.L.; methodology, A.L., E.-A.M., C.P. and I.L.; software, A.L., E.-A.M., C.P. and I.L.; validation, A.L., C.D., D.N. and I.L.; formal analysis, A.L., E.-A.M., C.P. and I.L.; investigation, S.S., A.L., E.-A.M., C.P. and I.L.; resources, C.D. and D.N.; data curation, S.S., A.L., E.-A.M., C.P. and I.L.; writing—original draft preparation, S.S. and I.L.; writing—review and editing, C.D., D.N. and I.L.; visualization, C.D., D.N. and I.L.; supervision, A.L. and I.L.; project administration, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Authors can provide raw data under request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kuhl, H. Pharmacology of estrogens and progestogens: Influence of different routes of administration. Climacteric 2005, 8 (Suppl. S1), 3–63. [Google Scholar] [CrossRef]
  2. Prossnitz, E.R.; Arterburn, J.B. International Union of Basic and Clinical Pharmacology. XCVII. G Protein-Coupled Estrogen Receptor and Its Pharmacologic Modulators. Pharmacol. Rev. 2015, 67, 505–540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Stanczyk, F.Z.; Archer, D.F.; Bhavnani, B.R. Ethinyl estradiol and 17β-estradiol in combined oral contraceptives: Pharmacokinetics, pharmacodynamics and risk assessment. Contraception 2013, 87, 706–727. [Google Scholar] [CrossRef]
  4. Alfredo, C.-A.; Noemí, C.-R.; Samuel, R.-L.; Daniel, O.-C.; Rodrigo, R.-B.; Paul, M.-T.; Mónica, E.-T.; José, C.-B.; Eleazar, L.-P.; Alfonso, A.-R.; et al. Effect of Norelgestromin and Ethinylestradiol in Transdermal Patches on the Clinical Outcomes and Biochemical Parameters of COVID-19 Patients: A Clinical Trial Pilot Study. Pharmaceuticals 2022, 15, 757. [Google Scholar] [CrossRef]
  5. Nappi, R.E.; Tiranini, L.; Sacco, S.; De Matteis, E.; De Icco, R.; Tassorelli, C. Role of Estrogens in Menstrual Migraine. Cells 2022, 11, 1355. [Google Scholar] [CrossRef] [PubMed]
  6. Ethinylestradiol—Drugbank Profile. Available online: https://go.drugbank.com/drugs/DB00977 (accessed on 18 June 2022).
  7. Alves, G.L.; Teixeira, F.V.; da Rocha, P.B.R.; Krawczyk-Santos, A.P.; Andrade, L.M.; Cunha-Filho, M.; Marreto, R.N.; Taveira, S.F. Preformulation and characterization of raloxifene-loaded lipid nanoparticles for transdermal administration. Drug Deliv. Transl. Res. 2021, 12, 526–537. [Google Scholar] [CrossRef] [PubMed]
  8. Calvino, M.M.; Lisuzzo, L.; Cavallaro, G.; Lazzara, G.; Milioto, S. Non-isothermal thermogravimetry as an accelerated tool for the shelf-life prediction of paracetamol formulations. Thermochim. Acta 2021, 700, 8940. [Google Scholar] [CrossRef]
  9. Asran, A.M.; Mohamed, M.A.; Khedr, G.E.; Eldin, G.M.G.; Yehia, A.M.; Mishra, R.K.; Allam, N.K. Investigation of the thermal stability of the antihypertensive drug nebivolol under different conditions: Experimental and computational analysis. J. Therm. Anal. Calorim. 2022, 147, 5779–5786. [Google Scholar] [CrossRef]
  10. Ahmed, A.A.M.; Asran, A.M.; Mohamed, M.A. Thermoanalytical and Kinetic Studies for the Thermal Stability of Nimesulide Under Different Heating Rates. Orient. J. Chem. 2022, 38, 343–347. [Google Scholar] [CrossRef]
  11. Fuliaş, A.; Vlase, G.; Şoica, C.; Bercean, V.; Vlase, T.; Ledeţi, I. Thermal behaviour of a modified encapsulation agent: Heptakis-6-iodo-6-deoxy-beta-cyclodextrin. J. Therm. Anal. Calorim. 2014, 118, 961–966. [Google Scholar] [CrossRef]
  12. Ledeti, A.; Olariu, T.; Caunii, A.; Vlase, G.; Circioban, D.; Baul, B.; Ledeti, I.; Vlase, T.; Murariu, M. Evaluation of thermal stability and kinetic of degradation for levodopa in non-isothermal conditions. J. Therm. Anal. Calorim. 2018, 131, 1881–1888. [Google Scholar] [CrossRef]
  13. Buda, V.; Andor, M.; Ledeti, A.; Ledeti, I.; Vlase, G.; Vlase, T.; Cristescu, C.; Voicu, M.; Suciu, L.; Tomescu, M.C. Comparative solid-state stability of perindopril active substance vs. pharmaceutical formulation. Int. J. Mol. Sci. 2017, 18, 164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Vyazovkin, S.; Achilias, D.; Fernandez-Francos, X.; Galukhin, A.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for analysis of thermal polymerization kinetics. Thermochim. Acta 2022, 714, 179243. [Google Scholar] [CrossRef]
  15. Muravyev, N.V.; Luciano, G.; Ornaghi, H.L.; Svoboda, R.; Vyazovkin, S. Artificial Neural Networks for Pyrolysis, Thermal Analysis, and Thermokinetic Studies: The Status Quo. Molecules 2021, 26, 3727. [Google Scholar] [CrossRef] [PubMed]
  16. Galukhin, A.; Nikolaev, I.; Nosov, R.; Vyazovkin, S. The Kinetics of Formation of Microporous Polytriazine in Diphenyl Sulfone. Molecules 2022, 27, 3605. [Google Scholar] [CrossRef]
  17. Vyazovkin, S. Determining Preexponential Factor in Model-Free Kinetic Methods: How and Why? Molecules 2021, 26, 3077. [Google Scholar] [CrossRef]
  18. Vyazovkin, S. Kissinger Method in Kinetics of Materials: Things to Beware and Be Aware of. Molecules 2020, 25, 2813. [Google Scholar] [CrossRef]
  19. Liavitskaya, T.; Vyazovkin, S. All You Need to Know about the Kinetics of Thermally Stimulated Reactions Occurring on Cooling. Molecules 2019, 24, 1918. [Google Scholar] [CrossRef] [Green Version]
  20. Kissinger, H.E. Variation of peak temperature with heating rate in differential thermal analysis. J. Res. Natl. Bur. Stand. 1956, 57, 217. [Google Scholar] [CrossRef]
  21. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  22. Sedighi, M.; Nasseri, S.; Ghotbi-Ravandi, A.A. Degradation of 17α-ethinylestradiol by Enterobacter tabaci Isolate and Kinetic Characterization. Environ. Process. 2019, 6, 741–755. [Google Scholar] [CrossRef]
  23. Cai, W.; Li, Y.; Niu, L.; Wang, Q.; Wu, Y. Kinetic study on the cometabolic degradation of 17β-estradiol and 17α-ethinylestradiol by an Acinetobacter sp. strain isolated from activated sludge. Desalin. Water Treat. 2016, 57, 22671–22681. [Google Scholar] [CrossRef]
  24. Nejedly, T.; Klimes, J. A model of natural degradation of 17-α-ethinylestradiol in surface water and identification of degradation products by GC-MS. Environ. Sci. Pollut. Res. 2017, 24, 23196–23206. [Google Scholar] [CrossRef]
  25. Ronderos-Lara, J.G.; Saldarriaga-Noreña, H.; Murillo-Tovar, M.A.; Alvarez, L.; Vergara-Sánchez, J.; Barba, V.; Guerrero-Alvarez, J.A. Distribution and Estrogenic Risk of Alkylphenolic Compounds, Hormones and Drugs Contained in Water and Natural Surface Sediments, Morelos, Mexico. Separations 2022, 9, 19. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Cao, J.; Ke, T.; Tao, Y.; Wu, W.; Wang, P.; Zhou, M.; Chen, L. Pollution Characteristics and Risk Prediction of Endocrine Disruptors in Lakes of Wuhan. Toxics 2022, 10, 93. [Google Scholar] [CrossRef] [PubMed]
  27. D’Amico, R.; Gugliandolo, E.; Cordaro, M.; Fusco, R.; Genovese, T.; Peritore, A.F.; Crupi, R.; Interdonato, L.; Di Paola, D.; Cuzzocrea, S.; et al. Toxic Effects of Endocrine Disruptor Exposure on Collagen-Induced Arthritis. Biomolecules 2022, 12, 564. [Google Scholar] [CrossRef] [PubMed]
  28. Martins, F.; Torrinha, Á.; Delerue-Matos, C.; Morais, S. Life Cycle Assessment and Life Cycle Cost of an Innovative Carbon Paper Sensor for 17α-Ethinylestradiol and Comparison with the Classical Chromatographic Method. Sustainability 2022, 14, 8896. [Google Scholar] [CrossRef]
  29. Luo, Z.; Li, H.; Yang, Y.; Lin, H.; Yang, Z. Adsorption of 17α-ethinylestradiol from aqueous solution onto a reduced graphene oxide-magnetic composite. J. Taiwan Inst. Chem. Eng. 2017, 80, 797–804. [Google Scholar] [CrossRef]
  30. Yasir, M.; Šopík, T.; Patwa, R.; Kimmer, D.; Sedlařík, V. Adsorption of estrogenic hormones in aqueous solution using electrospun nanofibers from waste cigarette butts: Kinetics, mechanism, and reusability. Express Polym. Lett. 2022, 16, 624–648. [Google Scholar] [CrossRef]
  31. Xu, M.; Huang, C.; Lu, J.; Wu, Z.; Zhu, X.; Li, H.; Xiao, L.; Luo, Z. Optimizing adsorption of 17α-ethinylestradiol from water by magnetic mxene using response surface methodology and adsorption kinetics, isotherm, and thermodynamics studies. Molecules 2021, 26, 3150. [Google Scholar] [CrossRef]
  32. Malsawmdawngzela, R.; Tiwari, D. 17α-Ethinylestradiol elimination using synthesized and dense nanocomposite materials: Mechanism and real matrix treatment. Korean J. Chem. Eng. 2022, 39, 646–654. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Chen, Z.; Tao, Y.; Wu, W.; Zeng, Y.; Liao, K.; Li, X.; Chen, L. Transcriptomic and Physiological Responses of Chlorella pyrenoidosa during Exposure to 17α-Ethinylestradiol. Int. J. Mol. Sci. 2022, 23, 3583. [Google Scholar] [CrossRef]
  34. Brown, M.E.; Maciejewski, M.; Vyazovkin, S.; Nomen, R.; Sempere, J.; Burnham, A.; Opfermann, J.; Strey, R.; Anderson, H.L.; Kemmler, A.; et al. Computational aspects of kinetic analysis Part A: The ICTAC Kinetics Project-data, methods and results. Thermochim. Acta 2000, 355, 125–143. [Google Scholar] [CrossRef]
  35. Roduit, B. Computational aspects of kinetic analysis. Part E: The ICTAC Kinetics Project—Numerical techniques and kinetics of solid state processes. Thermochim. Acta 2000, 355, 171–180. [Google Scholar] [CrossRef]
  36. Vyazovkin, S. Computational aspects of kinetic analysis. Part C. The ICTAC Kinetics Project—The light at the end of the tunnel? Thermochim. Acta 2000, 355, 155–163. [Google Scholar] [CrossRef]
  37. Burnham, A.K. Computational aspects of kinetic analysis.: Part D: The ICTAC kinetics project—multi-thermal–history model-fitting methods and their relation to isoconversional methods. Thermochim. Acta 2000, 355, 165–170. [Google Scholar] [CrossRef]
  38. Vyazovkin, S. Recent Advances, Techniques and Applications; Elsevier: Amsterdam, The Netherlands, 2008; Volume 5, ISBN 9780444531230. [Google Scholar]
  39. Budrugeac, P.; Cucos, A. Application of Kissinger, isoconversional and multivariate non-linear regression methods for evaluation of the mechanism and kinetic parameters of phase transitions of type I collagen. Thermochim. Acta 2013, 565, 241–252. [Google Scholar] [CrossRef]
  40. Budrugeac, P.; Criado, J.M.; Gotor, F.J.; Malek, J.; Pérez-Maqueda, L.A.; Segal, E. On the evaluation of the nonisothermal kinetic parameters of (GeS2)0.3(Sb2S3)0.7 crystallization using the IKP method. Int. J. Chem. Kinet. 2004, 36, 309–315. [Google Scholar] [CrossRef]
  41. Moreira, C.G.; Santos, H.G.; Bila, D.M.; da Fonseca, F.V. Assessment of fouling mechanisms on reverse osmosis (RO) membrane during permeation of 17α-ethinylestradiol (EE2) solutions. Environ. Technol. 2021, 1–13. [Google Scholar] [CrossRef]
  42. Minaeva, V.A.; Minaev, B.F.; Hovorun, D.M. Vibrational spectra of the steroid hormones, estradiol and estriol, calculated by density functional theory. The role of low-frequency vibrations. Ukr. Biokhimichnyi Zhurnal 2008, 80, 82–95. [Google Scholar]
  43. Ethinylestradiol—PubChem Profile. Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Ethinylestradiol (accessed on 18 June 2022).
  44. Pheasant, R. Polymorphism of 17-Ethinylestradiol. J. Am. Chem. Soc. 1950, 72, 4303–4304. [Google Scholar] [CrossRef]
  45. Friedman, H.L. Kinetics of thermal degradation of char-foaming plastics from thermogravimetry: Application to a phenolic resin. J. Polym. Sci. 1963, 6, 183–195. [Google Scholar]
  46. Ozawa, T. A New Method of Analyzing Thermogravimetric Data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  47. Koga, N. Ozawa’s kinetic method for analyzing thermoanalytical curves: History and theoretical fundamentals. J. Therm. Anal. Calorim. 2013, 113, 1527–1541. [Google Scholar] [CrossRef]
  48. Flynn, J.H.; Wall, L.A. A quick, direct method for the determination of activation energy from thermogravimetric data. J. Polym. Sci. Part B Polym. Lett. 1966, 4, 323–328. [Google Scholar] [CrossRef]
  49. Flynn, J.H. The isoconversional method for determination of energy of activation at constant heating rates—Corrections for the Doyle approximation. J. Therm. Anal. 1983, 27, 95–102. [Google Scholar] [CrossRef]
  50. Buda, V.; Baul, B.; Andor, M.; Man, D.E.; Ledeţi, A.; Vlase, G.; Vlase, T.; Danciu, C.; Matusz, P.; Peter, F.; et al. Solid State Stability and Kinetics of Degradation for Candesartan—Pure Compound and Pharmaceutical Formulation. Pharmaceutics 2020, 12, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Wang, Q.; Liu, S.-H.; Huang, A.-C.; Huang, C.-F.; Chuang, Y.-K.; Shu, C.-M. Effects of mixing malic acid and salicylic acid with metal oxides in medium- to low-temperature isothermal conditions, as determined using the thermal activity monitor IV. J. Therm. Anal. Calorim. 2018, 133, 779–784. [Google Scholar] [CrossRef]
Figure 1. Structural formula of EE.
Figure 1. Structural formula of EE.
Processes 10 01518 g001
Figure 2. FTIR spectrum of EE on spectral range 4000−400 cm−1.
Figure 2. FTIR spectrum of EE on spectral range 4000−400 cm−1.
Processes 10 01518 g002
Figure 3. Simultaneous−determined TG (thermogravimetric), DTG (derivative thermogravimetric), and DSC (differential scanning calorimetry) curves in air atmosphere at β = 2 °C/min for EE on temperature range 25–500 °C.
Figure 3. Simultaneous−determined TG (thermogravimetric), DTG (derivative thermogravimetric), and DSC (differential scanning calorimetry) curves in air atmosphere at β = 2 °C/min for EE on temperature range 25–500 °C.
Processes 10 01518 g003
Figure 4. Plotting of Kissinger kinetic method for EE.
Figure 4. Plotting of Kissinger kinetic method for EE.
Processes 10 01518 g004
Figure 5. The reaction progress vs. temperature dependence for thermal decomposition of EE: curve 1 corresponds to β = 2 °C/min, curve 2 to β = 4 °C/min, curve 3 to β = 6 °C/min, curve 4 to β = 8 °C/min and curve 5 to β = 10 °C/min, respectively.
Figure 5. The reaction progress vs. temperature dependence for thermal decomposition of EE: curve 1 corresponds to β = 2 °C/min, curve 2 to β = 4 °C/min, curve 3 to β = 6 °C/min, curve 4 to β = 8 °C/min and curve 5 to β = 10 °C/min, respectively.
Processes 10 01518 g005
Figure 6. The reaction rate vs. temperature dependence for thermal decomposition of EE; curve 1 corresponds to β = 2 °C/min, curve 2 to β = 4 °C/min, curve 3 to β = 6 °C/min, curve 4 to β = 8 °C/min and curve 5 to β = 10 °C/min, respectively.
Figure 6. The reaction rate vs. temperature dependence for thermal decomposition of EE; curve 1 corresponds to β = 2 °C/min, curve 2 to β = 4 °C/min, curve 3 to β = 6 °C/min, curve 4 to β = 8 °C/min and curve 5 to β = 10 °C/min, respectively.
Processes 10 01518 g006
Figure 7. Friedman’s linear plotting for decomposition of EE.
Figure 7. Friedman’s linear plotting for decomposition of EE.
Processes 10 01518 g007
Figure 8. Flynn–Wall–Ozawa linear plotting for decomposition of EE.
Figure 8. Flynn–Wall–Ozawa linear plotting for decomposition of EE.
Processes 10 01518 g008
Table 1. Evaluation of the decomposition Ea values vs. conversion degree obtained by the two isoconversional methods and the mean value of Ea.
Table 1. Evaluation of the decomposition Ea values vs. conversion degree obtained by the two isoconversional methods and the mean value of Ea.
Conversion Degree
α
Ea (kJ/mol) vs. α for EE
FrFWO
0.05101.92112.10
0.10108.76114.01
0.15103.54114.02
0.20102.06112.83
0.2599.95111.48
0.3099.61110.32
0.3599.02109.31
0.4098.52108.49
0.4597.81107.75
0.5097.08107.06
0.5596.17106.41
0.6094.96105.75
0.6593.33105.11
0.7091.06104.38
0.7588.28103.56
0.8085.31102.56
0.8582.16101.32
0.9079.3499.74
0.9573.7697.45
E a ¯ (kJ/mol)94.35 ± 9.00107.03 ± 4.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Simu, S.; Ledeţi, A.; Moacă, E.-A.; Păcurariu, C.; Dehelean, C.; Navolan, D.; Ledeţi, I. Thermal Degradation Process of Ethinylestradiol—Kinetic Study. Processes 2022, 10, 1518. https://doi.org/10.3390/pr10081518

AMA Style

Simu S, Ledeţi A, Moacă E-A, Păcurariu C, Dehelean C, Navolan D, Ledeţi I. Thermal Degradation Process of Ethinylestradiol—Kinetic Study. Processes. 2022; 10(8):1518. https://doi.org/10.3390/pr10081518

Chicago/Turabian Style

Simu, Sebastian, Adriana Ledeţi, Elena-Alina Moacă, Cornelia Păcurariu, Cristina Dehelean, Dan Navolan, and Ionuţ Ledeţi. 2022. "Thermal Degradation Process of Ethinylestradiol—Kinetic Study" Processes 10, no. 8: 1518. https://doi.org/10.3390/pr10081518

APA Style

Simu, S., Ledeţi, A., Moacă, E. -A., Păcurariu, C., Dehelean, C., Navolan, D., & Ledeţi, I. (2022). Thermal Degradation Process of Ethinylestradiol—Kinetic Study. Processes, 10(8), 1518. https://doi.org/10.3390/pr10081518

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop