Next Article in Journal
Cyclodextrin Inclusion Complexes for Improved Drug Bioavailability and Activity: Synthetic and Analytical Aspects
Next Article in Special Issue
Polysaccharides as a Hydrophilic Building Block of Amphiphilic Block Copolymers for the Conception of Nanocarriers
Previous Article in Journal
Discovery and Potential Utility of a Novel Non-Invasive Ocular Delivery Platform
Previous Article in Special Issue
Non-spherical Polymeric Nanocarriers for Therapeutics: The Effect of Shape on Biological Systems and Drug Delivery Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physical Ageing of Amorphous Poly(lactic acid)-Indapamide System Studied by Differential Scanning Calorimetry

1
Chair and Department of Pharmaceutical Technology, Poznan University of Medical Sciences, 60-780 Poznan, Poland
2
Department of Chemistry, Rzeszow University of Technology, 35-959 Rzeszow, Poland
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2023, 15(9), 2341; https://doi.org/10.3390/pharmaceutics15092341
Submission received: 2 August 2023 / Revised: 5 September 2023 / Accepted: 14 September 2023 / Published: 19 September 2023

Abstract

:
The process of isothermal and non-isothermal physical ageing of amorphous polylactide (PLA) with the active pharmaceutical ingredient, indapamide (IND), was investigated. A PLA–IND system with a 50/50 weight ratio was obtained and characterized using differential scanning calorimetry (DSC). In the 50/50 (w/w) mixture, two glass transitions were observed: the first at 64.1 ± 0.3 °C corresponding to the glass transition temperature (Tg) of PLA, and the second at 102.6 ± 1.1 °C corresponding to the Tg of IND, indicating a lack of molecular mixing between the two ingredients. The PLA–IND system was subjected to the isothermal physical ageing process at different ageing temperatures (Ta) for 2 h. It was observed that the highest effect of physical ageing (enthalpy relaxation change) on IND in the PLA–IND system occurred at Ta = 85 °C. Furthermore, the system was annealed for various ageing times at 85 °C. The relaxation enthalpies were estimated for each experiment and fitted to the Kohlrausch–Williams–Watts (KWW) equation. The KWW equation allowed for the estimation of the relaxation time and the parameter describing the distribution of relaxation times of the isothermal physical ageing process of IND in the PLA–IND system. The physical ageing of the PLA–IND mixture (50/50) was also discussed in the context of heat capacity. Moreover, the activation energy and fragility parameters were determined for the PLA–IND (50/50) system.

1. Introduction

Active pharmaceutical ingredients (APIs) may exist in a crystalline or amorphous form. The crystalline state is characterised by a regularly ordered lattice structure. In practical terms, the structures of these systems are generally thermodynamically stable and are relatively simple to study using techniques such as differential scanning calorimetry (DSC) or X-ray diffraction methods. On the other hand, there is no long-range order to amorphous forms and their “structures” are not easy to characterise by standard X-ray diffraction methods. A valuable technique for studying amorphous pharmaceuticals is DSC [1,2,3]. Amorphous drugs usually dissolve more readily and are more bioavailable than their crystalline counterparts [2]; however, amorphous APIs may recrystallise during the shelf-life of the formulation [4]. A common approach to stabilise amorphous drugs is polymeric amorphous solid dispersions (PASDs) [5]. In PASDs, the improved stability of an amorphous API is achieved by entrapping the drug in a high-energy glassy state between the polymer chains [6,7]. Although the excipients are often considered inert, it is known that they can interact with APIs, changing their stability, absorption and bioavailability [8,9]. Therefore, APIs must be investigated during preformulation studies at the early phase of the drug development process, in order to provide the necessary information to develop a stable formulation with increased bioavailability [10].
The amorphous forms of API may be desirable due to their improved apparent solubility and, as a consequence, their bioavailability in comparison with its crystalline counterparts. However, in contrast with crystals, glasses are not thermodynamically stable [11,12,13,14,15]. Thus, the stability of amorphous APIs is a primary issue associated with their use in the formulation. The amorphous forms of APIs, during storage below or above the glass transition temperature Tg, may revert to the crystalline form [6,11,16], losing their superior properties. Below the glass transition temperature and above a Kauzman temperature, glasses undergo a physical ageing process, i.e., structural relaxation towards thermodynamic equilibrium as a function of time and temperature [17]. In contrast with chemical or biological ageing, physical ageing is a reversible phenomenon involving the ordering of the amorphous phase, during which no breaking or forming of chemical bonding occurs. More significant structuring may cause, among other things, the deterioration of solubility [18], diffusivity and permeability [19], a decrease in physical stability [20], and a change of mechanical properties [21,22]. For instance, a decrease in solubility has been observed for physically aged cinnarizine-Soluplus solid dispersions [18]. Annealing may also positively affect amorphous APIs, for instance, it can increase chemical stability [23].
In this work, a polymer–drug system was obtained using amorphous poly(lactic acid) and an amorphous API—indapamide.
Poly(lactide) (PLA) belongs to the group of aliphatic polyesters. Depending on D-, L-isomers content, it can exist in a semi-crystalline or an amorphous state exhibiting different physicochemical properties [24,25,26]. It is widely used in the food, pharmaceutical and medical industries. In pharmaceutical applications, PLA is used as a drug carrier matrix in drug delivery systems as well as the bulk component of medical devices due to its biocompatibility and bioresorbability [27,28,29]. PLA is widely used in formulations in order to modify the dissolution profile of an API or to improve its stability [30,31]. The drug can be released from the polymer matrix in a controlled and prolonged manner. For example, Leroueil-Le Verger et al. successfully used solid dispersion with polylactide in an oral controlled release system for isradipine [30].
Indapamide (IND) is a thiazide-like diuretic drug used to treat hypertension [32]. IND is practically insoluble in water [33] and belongs to the biopharmaceutics classification system class II (low solubility, high permeability). Indapamide available on the market is formulated in a crystalline form, e.g., [34], however, it can also be obtained in an amorphous form [33,35,36].
The ability of a material to transform into its amorphous state is called glass-forming ability (GFA) [37]. APIs are categorised into three classes, I, II and III, based on their GFA [37,38]. The GFA of indapamide belongs to class III, where the material, after melting the crystalline form, does not recrystallise during the DSC cooling/heating cycle, and the sample remains amorphous.
As previously mentioned, amorphous materials, unlike crystals, are in a thermodynamic, non-equilibrium state, and, therefore, such materials undergo the physical ageing process. The process is a physical phenomenon related only to the amorphous solid phase. At a temperature lower than the glass transition temperature, the amorphous materials gradually evolve toward the thermodynamic equilibrium of the glassy state [17]. The aged amorphous material tends to be in a more equilibrium state, with less energy [39,40,41]. The physical ageing process can be studied, among other things, using differential scanning calorimetry [36,42]. On the scan obtained from the standard DSC measurement of an unaged amorphous material, a glass transition with a change in specific heat is observed in the heat flow or heat capacity versus temperature plot (Figure 1). When the amorphous sample is aged isothermally at a temperature lower than the glass transition temperature (called the ageing temperature) and is reheated through the glass transition region, an endothermic peak in the DSC curve is observed in addition to the change of heat capacity (Figure 1).
The physical ageing process can be described by measuring the enthalpy relaxation or enthalpy recovery, Δhr, which manifests an order in the amorphous structure and molecular mobility [17,43,44]. Enthalpy relaxation depends on the sample’s thermal history and is the additional energy burden required to reverse the ageing process—to destroy/reverse the ordering resulting from the evolution of the glass towards equilibrium. It can be calculated from the area difference between the heat flow or heat capacity curves of aged (A) and unaged (B) samples (ΔHr = A − B) (see Figure 1). The relaxation kinetic of the non-equilibrium glassy state can be determined using various models [43]. One of these models is the Kohlrausch–Williams–Watts model (KWW) [45,46].
One of the parameters describing the molecular mobility of the amorphous material is the fragility parameter (m). The fragility parameter has been used in the pharmaceutical industry as it can evaluate the tendency of amorphous API towards recrystallisation [14]. In accordance with the classification introduced by Angell, glasses can be referred to as “strong” or “fragile”, depending on the value of m [47]. “Strong” glasses are characterised by a low value of m, while “fragile” materials have a higher m value. Glasses classified as strong are expected to be more physically stable than fragile ones [48].
Heat capacity is an important quantity obtained by DSC for the characterisation of the thermal properties of amorphous systems. Quantitative thermal analysis of non-equilibrium processes and states needs reference baselines of solid (vibrational heat capacity) and liquid heat capacity [49,50,51]. These can be determined for complex systems using a linear combination of the vibrational heat capacity of each ingredient with their molar or weight fractions [49].
The aim of this work was to investigate the amorphous system of API (indapamide) with polymer (polylactide) in the 50:50 weight ratio. The system was not molecularly mixed—two glass transitions were detected. PLA–IND was comprehensively characterized by the differential scanning calorimetry. The process of isothermal physical ageing of IND in the binary PLA–IND system is shown for the different ageing temperatures and times. The kinetics of ageing of a two-component system was analysed by the Kohlrausch–Williams–Watts model. The data for aged and unaged samples of PLA–IND are presented in the frame of reference solid and liquid heat capacity. Furthermore, the activation energies of the structural relaxation at glass transition temperature and fragility parameters were determined for the system.

