Molecular modeling techniques were used to derive a predictive model for substrates of cytochrome P450 2D6, an isozyme known to metabolize only compounds with one or more basic nitrogen atoms. Sixteen substrates, accounting for 23 metabolic reactions, with a distance of either 5 A ("5-A substrates", e.g., debrisoquine) or 7 A ("7-A substrates", e.g., dextromethorphan) between oxidation site and basic nitrogen atom were fitted into one model by postulating an interaction of the basic nitrogen atom with a negatively charged carboxylate group on the protein. This acidic residue anchors and neutralizes the positively charged basic nitrogen atom of the substrates. In case of "5-A substrates" this interaction probably occurs with the carboxylic oxygen atom nearest to the oxidation site, whereas in the case of "7-A substrates" this interaction takes place at the other oxygen atom. Furthermore, all substrates exhibit a coplanar conformation near the oxidation site and have negative molecular electrostatic potentials (MEPs) in a part of this planar domain approximately 3 A away from the oxidation site. No common features were found in the neighbourhood of the basic nitrogen atom of the substrates studied so that this region of the active site can accommodate a variety of N-substituents. Therefore, the substrate specificity of P450 2D6 most likely is determined by the distance between oxidation site and basic nitrogen atom, by steric constraints near the oxidation site, and by the degree of complementarity between the MEPs of substrate and protein in the planar region adjacent to the oxidation site.(ABSTRACT TRUNCATED AT 250 WORDS)