Abstract
Background: Three-dimensional printing (3DP) is a promising technology for advancing pharmaceutical research by enabling the production of personalized drug products. However, its progress has been hindered by the conventional trial-and-error approach to excipient selection and optimization. Methods: In this study, the blend module was employed to determine the miscibility parameters-mixing energy (E(mix)) and Flory-Huggins interaction parameter (χ) to find the right excipients and drug-excipient ratio and examine the incorporation of plasticizers and lipids to enhance printability. Furthermore, molecular dynamics (MD) simulations were employed to calculate the cohesive energy density (CED) for predicting the dissolution behavior of 3DP formulations. Results: Data from 51 formulations were analyzed, enabling correlation and experimental validation of the in silico predictions. The predicted miscibility values demonstrated a strong correlation with experimental printability results. Furthermore, using a miscibility parameter, it was possible to accurately forecast minor changes in the drug-to-excipient ratio, plasticizer/lipid concentration, and hot-melt extrusion (HME) temperature that affect printability. Hydrophilic carriers exhibited lower CED values corresponding to faster drug release. In contrast, more hydrophobic carriers revealed high CED values, indicating stronger drug entrapment and sustained release. Conclusions: The miscibility parameters and MD-simulated CED values provide a practical framework for early-stage, high-throughput excipient screening. Overall, in silico prediction offers a viable strategy for modeling the entire 3DP workflow, minimizing the need for trial-and-error experimentation, and thereby accelerating the clinical translation of 3DP drug delivery systems.