Abstract

Differences in morphology, pH, and electric charge of chitosan (CS) based hydrogels prepared by complexation with carboxymethylcellulose (CMC), carboxymethylated starch (CMS), and alginic acid (AA) at different polymers ratios and changing the order of addition were studied. CMC/CS and AA/CS hydrogels were amorphous and porous three-dimensional networks, with smaller pores at higher anionic polymer/CS ratios. Gelation time increased the agglomeration in the case of CMC/CS and CMS/CS gels. CMC/CS gels showed negative zeta potential values around −372 mV to −51 mV and CMS/CS gels in the range of −526 mV and −158 mV.

1. Introduction

Within the applications of biotechnology are the developments of therapeutic agents having applications in tissue engineering and controlled release of drugs. Ten of the twenty main world sales of medical treatments are biological products because the biotechnology offers advantages over traditional medicine such as the reduction of side effects and specific treatments [1]. Nanotechnology provides fundamental tools to control the release of these biological products and to design three-dimensional tissue scaffolds to support cell growth due to their similarity to the native extracellular matrix [2]. That is how the ability to control the properties and functionalities of hydrogel broadens the understanding of possible cell-matrix interactions at the nanometric scale, and comprehensive understanding of the processes involved in the self-assembly of hydrogels is essential for the development of relevant biotechnologycal applications [1].

Hydrogels are three-dimensional networks capable of absorbing large quantities of water or biological fluids; therefore, they attract interest in medicine [3]. Hydrogels are formed by chemical or physical crosslinking. Chemical crosslinking occurs by covalent bonds and physical crosslinking by hydrophobic interactions, hydrogen bonds, stereo complexes, and electrostatic or charge interactions [3]. Physical hydrogels have the advantage of responding to changes in the surrounding environment such as pH, temperature, changes in ionic strength, or its state (swelling shrinkage) [4]. Noncovalent crosslinking allows the polymers to form a gel to be altered by the elimination of the weak links with changes of external factors such as pH, ionic strength, and temperature. These changes could be used as triggers to deliver a drug in a specific site and can also help keep the effect of the drug over a period of time, days, or even weeks [57].

Within the development of therapeutic biomaterials is the development of hydrogels derived from biopolymers that are biocompatible with the human organism. Polysaccharides have been studied during the last years as raw materials in the encapsulation of drugs forming a mesh structure by charge interaction between anions and cations. Polysaccharides have peculiar properties such as good solubility in aqueous environments, high stability, null toxicity, biocompatibility, and biodegradability [810].

Several procedures for the manufacture of hydrogels by noncovalent crosslinking with carboxymethylcellulose, carboxymethylated starch, or alginic acid as anionic polymers with chitosan as the cationic polymer have been reported; each couple (polyanion/polycation) generates a hydrogel with different morphology [8, 1113].

In hydrogel formation, both concentration of polymers and sequence of addition, among other factors, play important roles in the morphology and size of the hydrogel [1416]. The aim of this research is to investigate the morphology, state aggregation, pH, and electric charge of hydrogels obtained by polyanion (carboxylated starch and alginic acid) and polycation (chitosan) complexation, by varying the order of polymer addition and polymer-polymer ratios.

2. Materials and Methods

2.1. Materials

Carboxymethylcellulose (CMC) sodium salt, average molecular weight (Mw) ∼250,000 and 0.7 degree of substitution (CAS 9004-32-4), chitosan (CS), deacetylation 75–85%, medium molecular weight (CAS 9012-76-4), and alginic acid (AA) of medium viscosity (CAS 9005-38-3) all were the products of Sigma-Aldrich. Carboxymethylated starch (CMS) was obtained from DFE Pharma at Foxhol, The Netherlands. The polymers were used as supplied.

2.2. Sequence of Polymer Addition in the Ionotropic Gelation

In order to evaluate the effect of changing the sequence of addition and concentration of the polymers, the morphology of hydrogels was examined. The concentration of CMC was fixed at 0.75% w/v, while two CS preparations were employed (0.01 and 0.75% in acetic acid 1% v/v), varying also the sequence of addition. At higher concentrations, CMC solutions were very viscous and the polymer was not completely solubilized.

2.2.1. CMC on CS

