Acessibilidade / Reportar erro

Effect of Organic Solvent Composition on Dissociation Constants of Some Reversible Acetylcholinesterase Inhibitors

Abstract

Reversible acetylcholinesterase inhibitors are an important group of drug compounds that are used medicinally to treat Alzheimer's disease. In this study, dissociation constant values (pKa) of some reversible acetylcholinesterase inhibitors (donepezil, rivastigmine and galantamine) having different functional groups were determined in different percentages of acetonitrile (MeCN)-water and methanol (MeOH)-water binary mixtures using reversed phase liquid chromatography (RPLC) method. In this way, a complete description of the retention behavior of each solute in the space defined by the pH and organic modifier percentages variables was obtained. The linear relationships established between retention factors of the species and the polarity parameter of the mobile phase (ETN) was proved to predict accurately retention in liquid chromatography (LC) as a function of the acetonitrile and methanol content. In result, the estimated aqueous pKa values of studied basic compounds were assessed by comparing them with literature data.

Keywords:
reversible acetylcholinesterase inhibitors; dissociation constant; nonlinear regression; polarity parameter; RPLC


Introduction

Alzheimer's disease (AD) is a progressive neurological disorder, the most common form of dementia, characterized by memory loss and other intellectual abilities serious enough to interfere with daily life.11 Thompson, P. A.; Wright, D. E.; Counsell, C. E.; Zajicek J.;Int. Psychogeriatr. 2012 , 24, 689.Reversible acetylcholinesterase (AChE) inhibitors are drugs used for treatment of Alzheimer's disease. These inhibitors include compounds with different ionizable functional groups. Donepezil, rivastigmine and galantamine as reversible AChE inhibitors are basic and water-insoluble compounds.

The dissociation constant (pKa) of drugs is of main importance for their absorption, distribution, metabolism, and excretion (ADME). This physicochemical parameter is a main item in the biophysical characterization of a drug and may be helpful in predicting the behavior of a drug under in vivo conditions. Hence, it is important to calculate properly the drug dissociation constant value of a pharmaceutical compound.22 Gupta, V.; Jain, A. D. K.; Gill, N. S.; Gupta, K.; Int. Res. J. Pharm. Appl. Sci.2012 , 2, 17.

Various experimental modalities have been used for the determination of the dissociation constants, comprising potentiometric titration,33 Serjeant, E. P.; Potentiometry and Potentiometric Titrations, 3rd ed.; John Wiley & Sons: New York, 1984. spectrophotometry,44 Meloun, M.; Syrový, T.; Vrána, A.; Talanta2004 , 62, 511. high performance liquid chromatography (HPLC)55 Canbay, H. S.; Demiralay, E. C.; Alsancak, G.; Ozkan, S. A.;J. Chem. Eng. Data2011 , 56, 2071.

6 Yılmaz, H.; Demiralay, E. C.; J. Liq. Chromatogr. Relat. Technol.2015 , 38, 97.
-77 Gumustas, M.; Sanlı, S.; Sanlı, N.; Ozkan, S. A.;Talanta2010 , 82, 1528. and capillary electrophoresis (CE) methodologies.88 Ishihama, Y. M.; Nakamura, T.; Miwa, T.; Kajima, N.; Asakawa, A.;J. Pharm Sci.2002 , 91, 933. Among these techniques, most drugs can be analyzed by reversed phase liquid chromatography (RPLC) technique because of several advantages like rapidity, specificity, accuracy, precision and ease of automation in this method. Moreover, pKa values can also be predicted by computational methods on the basis of molecular structure.99 http://www.chemaxon.com/marvin/sketch/index.jsp accessed in October 2015.
http://www.chemaxon.com/marvin/sketch/in...

Analysis of the change in retention time of the analyte vs. the change in pH of the mobile phase gives an indirect measure of the dissociation constant. The theory for studying the pH dependence of chromatographic retention for ionizable compounds in liquid chromatography was proposed by Horváth et al.1010 Horváth, C.; Melander, W.; Molnár, I.; Anal. Chem.1977 , 49, 142. with the equation