2. Materials and Methods

2.1. Materials

Polylactide (PLA; Figure 2a) is a non-active pharmaceutical ingredient: a biodegradable and biocompatible polymer containing 16.4% of the D-isomer, making it a fully amorphous material [24,25,26,52]. PLA was obtained from Cargill Dow Nature-Works LLC (Blair, NE, USA).
The indapamide (IND; Figure 2b; M = 365.835 g·mol−1) used in this study was obtained from Polpharma S. A. (Starogard Gdański, Poland). Indapamide was received as a white crystalline powder in the form of a hemihydrate of pharmaceutical grade. The degradation of the studied material begins around 260 °C [35], thus the highest temperature used in this study was 190 °C.

Sample Preparation

An amorphous sample of pure IND was prepared immediately prior to experimental measurements in DSC pan in the instrument, by heating the sample to 190 °C, holding it for 3 min, and then cooling it at a 10 °C·min–1 rate to room temperature. The PLA–IND system (mass ratio of 50/50) was prepared by a solvent evaporation method. The mixture was prepared several times to estimate the repeatability of the results. Amounts of 10–20 mg of IND and PLA were used in equal proportions and dissolved in 50 mL of chloroform, followed by evaporation of the solvent under reduced pressure at about 55–60 °C in a rotary evaporator. The polymer–API system was dried at 40 °C for two days. The obtained samples were in the form of a brittle, thin film.

2.2. Methods

Calorimetric measurements were performed using differential scanning calorimeters (DSC): DSC Q1000 and DSC 2920 from TA Instruments USA (New Castle, DE, USA). The temperature and heat-flow rate were calibrated with indium (Tm = 156.65 °C and Δfh = 28.45 J·g–1) [41], and at least two tests were carried out on each sample. Sapphire (Al2O3) was used to calibrate the heat capacity [41].
Thin film samples of the PLA, PLA–IND system or powdered IND (2–10 mg) were placed into a standard aluminium pan with a pinhole. As a reference sample, an empty aluminium pan with a similar mass to that of the sample pan was used.
The standard DSC test included a heating run to 190 °C in order to erase thermal history, immediate cooling at a rate of 20 °C·min–1 to 5 °C followed by heating with a constant rate of 10 °C·min–1 to 190 °C.

2.2.1. Determination of the Enthalpy Relaxation

To examine the isothermal physical ageing process, the samples, after heating to 190 °C, were cooled down (20 °C·min–1) to the chosen ageing temperature (Ta) for ageing times between 10 min and 20 h. After ageing, the samples were immediately cooled to 5 °C at a 20 °C·min–1 rate, and the subsequent heating scan at a constant rate of 10 °C·min–1 provided the data for the aged sample.
Enthalpy relaxation was calculated from a heat-flow or heat capacity vs. temperature scan for aged and non-aged samples from the DSC measurement according to the Equation:
Δ h r = 1 q T 1 T 2 ( Φ aged Φ unaged ) d T = T 1 T 2 c p   aged c p   unaged d T
where q is the heating rate, T1 and T2 are the limits of the integration temperature (where T1 < Tg < T2), Φaged and Φunaged are heat flows originating from the aged and unaged samples, respectively, and cp aged, and cp unaged are heat capacities of the aged and unaged samples, respectively.
Equilibrium enthalpy relaxation, Δ h r inf , values were calculated by using the Equation:
Δ h r inf = T a T g Δ c p d T   Δ c p   T g T a
where Δcp is the change of heat capacity in solid–liquid transition at a glass temperature and Tg and Ta are the glass transition and ageing temperature, respectively.