Approximately 10 mL of CMC solution was dropped through a needle (21 G) from a plastic syringe into a beaker containing 10 mL of any CS solution (0.01 and 0.75%), under gentle stirring (600 rpm) for 5 minutes at room temperature. Then, the mix was stirred for another 25 minutes, and the hydrogel was kept 12 hours at room temperature.

2.2.2. CS on CMC

Approximately 10 mL of any CS solution (0.01 and 0.75%) was dropped through a needle (21 G) from a plastic syringe into a beaker containing 10 mL of CMC solution, under gentle stirring (600 rpm) for 5 minutes at room temperature. Then, the mix was stirred for another 25 minutes, and the hydrogel was kept 12 hours at room temperature.

All the samples were centrifuged at 3300 rpm for 30 minutes, and the precipitates were removed from the supernatants. The gel samples were then washed with water and separated again by centrifugation. The washing and centrifugation steps were repeated twice, and the samples were observed under a microscope (Leica DMREB, software LAS, version 7.1, Camera Leica DFC320).

2.3. Free Swelling Capacity

The free swelling capacity of hydrogels was determined as follows: the samples (gels) were dried until constant weight in an oven at 60°C; then they were weighed and rehydrated in distilled water for 24 hours, and the excess water was removed with a dry paper, depositing the gel on it a few seconds before weighting. The weight difference was reported, and the gel samples were dried again. This procedure was applied to gels formed with 0.75% of both (CMC and CS), varying the sequence of addition.

2.4. Variation of Polymer-Polymer Ratio

CS concentration was fixed at 0.01% w/v (dissolved in acetic acid 1% v/v) in order to avoid viscosity problems when mixed with anionic polymers at high concentrations. The concentration ranges chosen were as follows: CMC between 2% and 0.025% w/v; CMS between 3.2 and 0.4% w/v, and AA between 0.5 and 0.04% w/v. Six different concentrations were used within each range. Polymer-polymer ratios are shown in Table 1. The gels were prepared as described before. The addition order CMC/CS was chosen due to the higher production of hydrogel compared to the sequence CS/CMC.

The hydrogel samples obtained at different polymer-polymer ratios were dried in an oven at 60°C until constant weight to calculate the yield as follows [17]:where is the dry weight of the hydrogel and is the weight of products used for the gel preparation.

2.5. Gel Morphology

To characterize the morphology of the wet gels at different polymer-polymer ratios, a drop of the sample was put on a slide with a drop of water to disperse the gel before the observation. Then, a coverslip was put upon, and the sample was observed under an optical microscope (Leica DMREB, software LAS, version 7.1, Camera Leica DFC320).

The dried gel morphology was observed from gel samples frozen with liquid nitrogen and then lyophilized. Samples were analyzed by using a scanning electron microscope Hitachi TM-1000 (Hitachi High-Technologies Inc., Japan).

2.6. pH and Zeta Potential

The charge distribution (zeta potential) and pH of the gels were measured after a gelation time was established (between 20 and 24 hours after the preparation) to be associated with the polymer-polymer ratio used to prepare the hydrogels. The zeta potential was measured in a Zetasizer Nano ZS90 with DS170 cuvettes, using acetic acid as a solvent for chitosan and the hydrogels. Water was used as a solvent for the other polymer systems. The pH was measured using a benchtop pH meter.

2.7. FTIR Characterization of Gels

FTIR spectra of dried gel samples were recorded in a PerkinElmer FTIR spectrometer, model Spectrum GX with an attachment for attenuated total reflectance (ATR), and a diamond crystal. Spectra were taken between 4000 and 600 cm−1 with a resolution of 4.00 cm−1 and 16 scans.

2.8. Effects of Gelation Time on Hydrogel Morphology

In order to evaluate the effect produced by the time-to-reach the gelation on the morphology of hydrogels, the mix of polymers was agitated for 25 minutes and the complexation was stopped; at that time, the polymeric complex was centrifuged at 3300 rpm for 30 minutes and the precipitate was removed from the supernatants. The gel samples were washed with water and separated again by centrifugation. Hydrogel surfaces were analyzed by zeta potential in a Zetasizer Nano ZS90 with DS170 cuvettes, and changes in the morphology were evaluated by using optical microscopy (Leica DMREB, software LAS, version 7.1, Camera Leica DFC320).

3. Results and Discussion

3.1. Effect of the Order of Polymer Addition on Gel Formation

