Fabrication and Computational Study of pH-Responsive Chitosan/Poly(HEMA-co-2-HMBA) Microparticles for Controlled, Site-Specific Doxorubicin Delivery

TL;DR

Scientists created pH-sensitive microparticles from chitosan and a polymer to deliver doxorubicin, a cancer drug, in a controlled and targeted way.

Contribution

A novel pH-responsive chitosan-based drug delivery system was developed and computationally analyzed for site-specific doxorubicin delivery.

Findings

Microparticles achieved ~91% doxorubicin loading capacity due to carboxylic acid groups.

pH-dependent drug release was observed, with faster release at pH 5.5 (tumor-like) compared to pH 7.4.

DFT calculations showed hydrogen bonding between doxorubicin and microparticles influenced drug release behavior.

Abstract

As a chitosan (CTS)-based drug carrier (DC) for doxorubicin (DOX) delivery, poly(2-hydroxyethyl methacrylate-co-2-hydroxy-4-N-methacrylamidobenzoic acid) [poly(HEMA-co-2-HMBA)] (PHCH) was successfully grafted onto chitosan to fabricate DOX-loaded microparticles, and their in vitro release behavior was systematicaly investigated. Morphological characteristics were analyzed using scanning electron microscopy (SEM) and transmission electron microscopy (TEM), while DOX loading was validated through Fourier-transform infrared (FTIR) spectroscopy and thermogravimetric analysis (TGA), comparing pure and drug-loaded microparticles. The maximum loading capacity (~91%) was attributed to the presence of abundant carboxylic acid groups, which imparted pH responsiveness during in vitro DOX release. Furthermore, density functional theory (DFT) calculations revealed that hydrogen bonding interactions between DOX and the functional groups of the microparticles strongly influenced encapsulation efficiency (EE%), drug loading (DL%), and release behavior. The fabricated microparticles exhibited pH-dependent DOX release, with accelerated and more complete release at tumor microenvironment pH 5.5 compared to physiological pH 7.4. These results demonstrate that PHCH grafted CTS microparticles are promising candidates for controlled and site-specific anticancer drug delivery.

Figures

Full text

1. Introduction

Cancer is a complex genetic disease caused by specific changes in genes that regulate cell proliferation and homeostasis. Thus, cancer cells undergo unregulated growth, forming a mass or lump of tissue. Furthermore, these cells have the potential to spread to other locations in the body, where they can grow in an irregular manner to produce a malignant growth. Because chemotherapy is one of the best methods to eradicate cancerous tissues [1,2,3], efficient drug delivery systems (DDSs) should be developed to increase the bioavailability of drug molecules in target tissues.

Drug delivery efficiency is an important parameter in the treatment of diseases [4,5,6]. Therefore, the development of a high-efficiency DDS is a challenging task, as most pharmaceuticals face issues such as lower adsorption due to large molecular weight, slow onset of action due to systemic drug administration, and non-selectivity of drug delivery, leading to potential side effects in off-target tissues [7,8,9,10,11]. Nevertheless, the non-selectivity of these anti-cancer agents causes unintended damage to healthy normal tissues, resulting in undesirable side effects [12,13,14]. In addition, their poor bioavailability often requires higher dosages, increasing toxicity in normal cells. To overcome these limitations, the development of targeted chemotherapeutics through either passive or active mechanisms is essential to enhance efficacy while minimizing side effects. Among these, stimuli-responsive targeting systems offer a powerful platform for precision oncology. Stimuli-responsive targeting systems constitute a sophisticated class of drug delivery platforms distinguished by their exceptional precision, tunability, and therapeutic efficacy compared with conventional passive and active targeting strategies [15,16,17,18,19]. They are engineered to dynamically respond to specific endogenous stimuli such as pH, redox gradients, enzymatic activity, or reactive oxygen species and exogenous triggers like magnetic fields, light irradiation, heat, or ultrasound [20,21,22]. Therefore, smart stimuli-responsive systems enable precise spatial and progressive control over drug release within the tumor microenvironment, thereby optimizing therapeutic efficacy while significantly reducing systemic toxicity [23,24]. Drug delivery systems based on pH-responsive mechanisms are particularly effective for tumor treatment, as tumor tissues exhibit a relatively acidic microenvironment, with pH values ranging from approximately 5.6 to 6.8 [25,26,27,28,29,30,31].

