# The Thermodynamic Dissociation Constants of Clotrimazole, Terbinafine HCL, Acetylsalicylic Acid, Salicylic Acid, and Galanthamine by the Nonlinear Regression of Multiwavelength Spectrophotometric pH-Titration Data

Authors: Milan Meloun; Sylva Bordovská; Lubomír Galla
Journal: SRX Pharmacology
 (2010)

---

## Abstract

The mixed dissociation constants of five drugs—clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid, and galanthamine—at various ionic strengths I and at temperatures of 25°C and 37°C were determined with the use of multiwavelength and multivariate treatments of spectral data SPECFIT/32 nonlinear regression analysis and INDICES factor analysis. The factor analysis in the INDICES program predicts the number of components, when the data quality is high and the instrumental error is known. The thermodynamic dissociation constant pKaT was estimated by nonlinear regression of {pKa,I} data at 25°C and 37°C: for clotrimazole  pKa,1T = 4.38(1) and 4.16(3); for terbinafine HCl pKa,1T = 4.19(3) and 4.12(5); for acetylsalicylic acid pKa,1T = 3.49(25) and 3.41(15); for salicylic acid pKa,1T = 3.01(1) and 3.00(1) and for galanthamine pKa,1T = 8.21(1) and 7.99(2) where in brackets the standard deviation is in  the last significant digits. Goodness-of-fit tests for various regression diagnostics enabled the reliability of the parameter estimates found to be proven. Pharma Algorithms predicts pKa being based on the structural formulae of drug compounds.


---

## Body

## 1. Introduction

