Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate.

Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcoholß-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and ßC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the ßC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcoholß-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehydeß-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.

Identification and in silico prediction of metabolites of tebufenozide derivatives by major human cytochrome P450 isoforms.

Cytochrome P450 (CYP) enzymes constitute a superfamily of heme-containing monooxygenases. CYPs are involved in the metabolism of many chemicals such as drugs and agrochemicals. Therefore, examining the metabolic reactions by each CYP isoform is important to elucidate their substrate recognition mechanisms. The clarification of these mechanisms may be useful not only for the development of new drugs and agrochemicals, but also for risk assessment of chemicals. In our previous study, we identified the metabolites of tebufenozide, an insect growth regulator, formed by two human CYP isoforms: CYP3A4 and CYP2C19. The accessibility of each site of tebufenozide to the reaction center of CYP enzymes and the susceptibility of each hydrogen atom for metabolism by CYP enzymes were evaluated by a docking simulation and hydrogen atom abstraction energy estimation at the density functional theory level, respectively. In this study, the same in silico prediction method was applied to the metabolites of tebufenozide derivatives by major human CYPs (CYP1A2, 2C9, 2C19, 2D6, and 3A4). In addition, the production rate of the metabolites by CYP3A4 was quantitively analyzed by frequency based on docking simulation and hydrogen atom abstraction energy using the classical QSAR approach. Then, the obtained QSAR model was applied to predict the sites of metabolism and the metabolite production order by each CYP isoform.

Quantitative analysis of the relationship between structure and antioxidant activity of tripeptides.

Peptides from enzymatic hydrolysates of food proteins exhibit significant antioxidant activity. Several studies have attempted to determine the factors contributing to the antioxidant activity of peptides; however, the physicochemical properties and factors essential for the antioxidant activity of peptides are still unclear. In this study, in order to clarify the factors important for peptide antioxidant activity based on the properties of component amino acids, 55 tripeptides were synthesized from 20 natural amino acids and their antioxidant activity was measured using the Trolox equivalent antioxidant capacity (TEAC) assay system. The tripeptides were divided into two data sets: a training set comprising 50 compounds and a validated set comprising five compounds. The structure-activity relationship of the training set was then analyzed using classical quantitative structure-activity relationship (QSAR) analysis. The study findings demonstrate that the presence of a cysteine residue at any position, an aromatic amino acid at the C-terminus, higher hydrophobicity of the N-terminal residue, and smaller HOMO-LUMO energy gap of the middle residue can significantly enhance the antioxidant activity. The activities of the five validated compounds were predicted using the constructed QSAR model, and a good correlation between measured and predicted activities was observed. The information obtained from the QSAR model could be useful for effective production of antioxidant peptides from food proteins such as egg white proteins.

Microchip electrophoresis utilizing an in situ photopolymerized Phos-tag binding polyacrylamide gel for specific entrapment and analysis of phosphorylated compounds.

A method was developed for the specific entrapment and separation of phosphorylated compounds using a Phos-tag polyacrylamide gel fabricated at the channel crossing point of a microfluidic electrophoresis chip. The channel intersection of the poly(methyl methacrylate)-made microchip was filled with a solution comprising acrylamide, N,N-methylene-bis-acrylamide, Phos-tag acrylamide, and 2,2'-azobis[2-methyl-N-(2-hydroxyethyl)propionamide], which functioned as a photocatalytic initiator. In situ polymerization at the channel crossing point was performed by irradiation with a UV LED laser beam. The fabricated Phos-tag gel (100 × 100 × 30 µm) contains ca. 20 fmol of the Phos-tag group and therefore could entrap phosphorylated compounds at the femtomolar level. The electrophoretically trapped phosphorylated compounds were released from the gel by switching the voltage to deliver high concentrations of phosphate and EDTA in a background electrolyte. The broad sample band eluted from the gel was effectively reconcentrated at the boundary of a pH junction generated by sodium ions delivered from the outlet reservoir. The reconcentrated sample components were then separated and fluorometrically detected at the end of the separation channel. Under the optimized conditions, the phosphorylated compounds were concentrated by a factor of 100-fold, and the peak resolution was comparable to that obtained by pinched injection. This method was successfully utilized to preconcentrate and analyze phosphorylated peptides in a complex peptide mixture.

