Heat shock proteins (Hsps) 60 and 70 are postulated as a potential drug target for toxoplasmosis due to its importance in the developmental and survival of Toxoplasma gondii (T. gondii). As of today, there have been no reports on three-dimensional (3D) structure of Hsp60 and Hsp70 deposited in the Brookhaven Protein Data Bank. Hence, this study was conducted to predict 3D structures for Hsp60 and Hsp70 in T. gondii by homology modeling. Selection of the best predicted model was done based on multiple scoring functions. In addition, virtual screening was performed to short-list chemical compounds from the National Cancer Institute (NCI) Diversity Set III in search of potential inhibitor against Hsp60 and Hsp70 in T. gondii. Prior to virtual screening, binding sites of Hsp60 and Hsp70 were predicted using various servers and were used as the center in docking studies. The Hsps were docked against known natural ligands to validate the method used in estimating free energy of binding (FEB) and possible interactions between ligand and protein. Virtual screening was performed with a total of 1560 compounds from the NCI Diversity Set III. The compounds were ranked subsequently according to their FEB. Molecular basis of interactions of the top five ranked compounds was investigated using Ligplot. The major interactions exhibited were hydrogen bonding and hydrophobic interactions in binding to Hsp60 and Hsp70. The results obtained provided information and guidelines for the development of inhibitors for Hsp60 and Hsp70 in T. gondii.