More recently, considerable efforts have been focused on the localized delivery of DOX with macromolecular-based drug carriers, including natural and synthetic polymeric structures (micro/nano particles), polymeric micelles, and liposomes [32,33,34,35]. Among these carriers, carbohydrate-based polymeric structures such as chitosan, pectin, etc., have gained significant interest due to their interesting properties [36,37,38]. Owing to their many attractive properties, such as excellent biocompatibility, biodegradability, non-toxicity and non-immunogenicity, chitosan has been widely used as an effective material in drug delivery systems [37,39]. CTS is second most abundant naturally occurring polymer after cellulose, predominantly found in the outer shells of a wide range of invertebrates and in the cell walls of fungi [40,41]. It is a polycationic heteropolysaccharide composed of β-1,4-linked D-glucosamine and partially N-acetylated D-glucosamine, which is generally obtained by the de-acetylation of chitin, extracted from the shells of crustaceans [41].

Doxorubicin (DOX) is a highly potent, hydrophilic anthracycline-type anticancer agent derived from Streptomyces peucetius, widely used to treat a broad spectrum of cancers [42,43,44]. However, its clinical use is limited by dose-dependent side effects, including severe cardiotoxicity and nephrotoxicity [45,46]. Recently, Banendu et al. developed a dual-targeting pH-responsive nanocarrier composed of superparamagnetic Fe_3_O_4_ and oleoyl-modified chitosan conjugated with folic acid (FA) to encapsulate Fe_3_O_4_ and doxorubicin (DOX) [47]. The nanoparticles showed pH-sensitive behavior, with faster DOX release (71%) at endosomal pH (~5.5) than at physiological pH (45%), and enhanced uptake in U87 glioblastoma cells via folate receptor-mediated endocytosis. In another study, a nanocarrier system was developed for controlled drug release within the tumor microenvironment, enabling preferential release of chemotherapeutic agents under tumor-specific pH and redox conditions [21]. This stimuli-responsive behavior enhanced drug accumulation at the target site and improved intracellular delivery, resulting in superior therapeutic efficacy compared to non-responsive systems. A dual pH- and temperature-responsive carrier composed of a PNIPAM–DOX hydrogel inside pH-sensitive polymersomes enabled controlled DOX release, showing faster release at acidic pH (5.5–6.5) and slower release at 37 °C due to hydrogel gelation [18]. In C26 tumor-bearing mice, this dual-responsive system achieved strong tumor inhibition with reduced systemic toxicity and minimal body-weight loss compared to free DOX. pH-responsive gelatin nanoparticles were developed for controlled doxorubicin (DOX) release, exhibiting significantly higher release (~95%) under acidic conditions (pH 3) compared to physiological pH (24%) [17]. Ultimately, this system offers a smart, predictive platform that combines pH-sensitive carrier design with data-driven modeling, enabling optimized drug release profiles and minimizing experimental complexity.

In this work, chitosan/PHCH DC microparticles were produced by free radical polymerization followed by crosslinking with GA. The microparticles were then encapsulated with DOX and characterized to investigate their drug release properties and the mechanism of release. Furthermore, computational studies were conducted to examine the interactions between DOX and the various functionalities present in the microparticles.

2. Results and Discussion

2.1. Characterization

2.1.1. Fourier-Transform Infrared (FTIR) Spectroscopy

