Interactions of chitosan/genipin hydrogels during drug delivery: a QSPR approach

ARTIGO

Interactions of chitosan/genipin hydrogels during drug delivery: a QSPR approach

Nancy L. Delgadillo-Armendariz^I; Norma A. Rangel-VazquezI,^*; Edgar A. Marquez-Brazon^II,III; Blanca Rojas-De Gascue^II

^IDivisión de Estudios de Posgrado e Investigación, Instituto Tecnológico de Aguascalientes, Ave. López Mateos #1801 Ote. Col. Fracc. Bona Gens CP 20256, Aguascalientes, Ags. México

^IIUniversidad de Oriente, Instituto de Investigaciones en Biomedicina y Ciencias Aplicadas, “Dra. Susan Tai”, IIBCA-UDO, Apdo. Postal 245, Cumaná, Estado Sucre, Venezuela

^IIIUniversidad de Oriente, Departamento de Química, Núcleo de Sucre, Apdo. Postal 245, Cumaná, Estado Sucre, Venezuela

ABSTRACT

A hydrogel comprised of chitosan crosslinked using the low-toxicity crosslinker genipin was prepared, and the absorption of glibenclamide by the hydrogel was investigated. Optimized structures and their molecular electrostatic potentials were calculated using the AM1 method, and the results were used to evaluate the molecular interactions between the three compounds. The quantitative structure-property relationship model was also used to estimate the activity of the chemicals on the basis their molecular structures. In addition, theoretical Fourier transform infrared spectra were calculated to analyze the intermolecular interactions in the proposed system. Finally, the hydrophilicity of the hydrogel and its influence on the absorption process were also estimated.

Keywords: chitosan; genipin; glibenclamide; FTIR; QSPR modeling.

INTRODUCTION

The use of natural polymers (e.g., gelatin, chitosan, and silk proteins) for the development of hydrogels is gaining in importance owing to their inherent biocompatibility. ¹ Hydrogels are high-water content materials prepared from crosslinked polymers that can provide sustained, local delivery of a variety of therapeutic agents. Use of the natural polymer chitosan as a scaffold material for hydrogels has been intensively investigated because of the polymer's biocompatibility, low toxicity, and biodegradability. Therefore, the advanced development of chitosan hydrogels has led to new drug delivery systems that release their payloads in response to varying environmental stimuli. ^2,3

Chitosan (see Figure 1), the product of the N-deacetylation of chitin, has also received particular attention because it is produced as a waste product during crustacean (shrimp and crab) processing. ^4-8 In addition, it offers other advantages; chitosan can not only be used to control the release of active agents but also be prepared without the use of hazardous organic solvents because it is soluble in aqueous acidic solutions. Given the above-mentioned properties, chitosan is extensively used for the development of drug delivery systems.

Chitosan has been effectively used in drug delivery applications as hydrogel systems, drug conjugates, biodegradable release systems, and other components. ^9,10 Chitosan-based systems are used for the delivery of proteins/peptides, growth factors, anti-inflammatory drugs, and antibiotics, as well as, in gene therapy and bioimaging applications. ^9,11,12

The crosslinking of chitosan is considered a relevant approach for improving control of drug delivery. Among other reagents, glutaraldehyde and tripolyphosphate have been widely used to crosslink this biodegradable polymer. However, there are concerns over the toxicity of most of the investigated crosslinking agents, especially glutaraldehyde, which may reduce the biocompatibility of chitosan delivery systems. Therefore, there is a need for the development of crosslinking agents with low cytotoxicity. ^1,13

Genipin, which is a natural, water-soluble, and bifunctional crosslinker ¹⁴ (Figure 2) obtained from the fruits of Gardenia jasminoides¹, was selected for the study as a possible biocompatible crosslinker for chitosan. Genipin is an excellent natural crosslinker for proteins, collagen, gelatin, and chitosan. It has a low acute toxicity with an in vivo LD₅₀ of 382 mg/kg in mice and is 10,000 times less toxic than glutaraldehyde and degrades more slowly than glutaraldehyde and many other commonly used synthetic crosslinking reagents. Furthermore, genipin is also used as a regulating agent for drug delivery, as the raw material for gardenia blue pigment preparation, and as an intermediate for alkaloid synthesis. ^1,15-21

Glibenclamide (see Figure 3), also known as glyburide, is an antidiabetic drug in the class of medications known as sulfonylureas, which are closely related to sulfa drugs. Sulfonylureas bind to ATP-dependent K⁺ channels in beta cells of the pancreas, depolarizing them and stimulating the release of Ca²⁺, which in turn stimulates insulin production. ^22,23

