Quantitative Structure–Activity Relationship Study of Binding Affinity of Azole Compounds with CYP2B and CYP3A

The synergistic potential of azole fungicides does not directly correlate to the inhibition of cytochrome P450 activity in aquatic invertebrates

The ability of azole fungicides to inhibit cytochrome P450 dependent metabolism is proposed to be the main mechanism for their synergizing effect on pyrethroid insecticide toxicity in aquatic invertebrates. This study investigates the correlation between inhibition strength and synergistic potential of azole fungicides in the crustacean Daphnia magna and the insect larvae Chironomus riparius. Inhibition strength was measured in vivo toward the cytochrome P450 catalysed conversion of 7-ethoxycoumarin to 7-hydroxycoumarin (ECOD). Synergistic potentials were determined as the ratio between predicted and observed toxicity of mixtures based on the model of concentration addition (CA) and independent action (IA). Azoles (n = 9-11) enhanced the toxicity of α-cypermethrin in D. magna (Synergy ratios CA: 0.8 - 16; IA: 1.1 - 22) and inhibited cytochrome P450 activity by different degrees (IC50: 0.0023 - 36 μM for D. magna and 0.08 - 24 μM for C. riparius). Inhibition strengths were strongly correlated in the two organisms (r: 0.937 p: 0.019 for triazoles and r: 0.903 p: 0.097 for imidazoles). Lipophilicity governed the inhibition strength of triazoles in both species (r > 0.9, p < 0.05). No correlation was observed between inhibition strengths and synergistic potentials. Several reasons for the apparent lack of correlation were discussed.

Can the inhibition of cytochrome P450 in aquatic invertebrates due to azole fungicides be estimated with in silico and in vitro models and extrapolated between species?

Azole fungicides, designed to halt fungal growth by specific inhibition of fungal cytochrome P450 (CYP51), inhibit cytochrome P450s involved in the metabolism of xenobiotics in several non-target organisms thus raising environmental concern. The present study investigates the degree by which inhibition strengths of azoles toward cytochrome P450 in rat liver, the insect Chironomus riparius larvae and the snail Lymnaea stagnalis can be extrapolated from estimated in silico affinities. Azoles' affinities toward human cytochrome P450 isoforms involved in xenobiotic metabolism (CYP3A4, CYP2C9 and CYP2D6) as well as fungal CYP51 were estimated with a ligand-protein docking model based on the ChemScore scoring function. Estimated affinities toward the selected enzymatic structures correlated strongly with measured inhibition strengths in rat liver (ChemScore vs. logIC₅₀ among cytochrome P450 isoforms: -0.662 < r < -0.891, n = 17 azoles), while weaker correlations were found for C. riparius larvae (-0.167 < r < -0.733, n = 9) and L. stagnalis (-0.084 < r < -0.648, n = 8). Inhibition strengths toward C. riparius and rat liver activities were found to be highly correlated to each other (r: 0.857) while no significant relationship was found between either of the species and L. stagnalis. The inhibition of cytochrome P450 due to azole fungicides could be estimated in vitro and to a lesser extent in silico for C. riparius but not for L. stagnalis, possibly due to different enzymatic susceptibility toward azole inhibition among the species.

What causes the difference in synergistic potentials of propiconazole and prochloraz toward pyrethroids in Daphnia magna?

Azole fungicides (imidazoles and triazoles) are known to function synergistically with several compounds, especially with pyrethroid insecticides, most likely by inhibiting cytochrome P450. Different azole fungicides have been shown to differ in their synergistic potentials usually with the imidazoles being stronger synergists than the triazoles. This study investigated whether the toxicokinetic and toxicodynamic (TKTD) properties of the imidazole prochloraz and triazole propiconazole can explain their different synergistic potential toward the freshwater macroinvertebrate Daphnia magna. Pulse exposure to external concentrations of propiconazole (1.4μM) and prochloraz (1.7μM) for 18h resulted in internal concentrations of 22.7 and 53.5μmolkg(-1)w.w. for propiconazole and prochloraz, respectively. This 2-fold difference in bioaccumulation corresponded very well with the observed 2.7-fold lower external EC50-estimate (7 days) for prochloraz compared to propiconazole. The estimated IC50 for the in vivo inhibition of cytochrome P450 (ECOD) activity, however, measured as transformation of 7-ethoxycoumarin into 7-hydroxycoumarin, was almost 500-fold higher for prochloraz (IC50: 0.011±0.002μM) compared to propiconazole (IC50: 4.9±0.06μM). When indirectly measuring the binding strength of the two azoles, daphnids exposed to propiconazole recovered roughly 80% of their ECOD activity compared to the control shortly after being moved to azole-free medium, indicating that propiconazole causes reversible inhibition of cytochrome P450. In contrast, the ECOD-activity remained inhibited in the prochloraz-exposed daphnids for 12h following transfer to azole-free medium, which correlated with elimination of the measured internal prochloraz concentration (DT95≈13h). These results indicate that lethal toxicity of the azole fungicides is mainly driven by toxicokinetics through their hydrophobicities resulting in different internal concentrations. Their synergistic potential toward pyrethroid toxicity, on the other hand, is mainly governed by their toxicodynamic effects measured as the differences in IC50-values toward in vivo cytochrome P450 (ECOD) activity together with the proposed binding strength measured indirectly through the recovery of ECOD activity as a function of internal azole concentrations.