2.2.2. Heat Capacity Calculations

The heat capacity (cp) of the PLA–IND (50/50) system was calculated according to the Equation:
c p = Φ q
where Φ is the measured heat-flow rate and q is the heating rate.
In order to compare the experimental heat capacity of the aged and unaged PLA–IND mixtures, the vibration heat capacity and the liquid heat capacity of examined samples were established. The heat capacity of the solid and liquid was estimated from a sum of linear combinations of the weight fractions of the solid and liquid heat capacity of PLA and IND [35,53]. The vibration of the heat capacity in the solid state of the PLA–IND (wt% 50/50) system, cp(vibration), was estimated by using equation [50]:
cp(vibration) = xPLA·cpPLA(vibration) + xIND·cpIND(vibration)
where xPLA and xIND are the weight ratios, and cpPLA(vibration) and cpIND(vibration) are the vibrational heat capacities of PLA and IND in the solid state, respectively. The vibrational heat capacities of PLA and IND have already been published and were referenced in this study from [35,54].
The heat capacity of the liquid state for PLA–IND 50/50 (w/w) system was also calculated as a linear combination of two components as follows:
cp(liquid) = xPLA·cpPLA(liquid) + xIND·cpIND(liquid)
where xPLA and xIND are the weight ratios, cpPLA(liquid) = 1.67 + 1.05 × 10–3 T [J·K–1·g–1] and was referenced from [54], while cpIND (liquid) = 0.954 + 2.30 × 10–3 T [J·K–1·g–1] was referenced from [35].

2.2.3. Determination of Fragility Parameter

The energy activation (Ea) and fragility parameter (m) was obtained from the DSC scans using different cooling rates (0.2–10 °C·min–1) followed by heating with a constant heating rate of 10 °C·min–1. The measurements were performed within the temperature range of 25–190 °C. At 25 °C the sample was held for 5 min and at 190 °C for 3 min. The measurement from which the different cooling rates were used, the glass transition (during heating), is overlapped with the enthalpy relaxation peak caused by the physical ageing process. Due to this phenomenon, the glass transition temperature cannot be determined, so a fictive temperature is used to characterise a sample. In order to assess the fragility parameter, the value of Ea of the glass transition must be known. Ea was calculated based on the following Equation:
E a T g R = d ln q d T 1
where q is the heating/cooling rate, R is the gas constant and T is the fictive temperature (Tf) in kelvin.
Knowing the value of activation energy, the fragility parameter can be described by the following Equation:
m = E a T g 2.303 · R · T g
where Ea is the activation energy and Tg is the glass transition temperature in kelvin.

3. Results and Discussion

