Forecasting in vivo pharmacokinetics of metformin HCl floating beads using Gastroplus® PBPK

Materials

Na alginate was purchased from (Fine chem limited) and the Metformin hydrochloride was kindly gifted from Samara Drug Industry (Iraq). Calcium chloride (99.89%), Calcium carbonate and isopropyl alcohol were purchased from BDH, Ind., England (as listed in Table 1).

Methods

Formulation of floating beads

Different concentrations of sodium alginate solution were prepared for use in the MH floating beads formulation, as shown in Table 2. Each solution was made by using a magnetic stirrer to dissolve the appropriate amount of Na alginate in 100 ml of deionized water with gentle strirring. An alginate solution containing 1 gm MH and 1.5 mg calcium carbonate as a gas-forming agent was mixed well by stirring at a constant speed 200 rpm at room temperature. All remaining air bubbles in the dispersion were removed using a sonicator for 30 minutes. The resulting dispersion was introduced to a 100 ml solution of 5% (w/v) calcium chloride (fused) in isopropyl alcohol 5% v/v (isopropyl alcohol is used as dispersing agent and cross-linking agent, so might be these properties play an important role in uniform bead formation) at room temperature via a 23-gauge syringe needle.⁵ The mechanical strength of the beads was improved by leaving them in the curing solution for 30 minutes. The beads went through many rounds of filtration and washing with water, after which they were air-dried for 24 hours.⁴

Percentage of yield

The percentage of yield depends on the ratio of polymer and gas forming agent, which the following formula can calculate¹⁴:

Percentage of entrapment efficiency %

Fifty milligrams of MH beads were crushed and added to a 100 ml solution of 0.1N HCl pH 1.2 to be filtered after 24 hr of incubation with constant stirring; then MH was analyzed using a UV/VIS spectrophotometer (UV Spectrophotometer, 1601, Shimadzu, Japan) at 233 nm, the maximum wavelength and dilution of the filtrate with 0.1N HCl pH 1.2. The following equation was used to calculate the entrapment efficiency. The individual batch should be examined for drug content in triplicate.¹⁴

Floating properties (Buoyancy lag time and duration of buoyancy)

Buoyancy lag time starts from the beads’ introduction into the medium until their buoyancy to surface of the dissolution vessel, and the buoyancy duration means the time for which the beads constantly floated on the surface of the medium. Both tests were performed during the in vitro release study by visual observation. Buoyancy lag time was determined by weighing equivalent to 500 mg of MH beads and placing them into a dissolution vessel paddle type containing 900 ml of 0.1 N HCl, pH 1.2 at 37 ± 0.5°C. at 50 rpm.¹⁵ All the determinations were conducted in triplicate.

In vitro release study

In vitro, release study investigations were carried out using a USP Dissolution apparatus Type II. The dissolution medium was 900 ml of simulated gastric fluid 0.1 N HCl (pH 1.2) at 37±0.5 ^oC. A sufficient quantity of beads equivalent to 500 mg of MH was placed in the dissolution medium. The paddle speed was limited to 50 rpm, and 5 ml samples were taken and replaced with fresh medium hourly for up to 12 hr. The withdrawel samples were analyzed using UV/VIS spectrophotometer (UV Spectrophotometer, 1601, Shimadzu, Japan) at 233 nm.

Morphological examination

The selected beads were viewed under Field Emission Scanning Electron Microscope (FESEM) (Inspect^TM F50, FEI company, USA). The surface of the beads and their cross sections were coated with gold-palladium under an argon atmosphere using a gold sputter module in a high vacuum evaporator.

Fourier Transmittance Infrared (FTIR)

Shimadzu- 8300, Japan FT-IR spectroscopy was utilized for Na alginate, MH as powder, and the selected drug-loaded beads were milled and then put in a KBr press. The spectra were taken from 4000 to 400 cm⁻¹.

In silico modeling for metformin absorption