where Ka is defined as the dissociation constant of the protonated base (BH+) in the hydroorganic mixture used as mobile phase, and kBH+ and kBare the limiting retention factors of the protonated and undissociated base, respectively. aHm+ is the hydrogen ion activity in the mobile phase. γBHm+ is the activity coefficient of the dissociated base in the mobile phase that can be calculated by the classical Debye-Huckel equation.1111 Barbosa, J.; Berges, R.; Sanz-Nebot, V.; J. Chromatogr. A1998 , 823, 411.

For any chemical process occurring in solution, the polarity of the solvent plays a crucial role in determining the outcome.1212 Johnson, B. P.; Khaledi, M. G.; Dorsey, J. G.; Anal. Chem.1986 , 58, 2354. A widely employed measure of solvent polarity is Dimroth and Reichardt's ET(30).1313 Reichardt, C.; Solvents and Solvent Effects in Organic Chemistry, 2nd ed.; VCH Verlagsgesellschaft: Weinheim, 1988. ET(30) polarity of organic solvent-water binary mixtures mobile phases used in RPLC were measured and compared with chromatographic retention and selectivity. This polarity index is often used (also in this study) in its normalized, dimensionless form, the so-called empirical parameter of solvent polarity,ETN. Linear plots of log kvs. φ are also satisfactory for most solutes for narrow ranges of organic solvent. However, Johnson et al.1212 Johnson, B. P.; Khaledi, M. G.; Dorsey, J. G.; Anal. Chem.1986 , 58, 2354. demonstrated linear relationships in a wide range of organic solvents between the solute log k andETN.1414 Rosés, M.; Bosch, E.; Anal. Chim. Acta1993 , 274, 147.ETN polarity parameters have been used to predict the chromatographic behavior of ionizable compounds.

Variation of the retention factor of neutral form of a base (kB) and the retention factor of protonated form of a base (kBH+) with percentage of organic modifier in the mobile phase is represented by the normalized Dimroth and Reichardt polarity parameter (equations 2 and3).

Substituting equations 2 and 3 into equation 1 the theoretical expression describing the dependence of the retention factor for basic solutes as a combined function of pH andETN may be expressed as follows:

where CBH+ and CBare intercept values of the protonated and neutral species, respectively. Similarly,eBH+ andeB are slope values of these species.1515 Bergés, R.; Sanz-Nebot, V.; Barbosa, J.; J. Chromatogr. A2000 , 869, 27.

Aqueous dissociation constants are a relevant technological property to know stability and solubility of drug. A critical evaluation of different extrapolation approaches for compounds poorly soluble in water should be very useful in drug development.1616 Ruiz, R.; Ràfols, C.; Rosés, M.; Bosch, E.; J. Pharm. Sci.2003 , 92, 1473. Many drugs are poorly soluble in water and therefore literature proposes several different approaches for their aqueous pKa estimation. In this work, the chromatographic dissociation constant values of donepezil, rivastigmine and galantamine were determined in various MeCN-water and MeOH-water binary mixtures. pKa values of each compounds were measured in five different organic solvent-water mixtures. To obtain the best aqueous pKa value from pKa data, two different extrapolation methods have been tried. In the first approach, pKa values were plotted against MeCN and MeOH mole fraction (X of organic solvents).1717 Rosés, M.; Bosch, E.; J. Chromatogr. A2002 , 1, 982. This method was applied using the following equation:

where the intercepts of these linear equations obtained from this approach were the aqueous pKa values of studied compounds

In second approach, Yasuda-Shedlovsky equation1818 Yasuda, M.; Bull. Chem. Soc. Jpn.1959 , 32, 429.,1919 Shedlovsky, T. In Electrolytes, 3rd ed.; Peace, B., ed.; Pergamon Press: New York, 1962. was used in order to predict aqueous pKa values from the pKa values. This extrapolation method establishes a correlation with the dielectric constant. The following equation has been adopted:

where log[H2O] is the molar water concentration of the given solvent mixture, ε is the dielectric constant of the mixture and a and b are the slope and intercept, respectively. This method is the most widely used procedure in organic solvent-water binary mixtures.

