Identifying antimicrobials and their metabolites in wastewater and surface water with effect-directed analysis

This study aimed to identify antimicrobial contaminants in the aquatic environment with effect-directed analysis. Wastewater influent, effluent, and surface water (up- and downstream of the discharge location) were sampled at two study sites. The samples were enriched, subjected to high-resolution fractionation, and the resulting 80 fractions were tested in an antibiotics bioassay. The resulting bioactive fractions guided the suspect and nontargeted identification strategy in the high-resolution mass spectrometry data that was recorded in parallel. Chemical features were annotated with reference databases, assessed on annotation quality, and assigned identification confidence levels. To identify antibiotic metabolites, Phase I metabolites were predicted in silico for over 500 antibiotics and included as a suspect list. Predicted retention times and fragmentation patterns reduced the number of annotations to consider for confirmation testing. Overall, the bioactivity of three fractions could be explained by the identified antibiotics (clarithromycin and azithromycin) and an antibiotic metabolite (14-OH(R) clarithromycin), explaining 78% of the bioactivity measured at one study site. The applied identification strategy successfully identified antibiotic metabolites in the aquatic environment, emphasizing the need to include the toxic effects of bioactive metabolites in environmental risk assessments.

Benchtop NMR spectroscopy in the analysis of substandard and falsified medicines as well as illegal drugs

Substandard and falsified medical products may cause harm to patients and fail to treat the diseases or conditions for which they were intended. It is therefore required to have analytical methods available to assess medical product quality. Benchtop NMR spectroscopy provides a generic, inherently quantitative, analytical method capable of separating specific signals from those of a matrix. We have developed an analytical method for the analysis of active ingredients in pharmaceutical products and illegal drugs, based on benchtop NMR spectroscopy. Within its resolution limits, benchtop NMR spectroscopy is useful in determining the identity of the active ingredients in products containing acetaminophen, aspirin, caffeine, diclofenac, ibuprofen, naproxen, sildenafil, tadalafil and sibutramine, cocaine, and gamma hydroxybutyric acid, with a limit of detection of about 1 mg/mL. Furthermore, the content of the active ingredient can be determined with an error of 10%. Additionally, a chemometrics approach is shown to be useful to classify spectra in order to identify the active substances present in the sample, reducing the need for expert interpretation of the spectra acquired.

Fora fuelling the discovery of fortified dietary supplements - An exploratory study directed at monitoring the internet for contaminated food supplements based on the reported effects of their users

Dietary supplements are products that are widely used for instance as energisers or to lose weight. There have been cases reported where undeclared ingredients present in such supplements have caused adverse effects on the health of the user. As there are many different products to choose from, it seems impossible to predict which might contain harmful components and to ban them from the market. Nonetheless, the use of dietary supplements and the experiences of users are shared in online discussions. We describe the development of a search engine to retrieve products associated with certain effects. Upon application we were able to retrieve a list of dietary supplements that are repeatedly associated with excessive effects by users on public fora. The top of the list contains supplements that have previously been banned because they contained undeclared harmful components. The use of the search engine as described here is a powerful method for making a risk-based selection of dietary supplements which can then be analysed for the presence of illegal or other unwanted components.

A Focal Adhesion Kinase-Derived Peptide Binds the Src SH3 Domain in Two Orientations, As Demonstrated Using Paramagnetic Nuclear Magnetic Resonance

SH3 binding peptides contain polyproline helices and are classified according to their binding orientations as N-to-C-terminal or C-to-N-terminal. We have tested the hypothesis that such a peptide binds in both orientations but with different populations. A focal adhesion kinase (FAK)-derived peptide was tested for its binding orientation on the Src SH3 domain. Paramagnetic tags were introduced at several positions on the SH3 domain, and on the basis of the paramagnetic relaxation enhancements (PREs) of the amide protons, the positions of the paramagnetic centers were determined. Two peptides were synthesized with (13)C-enriched Ala or Pro, at the N-terminal or C-terminal side of the peptide, and the intermolecular PREs were measured. The results provide compelling evidence that the FAK-derived peptide binds the SH3 domain in two orientations. In the major state, the SH3 domain binds the peptide in the N-C orientation, whereas 20% of the time, the peptide binds in the C-N orientation. We conclude that the distinction between N-C and C-N orientations, which is based on crystal structures, might be artificial. The pseudosymmetric nature of the polyproline helix might allow for binding in both orientations in the solution state.