The Gastroplus^® software (version 9.8.2, SimulationPlus Inc., Lancaster, CA, USA ‘SLP Cloud Access /acadmic access’ helped construct an MH model to gain a prediction of in vivo pharmacokinetic parameters for in vitro multiple unit beads release data. Table 3 presents input parameters related to MH physicochemical and pharmacokinetic properties taken from the literature and/or in silico estimated. The constructed model relies on divided segments of the gastrointestinal tract according to Advanced Compartmental Absorption and Transient (ACAT) model composed of the following: stomach, duodenum, jejunum 1 and 2, ileum 1-3, caecum and ascending colon. The model construction depended on 500 mg of each intravenous bolus and immediate–release oral tablet of MH. The GetData Digitizer version 2.26.0.2 software was used to extract metformin 500 mg of intravenous and oral immediate-release tablet data of healthy individual aged 42 and weight 63.4 kg, which was the base of the PBPK.¹⁶ The clearance values shown in Table 3 were obtained by applying the PKPlus software module to 500 mg intravenous metformin and clarified a three-compartmental model.

As MH was classified III according to the Biopharmaceutical classification system (BCS), which is water soluble with low permeability;¹⁷ hence, the PBPK model was set as permeability limited.

MH was reported as a substrate for the organic cation transporters (OCTs), which are influx transporters starting with OCT1, primarily found in the human liver, and OCT2 transporters are located mainly in the kidney. Thus, these two transporters were included in the PBPK model. Nonetheless, the OCT3 were not included the simulation of the MH model as they showed low affinity to MH and its expression in the region of the human small intestine.¹⁸ Efflux transporters the multidrug and toxin extrusions (MATE1 and MATE2-K) were MH substrate as the MATE1 expression is mainly in the liver and kidney cells membrane. The kidney cell’s membrane is the prominent place of MATE2-K.¹⁹ In addition, Plasma Membrane Monoamine Transporter (PMAT) identified affinity for MH uptake and expressed in the human small intestine. All the km and Vmax values of the transporters are presented in Table 3.

In the prediction process after model building, many inputs were added, starting with beads in vitro release results were entered with a change of gastric physiological emptying time to the floating beads time, and the (controlled release) CR gastric was chosen as the input of the dosage form.

The predicted model verification was by using the error percentage equation^20,21:

Beads formulation

The ionotropic gelation method was effective to prepare the multiparticulate bead with all conditions that used such as stirring rate, temperature of preparation and the percentage of calcium chloride solution. All the formulations in Table 2 successfully and visually showed beads and collected to be stored for further investigations.

Percentage of yield

All beads were subjected to the yield percentage test to assess the effect of the bead’s content, primarily the Na alginate, and the results were shown in Table 4. The results found that the increase in Na alginate concentration (F6, F7, F1-F4) increased the percentage yield of produced beads. This finding was similar to Singh et al finding of floating microspheres of famotidine.²⁶ This means the increased Na alginate amount helped in better bead formulations. Except for F5 the percentage of yield was diminished, this may be related to the formation of a thicker solution that acts as a barrier to form beads easily. The same inference was observed by previous researcher Tønnesen HH.²⁷

Percentage of entrapment efficiency %

An important consideration when designing new bead formulations is how effectively they hold the drug. Table 4 shows that all formulations (F1-F7) showed acceptable entrapment efficiency percent %, referring that Na alginate polymer concentrations are satisfactory for enclosing the drug in bead formulations, similar finding was noticed in Abbas et al work of enalapril floating microspheres as a controlled release dosage form.²⁸ Despite the small differencs in entrapment efficiency values of all beads formulations, this might be attributed to the layers of Na algiante that were formed led to decrease the inboxed of MH.

Floating properties (buoyancy lag time and duration of buoyancy)

It was essential for the current study to investigate the floating property, which was one of the study goals. All the formulations of beads floated without delay except F6, which contains the lowest amount of the polymer (0.250 gm). There was a lag time of 30 min, however, collectively the formulas remained floated for the whole release study time. This finding is comparable with what was observed with floating alginate beads of curcumin.²⁹

In vitro drug release

The MH release from beads was required to compare the drug release profile from different bead formulations, and the results were illustrated in Figure 1. The MH release was higher with a lower Na alginate amount, as in F6 and F7. The gelation process is based on the formation of tight between the guluronic acid residues in Na alginate and CaCl_2, which might increase as the Na alginate concentration increases, thus resulting in prolonged drug release as was noticed with the rest of the prepared formulations, this observation was justified by Mandal et al either.³⁰ Also, it was clear from Figure 1 that the beads formulations F1, F2, F3, F4 and F5 exhibited very similar release profiles and around 80 % MH was released gradually after 12 hours, while F6 and F7 released the MH after 6 hr and 8 hr, respectively.

Figure 1. The in vitro release study in 0.1 N HCl pH 1.2 as each release profile represents an average triplicate ±SD for each bead formulation.

Morphological examination

The F4 (representative beads of F1, F2 and F3), F6 and F7 were selected for morphological examination and FTIR. The morphological investigation by FESEM gives an idea about the shape and the surfaces of the formed bead. Figure 2 reveals that the beads are rounded in shape, suggesting that the concentrations of the polymer were satisfactory for beads formation, the beads apparently exhibit smoother surfaces as the concentration of Na alginate increases as in F4 in comparison to F6 and F7; this may be due to the higher density of crosslinked polymer matrices at the surface. The inside of the beads as it appears in cross-sectional view didn’t show distinguished morphology, a similar finding was noticed by Wan-PingVoo and et al.³¹ Visually, beads diameter of F4, F6 and F7 was around 1mm, as exhibited in the first left column of Figure 2. The same value of the diameter (1mm) was obtained in a different study prepared alginate beads.³²

Figure 2. (A) Field emission scanning electron micrographs (FESEM) micrographs for surface view as the magnification bars are 1mm and 30 μm for F4, F6 and F7 bead formulations, (B) (FESEM) microimages for cross-sectional view, magnification bar is 5 μm for F4, F6 and F7 bead formulations.

Fourier transmittance infrared (FTIR)

To understand the interactions between molecules formulating beads, FTIR was applied. The FTIR spectrogram of MH, as presented in Figure 3(A), showed similar peaks to the previous study of MH spectrogram as the peaks at 3367 cm^-1, 3300 cm^-1 and 1580 cm^-1 associated with the distinct group of N-H asymmetric stretching, N-H symmetric stretching and N-H bending, respectively.³³ Also, the Na alginate spectrogram showed a peak at 1619 cm^-1, which is the position of the carbonyl group of Na alginate beads.³⁴ Figure 3(B) demonstrates beads F4, F6 and F7 spectrograms emphasizing a peak deviation associated with the carbonyl group to 1629 cm^-1 as this might refer to Na alginate molecules bindings with calcium ions that helped to arrange and build the beads.³⁵ The N-H asymmetric stretching, N-H symmetric at 3367 cm⁻¹ and 3300 cm⁻¹ overlaid with OH region of the hydroxyl group of sodium alginate as it hard to tell if physical interaction happened in the spectrograms of F4, F6 and F7 related to the mentioned groups; however, MH might physically interact at N-H group as there is a shift of the band at 1580 cm⁻¹ to a lower region. Although it is not inconclusive, it can indicate the hydrogen bondings.

Figure 3. The whole FTIR spectrograms in Figure A representing the full range spectrograms for F4, F6, F7, pure drug, Na alginate, and CaCO₃; Figure B indicates carbonyl group shifting in F4, F6, and F7.

Metformin absorption in silico model

The parameters as an outcome of the in silico simulation that assisted in building a physiological model relied on the physicochemical and pharmacokinetic information are presented in Table 5. Also, the constructed model parameters Cmax, Tmax, AUC 0-inf, and AUC 0-t were validated, as screened in Figure 4, depending on the error percentage of observed and calculated data. The acceptance of %PE is valid as the calculated values do not double the observed values or the fold error value is not doubled.³⁶ The model showed a very acceptable error in percentage, as shown clearly in Table 5.

Figure 4. 500 mg oral immediate tablet physiological model as the bold line represents observed and the dotted line represents the calculated values.

Moreover, Figure 5A demonstrates that absorption does not take place in the stomach, whereas the maximum MH absorption occurred in jejeunum 1, then followed by the duodenum and jejeunum 2, and then almost nowhere else. Similarly, as the prediction results revealed, MH multiple unit floating gastroretentive drug delivery system exhibited no stomach absorption site.¹¹ Also, the fraction absorbed (total availability) was 23.9, of Metformin IR of the model, as seen in Figure 5A, which was close to 27%, the predicted fraction absorbed that was found by the Dahan study.³⁷ Consequently, Figure 5B indicates that F1 to F5 beads formulas provided high stomach amount following their floating duration and showed gradual with slight differences decrease in MH amount within 12 hrs representing its floating duration, while F7 showed a drastic decrease in the stomach amount of MH within the first hour of simulation time. Additionally, F6 presented a gradual reduction in the amount of MH in the stomach within simulation time. These outcomes referred that these gastroretentive multiple unit floating beads play an essential role in restricting the MH release rate in the site-specific region, which in turn guarantees the release of MH into the appropriate absorption site, thus, may improve the bioavailability.

Figure 5. Predicting data, Figure A The regional absorption of the 500 mg MH selected gastroretentive multiple unit floating beads; Figure B the MH amount retained in the stomach for 12 h as IR represents the immediate release tablet and blue curve denotes F1, F2, F3, F4 and F5 in Figure B.

In the same figure, for comparison purposes, immediate-release kinetic was taken from an immediate-release tablet of 500 mg of MH that showed (rapid MH decline).¹⁷ The simulating curves of plasma concentration-time as in Figure 6 parameters reveal that F1-F5 exhibit low Cmax with non-declining curves, thus referring to the gastroretentive property of slow-release MH multiple unit floating beads. At the same time, F6 and F7 showed a higher Cmax, as F7 showed a decline in the curve; this might be attributed to the total released amount of MH in the stomach. However, the floating period was persistent during the in vitro release study. Interestingly, the US Food and Drug Administration revealed the pharmacokinetics parameters of a 500 mg extended-release tablet as its Cmax was 0.6 μg/ml, which was close to the Cmax of 0.449 μg/ml of F7.³⁸ Taking this into account, the MH simulation helped to decide the better formulation of floating beads that achieved the aim of this study.

Figure 6. Predicted plasma concentration-time curves of MH multiple unit beads.