When the pKa value is not directly available, it can often be estimated from the chemical structure of the compound. Marvin Sketch program99 http://www.chemaxon.com/marvin/sketch/index.jsp accessed in October 2015.
http://www.chemaxon.com/marvin/sketch/in...
was used for this aim. For a given structure, it estimates not only pKa values but also other physicochemical parameters.

Experimental

Chemicals and reagents

All chemicals were used directly without any further purification. Donepezil hydrochloride monohydrate and galantamine hydrobromide were purchased from Sigma-Aldrich (St. Louis, USA). Rivastigmine tartrate was kindly provided by Novartis Pharm. Inc. (Istanbul, Turkey). Orthophosphoric acid (min. 85%) was obtained from Riedel-de-Haen (Darmstadt, Germany); other chemicals (acetonitrile, methanol, sodium hydroxide, ammonia, ammonium bicarbonate, potassium hydrogen phthalate) employed were of analytical grade (Merck, Darmstadt, Germany). Double distilled water was used for the preparation of the mobile phase.

Stock solution of the investigated compounds (100 µg mL-1) were prepared in mobile phase, stored at 4 ºC, and then further diluted with mobile phase to the desired concentration. All solutions were protected from light and were used within 24 h to avoid decomposition.

Apparatus

The chromatographic hardware used comprised of a Shimadzu HPLC system (Shimadzu Technologies, Kyoto, Japan) equipped with a pump (LC-20AD), a UV-Vis detector (SPD-20A), a column oven (CTO-20A) and a degasser system (DGU-20A3). pH measurements of the mobile phase were performed with a Mettler Toledo MA 235 pH/ion analyzer (Schwerzenbach, Switzerland) using M-T assembly pH electrode. For the standardization of potentiometric system according to the IUPAC rules,2020 Rondinini, S.; Mussini, P. R.; Mussini, T.; Pure Appl. Chem.1987 , 59, 1549. potassium hydrogen phthalate solution (0.05 mol kg-1) was used.

Chromatographic procedure

In the present work, the effect of solvent composition in mobile phases was analyzed at five solvent levels. At each composition, different pH values were studied, spread over the pH range from 5.0 to 10.5. In this study, mobile phases used were different proportions of MeCN (25 to 55 vol%) and MeOH (35 to 65 vol%), respectively. o-Phosphoric acid and ammonium bicarbonate were used as buffer components because of their appropriate pKa values and the symmetrical peak shape of studied compounds in these buffer solutions. The concentration of the buffer systems was in all cases 30 mmol L-1 (concentration after mixing the organic solvent-water binary mixtures). Chromatographic pKa determination was carried out using a X Terra C18 column (250 mm × 4.6 mm i.d., 5 µm, Waters, Milford, Massachusetts, USA). Chromatographic measurements were done at 25 ºC with an eluent flow rate of 1 mL min-1. The volume of solution injected into the column was 20 µL for each run. The compounds studied had different optimal wavelengths (for donepezil and galantamine: 230 nm; for rivastigmine: 225 nm; for uracil: 254 nm).

Results and Discussion

In the present study, the retention factors for the investigated compounds were obtained at different percentages of acetonitrile/methanol in the mobile phase and at different mobile phase pH values. Selected organic modifier concentration in MeCN-water binary mixtures is much smaller than that in MeOH-water binary mixtures, because MeCN is a stronger solvent than MeOH. Therefore, concentrations of these two solvents are not the same.

The pKa values, as well as the retention factors of the studied compounds, were obtained by fitting the retention factors-pH data to equation 1, usually by non-linear regression. Nowadays, there are many commercial computer programs to perform these fittings, including ORIGIN Pro 8 software. The pKavalues and the retention factors of investigated compounds in different MeCN and MeOH ratios calculated by NLREG2121 Sherrod, P.H.; Non-Linear Regression Analysis Program, NLREG Version 5.0; 2000. Available at: http://www.nlreg.com/index.htm, accessed in October 2015.
http://www.nlreg.com/index.htm...
and ORIGIN2222 http://www.originlab.com/ accessed in October 2015.
http://www.originlab.com/...
programs are given in Tables 1 and 2, together with respective standard deviations. As expected, dissociation constant values decrease with the increase of MeCN and MeOH contents, according to the general behavior of weak basic compounds. The reason is that MeCN is less basic than methanol. Hence, bases are stronger in these solvents.