The functionality of the DC microparticles and the loading of DOX were characterized by FTIR spectroscopy. The FTIR spectra of DOX, DC microparticles, and DOX-loaded DC microparticles along with CTS and PHCH polymer are presented in Figure 1. The FTIR spectrum of CST exhibits characteristic absorption peaks corresponding to -O-H, -N-H, C-H, and amide functional groups. A broad band in the range of 3200−3500 cm^−1^ is attributed to -O-H and -N-H stretching vibrations, while peaks around 2800–2900 cm^−1^ correspond to C–H stretching. Additionally, distinct bands at approximately 1642 cm^−1^ and 1542 cm^−1^ are assigned to the amide I (C=O stretching) and amide II (-N-H bending) vibrations, respectively. In contrast, the FT-IR spectrum of PHCH displays absorption peaks at 1719, 1665, and 1510 cm^−1^, corresponding to ester, amide I, and amide II groups, respectively. In the case of DC microparticles, a characteristic broad peak is observed at 3371 cm^−1^ along with shoulder peaks at 3432 and 3298 cm^−1^ corresponding to -O-H and -N-H stretching vibrations of chitosan and 2-HMBA. The peaks appearing at 1657 and 1590 cm^−1^ are due to carbonyl of primary amide group stretching and -NH bending vibrations, respectively. The absorption bands at 1153, 1078, and 1033 cm^−1^ are skeletal vibrations involving C-O stretching [39]. In the FTIR spectrum of pure DOX, the characteristic peaks at 3525, 3332, 1730, 1617, 1581, 870, and 804 cm^−1^ are assigned for -N-H, -O-H, C=O, N-H bend and wag, respectively [10]. The FT-IR spectrum of DOX-loaded DC microparticles, shown in Figure 1, shows peculiar peaks at 1720 cm^−1^ (C=O in ester group), 1619, and 1583 cm^−1^ (C=O, -NH bending in amide), which were broadened and shifted to higher frequency, with the concurrent disappearance of -N-H wagging peaks at 870 and 805 cm^−1^ in pure DOX. Thus, the above spectral information confirms the successful fabrications of both microparticles and the loading of DOX into polymer microparticles.

2.1.2. Thermo-Gravimetric Analysis (TGA)

TGA is an analytical technique used to evaluate the thermal stability of a material by monitoring the weight loss of the sample with temperature changes, when heated at a constant rate in a controlled environment. As presented in Table S1 and Figure 2, three samples, such as DOX.HCl, as prepared DC microparticles and drug-incorporated microparticles were characterized with TGA studies to understand their thermal stabilities and the entrapment of DOX.

The thermogram of DOX.HCl (Figure 2a) displays a characteristic two-stage decomposition curve, and 57.7% of total weight loss was observed in the analysis, ranging from 30650 °C, which states that the DOX would not completely decompose to volatile gases but remain as a charred mass in the crucible [48]. The first decomposition was observed in the temperature range between 147 and 257 °C, with weight loss of 27.50%, while the second weight stage observed in the temperature range of 257364 °C accounted for up to 12.79% of loss of DOX weight. The TGA curve of CTS/PHCH (Figure 2b) exhibits a single-stage decomposition beginning at 156 °C and continuing up to 384 °C, corresponding to a weight loss of approximately 54.70%, after which the residual mass stabilizes at 15.45%. In the case of the TGA profile of DOX-incorporated microparticles (Figure 2c), approximately 70% of the material decomposed across the entire temperature range of 30–650 °C. This weight loss is attributed to the reduced thermal degradation of the DOX-loaded microparticles due to the formation of the DOX–PHCH complex. The increased thermal stability of the DC microparticles upon DOX incorporation suggests strong ionic interactions between the NH_3_^+^ ions of DOX and the carboxylate ions of the 2-HMBA side chains in the PHCH backbone, along with additional hydrogen bonding interactions, as illustrated in Figure 2. Similar observations have been previously reported by Manocha et al. [49].

2.1.3. Morphologies

Figure 3 illustrates the SEM and TEM images of the morphologies of DC microparticles, which shows irregular shapes with average diameters of 300 ± 50 nm. The SEM images of DC microparticles, captured at 12 KX and 30 KX magnifications, reveal irregularly shaped particles with smooth surfaces and variable dispersity due to attachment of particles. The TEM images further support the SEM observations, revealing that the DC microparticles possess irregular yet well-defined shapes with uniform contrast, indicative of a dense internal structure. The particle sizes determined from TEM analysis range between 150 and 400 nm, which aligns with the microscale dimensions desirable for efficient drug carrier systems. These observations confirm the successful fabrication of well-defined and stable DC microparticles via redox-mediated graft polymerization, followed by phase inversion and GA crosslinking.

2.2. Drug Entrapment and Drug Release Studies