CS and CMC form intermacromolecular complexes by the strong interaction between their functional groups, such as electrostatic interaction and hydrogen bonding [14]. Hydrogels were obtained for all studied sequences of addition, CMC on CS (CMC/CS) and CS on CMC (CS/CMC) as reported by several authors [14, 16, 18]. It was observed that the order of polyelectrolyte mixing affected the hydrogels where, for example, CS/CMC hydrogels were formed at pH higher than that of CMC/CS [14, 19, 20], and it produces differences in morphology and yield, similar to that reported by Fukuda [14].

At a CMC/CS ratio of 75 : 1 (low amount of CS), no precipitation or turbidity was observed when adding CS over CMC; this behavior may be due to an excess of anionic charge (COO) of CMC and a low pH that reduced the number of interactions with groups. By changing the sequence of addition and dropping the CMC solution onto a diluted CS solution, the contact between polymers was improved and produced a cloudy suspension (Figure 1) yielding more gel particles [21]. These conditions rendered gel particles in the form of sponges (Figure 1(b)) having a pore size around 3 micrometers.

In the case of CMC added on CS at a ratio of 75, agglomeration could be observed under the microscope, while the addition of CS over CMS showed only a few visible particles (Figure 2).

Moreover, at a CMC/CS ratio of 1/1, particles of the order of millimeters were observed. This gel showed more swelling than the hydrogel obtained by the addition of CS onto CMC (Figures 3(a) and 3(b)), thus indicating the importance of the order of addition of the polymers on gel morphology.

Complex and highly aggregated heterogeneous hydrogels were obtained by mixing polyelectrolytes of high or similar molecular weight. In this work, CMC and CS had an average molecular weight of 250 kDa. Molecular weight is a factor in the process of agglomeration. Among other factors (for instance, pH), the higher the molecular weight, the more the ionic sites that would be available to interact [22].

The CMS/CS dried hydrogel prepared with a ratio of 1/1 was a slim and translucent sponge that presented a network of particles of different sizes. The yield also changed with the sequence of addition. Thus, the hydrogel made adding CMC/CS had a yield of 70.85% with a deviation of 1.67%, and the inverse order of addition produced a yield close to 100% with a deviation of 1.86%.

Hydrogels prepared at a ratio of 1/1 had a swelling percent of 4325 ± 46% when CMC was added onto CS, while the one formed with the inverse order had only 856 ± 20% swelling. The swelling was calculated with the following equation [17]:where is the dry weight and is the weight of the swollen gel.

3.2. Influence of Varying Polymer Ratio on Hydrogel Yield

Table 2 shows the hydrogel yields for each different combination of polymers. The CMC/CS hydrogels had a yield around 20% at almost all ratios; this result differs from those reported by Fukuda [14], in which the yield of CMC/CS hydrogels increased with the concentration of CS. The yield for AA/CS and CMS/CS hydrogels had a tendency to decrease with more CS added.

According to Gåserød et al. [23], the core of the hydrogel is formed by the polymer that is dropped on a solution of the polyelectrolyte of opposite charge; in the present work, the anionic polymer (AA) was dropped over the cationic CS, thus forming a core of AA covered by the CS. It was observed that increasing the alginate concentration improved the yield of the gel. Likewise, with an increase in the amount of AA, the pH of the gel also increased [23]. Generally speaking, at a basic pH, CMC chains adopt a linear chain configuration and shrink as the pH decreases. AA may exhibit a loop formation at a pH around 3, and the CMS is in a state of pregel [24]. Therefore, in the case of CMC/CS, the yield of the gel improved as the pH decreased. The yield effect was probably due to a tighter structure of the CMC, which in turn could be more easily coated by the CS. These results are in agreement with those reported elsewhere [25].

The initial pH to polymer solutions were CMC 6.31, CMS 6.61, AA 7.76, CS 2.28, and the different interactions changed the pH between 2.6 and 3.8 for CMC/CS, 2.95 and 3.55 for CMS/CS, and 2.65 and 3.05 for AA/CS.

3.3. Effect of Polymer-Polymer Ratio on Hydrogel Morphology

CMC/CS hydrogels produced cloudy suspensions at higher ratios (ratio > 200), which did not precipitate quickly, whilst an agglomerated hydrogel was observed when the CMC/CS ratio decreased until a maximum agglomeration at a value ratio of 21.