[Not Available]

Online shopping is on the rise and this extends to purchasing of medicines. Patients can deliberately or unwittingly purchase medicines from one of the many illegal suppliers. On this illegal market, erectile dysfunction medicines play a prominent role. We estimate that at least 70% of sildenafil used in the Netherlands is purchased from illegal suppliers. Price cuts on the legal market with the introduction of generics have not made a difference. It is important to gain an insight into the scale of illegal use of medicines and the potential harmful effects we may encounter.

History full circle: 'Novel' sympathomimetics in supplements

Since the banning of ephedrine in over-the-counter nutritional supplements a decade ago, a plethora of untested and/or unsafe sympathomimetic stimulants have taken its place. This paper argues that these 'novel' stimulants in supplements recapitulate the work of synthetic chemists at commercial pharmaceutical firms during the 1930s and 1940s, all seeking substitutes for recently successful products based on ephedrine and amphetamine. Copyright © 2015 John Wiley & Sons, Ltd.

Chemical fingerprinting of silicone-based breast implants

With millions of women worldwide carrying them, silicone-based breast implants represent a large market. Even though silicone breast implants already have a history of use of more than 50 years, the discussion on their safety has not yet come to an end. To improve safety assessment, regulatory authorities should have the availability of a set of tests to be able to determine parameters of implant identity and quality. Therefore, the gels and envelopes of various brands and types of silicone-based breast implants have been subjected to infrared, Raman and NMR spectroscopy. We show that by using a combination of complementary spectroscopic techniques breast implants of various origins can be distinguished on typical chemical hallmarks. It was found that typical silicone-based implants display a surplus of vinyl signals in the gel, cyclosiloxane impurities are tolerable at low levels only and a barrier layer is present in the implant envelope. The techniques presented here and the results obtained offer a good starting point for market surveillance studies.

Small-molecule binding sites on proteins established by paramagnetic NMR spectroscopy

Determining the three-dimensional structure of a small molecule-protein complex with weak affinity can be a significant challenge. We present a paramagnetic NMR method to determine intermolecular structure restraints based on pseudocontact shifts (PCSs). Since the ligand must be in fast exchange between free and bound states and the fraction bound can be as low as a few percent, the method is ideal for ligands with high micromolar to millimolar dissociation constants. Paramagnetic tags are attached, one at a time, in a well-defined way via two arms at several sites on the protein surface. The ligand PCSs were measured from simple 1D (1)H spectra and used as docking restraints. An independent confirmation of the complex structure was carried out using intermolecular NOEs. The results show that structures derived from these two approaches are similar. The best results are obtained if the magnetic susceptibility tensors of the tags are known, but it is demonstrated that with two-armed probes, the magnetic susceptibility tensor can be predicted with sufficient accuracy to provide a low-resolution model of the ligand orientation and the location of the binding site in the absence of isotope-labeled protein. This approach can facilitate fragment-based drug discovery in obtaining structural information on the initial fragment hits.

PARAssign--paramagnetic NMR assignments of protein nuclei on the basis of pseudocontact shifts

The use of paramagnetic NMR data for the refinement of structures of proteins and protein complexes is widespread. However, the power of paramagnetism for protein assignment has not yet been fully exploited. PARAssign is software that uses pseudocontact shift data derived from several paramagnetic centers attached to the protein to obtain amide and methyl assignments. The ability of PARAssign to perform assignment when the positions of the paramagnetic centers are known and unknown is demonstrated. PARAssign has been tested using synthetic data for methyl assignment of a 47 kDa protein, and using both synthetic and experimental data for amide assignment of a 14 kDa protein. The complex fitting space involved in such an assignment procedure necessitates that good starting conditions are found, both regarding placement and strength of paramagnetic centers. These starting conditions are obtained through automated tensor placement and user-defined tensor parameters. The results presented herein demonstrate that PARAssign is able to successfully perform resonance assignment in large systems with a high degree of reliability. This software provides a method for obtaining the assignments of large systems, which may previously have been unassignable, by using 2D NMR spectral data and a known protein structure.

A pH-sensitive, colorful, lanthanide-chelating paramagnetic NMR probe

