The expression of cytochrome P450 (CYP) enzymes is highly variable and associated with factors, such as age, genotype, sex, and disease states. In this study, quantification of metronidazole metabolizing CYP isoforms (CYP2A6, CYP2E1, CYP3A4, CYP3A5, and CYP3A7) in human liver microsomes from 115 children and 35 adults was performed using a quantitative proteomics method. The data confirmed age-dependent increase in CYP2A6, CYP2E1, and CYP3A4 abundance, whereas, as expected, CYP3A7 abundance showed postnatal decrease with age. In particular, the fold difference (neonatal to adulthood levels) in the protein abundance of CYP2A6, CYP2E1, and CYP3A4 was 14, 11, and 20, respectively. In contrast, protein abundance of CYP3A7 was > 125-fold higher in the liver microsomes of neonates than of adults. The abundance of CYP2A6 and CYP3A5 was associated with genotypes, rs4803381 and rs776746, respectively. A proteomics-informed physiologically based pharmacokinetic (PBPK) model was developed to describe the pharmacokinetics of metronidazole and its primary metabolite, 2-hydroxymethylmetronidazole. The model revealed an increase in the metabolite-to-parent ratio with age and showed a strong correlation between CYP2A6 abundance and metabolite formation (r 2 = 0.75). Notably, the estimated contribution of CYP3A7 was ~ 75% in metronidazole clearance in neonates. These data suggest that variability in CYP2A6 and CYP3A7 in younger children poses the risk of variable pharmacokinetics of metronidazole and its active metabolite with a potential impact on drug efficacy and safety. No sex-dependent difference was observed in the protein abundance of the studied CYPs. The successful integration of hepatic CYP ontogeny data derived from a large liver bank into the pediatric PBPK model of metronidazole can be extended to other drugs metabolized by the studied CYPs.