Gels prepared from CMS also showed changes in appearance, in relation to the order of addition and polymer ratios. Hydrogels at the highest CMS/CS ratio (320) formed a viscous suspension; when the fraction of CMS was reduced, a cloudy suspension was produced. The increase in viscosity is due to a large amount of polymer in the suspension, which prevents the CS from dispersing quickly into the CMS; then by reducing the CMS/CS ratio, the viscosity of the mix is reduced, and CS can interact more easily with the polyanion. Since CMS is a pregel in the form of particles, a lower concentration favors the CMS-CS hydrogen interactions because there are no the impediments for hydrogen bonds to be created between the same starch chains due to the minor amounts of COO that are dissociated at lower pH [23, 26]. The decrease in viscosity is a demonstration that intermolecular complexation between starch chains does not occur at these conditions [27]. It might be assumed that, at the beginning of the addition, the CMS–COOH groups at the surface of the starch granules inside the aggregates formed strong hydrogen bonds that prevent interaction between starch chains [28].

On the contrary, for all the AA/CS ratios, cloudy suspensions were obtained. Other ratios outside the range not reported here also produced precipitates.

According to Devi et al. 2012, the optimal ratio is achieved when the supernatant has a lower viscosity [18]. For this reason, an insoluble complex is formed when the polyanion concentration decreased. Schatz et al. proposed that hydrogel formation leads to flocculation only at balanced charge ratios; however, if interaction may occur well at other ratios [29], it would prevent diffusion and good contact between polymer molecules.

3.4. Morphological Gel Characterization

Figure 4 shows an increase in the agglomeration of particles for CMC/CS hydrogels when the ratio decreased. Lower ratios gave a more structured hydrogel network with apparent longer segments.

Porous sponges were obtained when samples of CMC/CS hydrogels were freeze-dried. The SEM micrograph at a ratio of 200 (Figure 5(a)) shows a sponge with smaller pores than those exhibited when the ratio equals 75 (Figure 5(b)).

Devi et al. reported that the nature of gel formation for CMC/CS complexes depended on the relative quantity of CS until reaching a point where the interaction improved, leading to an increase in gel yield [18]. Thus, the rise in the agglomeration could be related to a higher interaction between the polymers.

A micrograph of a CMS/CS hydrogel sample (Figure 6) revealed oval particles similar to starch grains. When the polymer-polymer ratio decreased, the hydrogels exhibited a less dispersion of the particles.

The CMS/CS hydrogel is the only one that remained in powder form after freeze-drying. SEM micrographs (Figure 7) confirmed the presence of oval particles like those reported by Assaad et al. [30].

AA/CS formed hydrogels that looked like an agglomerated network (Figure 8), similar to what was reported elsewhere [31]. Few particles lower than 10 μm could also be observed. At this scale, morphological differences are not very noticeable between different polymer-polymer ratios.

The effect of the AA/CS ratio on hydrogel morphology could be observed under SEM microscopy (Figure 9). The mixture of alginic acid as polyanion and chitosan as polycation produced a three-dimensional crosslinked hydrogel with diverse porous structure; the gels were more homogeneous at lower ratios. This fact was also reported by Baysal et al., who observed that a similar polyelectrolyte complex led to a fibrous and porous hydrogel [32].

In all ratios studied for each polymer combination, a hydrogel was formed. It suggests that an effective interaction occurred between functional groups of polyanions and chitosan [19], making clear that the morphology of hydrogels depended on the preparation conditions [15].

3.5. pH and Zeta Potential

Table 3 reports the pH of the polymer solutions before complexation. The pH of the suspension changed with the polymer-polymer ratio, and some studies reported elsewhere have experimented controlling the pH with buffers, acids, or bases [14, 18]. In this work, the pH of the suspensions was low at lower ratios of the anionic polymer but increased at higher ratios [19]. The greater pH change was obtained for CMC/CS hydrogels (Figure 10).

The acidic pH in the polyelectrolyte complex is due to the presence of H+ ions (protonation) partially released by the acetic acid in which chitosan was solubilized, which also contributed to an increase in the conductivity. At acidic pH values, lower than 4-5, carboxyl groups are in the protonated form, which allows the formation of intermolecular hydrogen bonds between carboxyl groups of the polyanions and chitosan hydroxyl groups [33]. Therefore, the increase or decrease on hydrogel agglomeration for each combination of polymers used in the work that was mentioned before can be explained by the formation or absence of hydrogen bonds (Figure 10).