Paramagnetic lanthanides ions are broadly used in NMR spectroscopy. The effects of unpaired electrons on NMR spectral parameters provide a powerful tool for the characterization of macromolecular structures and dynamics. Here, a new lanthanide-chelating NMR probe, Caged Lanthanide NMR Probe-7 (CLaNP-7), is presented. It can be attached to protein surfaces via two disulfide bridges, yielding a probe that is rigid relative to the protein backbone. CLaNP-7 extends the application range of available probes. It has a yellow color, which is helpful for sample preparation. Its effects are comparable to those of CLaNP-5, but its charge is two units lower (+1) than that of CLaNP-5 (+3), reducing the change in surface potential after probe attachment. It also has a different magnetic susceptibility tensor, so by using both tags, two sets of structural restraints can be obtained per engineered cysteine pair. Moreover, it was found that the orientation of the magnetic susceptibility tensor is pH dependent (pK(a) ≈ 7) when a histidine residue is located in the neighborhood of the probe attachment site. The results show that the His imidazole group interacts with the CLaNP-7 tag. It is proposed that the histidine residue forms a hydrogen bond to a water/hydroxyl molecule that occupies the ninth coordination position on the lanthanide, thus breaking the two-fold symmetry of the CLaNP tag in a pH-dependent way.

Narrowing the conformational space sampled by two-domain proteins with paramagnetic probes in both domains

Calmodulin is a two-domain protein which in solution can adopt a variety of conformations upon reorientation of its domains. The maximum occurrence (MO) of a set of calmodulin conformations that are representative of the overall conformational space possibly sampled by the protein, has been calculated from the paramagnetism-based restraints. These restraints were measured after inclusion of a lanthanide binding tag in the C-terminal domain to supplement the data obtained by substitution of three paramagnetic lanthanide ions to the calcium ion in the second calcium binding loop of the N-terminal domain. The analysis shows that the availability of paramagnetic restraints arising from metal ions placed on both domains, reduces the MO of the conformations to different extents, thereby helping to identify those conformations that can be mostly sampled by the protein.

Validation of a lanthanide tag for the analysis of protein dynamics by paramagnetic NMR spectroscopy

Paramagnetic lanthanide tags potentially can enhance the effects of microsecond to millisecond dynamics in proteins on NMR signals and provide structural information on lowly populated states encoded in the pseudocontact shifts. We have investigated the microsecond to millisecond mobility of a two-point attached lanthanide tag, CLaNP-5, using paramagnetic (1)H CPMG relaxation dispersion methods. CLaNP-5 loaded with Lu(3+), Yb(3+), or Tm(3+) was attached to three sites on the surface of two proteins, pseudoazurin and cytochrome c. The paramagnetic center causes large relaxation dispersion effects for two attachment sites, suggesting that local dynamics of the protein at the attachment site causes mobility of the paramagnetic center. At one site the relaxation dispersions are small and limited to the immediate environment of the tag. It is concluded that paramagnetic relaxation dispersion could represent a sensitive method to probe protein dynamics. However, the selection of a rigid attachment site is of critical importance.

A solution model of the complex formed by adrenodoxin and adrenodoxin reductase determined by paramagnetic NMR spectroscopy

Lanthanide tags offer the opportunity to retrieve long-range distance information from NMR experiments that can be used to guide protein docking. To determine whether sufficient restraints can be retrieved for proteins with low solubility and availability, Ln tags were applied in the study of the 65 kDa membrane-associated protein complex formed by the electron carrier adrenodoxin and its electron donor, adrenodoxin reductase. The reductase is only monomeric at low concentration, and the paramagnetic iron-sulfur cluster of adrenodoxin broadens many of the resonances of nuclei in the interface. Guided by the paramagnetic restraints obtained using two Ln-tag attachment sites, protein docking yields a cluster of solutions with an rmsd of 3.2 A. The mean structure is close to the crystal structure of the cross-linked complex, with an rmsd of 4.0 A. It is concluded that with the application of Ln tags paramagnetic NMR restraints for structure determination can be retrieved even for difficult, low-concentration protein complexes.

Metabolism ofN-substituted 7-methoxy-4-(aminomethyl) -coumarins by cytochrome P450 2D6 mutants and the indication of additional substrate interaction points