After loading with free DOX, the DC microparticles were characterized using TGA, FTIR, and drug loading (% DL) prior to the assessment of drug release. The drug encapsulation efficiency (% EE) and DL (%) are the general terms in the drug delivery, which affect the drug release studies directly and are important to maximize the therapeutic efficacy of the drug delivery system. However, the drug entrapment to the polymeric gel and its release depends on several factors, such as (i) the initial concentration of polymer and its functionality, (ii) the ratio between polymer and drug, (iii) the type of interactions between drug and polymer, and (iv) the type of and amount of cross-linker, which controls the swelling ratio of polymer gels. Table 1 shows the encapsulation efficiencies of various DCFs based on CTS microparticles. The EE (%), varying from 30.25 ± 1.6 to 39.91 ± 1.2% depending on the different DCFs, increased with increasing 2-HMBA and cross-linker amount for DCFs A to D and F, C and G, respectively.

Similar results were also found in DL (%), ranging between 54.51 ± 2.5 and 91.77 ± 0.9. These results indicate that the increasing EE (%) is attributed to the increase of polar carboxylic acid groups and cross-linking density in the network. Here, the higher values of EE (%) and DL (%) for DCF-A can be attributed to the higher possibility of hydrogen bonds with DOX because of the highest amount of 2-HMBA, whereas the lowest values for DCF-G, which has the highest amount of cross-linking density, are brought on by the least porosity of the microparticles.

As the microparticles are fabricated with the help of CTS and 2-HMBA, they possess a pH-responsive property, owing to the protonation in response to the slightly acidic medium, which makes the switching of DOX release from the DCF. Thus, the microparticles can show the on-demand release of DOX in tumor tissues due to their relatively lower pH compared to normal tissues. Considering this divergence, the two sets of in vitro release studies were performed for all seven DCFs under two different pH conditions, such as 7.4, a physiological pH, and 5.5, an acidic tumor condition. The required amount of DOX-encapsulated microparticles was dispersed in aqueous buffers and continuously stirred in the shaking incubator. Then, the release medium was subjected to UV-visible spectroscopy at predetermined intervals of time to determine the release of DOX from the microparticles, which was further used for the calculation of cumulative release. Finally, plots were drawn for the cumulative DOX release as a function of time.

DOX release kinetics of DC microparticles were studied with phosphate buffer saline using a bath method. Figure 4 shows the plot of cumulative DOX release behaviors from microparticles at pH 7.4. The release kinetics of microparticles are different for the all seven DCFs and depend on the amount of 2-HMBA and cross-linker used in the preparation of microparticles. All the DCFs showed an initial burst release of DOX from microparticles in the first hour of the study due to the drug molecules adsorbed on the surface of microparticles, followed by slower release rates, and eventually reaching a maximum cumulative release at 36 h. In the case of DCF-A, a slower and lower drug release was initially anticipated; however, an opposite trend was observed. This unexpected behavior can be attributed to the enhanced interactions between the main polymer chain (CTS) and the graft homopolymer chains of 2-HMBA. The phenolic -OH group of 2-HMBA possesses a low radical propagation rate and mild terminating characteristics, resulting in the formation of shorter polymer chains grafted onto CTS [50]. Consequently, the strong repulsive interactions between the carboxylate ions (-COO^−^) of the 2-HMBA chains promote faster swelling of the polymer network. During this initial polymer relaxation period, the soluble DOX molecules diffuse rapidly via concentration gradient. Over time, the release rate decreases due to stronger electrostatic interactions between the carboxylate ions and DOX molecules, resulting in a slower, sustained release profile. However, as the 2-HMBA content decreases, such repulsive forces are reduced due to the dispersion of negative charges along the copolymer chains.

Figure 5 displays the cumulative release of DOX at pH 5.5, showing a higher release rate of DOX from microparticles than at pH 7.4. This effect can be ascribed to both the inherent swelling ability of the polymer matrix and the interactions occurring between the polymer network and drug molecules [47]. However, CTS possesses a pKₐ in the range of 6.3 to 6.5, depending on its degree of deacetylation, and exhibits protonation of its amino groups below this pH, leading to enhanced water uptake by CTS-based materials [51]. Thus, the resulting switchable delivery of DOX was achieved through the pH-responsive properties of CTS. However, the -NH_2_ groups in chitosan are present as neutral molecules at physiological pH and protonated in acidic tumor microenvironments at pH 5.5 to form their ionized -NH_3_^+^. All the DCFs show slower release of DOX at pH 7, which can be attributed to the strong ionic and hydrogen bonding interactions between network groups and drug molecules in comparatively more swollen gels, whereas the faster and higher drug release observed at pH 5.5 may be due to the strong ionic repulsions between -NH_3_^+^ groups of CTS and DOX at a sufficiently close distance in a compact gel network. As a result, the higher drug release was observed at pH 5.5 due to the pronation of amine groups in the CTS and DOX causing the electrostatic repulsion between chains as well as chains and DOX; hence, the swelling and solvent penetration increased. Further, the amount of cross-linking has a bigger effect, since the amount of network swelling largely controls the release of DOX in DC microparticles. The network’s porosity plays a significant role in the diffusion of solvent and medication molecules. As a result, the sequence DCF-F, DCF-C, and DCF-G showed a slower release rate as the cross-linking quantity increased.