Another reason for the agglomeration is the unbalanced COO/N ratio in the solutions [19]. Near to the transition from positive to a negative charge (near to neutrality), the aggregation of particles was observed.

Figure 11 shows the effect of the polymer-polymer ratio on the zeta potential of the final solution. The solution charge (mV) has a monotonic tendency to decrease as the polymer-polymer ratio increases, probably due to an increment of negative ions from the polyanions [19].

A high CMC/CS ratio (i.e., 200) resulted in a cloudy suspension with negative zeta potentials (−962 mV). Sæther et al. indicated that working with an excess of the major component stabilizes the particles and no aggregation is obtained [19]. For zeta potentials between −372 mV and 700 mV, the hydrogel tended to flocculate. An abrupt change from negative to a positive charge is also observed. Although the zeta potential of AA/CS systems was negative, the interaction between polymers produced a gel.

On the contrary, regarding the CMS/CS hydrogels, flocculation was stronger at lower ratios (i.e., positive zeta potential). However, no turbidity was observed at +850 mV. Micrographs of individual particles are shown in Figure 12.

With an excess of one of the polymers, the zeta potential had a large positive or negative value. Nonaggregating particles were formed due to an excess of polyanion, which stabilized the hydrogel and prevented the agglomeration [19]. In each case, higher polyanion concentrations did not produce aggregation. This may be due to adsorption of excess polyanion, imparting a negative gel surface charge that would prevent agglomeration.

3.6. FTIR Characterization of Gels

Figures 1315 show the IR spectra of CS, CMC, CMS, and AA and polyelectrolyte complex of CMC/CS, CMS/CS, and AA/CS.

The spectrum of chitosan had the characteristic bands at 1640 and 1560 cm−1 related to the vibration of carbonyl bonds (C=O) of the secondary amide and protonated amino groups in CS [34, 35] (Figure 13). Bands identified at 1420 and 1380 cm−1 belonged to methylene and amide II [36] groups bending, respectively. The band at 1310 cm−1 corresponded to -CH3 stretch in acetyl groups, while the band at 1260 cm−1 was related to C-O-H groups [34]. Bands at 1150, 1070, and 1030 cm−1 were related to the symmetric and asymmetric stretches of C-O-C bonds [3436]. The band at 890 cm−1 was attributed to C-H wagging of anomeric carbon in the saccharide structure [34].

The characteristic bands for the CMC are 1590, 1412, 1315, 1154, 1102, 1050, and 1020 cm−1 (Figure 13). The bands of interest are 1590 and 1412 cm−1, which correspond to asymmetrical and symmetrical stretching of -COO groups [36, 37]. The rest of the signals correspond to specific groups of carbohydrates as discussed previously.

The CMC/CS complex had the characteristic bands of both chitosan and carboxymethylcellulose (Figure 13). The ratio CMC/CS of 200 had all the band characteristics of CMC but also showed the band characteristic of chitosan, i.e., 1380 cm−1; the band at 1060 cm−1 (displacement of 1070 cm−1 of CS and 1050 cm−1 of CMC) was more pronounced than the band at 1030 cm−1. The ratio CMC/CS of 32 had a new band in 1725 cm−1 that began to show up due to C=O stretching [37]. The signal at 1640 cm−1 of the amide in CS diminished greatly due to the greater amount of CMC, and thus the CMC signal at 1590 cm−1 predominating, and it is the characteristic of intermacromolecular complex formation [36]. In the CMC/CS ratio 200, the signal at 1590 cm−1 (-COO) is the predominant band. Both bands at 1412 and 1315 cm−1 corresponded to CMC. Signals at 1380 cm−1 corresponding to -CH2 of CS appeared also as a weak peak in this complex system.

The FTIR spectrum of CMS (Figure 14) has the characteristic carboxylate absorption band at 1570 cm−1 [38]. The -CH2 symmetrical stretch band was found shifted to 1420 cm−1 [39]; previous bands, along with the 1320 cm−1 band, were the characteristics of CMS [40]. In addition, bands occurred at 1150, 1080, and 1000 cm−1, arising from C-O bond stretching vibration [40].

