Human farnesyl pyrophosphate synthase (hFPPS) is a well-settled therapeutic target and it is an enzyme of the mevalonate pathway which catalyzes the biosynthesis of the C-15 isoprenoid farnesyl pyrophosphate. QSAR studies by using Monte Carlo method for human farnesyl pyrophosphate synthase inhibitors has been carried out using balance of correlation technique with Index of ideality correlation. For construction of QSAR models, six random splits were prepared from the data of 73 phosphonates and hybrid optimal descriptors procured from graph (HFG) and SMILES based notations were employed. The developed QSAR models have robustness, good fitting ability, generalizability and internal predictive ability. The external predictive ability has been certified by testing various precedents. The values of R2, IIC, Q2 and ∆R2 m for the best model are 0.9304, 0.9614, 0.9061 and 0.0861 respectively. The developed QSAR models met with the specified standards given in OECD guideline and applicability domain. The structural feature promoters for the end point increase and promoters for end point decrease have been extracted. The predicted pIC50 for the new proposed compounds have also been reported.