SUMMARYAn integrated approach using targeted metabolite profiles and modest EST libraries each containing approximately 3500 unigenes was developed in order to discover and functionally characterize novel genes involved in plant-specialized metabolism. EST databases have been established for benzylisoquinoline alkaloid-producing cell cultures of Eschscholzia californica, Papaver bracteatum and Thalictrum flavum, and are a rich repository of alkaloid biosynthetic genes. ESI-FTICR-MS and ESI-MS/MS analyses facilitated unambiguous identification and relative quantification of the alkaloids in each system. Manual integration of known and candidate biosynthetic genes in each EST library with benzylisoquinoline alkaloid biosynthetic networks assembled from empirical metabolite profiles allowed identification and functional characterization of four N-methyltransferases (NMTs). One cDNA from T. flavum encoded pavine N-methyltransferase (TfPavNMT), which showed a unique preference for (AE)-pavine and represents the first isolated enzyme involved in the pavine alkaloid branch pathway. Correlation of the occurrence of specific alkaloids, the complement of ESTs encoding known benzylisoquinoline alkaloid biosynthetic genes and the differential substrate range of characterized NMTs demonstrated the feasibility of bilaterally predicting enzyme function and species-dependent specialized metabolite profiles.