The CMS/CS complex, ratio 320, had more characteristic bands of the carboxymethyl starch than those of chitosan (Figure 14). The spectrum shows the band 1570 cm−1 characteristic of COO groups of CMS [38].

The biggest difference between the two CMS/CS ratio 50 and CMS/CS ratio 320 complexes was in the wavenumber range 1570–1325 cm−1. The CMS/CS gel at ratio 50 did not show representative bands, while the CMS/CS gel at ratio 320 was more similar to the CMS spectrum, specifically the bands that corresponded to the vibration of carbonyl bonds (COO, 1570 cm−1) of the CMS [38].

The AA spectrum (Figure 15) shows an important band at 1600 cm−1 due to C=O stretch of the acid group [35, 41]. Furthermore, it presented a band at 1400 cm−1 probably due to C-OH deformation vibrations with the contribution of O-C-O symmetric stretching vibration of the carboxylate group [41]. Between 1290 and 1080 cm−1 appeared bands of C-O assigned to C-C-H and O-C-H deformation, C-O stretching, and C-O and C-C stretching. At 1020 cm−1, there was the characteristic absorption band of C-O stretching [35, 41]. The bands at 940 and 880 cm−1 were assigned to C-O of uronic acids and to C1-H deformation of β-mannuronic acid units, respectively. The band at 820 cm−1 possibly belonged to mannuronic residues [41].

The AA/CS complexes at a ratio of 5 and 50 had characteristic bands more related to AA than CS (Figure 15). The spectrum of the gel at ratio 5 had a band at 1730 cm−1 due to C=O stretch. The band at 1600 cm−1 appeared in both complexes, but more accentuated in AA/CS ratio 5. It was also noticed a shoulder at 1530 cm−1 probably due to N-H from CS. However, the signal at 1080 in AA/CS ratio 50 was better defined than in the gel at ratio 5. Also, the band at 1020 cm−1 appeared to be wider than the same band in AA. The band at 1150 cm−1 resembled more likely the one in CS.

Higher polymeric ratios decreased the agglomeration for hydrogels of the types CMC/CS and CMS/CS, while the hydrogel of AA/CS was porous and more homogeneous at smaller ratios. Agglomerated hydrogels (lower ratios) exhibited greater zeta potential values and smaller pHs.

3.7. Effect of Cure Time on Hydrogel Morphology

A more agglomerated hydrogel was obtained 12 hours after mixing for gelation. An extended aging time increases the apparent degree of crosslinking because more acid groups of CMC can react with amino groups of CS [42]. This process was accompanied by a decrease in the zeta potential of the system. CMC/CS hydrogel at ratio 32 had a charge of −276 mV before the gelation time, changing to −372 mV after gelling (Figure 16).

The CMS/CS had the same behavior as CMC/CS complexes. Before gelation time, CMS/CS at ratio 64 had a charge of −25 mV, while the gel with cure time had a potential of −293 mV (Figure 17).

Differences occurred with AA/CS at a ratio of 5 (alginic acid 0.28% w/v), which had a potential of −476 mV before gelation, keeping approximately the same z-potential after the gelation time (−444 mV). Also, there were no major changes in the morphology of the gel, except that more swelling occurred as the gelation time was reached (Figure 18).

4. Conclusions

The main objective of this work was to determine the effect of the polymer-polymer ratio and the order of addition on the morphology of hydrogels that could be promising materials for drug encapsulation. It can be concluded that the sequence of addition of polymer had indeed an effect on hydrogel formation, causing more aggregates by dropping the anionic polymer over the cationic polymer at lower concentrations.

The polymer-polymer ratio had a great influence on hydrogel morphology of CMC/CS, CMS/CS, and AA/CS. Higher polymeric ratios decreased the agglomeration for hydrogels of the type CMC/CS and CMS/CS, while the hydrogel of AA/CS was porous and more homogeneous at lower ratios. Agglomerated hydrogels at lower ratios exhibited greater zeta potential values and lower pH values.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest in the publication of these results.

Acknowledgments

The authors are grateful to the Alianza del Pacifico México for the financial support of the research stay at the Departamento de Madera, Celulosa y Papel of the Universidad de Guadalajara. The authors also acknowledge the support provided by Colciencias for financing the researcher and CIDI (UPB) for the financial support through internal project.