2.3. Drug Release Modelling

In order to determine the rate and mechanism of DOX release, the data obtained from the DOX release studies were fitted with three empirical kinetic models, such as Higuchi, zero-order, and first order. As presented in Table 2, for both pH conditions, the first order model was considered to be the best one, better than the others due to its higher regression (r^2^) and lowest sum of squares of residuals (SSR) values, which were used to decide the best-fit model. Furthermore, it clearly states that the DOX release behaviors from all DCFs were diffusion controlled.

Semi-empirical models such as the Ritger–Peppas and Peppas–Sahlin models were also applied to confirm the release type of diffusion that plays in the release of DOX from the microparticles. Table 3 reveals that the release of DOX at pH 5.5 was best with the Peppas–Sahlin model due to its highest r^2^ values of 0.9972 to 0.9733 and lowest SSR values of 0.0224 to 0.002, whereas at pH 7.4, the Ritger–Peppas model was best fitted because it had the highest r^2^ values of 0.9779 to 0.9601 and lowest SSR values of 0.021 to 0.0031. However, from the values of ‘n’, it can be known that the release of DOX followed a Fickian type due to the fact that all the ‘n’ values are less than 0.45. The fraction of drug release (F) calculated from Equation (9) indicates that DOX release at pH 7.4 was predominantly diffusion-controlled (Fickian behavior), and the diffusion is facilitated by electrostatic interactions. In contrast, at pH 5.5, drug release is influenced not only by diffusion but also by polymer chain relaxation and electrostatic interactions, resulting in a more complex, non-Fickian release mechanism. Furthermore, for DCF-C, DCF-F, and DCF-G, the porosity of the DC microparticles plays a critical role in drug release. As porosity decreases, the release profiles slow down, and the mechanism becomes more diffusion-controlled. In particular, DCF-G, which is expected to have the smallest pore size, exhibits a larger fraction of DOX release via Fickian diffusion at pH 5.5 due to stronger electrostatic repulsive forces. In contrast, at pH 7.4, the diffusion is more restricted by electrostatic attractions, leading to a slower, diffusion-controlled release.

2.4. Computational Study

All the density functional theory (DFT) simulations were performed using the Gaussian 16 program [52]. The geometry optimizations of DOX and DCFs A, C, and E of Table 1 were performed using a 6-31G(d) basis set and M062X functional [53]. The complexes formed from DOX and each of the DCFs were optimized at the M062X/3-21G level of theory, and the local minima of the optimized geometries were further confirmed with no imaginary frequencies using harmonic vibrational frequency analysis. Further, the corresponding complexation energies were obtained from single point calculation at the M062X/6-31G(d) level of theory. The selected hybrid meta-GGA-based Minnesota functional has been proven to be efficient at calculating the complexation energy [54]. The complexation energies obtained for DCFs A, C, and E with DOX were −114.28, −106.83, and −67.80 kJ/mole, respectively.

As shown in Figure 6, the high complexation energies can be attributed to the intermolecular hydrogen bonds (H-bonds) and π–π interactions formed between DOX and DCFs. It is interesting to observe that the number of H-bonds formed in each complex varies: six for DCF-A-DOX, four for DCF-C-DOX, and three for DCF-E-DOX. All H-bonds observed in complexes DCF-A-DOX and DCF-C-DOX are those formed between hydrogen and oxygen atoms, whereas in the case of complex DCF-E-DOX, two H-bonds were formed between hydrogen and oxygen atoms, and the remaining H-bond was formed between hydrogen and nitrogen atoms. In addition, as shown in Figure 7, the strong interaction between DOX and DCF is the supporting evidence of the charge transfer phenomenon observed in the frontier molecular orbitals. From the theoretical simulation results, it can be confirmed that DCF-A is the most efficient, over the other two DCFs.

