Aromatase is an important pharmacological target in the anti-cancer therapy as the intratumoral aromatase is the source of local estrogen production in breast cancer tissues. Suppression of estrogen biosynthesis by aromatase inhibition represents an effective approach for the treatment of hormone-sensitive breast cancer. Because of the membrane-bound character and heme-binding instability, no crystal structure of aromatase was reported for a long time, until recently when crystal structure of human placental aromatase cytochrome P450 in complex with androstenedione was deposited in PDB. The present study is towards understanding the structural and functional characteristics of aromatase to address unsolved mysteries about this enzyme and elucidate the exact mode of binding of aromatase inhibitors. We have performed molecular docking simulation with twelve different inhibitors (ligands), which includes four FDA approved drugs; two flavonoids; three herbal compounds and three compounds having biphenyl motif with known IC(50) values into the active site of the human aromatase enzyme. All ligands showed favorable interactions and most of them seemed to interact to hydrophobic amino acids Ile133, Phe134, Phe221, Trp224, Ala306, Val370, Val373, Met374 and Leu477 and hydrophilic Arg115 and neutral Thr310 residues. The elucidation of the actual structure-function relationship of aromatase and the exact binding mode described in this study will be of significant interest as its inhibitors have shown great promise in fighting breast cancer.