Figure 3 shows the heat-flow rate as a function of temperature for PLA, IND and PLA–IND systems with a weight ratio of 50/50, obtained from the second heating scan with a heating rate of 10 °C·min–1. The glass transition of pure PLA is observed at 56.7 ± 0.8 °C. The amorphous form of IND was obtained by a rapid cooling of the sample from the molten state to a temperature below Tg. The glass transition of the amorphous indapamide was observed at a temperature of 103.8 ± 1.1 °C and is in good agreement with the previous reports (102 °C) [35]. In the investigated PLA–IND 50/50 (w/w) system, two glass transitions were observed, the first at a temperature of 64.1 ± 0.3 °C, corresponding to the glass transition of PLA, and the second at 102.6 ± 1.1 °C, corresponding to the Tg of IND. The drug–polymer miscibility may be regarded as the “solubility” of the amorphous drug in the polymer. If the components are molecularly mixed (single amorphous phase), one glass transition is observed [55,56]. Two amorphous phases (two Tg’s) were observed, indicating a lack of molecular miscibility between the components. Two-phase amorphous–amorphous solid dispersion is obtained during quenching. Similar behaviour has been observed for polylactide–valsartan and indomethacin–glucose systems [57,58]. Crystallisation and the melting of samples were not observed in IND and PLA–IND thermograms, suggesting the stability of the materials in the amorphous state against crystallisation even when heated up. The change of the heat capacity at the glass transition temperature for PLA was found to be 0.56 J·g–1·°C–1, for IND 0.44 J·g–1·°C–1, and for the PLA–IND (50/50) system 0.35 J·g–1·°C–1 and 0.145 J·g–1·°C–1 for the polymer and the drug, respectively (see Table 1). The changes in the heat capacity in the glass transition regions for the components in the PLA–IND system are lower than for the pure ingredients, disproportionate to mass ratios, which may indicate intermolecular interactions in the system. Otherwise, Δcp at Tg of the components in the system with a mass ratio of 50/50 would be two times lower than the Δcp in Tg for pure ingredients.
Figure 4 shows a thermogram of PLA–IND 50/50 (w/w) for samples annealed for 2 h at different ageing temperatures below the Tg of IND, i.e., 30, 40, 50, 60, 70, 80, 85 and 90 °C. The annealing at 30 °C did not show any significant changes in enthalpy relaxation in the glass transition region of both the PLA and IND. The changes due to the physical ageing of IND in the PLA–IND system begin to be noticeable for ageing temperature (Ta) at 70 °C and are the highest for Ta = 85 °C. For the annealing at 90 °C, enthalpy relaxation was lower than for 85 °C. This may be due to the fact that, at 90 °C, the material is already at the beginning of the glass transition region. The highest value of the enthalpy relaxation of IND for the PLA–IND 50/50 (w/w) system was observed for the physical ageing process at 85 °C and, therefore, the physical ageing process of the system for longer timescales was carried out at this temperature. It can also be observed that, in the PLA–IND 50/50 (w/w) system, PLA is aged the most at Ta = 50 °C while IND, after 2 h of annealing at this temperature, does not show any change. On the other hand, PLA at 70, 80, 85 and 90 °C is already in its liquid state and thus the physical ageing of the PLA polymer in PLA–IND 50/50 system at this temperature does not occur.
Figure 5 shows DSC curves for the PLA–IND system annealed at 85 °C for ageing times from 10 min to 20 h. It can be seen that the glass transition of polylactide remains unchanged as polylactide is already in a liquid state and, therefore, does not undergo physical ageing. IND underwent a physical ageing process and, as the ageing time increases, the enthalpy relaxation becomes larger and the maximum of the peaks also shifts towards higher temperature values.
The enthalpy relaxation values of indapamide were calculated according to Equation (1), and the results were also compared with the results of the ageing of pure indapamide at the same ageing conditions as those in [36] (see Figure 6). The experimental enthalpy relaxations were fitted using the Kohlrausch–Williams–Watts (KWW) model [45,46], according to the relationship:
Δ h r = Δ h r inf · 1 exp t /   τ KWW β
where Δhr is measured enthalpy relaxation, Δ h r inf is the equilibrium enthalpy relaxation, t is the ageing time, τKWW is the relaxation time, and β is the coefficient describing the distribution of relaxation times.
Figure 6 shows the experimental and calculated enthalpy relaxation (from the best fit) for pure IND and IND from the PLA–IND (50/50) system aged at 85 °C for various ageing times. Presented experimental results are the average values with error bars for the standard deviation (maximum of ±0.3 J·g–1) for pure IND (maximum of ±0.2 J·g–1) and for IND from the PLA–IND system. Initially, the changes in the enthalpy relaxation (ta) are noticeable (relatively significant), and, after about 4 h, the differences between relaxation enthalpies for subsequent ageing times decrease, and the value of enthalpy relaxation approaches the equilibrium enthalpy relaxation. The solid lines in Figure 6 represent the calculated enthalpy relaxation, Δhr, from the best fit of experimental data to the KWW function (Equation (8)). From this best fit, the parameters of relaxation time τKWW and coefficient β were obtained. Both describe the kinetics of the physical ageing process of indapamide in a mixture of PLA and IND. The equilibrium enthalpy relaxation Δ h r inf used in the KWW model (Equation (8)) was estimated according to Equation (2) using the data of Δcp, Ta and Tg from Table 1. All of the calculated parameters from the KWW function are presented in Table 2.
Another parameter used for the description of the physical ageing kinetics is the recovery parameter, φ , which was obtained by transforming Equation (8) to the following form:
φ t = exp t   τ KWW β = 1   Δ h r / Δ h r inf
where all of the quantities in Equation (9) have the same meaning as in Equation (8).
The recovery parameter not only shows changes in enthalpy relaxation occurring during the isothermal physical ageing process but also shows how much more the system can change under given conditions and how fast the material reaches an equilibrium state. The recovery parameter can have values from 1 for non-aged material to 0 for material in which enthalpy relaxation corresponds to equilibrium enthalpy relaxation,   Δ h r inf .
Figure 7 shows the recovery parameter versus the ageing time (ta) for the results of enthalpy recovery presented in Figure 6. The solid lines in Figure 7 present the calculated recovery parameter according to Equation (9) with parameters of τKWW and β obtained from the best fit of the experimental data of Δhr to the KWW Equation (8).
According to the values of τKWW (see Table 2), the recovery parameter lowered much faster for IND in the PLA–IND system than for pure IND. As can be seen, IND in the mixture achieves equilibrium faster than pure API despite lower relaxation enthalpy values. Thus, the changes in enthalpy caused by physical ageing are lower and go down faster for the drug in the system with polylactide compared with the pure API. For IND from the mixture, Δhr and φ reach a saturation relatively quickly and, after 1200 min (20 h), further changes in the enthalpy relaxation are not as significant. Additionally, the recovery parameter after 1200 min of annealing at 85 °C reaches a value of around 0.1, which is close to the equilibrium value. The enthalpy relaxation for pure IND after 1200 min reaches 50% of the equilibrium enthalpy relaxation. It can therefore be assumed that adding the PLA to amorphous IND results in an improvement of the physical stability (in terms of the physical ageing) of the drug because, after 2400 min (40 h) at ageing temperature, further changes of the enthalpy relaxation are very small and therefore changes in physicochemical properties are unlikely to occur. The relaxation time of IND aged at 85 °C is nearly six times higher for the pure drug than for the drug in the PLA–IND (50/50) system. The β parameter may provide a measure of the physical stability of amorphous systems [59]. The glass formers with a broader distribution of structural relaxation times (low β value) may be more susceptible to nucleation resulting in reduced physical stability [60].
The non-isothermal physical ageing of PLA–IND (wt% 50/50) was characterised by activation energy (Ea) and fragility parameter (m).
Figure 8 presents the dependence of the fictive temperature derived by reheating the sample after the experiment with varying rates of cooling together with the linear regression line. The slopes in Figure 8 give the activation energy, Ea, according to Equation (6). Good linearity is observed between the temperatures and cooling rate. The slopes of the Arrhenius plots representing ln(q) as a function of 1/T were calculated using the least-square linear regression. The activation energy (Ea) and the value of kinetic fragility (m) were determined based on Equations (6) and (7), respectively (Table 3).
Relatively low values of the fragility parameter indicate that IND and its system with PLA belong to materials classified as ‘moderately fragile’. A value of m ≈ 200 denotes considerable fragile behaviour, while a value of m ≈ 16 and under indicates a strong glass. Most compounds of pharmaceutical interest have a value of fragility parameter between these two extremes, that is, 16 < m < 200, and are categorised as “moderately fragile” [47].
Figure 9 presents the heat capacity as the function of the temperature of the unaged and aged PLA–IND (50/50) system for 1200 min at 85 °C. The experimental data of cp(unaged) and cp(aged-1200 min) for heat capacity were compared with the reference cp for vibrational, solid heat capacity—cp(vibration)—and liquid heat capacity—cp(liquid)— of the PLA–IND 50/50 system. The heat capacity of the solid and liquid state of the tested mixture was calculated using Equations (4) and (5). Good agreement between the experimental data and the calculated heat capacity of the PLA–IND system was observed below and above glass transition regions.

4. Conclusions