Table 1
Retention factors of undissociated (kB) and protonated (kBH+) bases and the dissociation constant values of investigated compounds in various acetonitrile/water ratios calculated by NLREG and ORIGIN programs
Table 2
Retention factors of undissociated (kB) and protonated (kBH+) bases and the dissociation constant values of investigated compounds in various methanol/water ratios calculated by NLREG and ORIGIN programs

From the data given in Tables 1 and 2, it is immediately obvious that the nature of solvent plays a fundamental role in the base equilibria. As shown in Tables 1 and 2, all drugs studied have one pKa value corresponding to base functional group. pKa values obtained for donepezil are consistent with the dissociation equilibrium ofN-benzylpiperidine ring (Figure 1). The pKa of rivastigmine is similar to the value of donepezil and can be assigned to the dissociation of the secondary amine group. Galantamine has an azepine moiety and pKa value of this compound can be associated with nitrogen atom in azepine (Figure 1).

Figure 1
Dissociation behavior of (a) rivastigmine; (b) donepezil and (c) galantamine.

The retention factors were obtained over a pH range of 5.0-10.5, in order to determine pKa of these acetylcholinesterase inhibitors using RPLC method. In Figures 2 and 3, data pairs of k and pH for investigated compounds in different percentages of MeCN-water and MeOH-water binary mixtures are shown, together with the corresponding experimental and calculated retention factors. It can be concluded that plots of sigmoidal curves of retention factor (k) vs. pH of the mobile phase are related to the influence of an organic modifier on the dissociation of basic solute.

Figure 2
Plots of the calculated retention factors vs. the pH of the mobile phase for (a) galantamine; (b) rivastigmine and (c) donepezil in several MeCN-water mobile phases.
Figure 3
Plots of the calculated retention factors vs. the pH of the mobile phase for (a) galantamine; (b) rivastigmine and (c) donepezil in several MeOH-water mobile phases.

Variation of the mobile phase pH is a key parameter to enhance the chromatographic selectivity and retention time for ionization compounds. Figure 4 shows the typical chromatograms of drugs eluted by different percentages of MeCN/MeOH in the mobile phase at different pH environment. It provided that the retention times (tR) of donepezil, rivastigmine and galantamine increased as the pH value increased.

Figure 4
Typical HPLC chromatograms of studied compounds eluted by the different percentages of mobile phase at different pH environment. (a) Donepezil (55 vol% MeOH); (b) rivastigmine (55 vol% MeOH); (c) galantamine (35 vol% MeOH); (d) donepezil (45 vol% MeCN); (e) rivastigmine (45 vol% MeCN); and (f) galantamine (25 vol% MeCN).

In this work, equation 4 was evaluated statistically for these drugs at different conditions. In this way, a complete description of the retention behavior of each solute in the space defined by the pH of the mobile phase and ETN variables was obtained. The coefficients describing a retention behavior of these solutes usingequation 4 are listed in Table 3. The retention behavior of these compounds was accurately modeled as a function of pH and polarity of the mobile phase with the equations proposed. The results shown in Table 3 demonstrate a good performance of the liquid chromatographic method for determination of dissociation constants.

Table 3
Results of the application of the retention model proposed for all the experimental retention data available for donepezil, rivastigmine and galantamine (coefficients describe chromatographic retention in terms of equation 4 for basic compounds)

In the first described approach, equation 5, was used to estimate the aqueous pKavalues of studied reversible AChE inhibitors. There is actually a linear relationship between pKa of the rivastigmine, donepezil, galantamine used and mole fraction of methanol and acetonitrile in the binary mixtures. The results of linear equations for compounds used are provided in Table 4. All these straight lines show reasonably satisfactory correlation coefficients and standard deviations.

Table 4
Linear equations between experimental dissociation constant values and the mole fraction of acetonitrile/methanol in the binary mixtures

