3.4. Segmental DCM Dissolution Profiles
Standard drug dissolution/release profiles in the pharmaceutical industry are presented as a single profile with no spatial information, as in
Figure 3. However, the lumen of the DCM, much like the intestine in vivo, is not perfectly mixed. This means that information to describe how drug concentration varies with location is key to understanding how much API has been released from the tablet and how effective the motility pattern is at distributing the dissolved compound. Samples were withdrawn from five locations along the DCM to obtain the distribution of dissolved API along the DCM tube over 24 h.
Figure 4 presents the concentration profiles of dissolved theophylline measured at S1, S3, S5, S7 and S10, progressively further from the point of tablet insertion, whilst the respective inset figure shows the location of the sampled segment.
Besides from the DCM, the architecture of most advanced in vitro models of the colon takes the form of a single tube [
5]. The most advanced physiologically based pharmacokinetic (PBPK) models use oversimplified first order transit rate models that consider the colon to be a single homogeneously mixed compartment. It is well-documented that the colon contracts in a segmental fashion wherein there is likely to be regional differences in both mixing behaviour and thus drug concentration. This has been demonstrated using the DCM in this study; where mimicked in vivo motility involved frequent contractions (to different degrees of occlusion) whilst other segments remained stationary. The spatiotemporal information in
Figure 4 demonstrates clearly that a drug concentration profile is highly dependent on the segments of the DCM lumen from which measurements are taken (
Figure 4). If the DCM were a perfectly mixed system each segment should yield 2 mg mL
−1 of theophylline. This demonstrates the importance of a segmented lumen in models of the colon, thus, there is a need to update the PBPK models to include a segmented colon. Furthermore, the segmental and cumulative dissolution data from the DCM could be integrated into PBPK models of the colon to better predict overall exposure; particularly for erodible formulations where motility is critical for drug release. This is relevant as it was observed that operating motility patterns with a low MI
DCM can cause the contents of some segments to be somewhat stagnant. Operating patterns with such low motility may therefore encounter potential solubility effects on the rate of release measured in segments containing the tablet body, especially considering vehicles with a poorly soluble payload.
S1 was the point of tablet insertion.
Figure 4A presents the concentration profiles of theophylline measured at S1, representing measured dissolution in the fluid immediately surrounding the tablet, which did not move from its insertion point for the duration of the experiments. The highest concentration reached in S1 was under static (2.82 ± 0.50 mg mL
−1) and baseline (3.07 ± 0.21 mg mL
−1) conditions. When the DCM was static, drug release from the dosage form and distribution along the DCM was driven primarily by diffusion. This led to accumulation around the point of tablet insertion and high local drug concentration. Mass transport in the static model was not purely diffusive however, since the media sampling and replenishment process introduced advective mixing, although this was minimised by operating sampling/replenishment syringe pumps at low flow rates (6 mL min
−1). Under baseline conditions, there were single, isolated contractions in S3, S6, S5 and S8, in chronological order. Additionally, there was a propagating pressure wave (PPW) from S2 to S4. This baseline contractile activity had a low MI
DCM, involving low occlusion degrees (25%) and slow occlusion rates (1.5 mm s
−1) compared to the other patterns (refer to
Table 2). Under such lethargic wall motion, it is unlikely that sufficient energy was imparted to the fluid to mix the area of high drug accumulation in S1, with the segments of the DCM closer to the hepatic flexure via advective transport through the viscous dissolution media. This was proven by the very low concentrations reached in S5, S7 and S10 (
Figure 4C–E), where the baseline and static concentration profiles were very similar and significantly lower than in all other motility conditions.
Furthermore, it is possible that the baseline contractions promoted back-mixing rather than antegrade propagation of the contents. However, it is evident that the rate of increase of concentration in S1 was higher under baseline conditions compared to static. This suggests that the local mixing caused by proximate wall motion; the PPW from S2 to S4 and the isolated contraction at S3, was sufficient to either impart sufficient shear stress to erode the outer layers of the dosage form or mix the contents within S1. To elaborate on the latter, it was assumed that the concentration within a segment was homogeneous. However, the robustness of this assumption is proportional to the mixing of the system; in a perfectly mixed system it holds true whilst in a stagnant system, such as the static DCM, it is less valid. The baseline pattern may have been sufficient to disturb the fluid surrounding the tablet and mix the contents of S1 such that the sample withdrawn from SP1 gave a more accurate reflection of how much drug had dissolved from the tablet. It is likely that there was a steep concentration gradient radiating from the tablet body outwards in these patterns with low MIDCM, exacerbated by the viscosity of the dissolution media used. This could mean that the highest concentrations in S1 under the low MIDCM static and baseline conditions were not reflected by the concentration measurements taken from SP1.
Although the static and baseline patterns showed the highest concentrations of dissolved theophylline in S1 (
Figure 4), these data also exhibited the highest variability. This could be due to a limitation of the UV-VIS methodology used where samples with a higher concentration were diluted to a higher extent, which increases the scale of any initial error due to an increased number of measurements required using the pipette.
Theophylline concentrations in S1 were generally lowest under the CPPWs, reaching 1.31 ± 0.06 and 1.41 ± 0.08 mg mL
−1 at 24 h under the antegrade and retrograde waves, respectively.
Figure 4 clearly shows a constant increase in concentration at S3, S5, S7 and S10 for the antegrade and retrograde CPPWs. After only 1 h, concentration of theophylline measured at S5, S7 and S10 was significantly higher than for the baseline pattern for both CPPWs. This clearly demonstrates that these patterns were effective at mixing along the length of the tube, significantly more so than the baseline motility pattern. Concentrations at S1 after 24 h were not significantly different between the antegrade and retrograde CPPWs. Interestingly however, there were significant differences in the course taken by the concentration profiles. The rate of increase in theophylline concentration measured at S1 was much higher for the antegrade pattern until t = 10 h, where a plateau was reached whilst the concentration profile generated by the retrograde CPPW followed the same trajectory (at a lower rate of increase) as the stimulated patterns until t = 8 h. This suggests that the antegrade wave may not have been as effective as the retrograde wave in driving advective transport of dissolved drug along the DCM tube, away from the tablet lying in S1. This is reinforced by the consistently lower mean drug concentrations measured at S5, S7 and S10, however the differences were not statistically significant (
p < 0.05). Additionally, this hypothesis is supported by a previous MRI study of the DCM that showed that the highest magnitude velocities generated by an antegrade CPPW were backflows (retrograde) rather than in the direction of wave propagation [
11].
The rate of increase of drug concentration at S1 followed similar trajectory to the retrograde CPPW, lower than the static, baseline and antegrade patterns, until t = 8 h. At this point, the concentration continued to increase at an almost constant rate until the peak at 24 h. The concentration at S3, S5, S7 and S10 also constantly increases over the duration of the experiment, consistently reaching higher values than in the CPPWs. This infers that in addition to extensive mixing capabilities, the hydrodynamics caused by the higher MIDCM stimulant-driven patterns generated a higher erosive effect than the CPPWs, constantly releasing and dissolving drug from the tablet body and mixing it throughout the DCM.
At S1 after 24 h, the concentration of theophylline measured under application of the stimulated patterns lied in between the static and baseline patterns and the CPPWs, at 1.96 ± 0.11 and 2.07 ± 0.14 for PEG and maltose, respectively. These values are close to 2 mg mL−1 which would be expected for a homogenously mixed lumen and complete drug release. At S10 however, drug concentration was lower, reaching 1.43 ± 0.13 and 1.53 ± 0.21, therefore this situation was not realised. The highest release measured at S7 and S10 was from the Maltose stimulated pattern, closely followed by the PEG stimulated pattern. This demonstrates that the stimulated motility patterns were effective at distributing dissolved drug throughout the DCM. In these patterns with the highest MIDCM, wave propagation speed and frequency were higher, though wave propagation distance was shorter, travelling a maximum of 4 segments compared to 10 in the CPPWs. This proves that inside the DCM it is not necessary for a contractile wave to propagate the length of the system (such as in the CPPWs) to effectively transport and mix fluid contents across the entire length of the tube.
The concentration profiles generated by the stimulated motility patterns at S7 and S10 increased with time in a linear modality, as did that of the baseline and CPPW patterns. However, in the maltose stimulated motility pattern there was a discrepancy in the linear trend at t = 6 h. In the DCM, this type of artefact may be caused by dispersed tablet material that have become detached from the principal tablet body and were local to the sample point when samples were taken. The likelihood of this occurring is bound to increase with the eroding capabilities of the motility pattern. This type of discrepancy is likely to be more prevalent in the DCM than in typical USPII apparatus, due to a higher sample point to media volume ratio.
3.5. Cumulative Temporal Theophylline Release Profiles in the DCM
The cumulative dissolution profile generated under each motility pattern is shown in
Figure 5. This profile was built from the segmental concentration measurements, from which the distribution of the dissolved drug throughout the model could be estimated to evaluate the total mass of dissolved theophylline in the DCM lumen. The shaded regions represent the 95% confidence intervals for the power law models fitted to the experimental data, further explained in
Table 4.
In the early stages of dissolution where t < 3 h, all profiles were similar to the static and baseline patterns, suggesting that diffusive mass transport governed release from the tablet. Subsequently, the importance of erosion becomes evident as the dissolution profiles under mimic stimulated conditions deviated significantly from the static and baseline waves. This is likely to align with hydration time of the HEC gel layer. From this point onwards, the hydrodynamics becomes important and the potential of the DCM to discriminate drug release profiles based upon physiologically relevant motility patterns is highly pertinent. The ability to recreate motility conditions demonstrates the power of the DCM to replicate inter- and intrasubject variability in motility, including extremes, and how this may influence colonic drug delivery. However, the authors recognise that in vivo the tablet would be subjected to the conditions of the upper GI tract prior to reaching the colon and that these aspects are not recreated in this study. A future study that included a biorelevant mimic of the conditions including the hydrodynamics of the upper GI tract that enabled delivery of the hydrated tablet without disruption of the swollen gel layer would be a great extension of this work.
The baseline motility was able to release 56.74 ± 2.00% of theophylline form the formulation over 24 h, significantly higher than achieved by the static DCM 49.98 ± 4.14%. However, the spatial information provided by
Figure 4 showed that it was highly ineffective at transporting dissolved API along the DCM tube to segments beyond S3. When directly compared, the dissolution profiles generated from the antegrade and retrograde CPPWs were not statistically different (
f1 = 5.5%,
f2 = 76.4%). However, when compared to the profile measured under static conditions, the antegrade CPPW was not significantly different to the static profile (
f1 = 19.4%,
f2 = 54.3%) whilst the retrograde CPPW was different (
f1 = 31.2%,
f2 = 48.0%). This shows that the advective motion generated by a retrograde CPPW are sufficient to elevate the release and distribution of the drug inside the DCM beyond simple diffusion. However, the mean release achieved by the retrograde wave was significantly higher than the antegrade wave at t = 24 h. Although the wave parameters of the antegrade and retrograde contractions were equal, the propagation distance of the wave and the location of the wave front relative to the tablet position are drastically different. This is likely to be the leading contributing factor to the significantly higher mean release at t = 24 h under the retrograde CPPW. Additionally, O’Farrell et al. [
5] showed the nature of a CPPW with 40% occlusion degree to cause high velocities in the opposite direction to wave propagation, which may be attributed to the increased rate of dissolved theophylline transport along the length of the DCM. Furthermore, in one experiment using the retrograde motility pattern, a large fragment of tablet had been eroded from the bulk tablet body and was located at the hepatic flexure from t = 16–24 h. This suggests that motility patterns which include more frequent retrograde propagating pressure waves may be more effective at achieving higher release of a therapeutic from a dosage form located towards the early stages of the proximal colon.
Figure 5 shows a decrease in rate of increase of theophylline concentration from t = 3 h in both the antegrade and retrograde CPPWs, compared to the stimulated profiles which maintain a steady rate of increase. This suggests that insufficient shear was generated in the CPPWs to maintain the release rate observed in the profiles generated by the stimulated patterns. The same conclusion can be drawn from the baseline pattern. After t = 16 h, the shear stresses generated by the baseline and CPPWs appeared to be insufficient to cause breakdown and release of the drug that lies in the core of the tablet body, since the stimulated motility patterns produced significantly higher release than all others. The maltose stimulated pattern performed best in achieving highest release of theophylline.
Comparing the PEG and maltose stimulated patterns, there was no significant difference between the release profiles obtained (
f1 = 6.00%,
f2 = 70.70%). However, the mean concentration was consistently higher for the maltose-stimulated pattern at most time points. Although the PEG-stimulated pattern had a higher motility index, different features of the pattern may have caused the lower release. For example, S2 was the closest segment to the tablet insertion point in both patterns. However, in the PEG-stimulated pattern, S2 only exhibited isolated contractions, compared to the maltose-stimulated pattern, where S2 was involved in propagating wave contractions. This might suggest that the flows of the contents of the lumen generated by propagating waves are more effective at eroding a dosage form in the local vicinity. Furthermore, it has been shown that the location of the contents with respect to contractile activity is vital in determining the velocities imparted to that fluid and therefore the shear rates experienced [
11]. Propagating pressure waves in the PEG-stimulated motility pattern involved S6 rather than S2 in the maltose-stimulated pattern. This may have caused lower peaks in flow around the tablet located at S1, since contractions in S6 are further from the tablet and more fluid lies in between the contracting segment and the tablet to dampen the flow.
The findings from this study align with those from the in silico study by Schutt et al.-that the mimicked motor patterns based on stimulants were significantly more effective at releasing a water-soluble drug from an erodible vehicle. Stimulated motility patterns had a higher MI
DCM and thus exhibited faster occlusion rates and a higher frequency of waves and single contractions, in addition to higher occlusion degrees. This means that there would be a higher number of peaks in shear stress per 120 s cycle; Schutt et al. [
16] predicted that drug release from a computationally simulated dosage form was significantly impacted by the number of shear stress peaks. This is likely to also be a contributing mechanism in vitro. Higher occlusion rates have been shown to generate higher pressures inside the DCM [
12], which (through conservation of energy) would cause higher velocities of the contents and subject a dosage form to more intense agitation and higher shear rates, achieving higher release. Comparisons between the hydrodynamics experienced by the tablet in this study and the in silico study are limited though, since the solid dosage form modelled by Schutt et al. [
16] was neutrally buoyant whereas the UC tablets sunk in the DCM fluid and had less surface exposed to high shear rates in the fluid. Overall, the findings from this study suggest the in silico model may have value for testing whether a particular motility pattern, combinations of patterns or changes in motility parameters have a significant effect on drug release. This could be conducted as a first step to prioritise in vitro trials.
As previously discussed, the significant difference in release resulting from the change in rpm of the USPII apparatus demonstrated by
Figure 3 shows that this formulation is shear-sensitive. However, the dissolution profiles inside the DCM in
Figure 5 show that extended-release formulations that reach the colon may need to have greater sensitivity to lower shear rates. If a significant portion of the API remains in the vehicle when a prolonged release formulation reaches the colon, the polymer gel layer may be too resistant to shear to release sufficient API for effective treatment. It may be beneficial for formulators to use in vitro tools that model different regions of the GI tract for different stages of release, for example an inner core that is more sensitive to shear after 5.3 h, when a formulation typically reaches the ascending colon (median colon arrival time in healthy adult humans, 95% confidence interval 4.51–5.48 h [
26]). From these findings, it may be more applicable run the USPII apparatus with a pulsed agitation, stepping between 0–25 rpm with intermediate steps representing shear rates expected in vivo. This is similar to the stress test device developed by Garbacz et al. [
25] based on physiological pressures a dosage form may experience during GI transit. However, precise information regarding in vivo shear rates is not currently available [
17]. Future work with the DCM should use a combination of baseline and stimulated patterns in the same dissolution experiment in order to more accurately reflect the range of motility a healthy human subject may experience day-to-day after ingesting a solid oral dosage form. It may also be possible to evaluate the difference of DCM orientation on baseline dissolution, by rotating the DCM through 0–90° at predetermined intervals to mimic the influence of gravity and human physical activity. This would likely cause a change in positioning or orientation of the tablet, exposing the side of the tablet that constantly in contact with the lower DCM wall in this study.
The knowledge gained from overall release profiles can also be used to further understand the segmental concentration profiles.
Figure 4 showed plateauing behaviour of the S1 baseline concentration profile from t = 6 h, suggesting a decrease in the rate of theophylline release. However, it is clear from
Figure 5 and
Table 4 that total release was not close to 100%. Therefore, the plateau could be due to the introduction of solubility effects as local concentration gradients between the gel layer of the tablet and the surrounding fluid in S1 were diminished and the fluid approached the aqueous solubility of theophylline (5.5 mg mL
−1). This would retard the rate of mass transport by diffusion. Another possibility is that diffusion from the outer layers of the tablet was complete and the rate of diffusion from the inner core was slow and crystalline theophylline was essentially ‘locked up’ and inaccessible without erosion of the outer layers, due to a low MI
DCM.
A power law model was fitted to the overall release profiles obtained in the DCM under different motility patterns (
Figure 5), with R
2 > 0.999 for all profiles, as shown in
Table 4. Similar to in the Korsmeyer-Peppas model, the exponent, m, generally increased with MI
DCM and is likely to be related to the frequency and intensity of shear experienced in an environment that exhibits colon-like pulsatile, peristaltic hydrodynamics. There is some discrepancy in this trend as the PEG stimulated pattern had a higher MI
DCM than the maltose stimulated pattern but a lower m and overall release. This suggests a contributing factor of the locality of the contractile activity with respect to the tablet, as previously discussed.
Table 5 presents the release comparison at t = 4, 10 and 24 h for each DCM motility pattern versus the USPII at 25, 50 and 100 rpm. Release at 4 h was generally higher in the USPII and discrepancy between the USPII and DCM increased with rpm, with 38.98–62.79% greater release under 100 rpm compared to the DCM. Overall discrepancy between the DCM and USPII was most extreme at t = 10 h, with the DCM showing >30% less release than the USPII when all motility patterns were compared at agitation speeds of 50 and 100 rpm. This indicates that the release rate in the USPII under these conditions was considerably higher between t = 4 and t = 10 h. After 24 h, the DCM had had sufficient time to reduce the gap between release achieved in the systems. The static, baseline and CPPWs never achieved higher release than the USPII at any agitation speed. At t = 4, 10 and 24 h, release was higher under mimic stimulated conditions in the DCM compared with the USPII operated at 25 rpm, however it is unrealistic that the human proximal colon would display this level of stimulated activity for a continuous period of 24 h. These findings suggest that operation of the USPII at 25 rpm and above generates hydrodynamics conditions that may be too intensive to be representative of unstimulated human colon. However, the stimulated colon may generate hydrodynamic conditions that lie somewhere between 25 and 50 rpm. This data builds upon previous conclusions that the use of > 50 rpm in the USPII is not recommended for colon-targeted dosage forms [
10].
Table 4 highlights that the baseline and CPPW patterns did not exceed 70% release over 24 h, this suggests that the HEC gel layer in UC tablets may be overly resistant to the shear stresses that may be present in the unstimulated colon. Higher release was achieved in the USPII at 25, 50 and 100 rpm. Although conclusions cannot be drawn on the shear rates likely to have occurred from this comparison, since a constant shear is applied in the USPII and shear in the DCM has a high spatiotemporal dependence [
10]. Previous MRI studies of the DCM found velocities to vary from −2.16–0.78 cm s
−1 and shear rate fluctuating between 0–8 s
−1 when slower CPPWs were applied (0.4 and 0.8 cm s
−1) at the same volume and viscosity as in this study.
The variability in the overall release profile in the DCM was high compared to the USPII, likely due to the increased number of steps included in data processing. However, there was similarly high variability in direct measurements, e.g., the concentrations presented in
Figure 4. This may better reflect the level of variation in vivo due to position of the tablet in relation to contractions and potential to experience elevated shear from being squeezed between contracted haustra, amongst many other factors.
Regarding use of the
f2 approach to classify profiles as statistically indifferent/different from one another in this study, caution must be taken around the level of significance used. This method is typically used to compare dissolution profiles between products, rather than of the same product using different dissolution methods [
22]. However, since the same formulation and dissolution media are used in all experiments and the dosage form is intended for extended release, it would be highly unusual for the product to exhibit >10% difference at the earlier stages of dissolution. Additionally, the majority of official guidance from governing bodies worldwide regarding acceptance criteria for
f2 comparisons are for immediate release formulations rather than extended-release formulations where time is integral for matrix hydration, gel formation and erosion to take place.