Acid dissociation constants (pKa values) are key parameters to predict the extent of ionization of functional groups with respect to pH and also to explain chemical phenomena such as absorption, distribution, elimination of substances, biological activity, biological uptake, biological transport and environmental fate [1–3]. This information is important in drug discovery and development since the pharmacokinetic and pharmacodynamic properties of different protonation/dissociation forms of the drug molecules may vary considerably. There are several methods for the determination of dissociation constants. Traditionally, potentiometry and UV-VIS spectrometry have been the most useful techniques for the determination of equilibrium constants. Spectrophotometric methods in combination with suitable chemometric tools can be used for the determination of acid dissociation constants pKa even for barely soluble drugs provided that the compound under consideration possesses chromophore(s) in proximity to the ionization centre, but the absorptivity should also change significantly on (de)protonation, indicating that the ionization centre is a part of the chromophore. The information can be extracted if multivariate spectrophotometric data are analyzed by means of an appropriate multivariate data analysis method. Hard-modelling methods, for example, SQUAD(84) [4–11], include traditional least squares curve fitting approaches, based on a previous postulation of the chemical (here protonation) model. The postulations are a set of variously protonated species defined by their stoichiometric coefficients and protonation constants, which are refined by the least-squares minimization. However, a hard-modelling analysis cannot be applied if crucial information is missing. Soft-modelling or model-free approaches, for example, SPECFIT [12–16], are based on much more general prerequisites, such as positive molar absorbance, positive concentration of all species, unimodality of concentration profiles, and closure, that is, a concentration of all species is the same for all solutions. Naturally, if the strengths of hard-modelling and soft-modelling methodologies are combined, a much more powerful method of data analysis can be expected.In addition, relevant software has been developed for the rapid estimation of pKa values based on the chemical structure [17]. These techniques have certain advantages, for example, calculations can be performed on large virtual compound libraries. By entering the compound topological structure descriptors graphically, pKa values of drug compounds are predicted using approximately hundreds of Hammett and Taft equations and quantum chemistry calculus. Still, erroneous data are often predicted for complex and flexible drug compounds containing several functional groups. Furthermore, these calculations are based on parameters in databases containing experimental data from the literature. Hence, sufficient data for new types of compounds, to give accurate predictions, may be missing.An antifungal drug is a medication used to treat fungal infections such as athlete’s foot, ringworm, candidiasis (thrush), serious systemic infection such as cryptococcal meningitis, and other. Such drugs are usually obtained by a doctor’s prescription or purchased over-the-counter. Antifungals work by exploiting differences between mammalian and fungal cells to kill off the fungal organism without dangerous effects on the host. Unlike bacteria, both fungi and humans are eukaryotes. The basic structure of fungal cells and human cells is nearly identical. This means it is more difficult to find a target for an antifungal drug to attack that does not also exist in the infected organism. Consequently, there are often side-effects to some of these drugs. Some of these side-effects can be life-threatening if not used properly.Clotrimazole (CLO) (chemically 1-[(2-chlorophenyl)-diphenyl-methyl]imidazole, CAS no. 23593-75-1, molecular formula: C22H11ClN2, molecular weight: 344.837, description: a colourless, odourless, tasteless and crystalline solid, solubility: it is practically insoluble in water (<0.01 mg/mL), soluble in chloroform and methanol, in ethanol and in diethyl ether, it is freely soluble in acetone and methyl alcohol, is a weak base having a pKa value of 4.7) is of the following structure:Clotrimazole, an imidazole derivative with a broad spectrum of antimycotic activity, inhibits biosynthesis of the sterol ergostol, an important component of fungal cell membranes. Its action leads to increased membrane permeability and apparent disruption of enzyme systems bound to the membrane. Betamethasone and clotrimazole are used together to treat cutaneous tinea infections. It is a well-established drug used in dermatology and gynaecology, available in the form of solution formulations, tablet and cream [18–20].Terbinafine hydrochloride(TER-HCl) (chemically (E)-N-(6,6-dimethyl-2-hepten-4-ynyl)-N-methyl-1-naphthalene methanamine hydrochloride, CAS no. 78628-80-5, molecular formula: C21H25ClN, molecular weight: 327.90, description: white to off-white finely crystalline powder, solubility: soluble in ethanol (45 mg/mL), DMSO (30 mg/mL), water (3 mg/mL)) is of the following structure:Terbinafine hydrochloride is a new potent antifungal agent of the allylamine class that selectively inhibits fungal squalene epoxidase. The drug has broad-spectrum activity against yeast, fungi, moulds, and dermatophytes and is indicated for both oral and topical treatment of mycoses [21–23].The salicylates are a group of analgesics, or painkilling drugs, that are derivatives of salicylic acid. The best known is acetylsalicylic acid, or aspirin. Now often made synthetically, they were originally derived fromsalicin, the active ingredient in willow bark, used for centuries in the treatment of pain and fever. Salicylates also occur naturally in many plants used as foods (e.g., strawberries, almonds, tomatoes). Methyl salicylate is the main component of wintergreen, sweet birch, gautheria, and betula oils; the compound is used in rubbing liniments to soothe muscular aches and as a flavouring. Sodium salicylate, traditionally used in the treatment of arthritis, is also used in dyes and as a non-edible preservative. In general, salicylates, especially aspirin, are used medically to reduce fever and inflammation and to relieve headache, menstrual pain, and pain in nerves, muscles, and joints. Because of the effects of salicylates on blood platelets and clotting, aspirin is often prescribed prophylactically for those at risk of stroke or heart attack. Salicylates are useful, relatively safe drugs, but normal doses can cause gastrointestinal disturbances in sensitive patients and large doses can be toxic or fatal, especially to children.Acetylsalicylic Acid (ASA) (chemically 2-Acetoxybenzoic acid, CAS no. 50-78-2, description: colourless or a white crystalline powder, molecular formula: C9H8O4, molecular weight: 180.16, melting point: 138–140°C, solubility: 1 in 300 of water, 1 in 5–7 in alcohol, 1 in 17 of chloroform and 1 in 20 of ether; soluble in solutions of acetates and citrates and, with decomposition, in solutions of alkali hydroxides and carbonates, pKa is 3.49 at 25°C [24]) of the structureis arguably the world’s oldest and best known pharmaceutical product with accepted anti-inflammatory, anti-thrombotic, anti-pyretic, anti-oxidant and analgesic properties [25–27]. It has been suggested that daily use of aspirin can reduce the risk of some types of cancer [28], and even provide a means of extending human life [29]. Aspirin was the first-discovered member of the class of drugs known as nonsteroidal anti-inflammatory drugs (NSAIDs), not all of which are salicylates, although they all have similar effects and a similar mechanism of action.Salicylic Acid (SA) (chemically 2-Hydroxybenzoic acid, CAS no. 69-72-7, description: colourless or white crystalline powder, molecular formula: C7H6O3, molecular weight: 138.12, melting point: 160°C, solubility: soluble 1 in 460 to 550 of water, 1 in 15 of boiling water, 1 in 3 to 4 in alcohol, 1 in 3 in ether, and 1 in 45 in chloroform, pKa is 3.00 at 25°C) of the structureis a beta hydroxy acid (BHA). This colourless crystalline organic acid is widely used in organic synthesis and functions as a plant hormone. It is derived from the metabolism of salicin. It is probably best known as a compound that is chemically similar but not identical to the active component of aspirin. It is widely used as kerotalic, antimicrobial and antifungal agent and as an external therapeutical agent (keratolytic) in many pharmaceutical preparations [30, 31].Galanthamine(GAL) belongs to a class of acetylcholinesterase inhibitors approved for symptomatic treatment of Alzheimer’s disease. Galanthamine (chemically (4aS,6R,8aS)-4a,5,9,10,11,12-Hexahydro-3-methoxy-11-methyl-6H-benzofuro[3a,3,2-ef][2]benzazepin-6-ol, CAS no. 357-70-0, description: a white crystalline powder, molecular formula: C17H21NO3, molecular weight: 287.358, melting point: 127-128°C, solubility: poorly soluble in water, pKa is 8.2 at 25°C) of structureis a parasympathomimetic, specifically, a reversible cholinesterase inhibitor. It is indicated for the treatment of mild to moderate dementia of the Alzheimer's type. Galanthamine is postulated to exert its therapeutic effect by enhancing cholinergic function. This is accomplished by increasing the concentration of acetylcholine through reversible inhibition of its hydrolysis by acetylcholinesterase. If this proposed mechanism of action is correct, galanthamine's effect may lessen as the disease process advances and fewer cholinergic neurons remain functionally intact. There is no evidence that galanthamine alters the course of the underlying dementing process [32–36].The procedure for the determination of the mixed protonation/dissociation constants has been described previously [37–40]. The details of the computer data treatment are collected in the Supporting Information. The nonlinear estimation of the thermodynamic dissociation constant KaT=aH+aL-/aHL, is simply a problem of optimization in the parameter space in which pKa  and I are known and given values, while the parameters pKaT, a, and C are unknown variables to be estimated [41, 42]. The adequacy of a proposed regression model with experimental data and the reliability of parameter estimates pKa,i found, being denoted for the sake of simplicity as bj, and εij, j=1,…,m, may be examined by a goodness-of-fit test,cf. a previous tutorial [37–40]. Repeatability and reproducibility (R&R) are critical factors closely related to precision and accuracy. It helps to think of repeatability in terms of how capable the gage is of providing the same reading to a single user when measuring a specific sample. Repeatability is the variability of the measurements obtained by one person while measuring the same item repeatedly. This is also known as the inherent precision of the measurement equipment. Reproducibility is the variability of the measurement system caused by differences in operator behaviour. Mathematically, it is the variability of the average values obtained by several operators while measuring the same item.The aim of the study was to determine thepKaT values of clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid, and galanthamine using pH–spectrophotometric titrations. The pKa values were also calculated theoretically using computer program (Pharma Algorithms) making predictions based on the structural formula of drug compounds.


## 2. Experimental

### 2.1. Materials

Clotrimazolewas purchased from AMOLI ORGANICS Ltd., with a purity of 100%. Terbinafine HCl was purchased from CHEMAGIS, with a purity of 100%. Acetylsalicylic Acidwas purchased from Sigma-Aldrich Co., with a purity of  ≥99%. Salicylic Acid was purchased from Sigma-Aldrich Co., with a purity of   ≥99%. Galanthamine was purchased from Ivax-Pharmaceuticals s.r.o., with a purity of   ≥99%. Perchloric acid, 1 M, was prepared from conc. HClO4 (p.a., Lachema Brno) using redistilled water and standardized against HgO and NaI with reproducibility of less than 0.20%. Sodium hydroxide, 1 M, was prepared from pellets (p.a., Aldrich Chemical Company) with carbondioxidefree redistilled water and standardized against a solution of potassium hydrogen-phthalate using the Gran Method with a reproducibility of 0.1%. Mercuric oxide, sodium iodide, and sodium perchlorate (p.a. Lachema Brno) were not further purified. The preparation of other solutions from analytical reagent-grade chemicals has been described previously [37–40].


### 2.2. Apparatus and Procedure