Then, the following expression, equation 6, was used to estimate the aqueous pKa of investigated substances. The parameters of the Yasuda-Shedlovsky equations are also summarized in Table 5. It can be seen that basic functional groups have negative slopes. Different functional groups show remarkably different slopes. It is interesting to compare the slopes of a given molecule in MeCN-water and in MeOH-water binary mixtures. The linearity of the plots is characterized by the regression coefficients (r) values which indicate significant linear correlation for the molecules examined. The average of the r values is 0.999.

Table 5
Aqueous dissociation constant values by Yasuda-Shedlovsky extrapolation in different co-solvents

On the other hand, the Marvin Sketch program99 http://www.chemaxon.com/marvin/sketch/index.jsp accessed in October 2015.
http://www.chemaxon.com/marvin/sketch/in...
was used to compute the drugs' aqueous pKa. The pKa values reported in the literature and predicted by this program99 http://www.chemaxon.com/marvin/sketch/index.jsp accessed in October 2015.
http://www.chemaxon.com/marvin/sketch/in...
are presented in Table 6. In literature reports, different pKa values for the same compound are calculated. In the present study, calculated pKa values for galantamine are not in full consistency with that reported by Meloun et al.2323 Meloun, M.; Bordovská, S.; Galla, L.; SRX Pharmacol. 2010 , 2010, 1. In that study,2323 Meloun, M.; Bordovská, S.; Galla, L.; SRX Pharmacol. 2010 , 2010, 1. galantamine was determined in different ionic strengths and temperatures. For rivastigmine, pKavalues were predicted by Hsieh et al.,2424 Hsieh, Y. H.; Yang, Y. H.; Yeh, H. H.; Lin, P. C.; Chen, S. H.;Electrophoresis2009 , 30, 644. and Luan et al.2525 Luan, F.; Ma, W.; Zhang, H.; Zhang, X.; Liu, M.; Hu, Z.; Fan, B.;Pharm. Res.2005 , 22, 1454. In these studies, there are no experimental values. There is only one study in the literature for donepezil.2626 Ishihama, Y.; Nakamuro, M.; Miwa, T.; Kajima, T.; Asakawa, N.;J. Pharm. Sci.2002 , 91, 933. This study was carried out by capillary electrophoresis method. Therefore, calculated pKa values for donepezil are not in full consistency with that reported by Ishihama et al.2626 Ishihama, Y.; Nakamuro, M.; Miwa, T.; Kajima, T.; Asakawa, N.;J. Pharm. Sci.2002 , 91, 933. Table 6 shows that despite the use of different equipment and determination by different persons on different days using different starting drug concentrations, close values were obtained for dissociation constants of the determinations. A comparison with representative dissociation constants available in the literature shows that either method gives satisfactory values.

Table 6
Aqueous pKa values of some acetylcholinesterase inhibitors obtained from different methodologies

Conclusions

This work represents the first study dealing with the chromatographic determination of pKa values of some reversible acetylcholinesterase inhibitors at different proportions of MeCN-water and MeOH-water binary mixtures. The combined effect of the two factors (solvent percentage and pH of the mobile phase) on the retention behavior of donepezil, rivastigmine, and galantamine was investigated. Dissociation constant values and limiting retention factors of these compounds were calculated chromatographically using derived equations. Mole fraction and Yasuda-Shedlovsky extrapolation methods were proposed to obtain the aqueous pKa values. The results obtained agree very well with those obtained by calculation by means of the Marvin Sketch program, which is based on the analysis of the chemical structure relative to the dissociation of the selected basic group. The important data extracted from this work can be used for pharmacokinetic and pharmacological studies of these drugs. Moreover, the knowledge of the chromatographic behavior of the drug compounds in different water-organic solvent media is a useful parameter to optimize analytical procedures for the separation of ionizable compounds by liquid chromatography.

Acknowledgements

This work was supported by BAP (Project 3301-YL2-12) and TUBİTAK (Project 113Z025). The authors greatly acknowledge PhD Jose Luis Beltran from University of Barcelona for kindly providing NLREG program.