Previous studies have shown the critical roles residues F120 and F483 play in the oxidative metabolism of 7-methoxy-4-(aminomethyl)-coumarin (MAMC) by cytochrome P450 2D6 (CYP2D6). In the present study, a series of N-alkyl-7-methoxy-4-(aminomethyl)-coumarins (MAMC analogues) were used as substrates for the F120A and F483A mutants in order to probe the CYP2D6 active site. The F120A and F483A mutants of CYP2D6 displayed significant activity towards the MAMC analogues. Automated docking studies of the MAMC analogues in a CYP2D6 homology model suggested a distal hydrophobic active site binding cleft for the substrate N-alkyl chains, consisting of the residues L213 and V308.

Intermolecular dynamics studied by paramagnetic tagging

Yeast cytochrome c and bovine adrenodoxin form a dynamic electron transfer complex, which is a pure encounter complex. It is demonstrated that the dynamic nature of the interaction can readily be probed by using a rigid lanthanide tag attached to cytochrome c. The tag, Caged Lanthanide NMR Probe 5, induces pseudocontact shifts and residual dipolar couplings and does not perturb the binding interface. Due to the dynamics in the complex, residual dipolar couplings in adrenodoxin are very small. Simulation shows that cytochrome c needs to sample a large part of the surface of adrenodoxin to explain the small degree of alignment observed for adrenodoxin. The applied method provides a simple and straightforward way to observe dynamics in protein complexes or domain-domain mobility without the need for external alignment media.

Design, synthesis, and evaluation of a lanthanide chelating protein probe: CLaNP-5 yields predictable paramagnetic effects independent of environment

Immobilized lanthanide ions offer the opportunity to refine structures of proteins and the complexes they form by using restraints obtained from paramagnetic NMR experiments. We report the design, synthesis, and spectroscopic evaluation of the lanthanide chelator, Caged Lanthanide NMR Probe 5 (CLaNP-5) readily attachable to a protein surface via two cysteine residues. The probe causes tunable pseudocontact shifts, alignment, paramagnetic relaxation enhancement, and luminescence, by chelating it to the appropriate lanthanide ion. The observation of single shifts and the finding that the magnetic susceptibility tensors obtained from shifts and alignment analyses are highly similar strongly indicate that the probe is rigid with respect to the protein backbone. By placing the probe at various positions on a model protein it is demonstrated that the size and orientation of the magnetic susceptibility tensor of the probe are independent of the local protein environment. Consequently, the effects of the probe are readily predictable using a protein structure only. These findings designate CLaNP-5 as a protein probe to deliver unambiguous high quality structural restraints in studies on protein-protein and protein-ligand interactions.

Active-site structure, binding and redox activity of the heme-thiolate enzyme CYP2D6 immobilized on coated Ag electrodes: a surface-enhanced resonance Raman scattering study

Surface-enhance resonance Raman scattering spectra of the heme-thiolate enzyme cytochrome P450 2D6 (CYP2D6) adsorbed on Ag electrodes coated with 11-mercaptoundecanoic acid (MUA) were obtained in various experimental conditions. An analysis of these spectra, and a comparison between them and the RR spectra of CYP2D6 in solution, indicated that the enzyme's active site retained its nature of six-coordinated low-spin heme upon immobilization. Moreover, the spectral changes detected in the presence of dextromethorphan (a CYP2D6 substrate) and imidazole (an exogenous heme axial ligand) indicated that the immobilized enzyme also preserved its ability to reversibly bind a substrate and form a heme-imidazole complex. The reversibility of these processes could be easily verified by flowing alternately solutions of the various compounds and the buffer through a home-built spectroelectrochemical flow cell which contained a sample of immobilized protein, without the need to disassemble the cell between consecutive spectral data acquisitions. Despite immobilized CYP2D6 being effectively reduced by a sodium dithionite solution, electrochemical reduction via the Ag electrode was not able to completely reduce the enzyme, and led to its extensive inactivation. This behavior indicated that although the enzyme's ability to exchange electrons is not altered by immobilization per se, MUA-coated electrodes are not suited to perform direct electrochemistry of CYP2D6.

Free energies of binding of R- and S-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences.

Altered spin state equilibrium in the T309V mutant of cytochrome P450 2D6: a spectroscopic and computational study

Cytochrome P450 2D6 (CYP2D6) is one of the most important cytochromes P450 in humans. Resonance Raman data from the T309V mutant of CYP2D6 show that the substitution of the conserved I-helix threonine situated in the enzyme's active site perturbs the heme spin equilibrium in favor of the six-coordinated low-spin species. A mechanistic hypothesis is introduced to explain the experimental observations, and its compatibility with the available structural and spectroscopic data is tested using quantum-mechanical density functional theory calculations on active-site models for both the CYP2D6 wild type and the T309V mutant.