The apparatus used and the pH-spectrophotometric titration procedure described in [37–40] were applied. The microscale pH-titrimetric method generally decreases the amount of drug sample required but the duration of experiment is as long as that of a conventional titration. This method uses pH electrodes which requires the calibration of the pH electrodes cell every day. Calibration procedures normally take 20 minutes and the E0 and slope of the electrode remain valid for less than 1 day. The free hydrogen-ion concentration h=[H+] was measured  via emfon an OP-208/1 digital voltmeter (Radelkis, Budapest) with a precision of ±0.1 mV using a G202B glass electrode (Radiometer, Copenhagen) and an OP-8303P commercial SCE reference electrode (Radelkis, Budapest). Although spectrophotometric multiple-wavelength pH-titration has gained popularity in recent years, the method also requires precise pH measurement: an aqueous solution 20.00 cm3 containing 10-5 mol.dm-3 drug, 0.100 mol·dm-3 hydrochloric acid and 10 cm3 indifferent solution KCl for adjustment of ionic strength was titrated with standard 1.0 mol·dm-3 KOH at 298 K and from 20 to 50 absorption spectra were recorded. Titrations were performed in a thermostatic water-jacketed double-walled glass vessel of 100 mL, closed with a Teflon bung containing the electrodes, an argon inlet, a thermometer, a propeller stirrer and a capillary tip from a micro-burette. All pH measurements were carried out at 25.0°C ± 0.1° and 37.0°C ± 0.1°. As all methods involving the use of pH electrodes have the same limitation, namely, relatively extensive experimental time. During titration, the average pH meter responses time can be as long as 0.5 minute for each recorded point even with the best pH electrodes. Therefore, a typical titration process takes approximately 30 min to cover the whole pH range with 0.1 pH unit intervals necessary to provide the high quality data suitable for analysis. As the uniformity of the solution can be achieved quickly by vigorous stirring, and typical proton dissociation processes are less than 10-8 s, much time during pH titration experiments is wasted, waiting for the electrode to provide a stable response.When the drug was titrated, a stream of argon gas was bubbled through the solution both to stir and to maintain an inert atmosphere. The argon was passed through an aqueous ionic medium by prior passage through one or two vessels also containing the titrand medium before entering the corresponding titrand solution. The burettes used were syringe micro-burettes of 1250μL capacity (META, Brno) with a 2.50 cm micrometer screw, [43]. The polyethylene capillary tip of the micro-burette was immersed into the solution when adding reagent, but withdrawn after each addition in order to avoid leakage of the reagent during the pH read out. The micro-burette was calibrated by ten replicate determinations of the total volume of delivered water by weighing on a Sartorius 1712 MP8 balance with results evaluated statistically, leading to a precision of ±0.015% in added volume over the whole volume range. The solution was pumped into the cuvette and spectrophotometric measurement was performed with the use of a Cintra 40 spectrophotometer (GBC, Australia).The computation scheme for the determination of the protonation constants of the multicomponent system is taken from Meloun et al.,cf. [1, page 226] and the five steps are described elsewhere [37–40, 44–48]: (1) instrumental error of absorbance measurements, sinst(A), (2) experimental design, (3) number of light-absorbing species with a factor analysis, (4) choice of computational strategy of the nonlinear regression, (5) diagnostics indicating a correct chemical (protonation) model: when a minimization process terminates, some curve-fitting diagnostics are examined to determine whether the results should be accepted: the physical meaning of parametric estimates, the physical meaning of the species concentrations, the goodness-of-fit test of spectra calculated, and the deconvolution of the experimental spectra set.


### 2.3. Computations

Computations relating to the determination of dissociation constants were performed by regression analysis of the UV/VIS spectra using the SPECFIT/32 [12–16] programmes. Most of graphs were plotted using ORIGIN 7.5 [49] and S-Plus [50]. The thermodynamic dissociation constant pKaT was estimated with the MINOPT nonlinear regression programme in the ADSTAT statistical system [51]. A qualitative interpretation of the spectra with the use of the INDICES programme [52] aims to evaluate the quality of the dataset and remove spurious data, and to estimate the minimum number of factors, that is, contributing aqueous species, which are necessary to describe the experimental data and determine the number of dominant species present in the equilibrium mixture. Pharma Algorithms [17] is a programme for making predictions based on the structural formulae of drug compounds. Entering the compound topological structure descriptors graphically, pKa values of organic compound are predicted using approximately hundreds of Hammett and Taft equations and quantum chemistry calculus.


### 2.4. Supporting Information Available

Complete experimental and computational procedures [37–40, 44–48], input data specimens, and corresponding output in numerical and graphical form for the programmes, INDICES, SQUAD(84), and SPECFIT/32 are available free of charge on line http://meloun.upce.cz at and in the block DOWNLOAD and block DATA.


## 2.1. Materials

Clotrimazolewas purchased from AMOLI ORGANICS Ltd., with a purity of 100%. Terbinafine HCl was purchased from CHEMAGIS, with a purity of 100%. Acetylsalicylic Acidwas purchased from Sigma-Aldrich Co., with a purity of  ≥99%. Salicylic Acid was purchased from Sigma-Aldrich Co., with a purity of   ≥99%. Galanthamine was purchased from Ivax-Pharmaceuticals s.r.o., with a purity of   ≥99%. Perchloric acid, 1 M, was prepared from conc. HClO4 (p.a., Lachema Brno) using redistilled water and standardized against HgO and NaI with reproducibility of less than 0.20%. Sodium hydroxide, 1 M, was prepared from pellets (p.a., Aldrich Chemical Company) with carbondioxidefree redistilled water and standardized against a solution of potassium hydrogen-phthalate using the Gran Method with a reproducibility of 0.1%. Mercuric oxide, sodium iodide, and sodium perchlorate (p.a. Lachema Brno) were not further purified. The preparation of other solutions from analytical reagent-grade chemicals has been described previously [37–40].


## 2.2. Apparatus and Procedure