Obituary: Toshio Fujita, QSAR pioneer.

This is the obituary for Toshio Fujita, pioneer of the quantitative structure activity relationship (QSAR) paradigm.

In vitro and in vivo evaluations of the P-glycoprotein-mediated efflux of dibenzoylhydrazines.

P-glycoprotein (P-gp) is a member of the ATP-binding cassette transporter family. It actively transports a wide variety of compounds out of cells to protect humans from xenobiotics. Thus, determining whether chemicals are substrates and/or inhibitors of P-gp is important in risk assessments of pharmacokinetic interactions among chemicals because P-gp-mediated transport processes play a significant role in their absorption and disposition. We previously reported that dibenzoylhydrazines (DBHs) such as tebufenozide and methoxyfenozide (agrochemicals) stimulated P-gp ATPase activity. However, it currently remains unclear whether these derivatives are transport substrates of P-gp and inhibit transport of other chemicals by P-gp. In the present study, in order to evaluate the interactions of DBHs with other chemicals in humans, we determined whether DBHs are P-gp transport substrates using both the in vitro bidirectional transport assay and the in vivo study of rats. In the in vivo study, we investigated the influence of P-gp inhibitors on the brain to plasma ratio of methoxyfenozide in rats. We also examined the inhibitory effects of DBHs on quinidine (a P-gp substrate) transport by P-gp in order to ascertain whether these derivatives are inhibitors of P-gp. Based on the results, DBHs were concluded to be weak P-gp transport substrates and moderate P-gp inhibitors. However, the risk of DBHs caused by interaction with other chemicals including drugs was considered to be low by considering the DBHs' potential as the substrates and inhibitors of P-gp as well as their plasma concentrations as long as DBHs are properly used.

Structure-activity relationships of dibenzoylhydrazines for the inhibition of P-glycoprotein-mediated quinidine transport.

We previously demonstrated that dibenzoylhydrazines (DBHs) are not only P-glycoprotein (P-gp) substrates, but also inhibitors. In the present study, we evaluated the inhibition of P-gp-mediated quinidine transport by two series of DBHs and performed a classical QSAR analysis and docking simulation in order to investigate the mechanisms underlying P-gp substrate/inhibitor recognition. The results of the QSAR analysis identified the hydrophobic factor as the most important for inhibitory activities, while electronic and steric effects also influenced the activities. The different substituent effects observed in each series suggested the different binding modes of each series of DBHs, which was supported by the results of the docking simulation.

Docking model of the nicotinic acetylcholine receptor and nitromethylene neonicotinoid derivatives with a longer chiral substituent and their biological activities.

In the present study, nitromethylene neonicotinoid derivatives possessing substituents that contain a sulfur atom, oxygen atom or aromatic ring at position 5 on the imidazolidine ring were synthesized to evaluate their affinity for the nicotinic acetylcholine receptor (nAChR) and their insecticidal activity against adult female houseflies. Comparing the receptor affinity of the alkylated derivative with the receptor affinity of compounds possessing either ether or thioether groups revealed that conversion of the carbon atom to a sulfur atom did not influence the receptor affinity, whereas conversion to an oxygen atom was disadvantageous for the receptor affinity. The receptor affinity of compounds possessing a benzyl or phenyl group was lower than that of the unsubstituted compound. Analysis of the three-dimensional quantitative structure-activity relationship using comparative molecular field analysis demonstrated that steric hindrance of the receptor should exist around the C3 of an n-butyl group attached at position 5 on the imidazolidine ring. A docking study of the nAChR-ligand model suggested that the ligand-binding region expands as the length of the substituent increases by brushing against the amino acids that form the binding region. The insecticidal activity of the compounds was positively correlated with the receptor affinity by considering logP and the number of heteroatoms, including sulfur and oxygen atoms, in the substituents, suggesting that the insecticidal activity is influenced by the receptor affinity, hydrophobicity, and metabolic stability of the compounds.

Identification and in silico prediction of metabolites of the model compound, tebufenozide by human CYP3A4 and CYP2C19.