Surface-enhanced resonance Raman scattering of cytochrome P450-2D6 on coated silver hydrosols

Surface-enhanced resonance Raman scattering (SERRS) from dilute solutions (down to nanomolar concentrations) of human mono-oxygenase CYP2D6 is observed using aqueous dispersions of Ag nanoparticles (hydrosol) coated with self-assembled monolayers (SAMs) of mercaptoalkanoic acids of two different lengths. From a direct comparison with its resonance Raman spectrum in solution, CYP2D6 appears to fully retain its native structure upon adsorption on coated hydrosol through electrostatic interaction, while a structural change in the active site is observed when uncoated citrate-reduced hydrosol is used. Using SERRS on these biocompatible coated hydrosols, the effects of dextromethorphan on the enzyme's active site can be observed, demonstrating that CYP2D6 ability of binding substrates is preserved. Moreover, by tuning the wavelength of the exciting laser away from the main absorption band of the heme, the vibrational bands of the SAM coating are observed and analyzed to see how the presence of the protein affects the SAM structure.

Molecular Modeling-Guided Site-Directed Mutagenesis of Cytochrome P450 2D6

Cytochrome P450 (CYP) 2D6 is one of the most important drug metabolizing enzymes and the rationalization and prediction of potential CYP2D6 substrates is therefore advantageous in the discovery and development of new drugs. Experimentally, the active site of CYP2D6 can be probed by site directed mutagenesis studies. Such studies can be designed from structural models of enzyme-substrate complexes. Modeling approaches can subsequently be used to rationalize the observed effect of mutations on metabolism and inhibition. The current paper will present the construction, refinement and validation of the CYP2D6 homology model used in our laboratory for the prediction and rationalisation of CYP2D6 substrate metabolism and CYP2D6-ligand interactions. The model could explain reported site-directed mutagenesis data (for example, mutation of E216 and D301). Furthermore, based on the model, new CYP2D6 mutants were constructed and studied in our lab, and also for these mutants a rationalization of experimentally observed characteristics could be achieved (I106E, F120A, T309V, F483A). CYP2D6-substrate interaction fingerprint analysis of docked substrates in our homology model suggests that several other active site residues are probably interacting with ligands as well, opening the way for further mutagenesis studies. Our homology model was found to agree with most of the details of the recently solved substrate-free CYP2D6 crystal structure [Rowland et al. J. Biol. Chem. 2006, 281, 7614-7622]. Structural differences between the homology model and crystal structure were the same differences observed between substrate-free and substrate-bound structures of other CYPs, suggesting that these conformational changes are required upon substrate binding. The CYP2D6 crystal structure further validates our homology modeling approach and shows that computational chemistry is a useful and valuable tool to provide models for substrate-bound complexes of CYPs which give insight into CYP-ligand interactions. This information is essential for successful pre-experimental virtual screening, as well as accurate hypothesis generation for in vitro studies in drug discovery and development.

Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

Binding of bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine to wild-type and F120A mutant cytochrome P450 2D6 studied by resonance Raman spectroscopy

Cytochrome P450 2D6 (CYP2D6) is one of the most important drug-metabolizing enzymes in humans. Resonance Raman data, reported for the first time for CYP2D6, show that the CYP2D6 heme is found to be in a six-coordinated low-spin state in the absence of substrates, and it is perturbed to different extents by bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine (MDMA). Dextromethorphan and MDMA induce in CYP2D6 a significant amount of five-coordinated high-spin heme species and reduce the polarity of its heme-pocket, whereas bufuralol does not. Spectra of the F120A mutant CYP2D6 suggest that Phe120 is involved in substrate-binding of dextromethorphan and MDMA, being responsible for the spectral differences observed between these two compounds and bufuralol. These differences could be explained postulating a different substrate mobility for each compound in the CYP2D6 active site, consistently with the role previously suggested for Phe120 in binding dextromethorphan and MDMA.

Topological role of cytochrome P450 2D6 active site residues