In this study, Fourier transform infrared (FTIR) spectroscopy was used to analyze the intermolecular interactions in the proposed system. As one of the ways to verify the obtained results, we apply the theoretical calculations by means of the computational chemistry tools. ²⁴ Using the quantitative structure-property relationship (QSPR) model/method, we estimate the activity of the chemicals on the basis of only their molecular structure. Sorption of organic chemicals in soils or sediments is usually described by sorption coefficients. ²⁵

In addition, the log P value of the hydrogel was determined as a measure of its hydrophobicity. ²⁶ The log P value of a compound, which is the logarithm of its partition coefficient between n-octanol and water (c_octanol/c_water), is a well-established measure of the compound's hydrophilicity. Low hydrophilicities, and therefore high log P values, cause poor absorption or permeation. It has been shown that compounds must have a log P value less than 5 in order to have a reasonable probability of being well absorbed. ²⁷ Typically, the log P value of a substance at pH 7.4 is considered as an index of the compound behavior in plasma. ²⁸

Finally, the molecular electrostatic potential (MESP) was investigated using the Austin Method 1 (AM1) calculations. This method provides information about the regions in which intermolecular interactions occur between compounds. ²⁹ The electrostatic potential is the energy of the interaction of a point positive charge (an electrophile) with the nuclei and electrons of a molecule. Negative electrostatic potentials indicate areas that are prone to electrophilic attack. The electrostatic potential can be mapped onto the electron density by using color to represent the value of the potential. ³⁰

EXPERIMENTAL

The AM1 method was initially used to optimize the structures of the compounds investigated in the present study because it generates lower-energy structures, even when the initial structures are far away from the minimum structures. The Polak-Ribiere algorithm was used for mapping the energy barriers of the conformational transitions. For each structure, 5715 iterations, a level convergence of 0.001 kcal/mol/Å, and a line search of 0.1 were performed.

Structural parameters

The optimized structural parameters were used for the vibrational wavenumber calculations with AM1 method in order to characterize all the stationary points as minima. The structural parameters were calculated by selecting the Constrain bond and length options from the Build menu for the two methods of analysis.

FTIR

An infrared spectrum is commonly obtained by passing infrared electromagnetic radiation through a sample that possesses a permanent or induced dipole moment and determining what fractions of the incident radiation are absorbed at particular energies. The energy of each peak in an absorption spectrum corresponds to the frequency of the vibration of part of the molecule, thus allowing qualitative identification of certain bond types in the sample.

The FTIR spectra were obtained by first selecting from the Compute menu, vibrational and rotational spectrum options; after this step, using the vibrational spectrum option, FTIR spectrum pattern is obtained for the two methods of analysis. The results of the analyses for the optimized structures for chitosan, genipin, and glibenclamide obtained using the AM1 method are listed in Tables 1–3, respectively.

Electrostatic potential

After obtaining the Gibbs free energy for the optimized geometries using the AM1 method, two-dimensional contour diagrams of the electrostatic potentials surrounding the three molecules, their total electron densities and spin densities, their molecular orbitals, and the electron densities of the individual orbitals were plotted.

HyperChem software displays the electrostatic potential as a contour plot when the appropriate option in the Contour Plot dialog box is selected. Atomic charges indicate where large negative values (sites for electrophilic attack) are likely to occur. However, the largest negative value for the electrostatic potential is not necessarily adjacent to the atom with the largest negative charge. ³¹

RESULTS AND DISCUSSION

Structural parameters

The thermodynamic data obtained are listed in Table 4. The Gibbs free energy value (–947.998 kcal/mol) indicates that the crosslinking reaction involves crosslinking of two free amino groups in chitosan with one molecule of genipin (see Figure 4). ^17,18 Nucleophilic attack by the primary amine occurs at the carbon atom in genipin, while the secondary amine attacks the aldehyde group. ²¹

QSPR methodology is one of the most powerful tools for describing the relationship between the biological activity and physicochemical characteristics of a molecule. Diffusion is also influenced by the nature of the crosslinker. In addition, the porosity, water uptake, and swelling of hydrogels are influenced by the chemical nature and length of the bridging units of the crosslinker.

Swelling and diffusion in hydrogels are also influenced by the hydrophilicity of the crosslinked polymer. ³² Therefore, the Log P value for chitosan/genipin was calculated and found to be –22.15. This value represents the hydrophilic property of the substance and is considered to be an indicator of the ability of the hydrogel to absorb molecules. ³³ In hydrogels formed by chitosan crosslinked with itself, release is largely controlled by the crosslinking density; consequently, the higher the crosslinking density, the lower the release rate.