References

  • 1
    Thompson, P. A.; Wright, D. E.; Counsell, C. E.; Zajicek J.;Int. Psychogeriatr 2012 , 24, 689.
  • 2
    Gupta, V.; Jain, A. D. K.; Gill, N. S.; Gupta, K.; Int. Res. J. Pharm. Appl. Sci.2012 , 2, 17.
  • 3
    Serjeant, E. P.; Potentiometry and Potentiometric Titrations, 3rd ed.; John Wiley & Sons: New York, 1984.
  • 4
    Meloun, M.; Syrový, T.; Vrána, A.; Talanta2004 , 62, 511.
  • 5
    Canbay, H. S.; Demiralay, E. C.; Alsancak, G.; Ozkan, S. A.;J. Chem. Eng. Data2011 , 56, 2071.
  • 6
    Yılmaz, H.; Demiralay, E. C.; J. Liq. Chromatogr. Relat. Technol.2015 , 38, 97.
  • 7
    Gumustas, M.; Sanlı, S.; Sanlı, N.; Ozkan, S. A.;Talanta2010 , 82, 1528.
  • 8
    Ishihama, Y. M.; Nakamura, T.; Miwa, T.; Kajima, N.; Asakawa, A.;J. Pharm Sci.2002 , 91, 933.
  • 9
    http://www.chemaxon.com/marvin/sketch/index.jsp accessed in October 2015.
    » http://www.chemaxon.com/marvin/sketch/index.jsp
  • 10
    Horváth, C.; Melander, W.; Molnár, I.; Anal. Chem.1977 , 49, 142.
  • 11
    Barbosa, J.; Berges, R.; Sanz-Nebot, V.; J. Chromatogr. A1998 , 823, 411.
  • 12
    Johnson, B. P.; Khaledi, M. G.; Dorsey, J. G.; Anal. Chem.1986 , 58, 2354.
  • 13
    Reichardt, C.; Solvents and Solvent Effects in Organic Chemistry, 2nd ed.; VCH Verlagsgesellschaft: Weinheim, 1988.
  • 14
    Rosés, M.; Bosch, E.; Anal. Chim. Acta1993 , 274, 147.
  • 15
    Bergés, R.; Sanz-Nebot, V.; Barbosa, J.; J. Chromatogr. A2000 , 869, 27.
  • 16
    Ruiz, R.; Ràfols, C.; Rosés, M.; Bosch, E.; J. Pharm. Sci.2003 , 92, 1473.
  • 17
    Rosés, M.; Bosch, E.; J. Chromatogr. A2002 , 1, 982.
  • 18
    Yasuda, M.; Bull. Chem. Soc. Jpn.1959 , 32, 429.
  • 19
    Shedlovsky, T. In Electrolytes, 3rd ed.; Peace, B., ed.; Pergamon Press: New York, 1962.
  • 20
    Rondinini, S.; Mussini, P. R.; Mussini, T.; Pure Appl. Chem.1987 , 59, 1549.
  • 21
    Sherrod, P.H.; Non-Linear Regression Analysis Program, NLREG Version 5.0; 2000. Available at: http://www.nlreg.com/index.htm, accessed in October 2015.
    » http://www.nlreg.com/index.htm
  • 22
    http://www.originlab.com/ accessed in October 2015.
    » http://www.originlab.com/
  • 23
    Meloun, M.; Bordovská, S.; Galla, L.; SRX Pharmacol 2010 , 2010, 1.
  • 24
    Hsieh, Y. H.; Yang, Y. H.; Yeh, H. H.; Lin, P. C.; Chen, S. H.;Electrophoresis2009 , 30, 644.
  • 25
    Luan, F.; Ma, W.; Zhang, H.; Zhang, X.; Liu, M.; Hu, Z.; Fan, B.;Pharm. Res.2005 , 22, 1454.
  • 26
    Ishihama, Y.; Nakamuro, M.; Miwa, T.; Kajima, T.; Asakawa, N.;J. Pharm. Sci.2002 , 91, 933.

Publication Dates

  • Publication in this collection
    Mar 2016

History

  • Received
    29 July 2015
  • Accepted
    22 Oct 2015
Sociedade Brasileira de Química Instituto de Química - UNICAMP, Caixa Postal 6154, 13083-970 Campinas SP - Brazil, Tel./FAX.: +55 19 3521-3151 - São Paulo - SP - Brazil
E-mail: office@jbcs.sbq.org.br