The metabolites of tebufenozide, a model compound, formed by the yeast-expressed human CYP3A4 and CYP2C19 were identified to clarify the substrate recognition mechanism of the human cytochrome P450 (CYP) isozymes. We then determined whether tebufenozide metabolites may be predicted in silico. Hydrogen abstraction energies were calculated with the density functional theory method B3LYP/6-31G(∗). A docking simulation was performed using FRED software. Several alkyl sites of tebufenozide were hydroxylated by CYP3A4 whereas only one site was modified by CYP2C19. The accessibility of each site of tebufenozide to the reaction center of CYP enzymes and the susceptibility of each hydrogen atom for metabolism by CYP enzymes were evaluated by a docking simulation and hydrogen abstraction energy estimation, respectively.

Exon 3 splicing and mutagenesis identify residues influencing cell surface density of heterologously expressed silkworm (Bombyx mori) glutamate-gated chloride channels.

Glutamate-gated chloride channels (GluCls) mediate fast inhibitory neurotransmission in invertebrate nervous systems. Insect GluCls show alternative splicing, and, to determine its impact on channel function and pharmacology, we isolated GluCl cDNAs from larvae of the silkworm (Bombyx mori). We show that six B. mori glutamate-gated chloride channel variants are generated by splicing in exons 3 and 9 and that exons 3b and 3c are common in the brain and third thoracic ganglion. When expressed in Xenopus laevis oocytes, the three functional exon 3 variants (3a, b, c) all had similar EC50 values for l-glutamate and ivermectin (IVM); however, Imax (the maximum l-glutamate- and IVM-induced response of the channels at saturating concentrations) differed strikingly between variants, with the 3c variant showing the largest l-glutamate- and IVM-induced responses. By contrast, a partial deletion detected in exon 9 had a much smaller impact on l-glutamate and IVM actions. Binding assays using [(3)H]IVM indicate that diversity in IVM responses among the GluCl variants is mainly due to the impact on channel assembly, altering receptor cell surface numbers. GluCl variants expressed in HEK293 cells show that structural differences influenced Bmax but not Kd values of [(3)H]IVM. Domain swapping and site-directed mutagenesis identified four amino acids in exon 3c as hot spots determining the highest amplitude of the l-glutamate and IVM responses. Modeling the GluCl 3a and 3c variants suggested that three of the four amino acids contribute to intersubunit contacts, whereas the other interacts with the TM2-TM3 linker, influencing the receptor response.

A single amino acid polymorphism in the Drosophila melanogaster Dα1 (ALS) subunit enhances neonicotinoid efficacy at Dα1-chicken ß2 hybrid nicotinic acetylcholine receptor expressed in Xenopus laevis oocytes.

Polymorphisms are sometimes observed in native insect nicotinic acetylcholine receptor (nAChR) subunits, which are important insecticide targets, yet little is known of their impact on insecticide actions. Here we investigated the effects of a polymorphism involving the substitution of histidine108 by leucine in the Drosophila melanogaster Dα1 subunit on the agonist actions of the neurotransmitter acetylcholine (ACh) and two commercial neonicotinoid insecticides (imidacloprid and clothianidin). There was no significant impact of the H108L substitution on either the ACh EC50, the concentration leading to a half maximal ACh response, or the maximum current amplitude in response at 10 µM ACh, of the Dα1-chicken ß2 nAChR expressed in Xenopus laevis oocytes. However, the response amplitudes to imidacloprid and clothianidin were significantly enhanced, indicating a role of His108 in the selective interactions of Dα1 with these neonicotinoids.

The systematic structure-activity relationship to predict how flavones bind to human androgen receptor for their antagonistic activity.

Although flavones act as potent androgen receptor (AR) antagonists, it remains unclear how flavones interact with AR. The aim of this in silico study was to investigate the molecular recognition processes of newly synthesized 5,4'-difluoroflavone with the highest activity (IC50 value=0.19 µM) in the AR-ligand binding domain (AR-LBD). The results demonstrated that at its 4'-position of 5,4'-difluoroflavone the substituents may face Arg752 and that in AR-LBD, the submolecular bulk of substituents is unfavorable for AR antagonists and the negative electrostatic interaction site prefers the stronger hydrogen bond capability of substituents of AR antagonists. The prediction model is a valuable tool for designing a novel AR antagonist.

Synthesis of imidacloprid derivatives with a chiral alkylated imidazolidine ring and evaluation of their insecticidal activity and affinity to the nicotinic acetylcholine receptor.