However, other system parameters, such as drug concentration, often play a major role. To our knowledge, there are no examples of hydrogels formed by chitosan crosslinked with itself that exhibit pH-sensitive swelling.

Indeed, the numerous inter-chain interactions formed by crosslinking inhibit swelling because most of the amino groups of chitosan react with the crosslinker. Such systems do not present a release profile that can be further modulated after administration; for example, drug release cannot be targeted in the gastro-intestinal tract, which limits their range of application. ^5,32

Table 4 shows the thermodynamic results of glibenclamide absorption process; here the Gibbs free energy value –243.08 kcal/mol is spontaneous because of the formation of hydrogen bonds between the carbonyl groups of glibenclamide and NH groups of chitosan. Nucleophilic attack of the amino group of chitosan at the carbonyl of the genipin leads to formation of a stable bond (amide), and an oxygen atom in genipin is replaced by a nitrogen atom from chitosan (see Figure 5). ³⁴

Figure 5

The optimized structural parameters for chitosan, chitosan crosslinked with genipin, and the chitosan/genipin hydrogel with absorbed glibenclamide are listed in Tables 5–7, respectively, in accordance with the atom numbering scheme given in Figures 4 and 5. The title system (chitosan-genipin/glibenclamide) is not planar, and based on observations, deformation of the ring structure depends on the nature of the substituents (OH, NH₂, and CH₂-CH₂). Therefore, we could compare the calculation results given in Tables 5–7 with experimental data. As discussed in the previous literature, several authors have explained the changes in frequency or bond length of the C–H bond upon substitution due to a change in the charge distribution on the carbon atom of the ring. ³¹ The substituents may either be electron withdrawing (Cl, Br, F) or electron donating (CH₃, C₂H₅).

The carbon atoms are bonded to hydrogen atoms via σ bonds, and substitution of a hydrogen atom for a halogen atom on the benzene ring reduces the electronic density at the carbon atom owing to an induction effect. The ring carbon atoms exhibit a larger attraction for the valence electron cloud of the remaining hydrogen atoms, resulting in an increase in the C–H force constant and a decrease in the corresponding bond lengths. The reverse holds true for substitution with electron donating groups.

The actual change in the C–H bond lengths is influenced by the combined effects of the inductive–mesmeric interaction and the electric dipole field of the polar substituent. The calculated geometric parameters can be used as the foundation for calculating other parameters for the compound. ³⁵

FTIR

The calculated FTIR results for chitosan/genipin and chitosan (crosslinked with genipin)/glibenclamide are presented in Table 8. The second column indicates the absorption bands at 2320 and 1097 cm^–1, which clearly appeared after crosslinking with genipin. The peak at 1097 cm^–1 represents the C–N stretching vibration for the tertiary nitrogen atom formed from the reaction of the lysine with genipin, while the band at 2320 cm^–1 was assigned to C–C, C–H, C–O, and C–N bonds that formed following crosslinking. Moreover, the adsorption band that appeared at 1236 cm^–1 (amide III) represents a mixed CO–N/N–H vibration mode, and the peak at 813 cm^–1 is a characteristic absorption band for the C–H stretching vibration associated with the heterocyclic ring of the product.

An amine group in the chitosan macromolecule undergoes nucleophilic attack at the C–OH group of genipin, resulting in the formation of a new covalent bond between the aldehyde group and the secondary amine and formation of a double bond at the carbon in the ortho-position. ^16,17,36-38 One of the most important changes is evident in the reduction of the carbonyl group of the genipin, which reacts with the primary amine in chitosan to form an amide, which is indicated by the growth of the absorption band at 1630 cm^–1 in the FTIR spectrum of the hydrogel.

As a consequence of these two main reactions, the intensity of the absorption band for the amino groups at 3597 cm^–1 is drastically reduced in the spectrum of the hydrogel. ^36-40

The third column in Table 8 lists the absorption bands that appeared after the absorption of glibenclamide by the chitosan/genipin crosslinked hydrogel. These peaks include a characteristic amide band at 3407 cm^–1 and peaks for SO and SO₂ stretching vibrations at 2683 and 1369 cm^–1, respectively. The site of interaction on glibenclamide is likely the C=O group, which would also affect the NH vibration. ³⁸