The isothermal and non-isothermal physical ageing processes of IND in a PLA–IND system with a 50/50 mass ratio and pure IND were studied using differential scanning calorimetry. The study demonstrated that the IND and the PLA in the PLA–IND system are immiscible for a mass ratio of 50/50, exhibiting two distinct Tg’s. The presence of PLA in a 50/50 system impacts the kinetics of the physical ageing process of the API.
The isothermal physical ageing process of IND in the PLA–IND system was studied at various ageing temperatures (30–90 °C) for an ageing time of 2 h. It was shown that the physical ageing process of IND in the system at Ta = 85 °C is characterised by the highest enthalpy relaxation. Moreover, the results were compared with data for the physical ageing process at 85 °C for pure IND. The values of the enthalpy relaxation of IND from the PLA–IND (50/50) system were lower than those for pure IND. The experimental data of enthalpy relaxation were fitted to the KWW equation. Relaxation time for IND from the mixture is lower than that for pure IND, which may indicate a decrease in the stability of the drug. Furthermore, IND in the system with PLA, compared with pure IND, reaches a saturation of enthalpy relaxation relatively quickly and, after 20 h at 85 °C, further changes in the enthalpy relaxation due to physical ageing are not noticeable, while the recovery parameter has a value of around 0.1, which is close to the equilibrium.
Non-isothermal physical ageing experiments of IND and IND from the PLA–IND (50/50) system were also carried out. The results show similar values of the activation energy and fragility parameter for IND from the system with PLA compared with the pure IND.
The heat capacity of the PLA–IND system was determined based on the contribution of heat capacities of the components: polylactide and indapamide. The heat capacity of the solid and liquid was estimated from a sum of the linear combinations of the mass fractions of the solid and liquid heat capacity of PLA and IND, respectively. A good agreement between the experimental data and the calculated heat capacity of the PLA–IND system was observed.

Author Contributions

Conceptualisation, M.S. and M.P.; methodology, M.S. and M.P.; formal analysis A.D., M.S. and M.P.; investigation A.D. and M.S.; writing—original draft preparation, A.D., M.S. and M.P.; writing—review and editing, M.S., M.P. and J.L.; supervision J.L. and M.P. 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 required.

Informed Consent Statement

Not required.

Data Availability Statement

