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.

Differential scanning fluorimetric analysis of the amino-acid binding to taste receptor using a model receptor protein, the ligand-binding domain of fish T1r2a/T1r3.

Taste receptor type 1 (T1r) is responsible for the perception of essential nutrients, such as sugars and amino acids, and evoking sweet and umami (savory) taste sensations. T1r receptors recognize many of the taste substances at their extracellular ligand-binding domains (LBDs). In order to detect a wide array of taste substances in the environment, T1r receptors often possess broad ligand specificities. However, the entire ranges of chemical spaces and their binding characteristics to any T1rLBDs have not been extensively analyzed. In this study, we exploited the differential scanning fluorimetry (DSF) to medaka T1r2a/T1r3LBD, a current sole T1rLBD heterodimer amenable for recombinant preparation, and analyzed their thermal stabilization by adding various amino acids. The assay showed that the agonist amino acids induced thermal stabilization and shifted the melting temperatures (Tm) of the protein. An agreement between the DSF results and the previous biophysical assay was observed, suggesting that DSF can detect ligand binding at the orthosteric-binding site in T1r2a/T1r3LBD. The assay further demonstrated that most of the tested l-amino acids, but no d-amino acid, induced Tm shifts of T1r2a/T1r3LBD, indicating the broad l-amino acid specificities of the proteins probably with several different manners of recognition. The Tm shifts by each amino acid also showed a fair correlation with the responses exhibited by the full-length receptor, verifying the broad amino-acid binding profiles at the orthosteric site in LBD observed by DSF.

Obituary: Toshio Fujita, QSAR pioneer.

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

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.

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.

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.

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.

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.

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.

Functional characterisation of a nicotinic acetylcholine receptor α subunit from the brown dog tick, Rhipicephalus sanguineus.

Ticks and tick-borne diseases have a major impact on human and animal health worldwide. Current control strategies rely heavily on the use of chemical acaricides, most of which target the CNS and with increasing resistance, new drugs are urgently needed. Nicotinic acetylcholine receptors (nAChRs) are targets of highly successful insecticides. We isolated a full-length nAChR α subunit from a normalised cDNA library from the synganglion (brain) of the brown dog tick, Rhipicephalus sanguineus. Phylogenetic analysis has shown this R. sanguineus nAChR to be most similar to the insect α1 nAChR group and has been named Rsanα1. Rsanα1 is distributed in multiple tick tissues and is present across all life-stages. When expressed in Xenopus laevis oocytes Rsanα1 failed to function as a homomer, with and without the addition of either Caenorhabditis elegans resistance-to-cholinesterase (RIC)-3 or X. laevis RIC-3. When co-expressed with chicken ß2 nAChR, Rsanα1 evoked concentration-dependent, inward currents in response to acetylcholine (ACh) and showed sensitivity to nicotine (100 µM) and choline (100 µM). Rsanα1/ß2 was insensitive to both imidacloprid (100 µM) and spinosad (100 µM). The unreliable expression of Rsanα1 in vitro suggests that additional subunits or chaperone proteins may be required for more robust expression. This study enhances our understanding of nAChRs in arachnids and may provide a basis for further studies on the interaction of compounds with the tick nAChR as part of a discovery process for novel acaricides.

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.

Synthesis and structure-activity relationship study of antimicrotubule agents phenylahistin derivatives with a didehydropiperazine-2,5-dione structure.

Plinabulin (11, NPI-2358) is a potent microtubule-targeting agent derived from the natural diketopiperazine "phenylahistin" (1) with a colchicine-like tubulin depolymerization activity. Compound 11 was recently developed as VDA and is now under phase II clinical trials as an anticancer drug. To develop more potent antimicrotubule and cytotoxic derivatives based on the didehydro-DKP skeleton, we performed further modification on the tert-butyl or phenyl groups of 11, and evaluated their cytotoxic and tubulin-binding activities. In the SAR study, we developed more potent derivatives 33 with 2,5-difluorophenyl and 50 with a benzophenone in place of the phenyl group. The anti-HuVEC activity of 33 and 50 exhibited a lowest effective concentration of 2 and 1 nM for microtubule depolymerization, respectively. The values of 33 and 50 were 5 and 10 times more potent than that of CA-4, respectively. These derivatives could be a valuable second-generation derivative with both vascular disrupting and cytotoxic activities.

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.

Importance of physicochemical properties for the design of new pesticides.

The physicochemical properties of candidate compounds play important roles in the design of new pesticides. Pesticides must be absorbed by pests, be transported to the target site, and then interact with proteins. Hydrophobicity is very important for these processes. Log P, where P is the partition coefficient in the 1-octanol/water system, is commonly used as a hydrophobic descriptor and correlates with membrane permeation and transport. It was recently reported that permeability by the parallel artificial membrane permeation assay (PAMPA) could be used to predict human oral absorption of passively transported compounds. PAMPA, which is a rapid high-throughput screening system, may be useful to predict pesticide absorption because PAMPA permeability can be calculated using log P and other parameters. Electronic and structural properties as well as hydrophobicity are important factors for protein-ligand interaction. To show the importance of physicochemical properties, the classic QSAR and CoMFA of neonicotinoids and prediction of bioavailability of pesticides in terms of membrane permeability in comparison with drugs are described.

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).

In silico prediction of human oral absorption based on QSAR analyses of PAMPA permeability.

The parallel artificial membrane permeation assay (PAMPA) was developed as a model for the prediction of transcellular permeation in the process of drug absorption. Our research group has measured the PAMPA permeability of peptide-related compounds, diverse drugs, and agrochemicals. This work led to a classical quantitative structure-activity relationship (QSAR) equation for PAMPA permeability coefficients of structurally diverse compounds based on simple physicochemical parameters such as lipophilicity at a particular pH (log P(oct) and |pKa-pH|), H-bond acceptor ability (SA(HA)), and H-bond donor ability (SA(HD)). Since the PAMPA permeability of lipophilic compounds decreased with their apparent lipophilicity due to the unstirred water layer (UWL) barrier on membrane surfaces and to membrane retention, a bilinear QSAR model was introduced to explain the permeability of a broader set of compounds using the same physicochemical parameters as those used for the linear model. We also compared PAMPA and Caco-2 cell permeability coefficients of compounds transported by various absorption mechanisms. The compounds were classified according to their absorption pathway (passively transported compounds, actively transported compounds, and compounds excreted by efflux systems) in the plot of Caco-2 vs. PAMPA permeability. Finally, based on the QSAR analyses of PAMPA permeability, an in silico prediction model of human oral absorption for possibly transported compounds was proposed, and the usefulness of the model was examined.