Recent reports have identified Phe120, Asp301, Thr309, and Glu216 as important residues in cytochrome P450 2D6 (CYP2D6) substrate binding and catalysis. Complementary homology models have located these amino acids within the binding pocket of CYP2D6 and in the present study we have used aryldiazenes to test these models and gain further insight in the role these amino acids have in maintaining the integrity of the active site cavity. When Phe120 was replaced to alanine, there was a significant increase in probe migration to pyrrole nitrogens C and D, in agreement with homology models which have located the phenyl side-chain of Phe120 above these two pyrrole rings. No changes in topology were observed with the D301Q mutant, supporting claims that in this mutant the electrostatic interactions with the B/C-loop are largely maintained and the loop retains its native orientation. The T309V mutation resulted in significant topological alteration suggesting that, in addition to its potential role in dioxygen activation, Thr309 plays an important structural role within the active site crevice. Replacement of Ile106 with Glu, engineered to cause electrostatic repulsion with Glu216, had a profound topological effect in the higher region within the active site cavity and impaired the catalytic activity towards CYP2D6 probe substrates.

Role of the conserved threonine 309 in mechanism of oxidation by cytochrome P450 2D6

Based on sequence alignments and homology modeling, threonine 309 in cytochrome P450 2D6 (CYP2D6) is proposed to be the conserved I-helix threonine, which is supposed to be involved in dioxygen activation by CYPs. The T309V mutant of CYP2D6 displayed a strong shift from O-dealkylation to N-dealkylation reactions in oxidation of dextromethorphan and 3,4-methylenedioxymethylamphetamine. This may be explained by an elevated ratio of hydroperoxo-iron to oxenoid-iron of the oxygenating species. In consistence, using cumene hydroperoxide, which directly forms the oxenoid-iron, the T309V mutant again selectively catalyzed the O-dealkylation reactions. The changed ratio of oxygenating species can also explain the decreased activity and changed regioselectivity that were observed in 7-methoxy-4-(aminomethyl)-coumarin and bufuralol oxidation, respectively, by the T309V mutant. Interestingly, the T309V mutant always showed a significantly increased, up to 75-fold, higher activity compared to that of the wild-type when using cumene hydroperoxide. These results indicate that T309 in CYP2D6 is involved in maintaining the balance of multiple oxygenating species and thus influences substrate and regioselectivity.

The role of phenylalanine 483 in cytochrome P450 2D6 is strongly substrate dependent

The polymorphic cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 30% of the drugs currently prescribed, and is thus clinically relevant. Typical CYP2D6 substrates generally contain a basic nitrogen atom and an aromatic moiety adjacent to the site of metabolism. Recently, we demonstrated the importance of active site residue F120 in substrate binding and catalysis in CYP2D6. On the basis of protein homology models, it is claimed that another active site phenylalanine, F483, may also play an important role in the interaction with the aromatic moiety of CYP2D6 substrates. Experimental data to support this hypothesis, however, is not yet available. In fact, in the only study performed, mutation of F483 to isoleucine or tryptophan did not affect the 1'-hydroxylation of bufuralol at all [Smith G, Modi S, Pillai I, Lian LY, Sutcliffe MJ, Pritchard MP, et al., Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity. Biochem J 1998;331:783-92]. In the present study, the role of F483 in ligand binding and metabolism by CYP2D6 was examined experimentally using site-directed mutagenesis. Replacement of F483 by alanine resulted in a 30-fold lower V(max) for bufuralol 1'-hydroxylation, while the K(m) was hardly affected. The V(max) for 3,4-methylenedioxy-methylamphetamine O-demethylenation on the other hand decreased only two-fold, whereas the effect on the K(m) was much larger. For dextromethorphan, in addition to dextrorphan (O-demethylation) and 3-methoxymorphinan (N-demethylation), two other metabolites were formed that could not be detected for the wild-type. The substrate 7-methoxy-4-(aminomethyl)-coumarin was not metabolised at all by CYP2D6[F483A], a phenomenon that was reported also for CYP2D6[F120A]. The presented data show that next to F120, residue F483 plays a very important role in the metabolism of typical CYP2D6 substrates. The influence of F483 on metabolism was found to be strongly substrate-dependent.

Metabolic regio- and stereoselectivity of cytochrome P450 2D6 towards 3,4-methylenedioxy-N-alkylamphetamines: in silico predictions and experimental validation