Available on request from authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hancock, B.C.; Shamblin, S.L. Molecular mobility of amorphous pharmaceuticals determined using differential scanning calorimetry. Thermochim. Acta 2001, 380, 95–107. [Google Scholar] [CrossRef]
  2. Craig, D.Q.M.; Royall, P.G.; Kett, V.L.; Hopton, M.L. The relevance of the amorphous state to pharmaceutical dosage forms: Glassy drugs and freeze dried systems. Int. J. Pharm. 1999, 179, 179–207. [Google Scholar] [CrossRef] [PubMed]
  3. Duncan, Q.M.; Craig, M.R. (Eds.) Thermal Analysis of Pharmaceuticals; CRC Press: Boca Raton, FL, USA, 2006; ISBN 1420014897/9781420014891. [Google Scholar]
  4. Zhou, D.; Zhang, G.G.Z.; Law, D.; Grant, D.J.W.; Schmitt, E.A. Physical stability of amorphous pharmaceuticals: Importance of configurational thermodynamic quantities and molecular mobility. J. Pharm. Sci. 2002, 91, 1863–1872. [Google Scholar] [CrossRef] [PubMed]
  5. Bikiaris, D.N. Solid dispersions, Part I: Recent evolutions and future opportunities in manufacturing methods for dissolution rate enhancement of poorly water-soluble drugs. Expert Opin. Drug Deliv. 2011, 8, 1501–1519. [Google Scholar] [CrossRef]
  6. Hancock, B.C.; Shamblin, S.L.; Zografi, G. Molecular Mobility of Amorphous Pharmaceutical Solids Below Their Glass Transition Temperatures. Pharm. Sci. 1995, 12, 799–806. [Google Scholar] [CrossRef]
  7. Qian, F.; Huang, J.; Hussain, M.A. Drug-polymer solubility and miscibility: Stability consideration and practical challenges in amorphous solid dispersion development. J. Pharm. Sci. 2010, 99, 2941–2947. [Google Scholar] [CrossRef]
  8. Hancock, B.C.; Zografi, G. Characteristics and Significance of the Amorphous State in Pharmaceutical Systems. J. Pharm. Sci. 1997, 86, 1. [Google Scholar] [CrossRef]
  9. Guo, J.H.; Robertson, R.E.; Amidon, G.L. Influence of Physical Aging on Mechanical Properties of Polymer Free Films: The Prediction of Long-Term Aging Effects on the Water Permeability and Dissolution Rate of Polymer Film-Coated Tablets. Pharm. Res. 1991, 8, 1500–1504. [Google Scholar] [CrossRef]
  10. Jones, T.M. CHAPTER 1: Preformulation Studies. In RSC Drug Discovery Series; Royal Society of Chemistry: London, UK, 2018; Volume 2018-January, pp. 1–20. [Google Scholar]
  11. Yoshioka, M.; Hancock, B.C.; Zografi, G. Crystallization of indomethacin from the amorphous state below and above its glass transition temperature. J. Pharm. Sci. 1994, 83, 1700–1705. [Google Scholar] [CrossRef]
  12. Kissi, E.O.; Grohganz, H.; Löbmann, K.; Ruggiero, M.T.; Zeitler, J.A.; Rades, T. Glass-Transition Temperature of the β-Relaxation as the Major Predictive Parameter for Recrystallization of Neat Amorphous Drugs. J. Phys. Chem. B 2018, 122, 2803–2808. [Google Scholar] [CrossRef]
  13. Hancock, B.C.; Parks, M. What is the true solubility advantage for amorphous pharmaceuticals? Pharm. Res. 2000, 17, 397–404. [Google Scholar] [CrossRef] [PubMed]
  14. Graeser, K.A.; Patterson, J.E.; Zeitler, J.A.; Gordon, K.C.; Rades, T. Correlating thermodynamic and kinetic parameters with amorphous stability. Eur. J. Pharm. Sci. 2009, 37, 492–498. [Google Scholar] [CrossRef] [PubMed]
  15. Dinunzio, J.C.; Miller, D.A.; Yang, W.; McGinity, J.W.; Williams, R.O. Amorphous compositions using concentration enhancing polymers for improved bioavailability of itraconazole. Mol. Pharm. 2008, 5, 968–980. [Google Scholar] [CrossRef] [PubMed]
  16. Fukuoka, E.; Makita, M.; Yamamura, S. Some physicochemical properties of glassy indomethacin. Chem. Pharm. Bull. 1986, 34, 4314–4321. [Google Scholar] [CrossRef] [PubMed]
  17. Struik, L.C.E. Physical Aging in Amorphous Polymers and Other Materials; Elsevier Scientific Pub. Co.: Amsterdam, The Netherlands; New York, NY, USA, 1977; ISBN 0-444-41655-2. [Google Scholar]
  18. Tian, B.; Zhang, L.; Pan, Z.; Gou, J.; Zhang, Y.; Tang, X. A comparison of the effect of temperature and moisture on the solid dispersions: Aging and crystallization. Int. J. Pharm. 2014, 475, 385–392. [Google Scholar] [CrossRef]
  19. Nagai, K.; Sugawara, A.; Kazama, S.; Freeman, B.D. Effects of physical aging on solubility, diffusivity, and permeability of propane and n-butane in poly(4-methyl-2-pentyne). J. Polym. Sci. Part B Polym. Phys. 2004, 42, 2407–2418. [Google Scholar] [CrossRef]
  20. Qi, S.; Avalle, P.; Saklatvala, R.; Craig, D.Q.M. An investigation into the effects of thermal history on the crystallisation behaviour of amorphous paracetamol. Eur. J. Pharm. Biopharm. 2008, 69, 364–371. [Google Scholar] [CrossRef]
  21. Cui, L.; Imre, B.; Tátraaljai, D.; Pukánszky, B. Physical ageing of Poly(lactic acid): Factors and consequences for practice. Polymer 2020, 186, 122014. [Google Scholar] [CrossRef]
  22. Rawat, A.; Burgess, D.J. Effect of physical ageing on the performance of dexamethasone loaded PLGA microspheres. Int. J. Pharm. 2011, 415, 164–168. [Google Scholar] [CrossRef]
  23. Luthra, S.A.; Hodge, I.M.; Utz, M.; Pikal, M.J. Correlation of annealing with chemical stability in lyophilized pharmaceutical glasses. J. Pharm. Sci. 2008, 97, 5240–5251. [Google Scholar] [CrossRef]
  24. Pölöskei, K.; Csézi, G.; Hajba, S.; Tábi, T. Investigation of the thermoformability of various D-Lactide content poly(lactic acid) films by ball burst test. Polym. Eng. Sci. 2020, 60, 1266–1277. [Google Scholar] [CrossRef]
  25. Feng, L.; Zhang, B.; Bian, X.; Li, G.; Chen, Z.; Chen, X. Thermal Properties of Polylactides with Different Stereoisomers of Lactides Used as Comonomers. Macromolecules 2017, 50, 6064–6073. [Google Scholar] [CrossRef]
  26. Ahmed, J.; Varshney, S.K. Polylactides—Chemistry, Properties and Green Packaging Technology: A Review. Int. J. Food Prop. 2011, 14, 37–58. [Google Scholar] [CrossRef]
  27. Auras, R.; Lim, L.; Selke, S.; Tsuji, H. (Eds.) Poly(lactic acid): Synthesis, Structures, Properties, Processing, and Applications; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2011; ISBN 978-1-118-08813-5. [Google Scholar]
  28. Garlotta, D. A literature review of poly(lactic acid). J. Polym. Environ. 2001, 9, 63–84. [Google Scholar] [CrossRef]
  29. Narayanan, G.; Vernekar, V.N.; Kuyinu, E.L.; Laurencin, C.T. Poly (lactic acid)-based biomaterials for orthopaedic regenerative engineering. Adv. Drug Deliv. Rev. 2016, 107, 247–276. [Google Scholar] [CrossRef]
  30. Le Verger, M.L.; Fluckiger, L.; Kim, Y.I.; Hoffman, M.; Maincent, P. Preparation and characterization of nanoparticles containing an antihypertensive agent. Eur. J. Pharm. Biopharm. 1998, 46, 137–143. [Google Scholar] [CrossRef]
  31. Rahman, M.A.; Islam, S. Study of Metoprolol Tartrate delivery from biodegradable polymeric in situ implants for parenteral administration. Int. J. Pharm. Pharm. Sci. 2011, 3, 147–151. [Google Scholar]
  32. Caruso, F.S.; Szabadi, R.R.; Vukovich, R.A. Pharmacokinetics and clinical pharmacology of indapamide. Am. Heart J. 1983, 106, 212–220. [Google Scholar] [CrossRef]
  33. Wojnarowska, Z.; Grzybowska, K.; Hawelek, L.; Dulski, M.; Wrzalik, R.; Gruszka, I.; Paluch, M.; Pienkowska, K.; Sawicki, W.; Bujak, P.; et al. Molecular dynamics, physical stability and solubility advantage from amorphous indapamide drug. Mol. Pharm. 2013, 10, 3612–3627. [Google Scholar] [CrossRef]
  34. FDA Approved Drug Product Lozol: Indapamide. Available online: https://www.accessdata.fda.gov/drugsatfda_docs/label/2009/018538s028lbl.pdf (accessed on 13 September 2023).
  35. Skotnicki, M.; Drogoń, A.; Calvin, J.J.; Rosen, P.F.; Woodfield, B.F.; Pyda, M. Heat capacity and enthalpy of indapamide. Thermochim. Acta 2019, 674, 36–43. [Google Scholar] [CrossRef]
  36. Drogoń, A.; Skotnicki, M.; Skotnicka, A.; Pyda, M. Physical ageing of amorphous indapamide characterised by differential scanning calorimetry. Pharmaceutics 2020, 12, 800. [Google Scholar] [CrossRef] [PubMed]
  37. Blaabjerg, L.I.; Lindenberg, E.; Löbmann, K.; Grohganz, H.; Rades, T. Glass forming ability of amorphous drugs investigated by continuous cooling and isothermal transformation. Mol. Pharm. 2016, 13, 3318–3325. [Google Scholar] [CrossRef] [PubMed]
  38. Baird, J.A.; Van Eerdenbrugh, B.; Taylor, L.S. A classification system to assess the crystallization tendency of organic molecules from undercooled melts. J. Pharm. Sci. 2010, 99, 3787–3806. [Google Scholar] [CrossRef]
  39. McKenna, G.B. Glass Formation and Glassy Behavior. Compr. Polym. Sci. Suppl. 1989, 2, 311–362. [Google Scholar] [CrossRef]
  40. McKenna, G.B.; Simon, S.L. The Glass Transition: Its Measurement and Underlying Physics. 2002. Available online: https://www.nist.gov/publications/glass-transition-its-measurement-and-underlying-physics (accessed on 13 September 2023).
  41. Wunderlich, B. Thermal Analysis of Polymeric Materials; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2005; ISBN 3540236295. [Google Scholar]
  42. Monnier, X.; Saiter, A.; Dargent, E. Physical aging in PLA through standard DSC and fast scanning calorimetry investigations. Thermochim. Acta 2017, 648, 13–22. [Google Scholar] [CrossRef]
  43. Hutchinson, J.M. Physical aging of polymers. Prog. Polym. Sci. 1995, 20, 703–760. [Google Scholar] [CrossRef]
  44. Hodge, I.M. Enthalpy relaxation and recovery in amorphous materials. J. Non. Cryst. Solids 1994, 169, 211–266. [Google Scholar] [CrossRef]
  45. Kohlrausch, R. Theorie des elektrischen Rückstandes in der Leidener Flasche. Ann. Phys. 1854, 167, 179–214. [Google Scholar] [CrossRef]
  46. Williams, G.; Watts, D.C. Non-symmetrical dielectric relaxation behaviour arising from a simple empirical decay function. Trans. Faraday Soc. 1970, 66, 80–85. [Google Scholar] [CrossRef]
  47. Angell, C.A. Relaxation in liquids, polymers and plastic crystals—Strong/fragile patterns and problems. J. Non-Cryst. Solids 1991, 131–133, 13–31. [Google Scholar] [CrossRef]
  48. Tanaka, H. Relationship among glass-forming ability, fragility, and short-range bond ordering of liquids. J. Non. Cryst. Solids 2005, 351, 678–690. [Google Scholar] [CrossRef]
  49. Pyda, M.; Hu, X.; Cebe, P. Heat capacity of silk fibroin based on the vibrational motion of poly(amino acid)s in the presence and absence of water. Macromolecules 2008, 41, 4786–4793. [Google Scholar] [CrossRef]
  50. Pyda, M. Conformational contribution to the heat capacity of the starch and water system. J. Polym. Sci. Part B Polym. Phys. 2001, 39, 3038–3054. [Google Scholar] [CrossRef]
  51. Pyda, M.; Zawada, P.; Drogon, A.; Skotnicki, M.; Cebe, P. Vibrational heat capacity of collagen and collagen–water. J. Therm. Anal. Calorim. 2019, 138, 3389–3401. [Google Scholar] [CrossRef]
  52. Malmgren, T.; Mays, J.; Pyda, M. Characterization of poly(lactic acid) by size exclusion chromatography, differential refractometry, light scattering and thermal analysis. J. Therm. Anal. Calorim. 2006, 83, 35–40. [Google Scholar] [CrossRef]
  53. Pionteck, J.; Pyda, M. Polymers·Part 2: Thermodynamic Properties—pVT-Data and Thermal Properties; Arndt, K.-F., Lechner, M.D., Eds.; Landolt-Boernstein-Polymer; Springer: Berlin/Heidelberg, Germany, 2014; Volume 6A2. [Google Scholar]
  54. Pyda, M.; Bopp, R.C.; Wunderlich, B. Heat capacity of poly(lactic acid). J. Chem. Thermodyn. 2004, 36, 731–742. [Google Scholar] [CrossRef]
  55. Lu, Q.; Zografi, G. Phase behavior of binary and ternary amorphous mixtures containing indomethacin, citric acid, and PVP. Pharm. Res. 1998, 15, 1202–1206. [Google Scholar] [CrossRef]
  56. Rumondor, A.C.F.; Ivanisevic, I.; Bates, S.; Alonzo, D.E.; Taylor, L.S. Evaluation of drug-polymer miscibility in amorphous solid dispersion systems. Pharm. Res. 2009, 26, 2523–2534. [Google Scholar] [CrossRef]
  57. Martínez, L.M.; Videa, M.; López Silva, T.; Castro, S.; Caballero, A.; Lara-Díaz, V.J.; Castorena-Torres, F. Two-phase amorphous-amorphous solid drug dispersion with enhanced stability, solubility and bioavailability resulting from ultrasonic dispersion of an immiscible system. Eur. J. Pharm. Biopharm. 2017, 119, 243–252. [Google Scholar] [CrossRef]
  58. Drogoń, A.; Skotnicki, M.; Pyda, M. Physical aging of polylactide-valsartan system investigated by differential scanning calorimetry. Polimery 2020, 65, 533–541. [Google Scholar] [CrossRef]
  59. Shamblin, S.L.; Hancock, B.C.; Dupuis, Y.; Pikal, M.J. Interpretation of relaxation time constants for amorphous pharmaceutical systems. J. Pharm. Sci. 2000, 89, 417–427. [Google Scholar] [CrossRef]
  60. Tombari, E.; Ferrari, C.; Johari, G.P.; Shanker, R.M. Calorimetric relaxation in pharmaceutical molecular glasses and its utility in understanding their stability against crystallization. J. Phys. Chem. B 2008, 112, 10806–10814. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Scheme for determination of enthalpy relaxation from DSC plot of heat capacity vs. temperature.