3. Materials and Methods

3.1. Materials

Medium molecular weight CTS, 4-amino salicylic acid (4-ASA), HEMA, methacrylic anhydride (MA), ceric ammonium nitrate (CAN), and DOX were procured from Sigma-Aldrich (Milwaukee, WI, USA). The NaNO_2_, NaOH, methanol, and acetic acid were analytical grade, which were purchased from Tokyo Chemical Industry (TCI) (Tokyo, Japan) and used without further purification.

3.2. Synthesis of Monomer

The 2-HMBA monomer was synthesized (Scheme S1) according to the previous literature [55] by selective amidation of 4-aminosalicylic acid (4-ASA) with methacrylic anhydride (MA) at 0 °C using dry acetone as a solvent. To 50 mL of 10% 4-ASA solution of dry acetone under constant stirring, 30 mL of acetone solution containing the required amount of MA was added dropwise under nitrogen atmosphere for 15 min and allowed to react for an additional 2 h with continuous stirring. After completion of reaction, the solvent was removed with a rotary evaporator to obtain a crude product. The remaining solid residue was washed several times with mixed solvent of methanol/water (4/1, v/v). The crystals were recrystallized from the mixed solvent of methanol/water (2/1. v/v) and dried under vacuum to obtain pure 2-HMBA (yield, 60%).

3.3. Depolymerization of Chitosan

Medium molecular weight CTS was depolymerized by a chemical oxidation reaction to achieve the desired molecular weights, according to the methods described by Mao et al. with a minor modification [56]. Briefly, five milliliters of NaNO_2_ solution (14.0 g/L) in deionized water was added to the solution of CTS (2 gms of CTS in 100 mL of 6% (v/v) acetic acid) at room temperature under stirring for 1 h. The reduced-molecular-weight CTS formed was precipitated by adjusting pH to 9.0 with 4 molar NaOH solution. The white-to-yellowish precipitate obtained was filtered and washed 2~3 times with acetone. The precipitate was again dissolved in a minimum volume of acetic acid (0.1 M) and purified by dialysis (Sigma dialysis tubes MW cutoff 12 kDa) against deionized water with periodic bath changes thrice a day. The dialyzed product was freeze-dried and stored at 4 °C.

3.4. Fabrication of Microparticles

Microparticles with different drug carrier formulations (DCFs) were fabricated via redox-mediated graft copolymerization onto the linear chains of CTS, as illustrated in Scheme 1. The required amount of CTS was dissolved in 50 mL of double-distilled water and purged with nitrogen for 10 min. When the temperature reached 50 °C, the required amount of CAN was added for radical initiation on the CTS chains. Then, grafting on CTS was achieved by the addition of 10 mL of pre-dissolved monomer mixture of 2-HMBA and HEMA and continued to react for 2 h. After the polymerization was completed, the reaction solution was cooled to room temperature, and the pH of the solution was increased to 7 by adding dilute (0.5 N) sodium hydroxide solution. When the solution became opaque, a certain amount of glutaraldehyde solution was added to complete the formation of microparticles by cross-linking.

3.5. Characterizations

The DC microparticles were characterized by Fourier-transform infrared (FTIR) spectroscopy, differential scanning calorimetry (DSC), and electron microscopies such as SEM and TEM. The FTIR spectra of DC microparticles and their drug-loaded entities were recorded in an Infrared Spectrometer (JASCO spectrophotometer, JASCO Inc., Easton, MD, USA) in KBr discs. The morphology of the DC microparticles was evaluated by scanning electron microscopy (SEM) observations. The microparticles dispersed on carbon tapes were sputter-coated with gold prior to the acquisition of SEM images, using a scanning electron microscope (Carl Zeiss AG- Supra 25, Oberkochen, Germany). TEM analysis was performed with a transmission electron microscope (Hitachi H7600 TEM, Tokyo, Japan) after the deposition of samples on the copper grids covered with a perforated carbon film with the aid of water-dispersed samples. The TGA results were acquired using a thermal analyzer (Netzsch STA 449 F3 Jupiter, Forschungszentrum Jülich, Germany). The powdered samples (8~10 mgs) were crimped in an Al pan and heated at a 10 °C per minute temperature ramp over a temperature range of 35 to 650 °C under N_2_ atmosphere.

