Current formulation development strongly relies on trial-and-error experiments in the laboratory by pharmaceutical scientists, which is time-consuming, high cost and waste materials. This research aims to integrate various computational tools, including machine learning, molecular dynamic simulation and physiologically based absorption modeling (PBAM), to enhance andrographolide (AG) /cyclodextrins (CDs) formulation design. The lightGBM prediction model we built before was utilized to predict AG/CDs inclusion's binding free energy. AG/γ-CD inclusion complexes showed the strongest binding affinity, which was experimentally validated by the phase solubility study. The molecular dynamic simulation was used to investigate the inclusion mechanism between AG and γ-CD, which was experimentally characterized by DSC, FTIR and NMR techniques. PBAM was applied to simulate the in vivo behavior of the formulations, which were validated by cell and animal experiments. Cell experiments revealed that the presence of D-α-Tocopherol polyethylene glycol succinate (TPGS) significantly increased the intracellular uptake of AG in MDCK-MDR1 cells and the absorptive transport of AG in MDCK-MDR1 monolayers. The relative bioavailability of the AG-CD-TPGS ternary system in rats was increased to 2.6-fold and 1.59-fold compared with crude AG and commercial dropping pills, respectively. In conclusion, this is the first time to integrate various computational tools to develop a new AG-CD-TPGS ternary formulation with significant improvement of aqueous solubility, dissolution rate and bioavailability. The integrated computational tool is a novel and robust methodology to facilitate pharmaceutical formulation design.