Functional Analysis of the Unique Cytochrome P450 of the Liver Fluke Opisthorchis felineus

The basic metabolic cytochrome P450 (CYP) system is essential for biotransformation of sterols and xenobiotics including drugs, for synthesis and degradation of signaling molecules in all living organisms. Most eukaryotes including free-living flatworms have numerous paralogues of the CYP gene encoding heme monooxygenases with specific substrate range. Notably, by contrast, the parasitic flatworms have only one CYP gene. The role of this enzyme in the physiology and biochemistry of helminths is not known. The flukes and tapeworms are the etiologic agents of major neglected tropical diseases of humanity. Three helminth infections (Opisthorchis viverrini, Clonorchis sinensis and Schistosoma haematobium) are considered by the International Agency for Research on Cancer (IARC) as definite causes of cancer. We focused our research on the human liver fluke Opisthorchis felineus, an emerging source of biliary tract disease including bile duct cancer in Russia and central Europe. The aims of this study were (i) to determine the significance of the CYP activity for the morphology and survival of the liver fluke, (ii) to assess CYP ability to metabolize xenobiotics, and (iii) to localize the CYP activity in O. felineus tissues. We observed high constitutive expression of CYP mRNA (Real-time PCR) in O. felineus. This enzyme metabolized xenobiotics selective for mammalian CYP2E1, CYP2B, CYP3A, but not CYP1A, as determined by liquid chromatography and imaging analyses. Tissue localization studies revealed the CYP activity in excretory channels, while suppression of CYP mRNA by RNA interference was accompanied by morphological changes of the excretory system and increased mortality rates of the worms. These results suggest that the CYP function is linked to worm metabolism and detoxification. The findings also suggest that the CYP enzyme is involved in vitally important processes in the organism of parasites and is a potential drug target.

The basic metabolic system CYP (cytochrome P450) is essential for biotransformation of sterols and xenobiotics, for synthesis and degradation of signaling molecules in all living organisms. Most eukaryotes including free-living flatworms evolved numerous paralogues of the CYP gene. Notably, by contrast, flukes and tapeworms–the etiologic agents of major neglected tropical diseases of humanity, have only one gene. However, the role of P450 in the physiology and biochemistry of helminths is not known. This report presents the first functional study of the CYP enzyme of any of the parasitic flatworms. We focused our research on the food-borne human liver fluke, Opisthorchis felineus, an emerging source of biliary tract diseases in Russia, Kazakhstan and central Europe. Here we report that this liver fluke has evolved a highly expressed functional monooxygenase system with broad substrate specificity. Tissue localization studies and suppression of CYP mRNA by RNA interference revealed the CYP function is linked to the excretory system and possibly to metabolism and detoxification. The fluke’s monooxygenase likely is a model for orthologues of the singular CYP of parasitic flatworms at large, where it plays a critical role in the pathogen’s metabolism that contributes to worm survival and drug resistance.

QSAR of cytochrome inhibitors

Cytochrome P450 (CYP450) enzymes are predominantly involved in the Phase I metabolism of xenobiotics. Metabolic inhibition and induction can give rise to clinically important drug-drug interactions. Metabolic stability is a prerequisite for sustaining the therapeutically relevant concentrations, and very often drug candidates are sacrificed due to poor metabolic profiles. Computational tools such as quantitative structure-activity relationships are widely used to study different metabolic end points successfully to accelerate the drug discovery process. There are a lot of computational studies on clinically important CYPs already reported in recent years. But other clinically significant families are to yet be explored computationally. Powerfulness of quantitative structure-activity relationship will drive computational chemists to develop new potent and selective inhibitors of different classes of CYPs for the treatment of different diseases with least drug-drug interactions. Furthermore, there is a need to enhance the accuracy, interpretability and confidence in the computational models in accelerating the drug discovery pathways.

Toward basic understanding of the partition coefficient log P and its application in QSAR

The log P value has been the first choice for the molecular hydrophobicity descriptor in QSAR studies. However, it is still now difficult to understand the partitioning phenomenon in terms of physical chemistry. First, an attempt to understand and predict log P is addressed. We formulated a simple model that expressed by the solvent accessible surface area and the solvation energy difference between aqueous and solvent phases. Next, an application of log P in QSAR analyses of ligand-CYP (Cytochrome P450) interaction was described. Azole compounds are widely used as antifungal agents. We showed that the binding affinity of 18 azole compounds with CYP2B and CYP3A were nicely expressed by the bilinear model of log P. These results suggest that molecular hydrophobicity plays a major role in the binding affinity. The binding mechanism was discussed based on the correlation equations.