A series of 3,4-methylenedioxy-N-alkylamphetamines (MDAAs) were automatically docked and subjected to molecular dynamics (MD) simulations in a cytochrome P450 2D6 (CYP2D6) protein model. The predicted substrate binding orientations, sites of oxidation, and relative reactivities were compared to the experimental data of wild-type and Phe120Ala mutant CYP2D6. Automated docking results were not sufficient to accurately rationalize experimental binding orientations of 3,4-methylenedioxy-N-methylamphetamine (MDMA) in the two enzymes as measured with spin lattice relaxation NMR. Nevertheless, the docking results could be used as starting structures for MD simulations. Predicted binding orientations of MDMA and sites of oxidation of the MDAAs derived from MD simulations matched well with the experimental data. It appeared the experimental results were best described in MD simulations considering the nitrogen atoms of the MDAAs in neutral form. Differences in regioselectivity and stereoselectivity in the oxidative metabolism of the MDAAs by the Phe120Ala mutant CYP2D6 were correctly predicted, and the effects of the Phe120Ala mutation could be rationalized as well.

Metabolic Regio- and Stereoselectivity of Cytochrome P450 2D6 towards 3,4-Methylenedioxy-N-alkylamphetamines: in Silico Predictions and Experimental Validation

A series of 3,4-methylenedioxy-N-alkylamphetamines (MDAAs) were automatically docked and subjected to molecular dynamics (MD) simulations in a cytochrome P450 2D6 (CYP2D6) protein model. The predicted substrate binding orientations, sites of oxidation, and relative reactivities were compared to the experimental data of wild-type and Phe120Ala mutant CYP2D6. Automated docking results were not sufficient to accurately rationalize experimental binding orientations of 3,4-methylenedioxy-N-methylamphetamine (MDMA) in the two enzymes as measured with spin lattice relaxation NMR. Nevertheless, the docking results could be used as starting structures for MD simulations. Predicted binding orientations of MDMA and sites of oxidation of the MDAAs derived from MD simulations matched well with the experimental data. It appeared the experimental results were best described in MD simulations considering the nitrogen atoms of the MDAAs in neutral form. Differences in regioselectivity and stereoselectivity in the oxidative metabolism of the MDAAs by the Phe120Ala mutant CYP2D6 were correctly predicted, and the effects of the Phe120Ala mutation could be rationalized as well.

Influence of phenylalanine 120 on cytochrome P450 2D6 catalytic selectivity and regiospecificity: crucial role in 7-methoxy-4-(aminomethyl)-coumarin metabolism

The polymorphic human debrisoquine hydroxylase, cytochrome P450 2D6 (CYP2D6), is one of the most important phase I drug metabolising enzymes. It is responsible for metabolising a large number of compounds that mostly share similarity in having a basic N-atom and an aromatic moiety. In homology modelling studies, it has been suggested that in fixation of this aromatic moiety, there may be an important role for phenylalanine 120 (Phe(120)). In this study, the role of Phe(120) in ligand binding and catalysis was experimentally examined by mutating it into an alanine. Strikingly, this substitution led to a completely abolished 7-methoxy-4-(aminomethyl)-coumarin (MAMC) O-demethylating activity of CYP2D6. On the other hand, bufuralol metabolism was hardly affected (K(m) of 1-hydroxylation mutant: 1.2 microM, wild-type: 2.9 microM, 4-hydroxylation mutant: 1.5 microM, and wild-type: 3.2 microM) and neither was affected dextromethorphan O-demethylation (K(m) mutant: 1.2 microM, wild-type: 2 microM, k(cat) mutant: 4.5 min(-1), and wild-type: 3.3 min(-1)). However, the Phe(120)Ala mutant also formed 3-hydroxymorphinan, the double demethylated form of dextromethorphan, which was not detected using wild-type CYP2D6. 3,4-Methylenedioxymethamphetamine (MDMA) was demethylenated by both mutant and wild-type CYP2D6 to 3,4-dihydroxymethamphetamine (3,4-OH-MA K(m) of mutant: 55 microM and wild-type: 2 microM). In addition, the mutant formed two additional metabolites; 3,4-methylenedioxyamphetamine (MDA) and N-hydroxy-3,4-methylenedioxymethamphetamine (N-OH-MDMA). Inhibition experiments of dextromethorphan O-demethylation showed a decreased affinity of the Phe(120)Ala mutant for quinidine (IC(50) mutant: 240 nM and wild-type, 40 nM), while IC(50)s for quinine were equal (1 microM). These data indicate the importance of Phe(120) in the selectivity and regiospecificity in substrate binding and catalysis by CYP2D6.