The present invention provides a set of in vitro and in silico methodologies for predicting in vivo pharmacokinetics and pharmacodynamics of multiple components the methodologies comprise mathematical models for solving multiple unknowns which are linearly independent and/or interacting with each other. The present invention can be applied to develop phytomedicines which contain multiple active ingredients without prior identification, isolation and purification of these components.