Figure 1. Scheme for determination of enthalpy relaxation from DSC plot of heat capacity vs. temperature.
Pharmaceutics 15 02341 g001
Figure 2. Chemical structure of (a) polylactide and (b) indapamide.
Figure 2. Chemical structure of (a) polylactide and (b) indapamide.
Pharmaceutics 15 02341 g002
Figure 3. Heat-flow rate as a function of temperature for polylactide (PLA; blue curve), indapamide (IND; red curve) and the PLA–IND 50/50 (w/w) system (purple curve). Two glass transitions (Tg) are observed for the PLA–IND system, indicating the immiscibility of the components.
Figure 3. Heat-flow rate as a function of temperature for polylactide (PLA; blue curve), indapamide (IND; red curve) and the PLA–IND 50/50 (w/w) system (purple curve). Two glass transitions (Tg) are observed for the PLA–IND system, indicating the immiscibility of the components.
Pharmaceutics 15 02341 g003
Figure 4. Heat-flow rate vs. temperature for the isothermal physical ageing showing enthalpy relaxation changes of the PLA–IND 50/50 (w/w) system after annealing for 2 h at different ageing temperatures, Ta = (30, 40, 50, 60, 70, 80, 85, 90) °C. The curve for the sample aged at 30 °C was shifted for clarity—no significant changes in enthalpy relaxation occurred.
Figure 4. Heat-flow rate vs. temperature for the isothermal physical ageing showing enthalpy relaxation changes of the PLA–IND 50/50 (w/w) system after annealing for 2 h at different ageing temperatures, Ta = (30, 40, 50, 60, 70, 80, 85, 90) °C. The curve for the sample aged at 30 °C was shifted for clarity—no significant changes in enthalpy relaxation occurred.
Pharmaceutics 15 02341 g004
Figure 5. Heat-flow rate as a function of the temperature of the PLA–IND (50/50) system after annealing at 85 °C for various ageing times.
Figure 5. Heat-flow rate as a function of the temperature of the PLA–IND (50/50) system after annealing at 85 °C for various ageing times.
Pharmaceutics 15 02341 g005
Figure 6. Experimental enthalpy relaxation for IND in the PLA–IND (50/50) system (dots) and for pure IND (squares) aged at 85 °C as a function of ageing time. The solid lines present the calculated enthalpy relaxation from the best fit of the experimental data to the KWW equation.
Figure 6. Experimental enthalpy relaxation for IND in the PLA–IND (50/50) system (dots) and for pure IND (squares) aged at 85 °C as a function of ageing time. The solid lines present the calculated enthalpy relaxation from the best fit of the experimental data to the KWW equation.
Pharmaceutics 15 02341 g006
Figure 7. The recovery parameter, φ , values for IND in the PLA–IND (50/50) system (dots) and for pure IND (squares) aged at 85 °C as a function of ageing time. Solid lines represent the fit of the calculated data to Equation (9).
Figure 7. The recovery parameter, φ , values for IND in the PLA–IND (50/50) system (dots) and for pure IND (squares) aged at 85 °C as a function of ageing time. Solid lines represent the fit of the calculated data to Equation (9).
Pharmaceutics 15 02341 g007
Figure 8. Arrhenius plots of temperature dependence on the cooling rates of the IND in the PLA–IND (50/50) system and of pure IND.
Figure 8. Arrhenius plots of temperature dependence on the cooling rates of the IND in the PLA–IND (50/50) system and of pure IND.
Pharmaceutics 15 02341 g008
Figure 9. Heat capacity of the PLA–IND (50/50) system studied by DSC, where cp(aged-1200 min) and cp(unaged) represent the experimental heat capacity of the aged and unaged sample, cp(vibration) represents the vibrational heat capacity of the PLA–IND (50/50) system and cp(liquid) represents the heat capacity of the liquid state of PLA–IND (50/50).
Figure 9. Heat capacity of the PLA–IND (50/50) system studied by DSC, where cp(aged-1200 min) and cp(unaged) represent the experimental heat capacity of the aged and unaged sample, cp(vibration) represents the vibrational heat capacity of the PLA–IND (50/50) system and cp(liquid) represents the heat capacity of the liquid state of PLA–IND (50/50).
Pharmaceutics 15 02341 g009
Table 1. The glass transition temperature (Tg) and changes of the heat capacity (Δcp) at Tg for PLA, IND and the PLA–IND 50/50 (w/w) system.
Table 1. The glass transition temperature (Tg) and changes of the heat capacity (Δcp) at Tg for PLA, IND and the PLA–IND 50/50 (w/w) system.
Sample (w%/w%)Tg ± SD/°CΔcp at Tg ± SD/J·g–1·°C–1
PLA (100/0)56.7 ± 0.80.56 ± 0.01
PLA–IND (50/50)64.1 ± 0.3/102.6 ± 1.10.35 ± 0.01/0.145 ± 0.01
IND (0/100)103.9 ± 1.10.44 ± 0.01
Table 2. KWW parameters for pure IND and IND from the PLA–IND 50/50 (w/w) system aged at 85 °C.
Table 2. KWW parameters for pure IND and IND from the PLA–IND 50/50 (w/w) system aged at 85 °C.
Sample Δ h r i n f   /J·g−1 β /– τ K W W   /h
IND *7.920.3328
IND in the PLA–IND (50/50) system2.560.464.75
* taken from Ref. [36].
Table 3. Energy activation and fragility parameters of the IND from the PLA–IND (50/50) system and of pure IND obtained from the dependence of Tf on heating rate.
Table 3. Energy activation and fragility parameters of the IND from the PLA–IND (50/50) system and of pure IND obtained from the dependence of Tf on heating rate.
SampleEnergy Activation,
ΔEa ± SD/kJ/mol
Fragility Parameter,
m ± SD/–
IND in the PLA–IND (50/50) system452.5 ± 18.563 ± 3
IND473 ± 5166 ± 7
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

