Modelling and Simulation of the Drug Release from a Solid Dosage Form in the Human Ascending Colon: The Influence of Different Motility Patterns and Fluid Viscosities

Model-Informed Drug Development: In Silico Assessment of Drug Bioperformance following Oral and Percutaneous Administration

The pharmaceutical industry has faced significant changes in recent years, primarily influenced by regulatory standards, market competition, and the need to accelerate drug development. Model-informed drug development (MIDD) leverages quantitative computational models to facilitate decision-making processes. This approach sheds light on the complex interplay between the influence of a drug's performance and the resulting clinical outcomes. This comprehensive review aims to explain the mechanisms that control the dissolution and/or release of drugs and their subsequent permeation through biological membranes. Furthermore, the importance of simulating these processes through a variety of in silico models is emphasized. Advanced compartmental absorption models provide an analytical framework to understand the kinetics of transit, dissolution, and absorption associated with orally administered drugs. In contrast, for topical and transdermal drug delivery systems, the prediction of drug permeation is predominantly based on quantitative structure-permeation relationships and molecular dynamics simulations. This review describes a variety of modeling strategies, ranging from mechanistic to empirical equations, and highlights the growing importance of state-of-the-art tools such as artificial intelligence, as well as advanced imaging and spectroscopic techniques.

Modelling the Evolution of Pore Structure during the Disintegration of Pharmaceutical Tablets

Pharmaceutical tablet disintegration is a critical process for dissolving and enabling the absorption of the drug substance into the blood stream. The tablet disintegration process consists of multiple connected and interdependent mechanisms: liquid penetration, swelling, dissolution, and break-up. One key dependence is the dynamic change of the pore space in a tablet caused by the swelling of particles while the tablet takes up liquid. This study analysed the changes in the pore structure during disintegration by coupling the discrete element method (DEM) with a single-particle swelling model and experimental liquid penetration data from terahertz-pulsed imaging (TPI). The coupled model is demonstrated and validated for pure microcrystalline cellulose (MCC) tablets across three porosities (10, 15, and 22%) and MCC with three different concentrations of croscarmellose sodium (CCS) (2, 5, and 8% w/w). The model was validated using experimental tablet swelling from TPI. The model captured the difference in the swelling behaviour of tablets with different porosities and formulations well. Both the experimental and modelling results showed that the swelling was lowest (i.e., time to reach the maximum normalised swelling capacity) for tablets with the highest CCS concentration, cCCS = 8%. The simulations revealed that this was caused by the closure of the pores in both the wetted volume and dry volume of the tablet. The closure of the pores hinders the liquid from accessing other particles and slows down the overall swelling process. This study provides new insights into the changes in the pore space during disintegration, which is crucial to better understand the impact of porosity and formulations on the performance of tablets.

Development of a digital twin of a tablet that mimics a real solid dosage form: Differences in the dissolution profile in conventional mini-USP II and a biorelevant colon model

The performance of colon-targeted solid dosage forms is commonly assessed using standardised pharmacopeial dissolution apparatuses like the USP II or the miniaturised replica, the mini-USP II. However, these fail to replicate the hydrodynamics and shear stresses in the colonic environment, which is crucial for the tablet's drug release process. In this work, computer simulations are used to create a digital twin of a dissolution apparatus and to develop a method to create a digital twin of a tablet that behaves realistically. These models are used to investigate the drug release profiles and shear rates acting on a tablet at different paddle speeds in the mini-USP II and biorelevant colon models to understand how the mini-USP II can be operated to achieve more realistic (i.e., in vivo) hydrodynamic conditions. The behaviour of the tablet and the motility patterns used in the simulations are derived from experimental and in vivo data, respectively, to obtain profound insights into the tablet's disintegration/drug release processes. We recommend an "on-off" operating mode in the mini-USP II to generate shear rate peaks, which would better reflect the in vivo conditions of the human colon instead of constant paddle speed.

The Effect of Biorelevant Hydrodynamic Conditions on Drug Dissolution from Extended-Release Tablets in the Dynamic Colon Model