A series of imidacloprid (IMI) derivatives with an alkylated imidazolidine ring were asymmetrically synthesized to evaluate their insecticidal activity against adult female housefly, Musca domestica, and affinity to the nicotinic acetylcholine receptor of the flies. The bulkier the alkyl group, the lower was the receptor affinity, but the derivatives methylated and ethylated at the R-5-position of the imidazolidine ring were equipotent to the unsubstituted compound. Quantitative structure-activity relationship (QSAR) analysis of the receptor affinity demonstrated that the introduction of a substituent into the imidazolidine ring was fundamentally disadvantageous, but the introduction of a substituent at the R-5-position was permissible in the case of its small size. The binding model of the synthesized derivatives with the receptor supported the QSAR analysis, indicating the existence of space for a short alkyl group around the R-5-position in the ligand-binding site. In addition, positive correlation was observed between the insecticidal activity and receptor affinity, suggesting that the receptor affinity was the primary factor in influencing the insecticidal activity even if the imidazolidine ring was modified.

Larvicidal activity of (-)-dihydroguaiaretic acid derivatives against Culex pipiens.

The larvicidal activity against Culex pipiens of all stereoisomers of dihydroguaiaretic acid (DGA) and secoisolariciresinol was measured, and these DGAs were found to be potent. Sixteen (-)-DGA derivatives were then newly synthesized to analyze their structure-activity relationship. Two derivatives monohydroxylated at the 3- or 4-position of the 7-phenyl group of DGA induced acute paralytic activity in the mosquitoes. Derivatives with several hydroxyl groups had lower activity than the natural compound, suggesting that hydrophobicity was probably an important factor for their insecticidal activity.

Effect of essential oils, such as raspberry ketone and its derivatives, on antiandrogenic activity based on in vitro reporter gene assay.

The effect of essential oils, such as raspberry ketone, on androgen (AR) receptor was investigated using a MDA-kb2 human breast cancer cell line for predicting potential AR activity. Among them, eugenol had the highest AR antagonistic activity with its IC(50) value of 19 microM. Raspberry ketone, which has threefold higher anti-obese activity than that of capsaicin, also had AR antagonist activity with its IC(50) value of 252 microM. Based on these findings, a more precise CoMFA model was proposed as follows: pIC(50) [log (1/IC(50))]=3.77+[CoMFA field terms] (n=39, s=0.249, r(2)=0.834, s(cv)=0.507, q(2)=0.311 (three components).

Diverse actions and target-site selectivity of neonicotinoids: structural insights.

The nicotinic acetylcholine receptors (nAChRs) are targets for human and veterinary medicines as well as insecticides. Subtype-selectivity among the diverse nAChR family members is important for medicines targeting particular disorders, and pest-insect selectivity is essential for the development of safer, environmentally acceptable insecticides. Neonicotinoid insecticides selectively targeting insect nAChRs have important applications in crop protection and animal health. Members of this class exhibit strikingly diverse actions on their nAChR targets. Here we review the chemistry and diverse actions of neonicotinoids on insect and mammalian nAChRs. Electrophysiological studies on native nAChRs and on wild-type and mutagenized recombinant nAChRs have shown that basic residues particular to loop D of insect nAChRs are likely to interact electrostatically with the nitro group of neonicotinoids. In 2008, the crystal structures were published showing neonicotinoids docking into the acetylcholine binding site of molluscan acetylcholine binding proteins with homology to the ligand binding domain (LBD) of nAChRs. The crystal structures showed that 1) glutamine in loop D, corresponding to the basic residues of insect nAChRs, hydrogen bonds with the NO(2) group of imidacloprid and 2) neonicotinoid-unique stacking and CH-pi bonds at the LBD. A neonicotinoid-resistant strain obtained by laboratory-screening has been found to result from target site mutations, and possible reasons for this are also suggested by the crystal structures. The prospects of designing neonicotinoids that are safe not only for mammals but also for beneficial insects such as honey bees (Apis mellifera) are discussed in terms of interactions with non-alpha nAChR subunits.

Combined roles of loops C and D in the interactions of a neonicotinoid insecticide imidacloprid with the alpha4beta2 nicotinic acetylcholine receptor.