Skotnicki, M.; Drogoń, A.; Lulek, J.; Pyda, M. Physical Ageing of Amorphous Poly(lactic acid)-Indapamide System Studied by Differential Scanning Calorimetry. Pharmaceutics 2023, 15, 2341. https://doi.org/10.3390/pharmaceutics15092341

AMA Style

Skotnicki M, Drogoń A, Lulek J, Pyda M. Physical Ageing of Amorphous Poly(lactic acid)-Indapamide System Studied by Differential Scanning Calorimetry. Pharmaceutics. 2023; 15(9):2341. https://doi.org/10.3390/pharmaceutics15092341

Chicago/Turabian Style

Skotnicki, Marcin, Agata Drogoń, Janina Lulek, and Marek Pyda. 2023. "Physical Ageing of Amorphous Poly(lactic acid)-Indapamide System Studied by Differential Scanning Calorimetry" Pharmaceutics 15, no. 9: 2341. https://doi.org/10.3390/pharmaceutics15092341

APA Style

Skotnicki, M., Drogoń, A., Lulek, J., & Pyda, M. (2023). Physical Ageing of Amorphous Poly(lactic acid)-Indapamide System Studied by Differential Scanning Calorimetry. Pharmaceutics, 15(9), 2341. https://doi.org/10.3390/pharmaceutics15092341

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