The apparatus used and the pH-spectrophotometric titration procedure described in [37–40] were applied. The microscale pH-titrimetric method generally decreases the amount of drug sample required but the duration of experiment is as long as that of a conventional titration. This method uses pH electrodes which requires the calibration of the pH electrodes cell every day. Calibration procedures normally take 20 minutes and the E0 and slope of the electrode remain valid for less than 1 day. The free hydrogen-ion concentration h=[H+] was measured  via emfon an OP-208/1 digital voltmeter (Radelkis, Budapest) with a precision of ±0.1 mV using a G202B glass electrode (Radiometer, Copenhagen) and an OP-8303P commercial SCE reference electrode (Radelkis, Budapest). Although spectrophotometric multiple-wavelength pH-titration has gained popularity in recent years, the method also requires precise pH measurement: an aqueous solution 20.00 cm3 containing 10-5 mol.dm-3 drug, 0.100 mol·dm-3 hydrochloric acid and 10 cm3 indifferent solution KCl for adjustment of ionic strength was titrated with standard 1.0 mol·dm-3 KOH at 298 K and from 20 to 50 absorption spectra were recorded. Titrations were performed in a thermostatic water-jacketed double-walled glass vessel of 100 mL, closed with a Teflon bung containing the electrodes, an argon inlet, a thermometer, a propeller stirrer and a capillary tip from a micro-burette. All pH measurements were carried out at 25.0°C ± 0.1° and 37.0°C ± 0.1°. As all methods involving the use of pH electrodes have the same limitation, namely, relatively extensive experimental time. During titration, the average pH meter responses time can be as long as 0.5 minute for each recorded point even with the best pH electrodes. Therefore, a typical titration process takes approximately 30 min to cover the whole pH range with 0.1 pH unit intervals necessary to provide the high quality data suitable for analysis. As the uniformity of the solution can be achieved quickly by vigorous stirring, and typical proton dissociation processes are less than 10-8 s, much time during pH titration experiments is wasted, waiting for the electrode to provide a stable response.When the drug was titrated, a stream of argon gas was bubbled through the solution both to stir and to maintain an inert atmosphere. The argon was passed through an aqueous ionic medium by prior passage through one or two vessels also containing the titrand medium before entering the corresponding titrand solution. The burettes used were syringe micro-burettes of 1250μL capacity (META, Brno) with a 2.50 cm micrometer screw, [43]. The polyethylene capillary tip of the micro-burette was immersed into the solution when adding reagent, but withdrawn after each addition in order to avoid leakage of the reagent during the pH read out. The micro-burette was calibrated by ten replicate determinations of the total volume of delivered water by weighing on a Sartorius 1712 MP8 balance with results evaluated statistically, leading to a precision of ±0.015% in added volume over the whole volume range. The solution was pumped into the cuvette and spectrophotometric measurement was performed with the use of a Cintra 40 spectrophotometer (GBC, Australia).The computation scheme for the determination of the protonation constants of the multicomponent system is taken from Meloun et al.,cf. [1, page 226] and the five steps are described elsewhere [37–40, 44–48]: (1) instrumental error of absorbance measurements, sinst(A), (2) experimental design, (3) number of light-absorbing species with a factor analysis, (4) choice of computational strategy of the nonlinear regression, (5) diagnostics indicating a correct chemical (protonation) model: when a minimization process terminates, some curve-fitting diagnostics are examined to determine whether the results should be accepted: the physical meaning of parametric estimates, the physical meaning of the species concentrations, the goodness-of-fit test of spectra calculated, and the deconvolution of the experimental spectra set.


## 2.3. Computations

Computations relating to the determination of dissociation constants were performed by regression analysis of the UV/VIS spectra using the SPECFIT/32 [12–16] programmes. Most of graphs were plotted using ORIGIN 7.5 [49] and S-Plus [50]. The thermodynamic dissociation constant pKaT was estimated with the MINOPT nonlinear regression programme in the ADSTAT statistical system [51]. A qualitative interpretation of the spectra with the use of the INDICES programme [52] aims to evaluate the quality of the dataset and remove spurious data, and to estimate the minimum number of factors, that is, contributing aqueous species, which are necessary to describe the experimental data and determine the number of dominant species present in the equilibrium mixture. Pharma Algorithms [17] is a programme for making predictions based on the structural formulae of drug compounds. Entering the compound topological structure descriptors graphically, pKa values of organic compound are predicted using approximately hundreds of Hammett and Taft equations and quantum chemistry calculus.


## 2.4. Supporting Information Available

Complete experimental and computational procedures [37–40, 44–48], input data specimens, and corresponding output in numerical and graphical form for the programmes, INDICES, SQUAD(84), and SPECFIT/32 are available free of charge on line http://meloun.upce.cz at and in the block DOWNLOAD and block DATA.


## 3. Results and Discussion

### 3.1. Estimation of Dissociation Constants

#### 3.1.1. Clotrimazole