3.6. Drug Loading and Drug Entrapment Efficiency

The drug DOX was encapsulated into DC microparticles by stirring overnight with a DOX solution at a 50 μg mL^−1^ concentration. After being centrifuged, the microparticles were washed three times to remove free DOX and dried for further use. The supernatant DOX solution was assayed by the UV-Vis absorption signals at 485 nm to know the encapsulation efficiency (% EE). Then, the DC microparticles (5 mg) loaded with DOX were dispersed into 10 mL of distilled water and stirred continuously at 100 rpm in a vibrating incubator at 37 °C for about one day. Finally, the microparticle suspension was centrifuged, and the supernatant solution was assayed by the UV-Vis absorption signals at 485 nm, then the amount of DOX encapsulated into the microparticles and the % of drug loading (% DL) using the absorption intensity was calculated. The percentage of EE (%) and the DL (%) were calculated using the following equations:

$$[eqn]$$ $$[eqn]$$

3.7. Drug Release and Model Fitting

The release of DOX from the pH-sensitive microparticles was studied through the bath method, with the aqueous dispersions of DOX-encapsulated microparticles in the buffer solutions [57]. Typically, 2 mg of DOX-encapsulated microparticles were dispersed in 2 mL of release medium and stirred continuously in a shaking incubator with a shaking speed of 50 rpm at 37 °C (±0.5 °C). At predetermined time intervals, the samples were centrifuged, and the supernatant solution was subjected to UV-Vis, and the same amount of fresh release medium was added to the withdrawn samples.

The cumulated release rates were calculated by using Equation (3):

$$[eqn]$$

Each experiment was conducted in triplicate, and the data collected in this study were expressed as the mean value ± standard deviation (SD).

Then, the drug release data was fitted with different mathematical empirical models, like Higuchi’s square root of time equation, first-order, and zero-order, and semi-empirical models like the power law equation of Ritger–Peppas and weighted sum of non-linear time function of Peppas–Sahlin were fitted as follows:

Higuchi model:

$$[eqn]$$

First-order:

$$[eqn]$$

Zero-order:

$$[eqn]$$

Ritger–Peppas:

$$[eqn]$$

Peppas–Sahlin:

$$[eqn]$$

where Q_0_ is initial amount of drug present in the microparticles; Q_t_ is percentage of drug released at time t; and k_0_, k_1_, k_H_, and k_C_ are the rate constants of the respective equations. Here the empirical models describe the information about the types of drug release, while the semi-empirical Ritger–Peppas model derives the actual mechanism of drug release, i.e., whether the drug release belongs to Fickian, Non-Fickian, or case II (anomalous) transport. However, the Peppas–Sahlin models gives the fraction of drug release due to the diffusion at a given instant of time t using the following equation:

$$[eqn]$$

3.8. Computational Study

The calculations of density functional theory with a hybrid functional B3LYP (Becke’s three-parameter hybrid functional) used Gaussian 16 software. The optimized molecular structures of DOX and DCFs (4 repeating units of polymer and CTS rings) and corresponding vibrational assignments of DOX and DC microparticles were investigated with a 6-31G(d) basis set and M062X functional.

4. Conclusions

This work has described the design, synthesis, physicochemical characterization, and drug release evaluation of chitosan-based DC microparticles. By varying the comonomer (2-HMBA, HEMA) ratio, different DCFs of DC microparticles were produced that exhibited a wide range of physicochemical and drug release properties. FTIR and TGA studies confirm the uniform distribution of DOX in the DC microparticles, and DFT calculations indicate that the DOX molecules can form hydrogen bond interactions with DC microparticles. All the DCFs have optimal drug loading capacity, and the DOX release profiles are pH-dependent and controlled in manner. The higher loading capacity and faster in vitro DOX release properties of microparticles at pH 5.5 were due to the large number of polar groups present in the microparticles. Further, the same was confirmed from computational studies. The experimental results of the drug cumulative release were as expected according to the extent of protonation predicted by the Ritger–Peppas equation. Finally, this data suggests that the prepared microparticles are suitable to release DOX at cancer tissue to maximize the drug’s efficacy and reduce the side effects.