Correlation between compactibility values and excipient cluster size using an in silico approach

In silico simulation and percolation theory are important tools in the study of physical and mechanical behavior of pharmaceutical compacts. The aim was to generate a new in silico simulation program that describes the mechanical structure of binary compacts formed from an excipient with excellent compactibility and a drug with null compactibility.

Paracetamol and microcrystalline cellulose powders were compressed under different pressures. Values for the indentation hardness and tensile strength were measured and fitted to the Leuenberger's model. On the other hand, compacts with different composition were in silico simulated. In each system, the biggest excipient cluster was identified and quantified using the Hoshen-Kopelman algorithm. Then, the size of the biggest in silico cluster was correlated with experimental compactibility values.

The Leuenberger's model resulted in good fit to the experimental data for all formulations over 40% of excipient load. Formulations with high drug load (≥0.8) had reduced range for forming compacts and gave low compactibility values. The excipient percolation threshold for the simulated system was 0.3395, indicating that over this excipient fraction, a compact with defined mechanical properties will be formed. The compactibility values presented a change in the range of 0.3-0.4 of excipient fraction load, just where the in silico excipient percolation threshold was found.

Physical measurements of the binary compacts showed good agreement with computational measurements. Subsequently, this in silico approach may be used for the optimization of pharmaceutical powder formulations used in tablet compression.