A suggested strategy for effective experimentation in dissociation constants determination ensued from spectral data treatment is presented on the protonation equilibria of clotrimazole. The experimental spectra are acquired for the titration of an acid4.29×10-5 M clotrimazole solution by a standard solution of 1 M NaOH to adjust pH value. As all variously protonated anions exhibit quite similar absorption bands, a part of the spectrum from 240–320 nm was selected as the most appropriate for an estimation of protonation constants. pH-spectrophotometric titration allows absorbance-response-surface data (Figure 1(a)) to be acquired for analysis with non-linear regression and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. In the first step of the spectra analysis, the number of light-absorbing species was estimated using the INDICES algorithm (Figure 1(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.37 which gives the value sk*(A)=0.43 mAU and which also represents the actual instrumental error sinst(A) of the spectrophotometer used. Very low values of sinst(A) prove that quite reliable spectrophotometer and experimental techniques were used. In the first run, the dissociation constant and two molar absorptivities of clotrimazole by SPECFIT/32 are estimated. The reliability of the regression parameter estimates may be tested using the following diagnostic which have been described previously [37–40]. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 1(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 1(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 1 and Figure 6.Table 1
Dependence of the mixed dissociation constants pKa of clotrimazole on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0250.0400.0480.0550.0860.116pKa4.360(32)4.325(24)4.337(24)4.318(13)4.338(23)4.333(24)4.337(17)s(A) [mAU]0.850.800.510.650.770.950.94Estimated pKa at 37°CSPECFITIonic strength0.0170.0250.0330.0400.0630.0710.116pKa4.152(27)4.072(26)4.037(18)4.073(11)4.030(18)4.066(21)4.037(17)s(A) [mAU]0.740.590.800.570.970.850.84The nonlinear regression analysis of the protonation equilibria model and factor analysis of clotrimazole: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.43 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of clotrimazole in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


#### 3.1.2. Terbinafine HCl

The experimental spectra are acquired for the titration of an acid4.10×10-5 M terbinafine HCl solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 2(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in absorbance spectra are significant within deprotonation, both of the variously protonated species L- and LH exhibit sufficiently different absorption bands. The best region of the spectrum seems to be 240–310 nm and pKa=4.19 (s=0.03).The nonlinear regression analysis of the protonation equilibria model and factor analysis of terbinafine hydrochloride: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.56 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of terbinafine hydrochloride in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)In the first step of the regression spectra analysis, the number of light-absorbing species is estimated by the INDICES algorithm (Figure2(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with corresponding coordinate log (sk*(A))=-3.25 which gives the value sk*(A)=0.56 mAU. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 2(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 2(d) the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 2 and Figure 6.Table 2
Dependence of the mixed dissociation constantspKa of terbinafine hydrochloride on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0170.0320.0400.0480.0550.070pKa4.212(2)4.347(4)4.527(1)4.617(2)4.599(1)4.768(2)4.900(3)s(A) [mAU]0.931.051.000.811.040.551.40Estimated pKa at 37°CSPECFITIonic strength0.0020.0320.0480.0630.078pKa4.136(2)4.276(1)4.282(2)4.470(1)4.474(2)s(A)  [mAU]0.820.651.000.500.98


#### 3.1.3. Acetylsalicylic Acid

The estimation of thermodynamic dissociation constantpKaTof acetylsalicylic acid (ASA) is quite a major problem since during preparation it may contain its main degradation product salicylic acid (SA). So during the determination of ASA,  the presence of SA must be considered. The experimental spectra are acquired for the fast titration of an alkaline 5.96×10-5 M acetylsalicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 3(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. The best region of the spectrum seems to be 220–350 nm and pKa=3.49(25). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 3(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.10 which gives the value sk*(A)=0.79 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 3(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 3(d) shows the distribution diagram of L- and LH species in dependence on pH. The repeatability of the estimated dissociation constant pKa and its mean value at two temperatures 25°C and 37°C at one ionic strength I=0.001 is in Table 3.Table 3
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No.Estimated pKa at 25°Cs(A)  [mAU]No.Estimated pKa at 37°Cs(A)  [mAU]13.514(7)0.8313.440(15)0.8323.481(10)0.9523.413(11)0.7533.447(11)0.9633.400(13)1.1943.510(6)1.0043.422(14)0.9653.505(13)0.9853.403(7)0.8363.498(11)0.8863.405(11)0.95Sample mean3.493(25)3.414(15)The nonlinear regression analysis of the protonation equilibria model and factor analysis of acetylsalicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.79 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of acetylsalicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


#### 3.1.4. Salicylic Acid

The experimental spectra are acquired for the titration of an alkaline4.18×10-5 M salicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 4(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. In cases of small changes in spectra, a precise measurement of absorbance is necessary for a reliable detection of the deprotonation equilibrium studied. The best region of the spectrum seems to be 230–350 nm and pKa=3.01 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 4(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.15 which gives the value sk*(A)=0.71 mAU and which may also be taken here as the actual instrumental error sinst(A) of  the spectrophotometer used. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 4(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 4(d) shows the distribution diagram of L- and LH species on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C on an ionic strength I is in Table 4 and Figure 6. The repeatability of the estimated dissociation constant pKa and its mean value at the temperature 25°C at one ionic strength I=0.001 are in Table 5.Table 4
Dependence of the mixed dissociation constantspKa of salicylic acid on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0010.0170.0470.0540.0620.0690.077pKa2.997(20)2.977(12)2.976(15)2.970(17)2.971(12)2.971(16)2.970(16)s(A) [mAU]1.020.920.800.860.720.830.79Estimated pKa at 37°CSPECFITIonic strength0.0010.0090.0170.0320.0390.0540.0620.077pKa2.995(13)2.939(15)2.964(11)2.917(15)2.935(14)2.929(15)2.930(17)2.917(13)s(A)  [mAU]0.800.620.911.000.770.950.980.98Table 5
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No..Estimated dissociation constants pKa at 25°Cs(A)  [mAU]12.914(13)0.8422.959(14)0.9132.965(15)1.0042.940(15)0.8952.935(18)0.8762.970(16)0.7972.925(21)0.8782.935(12)0.87Mean from the propagation-of-errors2.943(6)Sample means2.943(20)The nonlinear regression analysis of the protonation equilibria model and factor analysis of salicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.71 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


#### 3.1.5. Galanthamine

The experimental spectra are acquired for the titration of an acid4.71×10-5 M galanthamine solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 5(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. The best region of the spectrum seems to be 209–290 nm and pKa=8.21 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 5(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.00 which gives the value sk*(A)=1.00 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 5(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 5(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 6 and Figure 6.Table 6
Dependence of the mixed dissociation constantspKa of galanthamine on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKain the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0700.0960.1360.1360.1630.2550.3220.3480.4150.507pKa8.224(5)8.241(6)8.263(6)8.291(4)8.281(9)8.309(4)8.357(4)8.330(6)8.391(8)8.409(4)s(A) [mAU]0.901.071.120.871.740.930.841.031.430.83Estimated pKa at 37°CSPECFITIonic strength0.0170.0300.0830.1100.1760.2290.2690.362pKa7.989(11)8.027(6)8.049(5)8.015(11)8.109(7)8.091(5)8.147(5)8.183(6)s(A) [mAU]2.391.381.001.391.631.150.981.16The nonlinear regression analysis of the protonation equilibria model and factor analysis of galanthamine: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=1.00 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)Dependence of the mixed dissociation constant pKa of clotrimazole, terbinafine HCl, salicylic acid and galanthamine on the square root of the ionic strength, leading to the thermodynamic dissociation constant pKaT at 25°C, (ADSTAT, ORIGIN).
(a)(b)(c)(d)


### 3.2. Thermodynamic Dissociation Constants

The thermodynamic dissociation constants as the unknown parameterpKaT were estimated by applying a Debye-Hückel equation to the data in Tables 1, 2, 4, 6, and Figure 6 according to the regression criterion; Table 7 shows point estimates of the thermodynamic dissociation constants of the five drugs (clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid, and galanthamine) at two temperatures. Because of the narrow range of ionic strengths the ion-size parameter and the salting-out coefficient C could not be estimated.Table 7
Thermodynamic dissociation constantspKaT for clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid and galanthamine at two temperatures 25°C and 37°C. The standard deviations in the last valid digits are in brackets.
pKaT at 25°CpKaT at 37°CpKaT (Pharma)Clotrimazole4.38(1)4.16(3)7.50± 0.50Terbinafine HCl4.19(3)4.12(5)7.00± 0.50Acetylsalicylic Acid3.49(25)3.41(15)3.50± 0.50Salicylic Acid3.01(1)3.00(1)3.00± 0.50Galanthamine8.21(1)7.99(2)8.60± 0.50


## 3.1. Estimation of Dissociation Constants

### 3.1.1. Clotrimazole

A suggested strategy for effective experimentation in dissociation constants determination ensued from spectral data treatment is presented on the protonation equilibria of clotrimazole. The experimental spectra are acquired for the titration of an acid4.29×10-5 M clotrimazole solution by a standard solution of 1 M NaOH to adjust pH value. As all variously protonated anions exhibit quite similar absorption bands, a part of the spectrum from 240–320 nm was selected as the most appropriate for an estimation of protonation constants. pH-spectrophotometric titration allows absorbance-response-surface data (Figure 1(a)) to be acquired for analysis with non-linear regression and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. In the first step of the spectra analysis, the number of light-absorbing species was estimated using the INDICES algorithm (Figure 1(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.37 which gives the value sk*(A)=0.43 mAU and which also represents the actual instrumental error sinst(A) of the spectrophotometer used. Very low values of sinst(A) prove that quite reliable spectrophotometer and experimental techniques were used. In the first run, the dissociation constant and two molar absorptivities of clotrimazole by SPECFIT/32 are estimated. The reliability of the regression parameter estimates may be tested using the following diagnostic which have been described previously [37–40]. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 1(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 1(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 1 and Figure 6.Table 1
Dependence of the mixed dissociation constants pKa of clotrimazole on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0250.0400.0480.0550.0860.116pKa4.360(32)4.325(24)4.337(24)4.318(13)4.338(23)4.333(24)4.337(17)s(A) [mAU]0.850.800.510.650.770.950.94Estimated pKa at 37°CSPECFITIonic strength0.0170.0250.0330.0400.0630.0710.116pKa4.152(27)4.072(26)4.037(18)4.073(11)4.030(18)4.066(21)4.037(17)s(A) [mAU]0.740.590.800.570.970.850.84The nonlinear regression analysis of the protonation equilibria model and factor analysis of clotrimazole: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.43 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of clotrimazole in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


### 3.1.2. Terbinafine HCl

The experimental spectra are acquired for the titration of an acid4.10×10-5 M terbinafine HCl solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 2(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in absorbance spectra are significant within deprotonation, both of the variously protonated species L- and LH exhibit sufficiently different absorption bands. The best region of the spectrum seems to be 240–310 nm and pKa=4.19 (s=0.03).The nonlinear regression analysis of the protonation equilibria model and factor analysis of terbinafine hydrochloride: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.56 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of terbinafine hydrochloride in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)In the first step of the regression spectra analysis, the number of light-absorbing species is estimated by the INDICES algorithm (Figure2(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with corresponding coordinate log (sk*(A))=-3.25 which gives the value sk*(A)=0.56 mAU. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 2(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 2(d) the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 2 and Figure 6.Table 2
Dependence of the mixed dissociation constantspKa of terbinafine hydrochloride on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0170.0320.0400.0480.0550.070pKa4.212(2)4.347(4)4.527(1)4.617(2)4.599(1)4.768(2)4.900(3)s(A) [mAU]0.931.051.000.811.040.551.40Estimated pKa at 37°CSPECFITIonic strength0.0020.0320.0480.0630.078pKa4.136(2)4.276(1)4.282(2)4.470(1)4.474(2)s(A)  [mAU]0.820.651.000.500.98


### 3.1.3. Acetylsalicylic Acid

The estimation of thermodynamic dissociation constantpKaTof acetylsalicylic acid (ASA) is quite a major problem since during preparation it may contain its main degradation product salicylic acid (SA). So during the determination of ASA,  the presence of SA must be considered. The experimental spectra are acquired for the fast titration of an alkaline 5.96×10-5 M acetylsalicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 3(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. The best region of the spectrum seems to be 220–350 nm and pKa=3.49(25). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 3(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.10 which gives the value sk*(A)=0.79 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 3(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 3(d) shows the distribution diagram of L- and LH species in dependence on pH. The repeatability of the estimated dissociation constant pKa and its mean value at two temperatures 25°C and 37°C at one ionic strength I=0.001 is in Table 3.Table 3
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No.Estimated pKa at 25°Cs(A)  [mAU]No.Estimated pKa at 37°Cs(A)  [mAU]13.514(7)0.8313.440(15)0.8323.481(10)0.9523.413(11)0.7533.447(11)0.9633.400(13)1.1943.510(6)1.0043.422(14)0.9653.505(13)0.9853.403(7)0.8363.498(11)0.8863.405(11)0.95Sample mean3.493(25)3.414(15)The nonlinear regression analysis of the protonation equilibria model and factor analysis of acetylsalicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.79 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of acetylsalicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


### 3.1.4. Salicylic Acid

The experimental spectra are acquired for the titration of an alkaline4.18×10-5 M salicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 4(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. In cases of small changes in spectra, a precise measurement of absorbance is necessary for a reliable detection of the deprotonation equilibrium studied. The best region of the spectrum seems to be 230–350 nm and pKa=3.01 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 4(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.15 which gives the value sk*(A)=0.71 mAU and which may also be taken here as the actual instrumental error sinst(A) of  the spectrophotometer used. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 4(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 4(d) shows the distribution diagram of L- and LH species on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C on an ionic strength I is in Table 4 and Figure 6. The repeatability of the estimated dissociation constant pKa and its mean value at the temperature 25°C at one ionic strength I=0.001 are in Table 5.Table 4
Dependence of the mixed dissociation constantspKa of salicylic acid on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0010.0170.0470.0540.0620.0690.077pKa2.997(20)2.977(12)2.976(15)2.970(17)2.971(12)2.971(16)2.970(16)s(A) [mAU]1.020.920.800.860.720.830.79Estimated pKa at 37°CSPECFITIonic strength0.0010.0090.0170.0320.0390.0540.0620.077pKa2.995(13)2.939(15)2.964(11)2.917(15)2.935(14)2.929(15)2.930(17)2.917(13)s(A)  [mAU]0.800.620.911.000.770.950.980.98Table 5
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No..Estimated dissociation constants pKa at 25°Cs(A)  [mAU]12.914(13)0.8422.959(14)0.9132.965(15)1.0042.940(15)0.8952.935(18)0.8762.970(16)0.7972.925(21)0.8782.935(12)0.87Mean from the propagation-of-errors2.943(6)Sample means2.943(20)The nonlinear regression analysis of the protonation equilibria model and factor analysis of salicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.71 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


### 3.1.5. Galanthamine

The experimental spectra are acquired for the titration of an acid4.71×10-5 M galanthamine solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 5(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. The best region of the spectrum seems to be 209–290 nm and pKa=8.21 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 5(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.00 which gives the value sk*(A)=1.00 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 5(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 5(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 6 and Figure 6.Table 6
Dependence of the mixed dissociation constantspKa of galanthamine on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKain the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0700.0960.1360.1360.1630.2550.3220.3480.4150.507pKa8.224(5)8.241(6)8.263(6)8.291(4)8.281(9)8.309(4)8.357(4)8.330(6)8.391(8)8.409(4)s(A) [mAU]0.901.071.120.871.740.930.841.031.430.83Estimated pKa at 37°CSPECFITIonic strength0.0170.0300.0830.1100.1760.2290.2690.362pKa7.989(11)8.027(6)8.049(5)8.015(11)8.109(7)8.091(5)8.147(5)8.183(6)s(A) [mAU]2.391.381.001.391.631.150.981.16The nonlinear regression analysis of the protonation equilibria model and factor analysis of galanthamine: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=1.00 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)Dependence of the mixed dissociation constant pKa of clotrimazole, terbinafine HCl, salicylic acid and galanthamine on the square root of the ionic strength, leading to the thermodynamic dissociation constant pKaT at 25°C, (ADSTAT, ORIGIN).
(a)(b)(c)(d)


## 3.1.1. Clotrimazole

A suggested strategy for effective experimentation in dissociation constants determination ensued from spectral data treatment is presented on the protonation equilibria of clotrimazole. The experimental spectra are acquired for the titration of an acid4.29×10-5 M clotrimazole solution by a standard solution of 1 M NaOH to adjust pH value. As all variously protonated anions exhibit quite similar absorption bands, a part of the spectrum from 240–320 nm was selected as the most appropriate for an estimation of protonation constants. pH-spectrophotometric titration allows absorbance-response-surface data (Figure 1(a)) to be acquired for analysis with non-linear regression and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. In the first step of the spectra analysis, the number of light-absorbing species was estimated using the INDICES algorithm (Figure 1(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.37 which gives the value sk*(A)=0.43 mAU and which also represents the actual instrumental error sinst(A) of the spectrophotometer used. Very low values of sinst(A) prove that quite reliable spectrophotometer and experimental techniques were used. In the first run, the dissociation constant and two molar absorptivities of clotrimazole by SPECFIT/32 are estimated. The reliability of the regression parameter estimates may be tested using the following diagnostic which have been described previously [37–40]. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 1(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 1(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 1 and Figure 6.Table 1
Dependence of the mixed dissociation constants pKa of clotrimazole on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0250.0400.0480.0550.0860.116pKa4.360(32)4.325(24)4.337(24)4.318(13)4.338(23)4.333(24)4.337(17)s(A) [mAU]0.850.800.510.650.770.950.94Estimated pKa at 37°CSPECFITIonic strength0.0170.0250.0330.0400.0630.0710.116pKa4.152(27)4.072(26)4.037(18)4.073(11)4.030(18)4.066(21)4.037(17)s(A) [mAU]0.740.590.800.570.970.850.84The nonlinear regression analysis of the protonation equilibria model and factor analysis of clotrimazole: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.43 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of clotrimazole in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


## 3.1.2. Terbinafine HCl

The experimental spectra are acquired for the titration of an acid4.10×10-5 M terbinafine HCl solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 2(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in absorbance spectra are significant within deprotonation, both of the variously protonated species L- and LH exhibit sufficiently different absorption bands. The best region of the spectrum seems to be 240–310 nm and pKa=4.19 (s=0.03).The nonlinear regression analysis of the protonation equilibria model and factor analysis of terbinafine hydrochloride: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.56 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of terbinafine hydrochloride in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)In the first step of the regression spectra analysis, the number of light-absorbing species is estimated by the INDICES algorithm (Figure2(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with corresponding coordinate log (sk*(A))=-3.25 which gives the value sk*(A)=0.56 mAU. The goodness-of-fit establishes sufficiently reliable estimates of the dissociation constant and molar absorption coefficient. Figure 2(c) shows the curves of molar absorption coefficients of L- and LH species in dependence on wavelength and Figure 2(d) the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 2 and Figure 6.Table 2
Dependence of the mixed dissociation constantspKa of terbinafine hydrochloride on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0020.0170.0320.0400.0480.0550.070pKa4.212(2)4.347(4)4.527(1)4.617(2)4.599(1)4.768(2)4.900(3)s(A) [mAU]0.931.051.000.811.040.551.40Estimated pKa at 37°CSPECFITIonic strength0.0020.0320.0480.0630.078pKa4.136(2)4.276(1)4.282(2)4.470(1)4.474(2)s(A)  [mAU]0.820.651.000.500.98


## 3.1.3. Acetylsalicylic Acid

The estimation of thermodynamic dissociation constantpKaTof acetylsalicylic acid (ASA) is quite a major problem since during preparation it may contain its main degradation product salicylic acid (SA). So during the determination of ASA,  the presence of SA must be considered. The experimental spectra are acquired for the fast titration of an alkaline 5.96×10-5 M acetylsalicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 3(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. The best region of the spectrum seems to be 220–350 nm and pKa=3.49(25). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 3(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.10 which gives the value sk*(A)=0.79 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 3(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 3(d) shows the distribution diagram of L- and LH species in dependence on pH. The repeatability of the estimated dissociation constant pKa and its mean value at two temperatures 25°C and 37°C at one ionic strength I=0.001 is in Table 3.Table 3
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No.Estimated pKa at 25°Cs(A)  [mAU]No.Estimated pKa at 37°Cs(A)  [mAU]13.514(7)0.8313.440(15)0.8323.481(10)0.9523.413(11)0.7533.447(11)0.9633.400(13)1.1943.510(6)1.0043.422(14)0.9653.505(13)0.9853.403(7)0.8363.498(11)0.8863.405(11)0.95Sample mean3.493(25)3.414(15)The nonlinear regression analysis of the protonation equilibria model and factor analysis of acetylsalicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.79 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of acetylsalicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


## 3.1.4. Salicylic Acid

The experimental spectra are acquired for the titration of an alkaline4.18×10-5 M salicylic acid solution by a standard solution of 1 M HCl (or HClO4) to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 4(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. In cases of small changes in spectra, a precise measurement of absorbance is necessary for a reliable detection of the deprotonation equilibrium studied. The best region of the spectrum seems to be 230–350 nm and pKa=3.01 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 4(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.15 which gives the value sk*(A)=0.71 mAU and which may also be taken here as the actual instrumental error sinst(A) of  the spectrophotometer used. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 4(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 4(d) shows the distribution diagram of L- and LH species on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C on an ionic strength I is in Table 4 and Figure 6. The repeatability of the estimated dissociation constant pKa and its mean value at the temperature 25°C at one ionic strength I=0.001 are in Table 5.Table 4
Dependence of the mixed dissociation constantspKa of salicylic acid on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKa in the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0010.0170.0470.0540.0620.0690.077pKa2.997(20)2.977(12)2.976(15)2.970(17)2.971(12)2.971(16)2.970(16)s(A) [mAU]1.020.920.800.860.720.830.79Estimated pKa at 37°CSPECFITIonic strength0.0010.0090.0170.0320.0390.0540.0620.077pKa2.995(13)2.939(15)2.964(11)2.917(15)2.935(14)2.929(15)2.930(17)2.917(13)s(A)  [mAU]0.800.620.911.000.770.950.980.98Table 5
Repeatability of measurements. The standard deviations of the parameterpKa in the last valid digits are in brackets.
No..Estimated dissociation constants pKa at 25°Cs(A)  [mAU]12.914(13)0.8422.959(14)0.9132.965(15)1.0042.940(15)0.8952.935(18)0.8762.970(16)0.7972.925(21)0.8782.935(12)0.87Mean from the propagation-of-errors2.943(6)Sample means2.943(20)The nonlinear regression analysis of the protonation equilibria model and factor analysis of salicylic acid: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=0.71 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)


## 3.1.5. Galanthamine

The experimental spectra are acquired for the titration of an acid4.71×10-5 M galanthamine solution by a standard solution of 1 M NaOH to adjust pH value. pH spectrophotometric titration allows absorbance-response data (Figure 5(a)) to be acquired for analysis by nonlinear regression, and the reliability of parameter estimates (pKs′ and εs′) can be assessed on the basis of the goodness-of-fit test of residuals. As the changes in spectra are quite small within deprotonation, however, both of the variously protonated species L- and LH exhibit nearly similar absorption bands. The best region of the spectrum seems to be 209–290 nm and pKa=8.21 (s=0.01). Most of the selected methods of factor analysis by the INDICES algorithm lead to two light-absorbing components in the equilibrium mixture (Figure 5(b)). The position of the break point on the sk(A)=f(k) curve in the factor analysis scree plot is calculated and gives k*=2 with the corresponding coordinate log (sk*(A))=-3.00 which gives the value sk*(A)=1.00 mAU. The goodness-of-fit establishes sufficiently reliable estimate of the dissociation constant and molar absorption coefficient. Figure 5(c) shows curves of molar absorption coefficients of L- and LH species in dependence on wavelengths and Figure 5(d) shows the distribution diagram of L- and LH species in dependence on pH. The estimated dissociation constant pKa at two temperatures 25°C and 37°C in dependence on an ionic strength I is in Table 6 and Figure 6.Table 6
Dependence of the mixed dissociation constantspKa of galanthamine on ionic strength using regression analysis of pH-spectrophotometric data with SPECFIT at 25°C and 37°C. The standard deviations of the parameter pKain the last valid digits are in brackets.
Estimated pKa at 25°CSPECFITIonic strength0.0700.0960.1360.1360.1630.2550.3220.3480.4150.507pKa8.224(5)8.241(6)8.263(6)8.291(4)8.281(9)8.309(4)8.357(4)8.330(6)8.391(8)8.409(4)s(A) [mAU]0.901.071.120.871.740.930.841.031.430.83Estimated pKa at 37°CSPECFITIonic strength0.0170.0300.0830.1100.1760.2290.2690.362pKa7.989(11)8.027(6)8.049(5)8.015(11)8.109(7)8.091(5)8.147(5)8.183(6)s(A) [mAU]2.391.381.001.391.631.150.981.16The nonlinear regression analysis of the protonation equilibria model and factor analysis of galanthamine: (a) Absorption spectra in dependence on pH at 25°C, (b) Cattel’s scree plot of the Wernimont-Kankare procedure for determination of the number of light-absorbing species in the mixture k*=2 leads to the actual instrumental error of the spectrophotometer used sinst(A)=1.00 mAU (INDICES in S-Plus), (c) Pure spectra profiles of molar absorptivities versus wavelengths for species L and LH, (d) Distribution diagram of the relative concentrations of species L and LH of salicylic acid in dependence on pH at 25°C. The charges of species are omitted for the sake of simplicity. (SPECFIT, ORIGIN).
(a)(b)(c)(d)Dependence of the mixed dissociation constant pKa of clotrimazole, terbinafine HCl, salicylic acid and galanthamine on the square root of the ionic strength, leading to the thermodynamic dissociation constant pKaT at 25°C, (ADSTAT, ORIGIN).
(a)(b)(c)(d)


## 3.2. Thermodynamic Dissociation Constants

The thermodynamic dissociation constants as the unknown parameterpKaT were estimated by applying a Debye-Hückel equation to the data in Tables 1, 2, 4, 6, and Figure 6 according to the regression criterion; Table 7 shows point estimates of the thermodynamic dissociation constants of the five drugs (clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid, and galanthamine) at two temperatures. Because of the narrow range of ionic strengths the ion-size parameter and the salting-out coefficient C could not be estimated.Table 7
Thermodynamic dissociation constantspKaT for clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid and galanthamine at two temperatures 25°C and 37°C. The standard deviations in the last valid digits are in brackets.
pKaT at 25°CpKaT at 37°CpKaT (Pharma)Clotrimazole4.38(1)4.16(3)7.50± 0.50Terbinafine HCl4.19(3)4.12(5)7.00± 0.50Acetylsalicylic Acid3.49(25)3.41(15)3.50± 0.50Salicylic Acid3.01(1)3.00(1)3.00± 0.50Galanthamine8.21(1)7.99(2)8.60± 0.50


## 4. Conclusions

When drugs are very poorly soluble, pH-spectrophotometric titration may be used with non-linear regression of the absorbance-response-surface data instead of performing a potentiometric determination of the dissociation constants. The reliability of the dissociation constants of the five drugs (clotrimazole, terbinafine HCl, acetylsalicylic acid, salicylic acid, and galanthamine) may be proven with goodness-of-fit tests of the absorption spectra measured at various pH. The dissociation constant pKa was estimated by non-linear regression of {pKa,I} data at 25°C: for clotrimazole pKa,1T=4.38(1), for terbinafine HCl pKa,1T=4.19(3), for acetylsalicylic acid pKa,1T=3.49(25), for salicylic acid pKa,1T=3.01(1), and for galanthamine pKa,1T=8.21(1), where in brackets is the standard deviation in last significant digits. Goodness-of-fit tests for various regression diagnostics enable the reliability of the parameter estimates to be determined.



---

*Source: 527013-2010-04-19.xml*