Fragment molecular orbital (FMO) method provides a novel tool for ab initio calculations of large biomolecules. This method overcomes the size limitation difficulties in conventional molecular orbital methods and has several advantages compared to classical force field approaches. While there are many features in this method, we here focus on explaining the issues related to protein-ligand binding: FMO method provides useful interaction-analysis tools such as IFIE, CAFI and FILM. FMO calculations can provide not only binding energies, which are well correlated with experimental binding affinity, but also QSAR descriptors. In addition, FMO-derived charges improve the descriptions of electrostatic properties and the correlations between docking scores and experimental binding affinities. These calculations can be performed by the ABINIT-MPX program and the calculation results can be visualized by its proper BioStation Viewer. The acceleration of FMO calculations on various computer facilities is ongoing, and we are also developing methods to deal with cytochrome P450, which belongs to the family of drug metabolic enzymes.