The absorption peaks observed at 2412 and 793 cm^–1 are associated with C–C, CH, C–O, and OH bending vibrations, while the absorption bands at 1955, 1334, 1090, and 924 cm^–1 are due to the C–C and C–O bonds in the chitosan crosslinked with genipin. ⁴¹ The new broad peak that appeared near 1415 cm^–1 after crosslinking with genipin is associated with the ring-stretching vibrations of the heterocyclic amine. In addition, the C–N stretching of amide III at 1233 cm^–1 shifted to 1260 cm^–1 after crosslinking with genipin, while the peak at 1740 cm^–1, which corresponds to the carboxylic groups in the electrospun chitosan fibers, disappeared in the genipin-crosslinked chitosan. ³⁹

Electrostatic potential

The MESP values (Figure 6) ranged from –0.075 to 0.607 eV and 0.065 to 0.718 eV for the crosslinking of chitosan with genipin and glibenclamide absorption by the crosslinked hydrogel, respectively. The negative regions appeared near the OH groups (C–OH bonds) in the crosslinked chitosan.

The absorption of glibenclamide by the chitosan crosslinked with genipin mainly involved the formation of hydrogen bonds between the CH₂OH and CH groups. Figure 7 shows the structure of the genipin-crosslinked chitosan with absorbed glibenclamide.

CONCLUSIONS

Genipin, a natural crosslinking agent, reacts with compounds containing primary amine groups, such as chitosan, to form covalently crosslinked networks. In the case of chitosan, genipin reacts with the free amino groups present in the glucosamine units. The results of this study support the relevance of genipin as a valuable natural, nontoxic, crosslinking agent for controlled drug release in drug delivery systems based on chitosan.

The negative log P values for the crosslinked chitosan indicated that absorption properties of the hydrogel are dominated by its hydrophilic character and that absorption plays an important role in both glibenclamide absorption partitioning and receptor-binding. It was also shown that genipin changes the hydrophilicity of chitosan. FTIR results revealed that secondary amide linkages are formed via the reaction of genipin ester groups with chitosan amino groups, yielding a polymeric network structure. The MESP values indicated that the nucleophilic and electrophilic regions mainly involved the NH/C–OH and C=O bonds, respectively.

Recebido em 10/04/2014

aceito em 27/06/2014

publicado na web em 05/09/2014

• 1. Thakur, G.; Mitra, A.; Rousseau, D.; Basak, A.; Sarkar, S.; Pal, K.; J. Mater. Sci. : Mater. Med 2011, 22, 115.
• 2. Bhattarai, N.; Gunn, J.; Zhang, M.; Adv. Drug Delivery Rev. 2010, 62, 83.
• 3. Muzzarelli, R. A.; Carbohydr. Polym 2009, 77, 1.
• 4. Wu, C. S.; Polymer 2005, 46, 147
• 5. Berger, J.; Reist, M.; Mayer, J. M.; Felt, O.; Peppas, N. A.; Gurny, R.; Eur. J. Pharm. Biopharm 2004, 57, 19.
• 6. Diaconu, M.; Tache, A.; Victoria, S. A.; Gatea, F.; Litescu, S.; Radu, G.; UPB Scientific Bulletin, Series B: Chemistry and Materials Science 2010, 72, 115.
• 7. Patel, M.; Ravi, R.; Patel, J.; Patel, K.; J. Pharm. Pharm. Sci. 2010, 13, 536.
• 8. Sailaja, A. K.; Amareshwar, P.; Chakravarty, P.; J Pharm. Biol. Chem. Sci 2010, 1, 474.
• 9. Dash, M.; Chiellini, F.; Ottenbrite, R. M.; Chiellini, E.; Prog. Polym. Sci. 2011, 36, 981.
• 10. Grolik, M.; Szczubiałka, K.; Wowra, B.; Dobrowolski, D.; Orzechowska-Wylęgała, B. : Wylęgała, E.; Nowakowska, M.; J. Mater. Sci. : Mater. Med. 2012, 23, 1991.
• 11. Bansal, V.; Sharma, P. K.; Sharma, N.; Pal, O. P.; Malviya, R.; Adv. Biol. Res. 2011, 5, 28.
• 12. Bhattarai, N.; Gunn, J.; Zhang, M.; Adv Drug Deliver Rev 2010, 62, 83
• 13. Harris, R.; Lecumberri, E.; Heras, A.; Mar. Drugs 2010, 8, 1750.
• 14. Karnchanajindanun, J.; Srisa-Ard, M.; Srihanam, P.; Baimark, Y.; Nat. Sci. 2010, 2, 1061.
• 15. Zhang, C. Y.; Parton, L. E.; Ye, C. P.; Krauss, S.; Shen, R.; Lin, C. T.; Porco, J. A.; Cell Metab 2006, 3, 417.
• 16. Park, J. E.; Lee, Y.; Kim, H. G.; Hahn, T. R.; Paik, Y. S.; J. Agric. Food Chem 2002, 50, 6511.
• 17. Wang, L.; Wang, Y.; Qu, J.; Hu, Y.; You, R.; Li, M.; J. Biomater. Nanobiotechnol. 2013, 4, 213.
• 18. Imsombut, T.; Srisuwan, Y.; Srihanam, P.; Baimark, Y.; Powder Technol 2010, 203, 603.
• 19. Silva, S. S.; Maniglio, D.; Motta, A.; Mano, J. F.; Reis, R. L.; Migliaresi.; Macromol. Biosci 2008, 8, 1.
• 20. Taylor, M.; Ding, K.; Brown, E.; J. Am. Leather Chem. Assoc. 2009, 1.
• 21. Yoo, J. S.; Kim, Y. J.; Kim, S. H.; Choi, S. H.; The Korean Journal of Thoracic and Cardiovascular Surgery 2011, 44, 197.
• 22. Serrano-Martín, X.; Payares, G.; Mendoza-León, A.; Antimicrob. Agents Chemother. 2006, 50, 4214.
• ²³
  http://www.scbt.com/es/datasheet-200982-glyburide-glibenclamide.html, accessed July, 2014.
  » link