Neonicotinoid insecticides are widely used for crop protection based on their selective actions on insect nicotinic acetylcholine receptors (insect nAChRs). Loops C and D in insect nAChRs have been shown to possess structural features favorable for neonicotinoid-nAChR interactions. However, it remains to be resolved whether such features serve either co-operatively, or independently, to enhance neonicotinoid sensitivity of nAChRs. We therefore examined using voltage-clamp electrophysiology the effects on the response to imidacloprid of combinatorial substitutions of residues in loops C and D of the chicken alpha4beta2 nAChR by those present in insect nAChRs. The E219P mutation in loop C of the alpha4 subunit resulted in enhanced responses to imidacloprid of alpha4beta2, whereas E219S and E219T mutations barely influenced its actions. On the other hand, mutations in loop D (T77R; E79V and T77N; E79R) alone shifted the imidacloprid concentration-response curve to the left (lower concentrations). Interestingly, all three mutations did, however, further enhance the agonist efficacy of imidacloprid when combined with the mutations in loop D. Such synergistic effects of the two loops on the interactions with imidaclprid were observed irrespective of subunit stoichiometry. Computational modeling of the ligand binding domain of the wild-type and mutant alpha4beta2 nAChRs using the crystal structure of the acetylcholine binding protein from Lymnaea stagnalis also indicated that interactions with loop F of loops C and D may contribute to determining the response to imidacloprid.

Effect of flavonoids on androgen and glucocorticoid receptors based on in vitro reporter gene assay.

The effect of 32 flavonoids on androgen (AR) and glucocorticoid receptors (GR) was investigated using an MDA-kb2 human breast cancer cell line to predict potential AR and GR activities. Among them, 5-hydroxyflavone (7) had the highest AR antagonistic activity with an IC(50) value of 0.3 microM, whereas 6-methoxyflavone (11) had the highest induced luciferase activity with an EC(150) value of 0.7 microM. Genistein (2) and daizein (1) showed a sufficient increase of luciferase activities as their concentrations increased with EC(150) values of 4.4 and 10.1 microM, respectively. These findings provide evidence of a fundamental property of their structure-activity relationship with AR and/or GR.

QSAR application for the prediction of compound permeability with in silico descriptors in practical use.

The parallel artificial membrane permeation assay (PAMPA) was developed as a model for the prediction of transcellular permeation in the process of drug absorption. In our previous report, it was revealed that PAMPA permeability is governed by log P, pK (a), and the hydrogen-bonding ability of compounds. In order to construct a new filtering method for selecting informative compounds from the whole combinatorial library, this study tried to predict PAMPA permeability with in silico descriptors. Log P, pK(a), and polar surface areas (PSA) as a hydrogen-bonding descriptor were calculated by commercially available or free-accessible web programs. Five-fold cross-validations and conventional regression analyses were examined with the training set for the entire 81 combinations with nine log P, three pK(a) and three PSA descriptors. By comparison of statistical indices, four equations were selected and then the model with the best combination of in silico descriptors was determined based on the external validation. The PAMPA prediction equation obtained in this report could be applied for the prediction of both Caco-2 cell permeability and human intestinal absorption of mainly passively-transported drugs.

Evaluation of hydrogen bonds of ecdysteroids in the ligand-receptor interactions using a protein modeling system.

The insect molting hormone, 20-hydroxyecdysone (20E) and its analogs (ecdysteroids) specifically bind to the ecdysone receptor. Previously, we synthesized various ecdysteroids containing the side chain moiety of ponasterone A (PonA), and measured the binding activity against Drosophila Kc cells to study the structure-activity relationship. Here we quantitatively analyzed the structure-activity relationship for the ligand binding of ecdysteroids including 20E and PonA. Since the hydrogen bonding (HB) is one of the important physicochemical properties for ligand binding to the ecdysteroid receptor, the number of possible HBs between the ligand molecule and the receptor was manually counted in the modeled ligand-receptor complex for all compounds. The construction of the ligand-receptor model was executed by the full-automatic modeling system (FAMS) in which calculation was done by simulated annealing. The binding potency of 15 ecdysteroids to Kc-cells were linearly correlated (r(2)=0.63) with the number of HBs which are observed between ligand and receptor molecule. Contribution of steric and electrostatic effects on the ligand-receptor binding was also examined using a three-dimensional quantitative structure-activity relationship (3-D QSAR), comparative molecular field analysis (CoMFA).