The in vitro release of theophylline from an extended-release dosage form was studied under different hydrodynamic conditions in a United States Pharmacopoeial (USP) dissolution system II and a bespoke in vitro tubular model of the human colon, the Dynamic Colon Model (DCM). Five biorelevant motility patterns extracted from in vivo data were applied to the DCM, mimicking the human proximal colon under baseline conditions and following stimulation using polyethylene glycol or maltose; these represent the lower and upper bounds of motility normally expected in vivo. In the USPII, tablet dissolution was affected by changing hydrodynamic conditions at different agitation speeds of 25, 50 and 100 rpm. Applying different motility patterns in the DCM affected the dissolution profiles produced, with theophylline release at 24 h ranging from 56.74 ± 2.00% (baseline) to 96.74 ± 9.63% (maltose-stimulated). The concentration profiles of theophylline were markedly localized when measured at different segments of the DCM tube, highlighting the importance of a segmented lumen in intestine models and in generating spatial information to support simple temporal dissolution profiles. The results suggested that the shear stresses invoked by the unstimulated, healthy adult human colon may be lower than those in the USPII at 25 rpm and thus insufficient to achieve total release of a therapeutic compound from a hydroxyethyl cellulose matrix. When operated under stimulated conditions, drug release in the DCM was between that achieved at 25 and 50 rpm in the USPII.

Agglomeration Regimes of Particles under a Linear Laminar Flow: A Numerical Study

In this work, a combined smoothed particle hydrodynamics and discrete element method (SPH-DEM) model was proposed to model particle agglomeration in a shear flow. The fluid was modeled with the SPH method and the solid particles with DEM. The system was governed by three fundamental dimensionless groups: the Reynolds number Re (1.5 ~ 150), which measured the effect of the hydrodynamics; the adhesion number Ad (6 × 10−5 ~ 6 × 10−3), which measured the inter-particle attraction; and the solid fraction α, which measured the concentration of particles. Based on these three dimensionless groups, several agglomeration regimes were found. Within these regimes, the aggregates could have different sizes and shapes that went from long thread-like structures to compact spheroids. The effect of the particle–particle interaction model was also investigated. The results were combined into ‘agglomeration maps’ that allowed for a quick determination of the agglomerate type once α, Re, Ad were known.

Simulating the Hydrodynamic Conditions of the Human Ascending Colon: A Digital Twin of the Dynamic Colon Model

The performance of solid oral dosage forms targeting the colon is typically evaluated using standardised pharmacopeial dissolution apparatuses. However, these fail to replicate colonic hydrodynamics. This study develops a digital twin of the Dynamic Colon Model; a physiologically representative in vitro model of the human proximal colon. Magnetic resonance imaging of the Dynamic Colon Model verified that the digital twin robustly replicated flow patterns under different physiological conditions (media viscosity, volume, and peristaltic wave speed). During local contractile activity, antegrade flows of 0.06-0.78 cm s^-1 and backflows of -2.16--0.21 cm s^-1 were measured. Mean wall shear rates were strongly time and viscosity dependent although peaks were measured between 3.05-10.12 s^-1 and 5.11-20.34 s^-1 in the Dynamic Colon Model and its digital twin respectively, comparable to previous estimates of the USPII with paddle speeds of 25 and 50 rpm. It is recommended that viscosity and shear rates are considered when designing future dissolution test methodologies for colon-targeted formulations. In the USPII, paddle speeds >50 rpm may not recreate physiologically relevant shear rates. These findings demonstrate how the combination of biorelevant in vitro and in silico models can provide new insights for dissolution testing beyond established pharmacopeial methods.

Luminal Fluid Motion Inside an In Vitro Dissolution Model of the Human Ascending Colon Assessed Using Magnetic Resonance Imaging

Knowledge of luminal flow inside the human colon remains elusive, despite its importance for the design of new colon-targeted drug delivery systems and physiologically relevant in silico models of dissolution mechanics within the colon. This study uses magnetic resonance imaging (MRI) techniques to visualise, measure and differentiate between different motility patterns within an anatomically representative in vitro dissolution model of the human ascending colon: the dynamic colon model (DCM). The segmented architecture and peristalsis-like contractile activity of the DCM generated flow profiles that were distinct from compendial dissolution apparatuses. MRI enabled different motility patterns to be classified by the degree of mixing-related motion using a new tagging method. Different media viscosities could also be differentiated, which is important for an understanding of colonic pathophysiology, the conditions that a colon-targeted dosage form may be subjected to and the effectiveness of treatments. The tagged MRI data showed that the DCM effectively mimicked wall motion, luminal flow patterns and the velocities of the contents of the human ascending colon. Accurate reproduction of in vivo hydrodynamics is an essential capability for a biorelevant mechanical model of the colon to make it suitable for in vitro data generation for in vitro in vivo evaluation (IVIVE) or in vitro in vivo correlation (IVIVC). This work illustrates how the DCM provides new insight into how motion of the colonic walls may control luminal hydrodynamics, driving erosion of a dosage form and subsequent drug release, compared to traditional pharmacopeial methods.