• 24. Mon, J.; Flury, M.; Harsh, J. B.; J. Hydrol. 2006, 316, 84.
• 25. Zhang, X. H.; Zhan, Y. H.; Chen, D.; Wang, F.; Wang, L. Y.; Dyes Pigm 2012, 93, 1408.
• ²⁶
  http://www.molinspiration.com/services/logp.html, accessed July, 2014.
  » link
• ²⁷
  http://www.organic-chemistry.org/prog/peo/cLogP.html, accessed July, 2014.
  » link
• ²⁸
  http://www.admescope.com/learn-adme/physicochemical-properties/lipophilicity-logd-logp.html, accessed July, 2014.
  » link
• 29. Glish, L.; Hanks, T. W.; J. Chem. Educ. 2001, 84, 2001.
• ³⁰
  http://www.prudencepharma.com/Glibenclamide.html, accessed July, 2014.
  » link
• 31. Dash, S. K.; Khan, A. S.; Das, S. R.; Padhan, A.; Rout, D.; Behera, B. C.; Int. J. Pharm. Sci. Res. 2012, 3, 1433.
• 32. Rangel-Vázquez, N. A.; Rodríguez-Felix, F.; Computational Chemistry Applied in the Analyses of Chitosan/Polyvinylpyrrolidone/Mimosa Tenuiflore, Science Publishing Group: Hong Kong, 2013, cap. 2, pp. 25-26.
• 33. Hwang, Y. J.; Krasieva, T.; Lyubovitsky, J.; ACS Appl. Mater. Interfaces 2011, 3, 2579
• 34. Santoni, N.; Matos, M.; Müller-Karger, C.; Nicola, H.; Sabino, M.; Müller, A.; Revista Iberoamericana de Polímeros 2008, 9, 326.
• 35. Gunasekaran, S.; Rajalakshmi, K.; Kumaresan, S.; Spectrochim. Acta, Part A 2013, 112, 351.
• 36. Lertsutthiwong, P.; Noomun, K.; Khunthon, S.; Limpanart, S.; Prog. Nat. Sci. : Mater. Int. 2012, 22, 502.
• 37. Yao, C. H.; Liu, B. S.; Chang, C. J.; Hsu, S. H.; Chen, Y. S.; Mater. Chem. Phys 2004, 83, 204.
• 38. Levinton-Shamuilov, G.; Cohen, Y.; Azoury, M.; Chaikovsky, A.; J. Forensic Sci 2005, 50, 1367.
• 39. Bahri-Najafi, R.; Tavakoli, N.; Senemar, M.; Peikanpour, M.; Res. Pharm. Sci. 2014, 9, 213.
• 40. Rodríguez, R.; Tesis de Licenciatura, Universidad de Oriente, Cumaná, Venezuela, 2010.
• 41. Mirzaei, E.; Faridi-Majidi, R.; Shokrgozar, M. A.; Paskiabi, F. A.; Nanomed. J. 2014, 1,137.

Publication Dates

• Publication in this collection
  24 Oct 2014
• Date of issue
  2014

History

• Accepted
  27 June 2014
• Received
  10 Apr 2014