Human cytochrome P450 enzymes are membrane-embedded monooxygenases responsible for xenobiotic metabolism, steroidogenesis, fatty acid metabolism, and vitamin metabolism. Their active sites can accommodate diverse small molecules and understanding these interactions is key to decoding enzymatic functionality and designing drugs. The most common method for characterizing small molecule binding is quantifying absorbance changes that typically occur when ligands enter the active site near the heme iron. Traditionally, such titrations are monitored by a spectrophotometer, requiring significant manual time, protein, and increasing solvents. This assay was adapted for semi-automated high throughput screening, increasing throughput 50-fold while requiring less protein and keeping solvent concentrations constant. This 384-well assay was validated for both type I and II shifts typically observed for substrates and heme-coordinating inhibitors, respectively. This assay was used to screen a library of ∼100 diverse imidazole-containing compounds which can coordinate with the heme iron if compatible with the overall active site. Three human cytochrome P450 enzymes were screened: drug-metabolizing CYP2A6 and CYP2D6 and sterol-metabolizing CYP8B1. Each bound different sets of imidazole compounds with varying K values, providing a unique binding fingerprint. As a final validation, the K values were used to generate pharmacophores to compare to experimental X-ray structures. Applications for the high-throughput assay include the following: 1) facilitating generation of pharmacophores for enzymes where structures are not available, 2) screening to identify ligands for P450 orphans, 3) screening for inhibitors of P450s drug targets, 4) screening potential new drugs to avoid and/or control P450 metabolism, and 5) efficient validation of computational ligand binding predictions.