Exposition and reactivity optimization to predict sites of metabolism in chemicals

New Insights into Bitopic Orthosteric/Allosteric Ligands of Cannabinoid Receptor Type 2

Very recently, we have developed a new generation of ligands targeting the cannabinoid receptor type 2 (CB2R), namely JR compounds, which combine the pharmacophoric portion of the CB2R positive allosteric modulator (PAM), EC21a, with that of the CB2R selective orthosteric agonist LV62, both synthesized in our laboratories. The functional examination enabled us to identify JR14a, JR22a, and JR64a as the most promising compounds of the series. In the current study, we focused on the assessment of the bitopic (dualsteric) nature of these three compounds. Experiments in cAMP assays highlighted that only JR22a behaves as a CB2R bitopic (dualsteric) ligand. In parallel, computational studies helped us to clarify the binding mode of these three compounds at CB2R, confirming the bitopic (dualsteric) nature of JR22a. Finally, the potential of JR22a to prevent neuroinflammation was investigated on a human microglial cell inflammatory model.

In silico and in vitro ADME-Tox analysis and in vivo pharmacokinetic study of representative pan-PDE inhibitors from the group of 7,8-disubstituted derivatives of 1,3-dimethyl-7H-purine-2,6-dione

Phosphodiesterase (PDE) inhibitors represent a wide class of chemically different compounds that have been extensively studied in recent years. Their anti-inflammatory and anti-fibrotic effects are particularly desirable in the treatment of chronic respiratory diseases, including asthma and chronic obstructive pulmonary disease (COPD). Due to diversified expression of individual PDEs within cells and/or tissues as well as PDE signaling compartmentalization, pan-PDE inhibitors (compounds capable of simultaneously blocking various PDE subtypes) are of particular interest. Recently, a large group of 7,8-disubstituted derivatives of 1,3-dimethyl-7H-purine-2,6-dione (theophylline) was designed and synthesized. These compounds were characterized as potent pan-PDE inhibitors and their prominent anti-inflammatory and anti-fibrotic activity in vitro has been proved. Herein, we investigated a general in vitro safety profile and pharmacokinetic characteristics of two leading compounds from this group: a representative compound with N'-benzylidenebutanehydrazide moiety (38) and a representative derivative containing N-phenylbutanamide fragment (145). Both tested pan-PDE inhibitors revealed no cytotoxic, mutagenic, and genotoxic activity in vitro, showed moderate metabolic stability in mouse and human liver microsomes, as well as fell into the low or medium permeation category. Additionally, 38 and 145 revealed a lack of interaction with adenosine receptors, including A1, A2A, and A2B. Pharmacokinetic analysis revealed that both tested 7,8-disubstituted derivatives of 1,3-dimethyl-7H-purine-2,6-dione were effectively absorbed from the peritoneal cavity. Simultaneously, they were extensively distributed to mouse lungs and after intraperitoneal (i.p.) administration were detected in bronchoalveolar lavage fluid. These findings provide evidence that investigated compounds represent a new drug candidates with a favorable in vitro safety profile and satisfactory pharmacokinetic properties after a single i.p. administration. As the next step, further pharmacokinetic studies after multiple i.p. and p.o. doses will be conducted to ensure effective 38 and 145 serum and lung concentrations for a longer period of time. In summary, 7,8-disubstituted derivatives of 1,3-dimethyl-7H-purine-2,6-dione represent a promising compounds worth testing in animal models of chronic respiratory diseases, the etiology of which involves various PDE subtypes.

Modeling the Bioactivation and Subsequent Reactivity of Drugs

Electrophilically reactive drug metabolites are implicated in many adverse drug reactions. In this mechanism-termed bioactivation-metabolic enzymes convert drugs into reactive metabolites that often conjugate to nucleophilic sites within biological macromolecules like proteins. Toxic metabolite-product adducts induce severe immune responses that can cause sometimes fatal disorders, most commonly in the form of liver injury, blood dyscrasia, or the dermatologic conditions toxic epidermal necrolysis and Stevens-Johnson syndrome. This study models four of the most common metabolic transformations that result in bioactivation: quinone formation, epoxidation, thiophene sulfur-oxidation, and nitroaromatic reduction, by synthesizing models of metabolism and reactivity. First, the metabolism models predict the formation probabilities of all possible metabolites among the pathways studied. Second, the exact structures of these metabolites are enumerated. Third, using these structures, the reactivity model predicts the reactivity of each metabolite. Finally, a feedfoward neural network converts the metabolism and reactivity predictions to a bioactivation prediction for each possible metabolite. These bioactivation predictions represent the joint probability that a metabolite forms and that this metabolite subsequently conjugates to protein or glutathione. Among molecules bioactivated by these pathways, we predicted the correct pathway with an AUC accuracy of 89.98%. Furthermore, the model predicts whether molecules will be bioactivated, distinguishing bioactivated and nonbioactivated molecules with 81.06% AUC. We applied this algorithm to withdrawn drugs. The known bioactivation pathways of alclofenac and benzbromarone were identified by the algorithm, and high probability bioactivation pathways not yet confirmed were identified for safrazine, zimelidine, and astemizole. This bioactivation model-the first of its kind that jointly considers both metabolism and reactivity-enables drug candidates to be quickly evaluated for a toxicity risk that often evades detection during preclinical trials. The XenoSite bioactivation model is available at http://swami.wustl.edu/xenosite/p/bioactivation.

Understanding the Metabolism of Proteolysis Targeting Chimeras (PROTACs): The Next Step toward Pharmaceutical Applications

Hetero-bifunctional PROteolysis TArgeting Chimeras (PROTACs) represent a new emerging class of small molecules designed to induce polyubiquitylation and proteasomal-dependent degradation of a target protein. Despite the increasing number of publications about the synthesis, biological evaluation, and mechanism of action of PROTACs, the characterization of the pharmacokinetic properties of this class of compounds is still minimal. Here, we report a study on the metabolism of a series of 40 PROTACs in cryopreserved human hepatocytes at multiple time points. Our results indicated that the metabolism of PROTACs could not be predicted from that of their constituent ligands. Their linkers' chemical nature and length resulted in playing a major role in the PROTACs' liability. A subset of compounds was also tested for metabolism by human cytochrome P450 3A4 (CYP3A4) and human aldehyde oxidase (hAOX) for more in-depth data interpretation, and both enzymes resulted in active PROTAC metabolism.

Computer-aided Drug Design Investigations for Benzothiazinone Derivatives Against Tuberculosis

BACKGROUND

Tuberculosis (Mycobacterium tuberculosis) is an infectious bacterial disease with the highest levels of mortality worldwide, presenting numerous cases of resistance. In silico studies, which elaborate chemical and biological models in computational tools and make it possible to interpret molecular characteristics, are among the methods used in the search for new drugs.

OBJECTIVE

In this perspective, our aim was to use QSAR and molecular modeling to propose possible pharmacophores from benzothiazinone derivatives.

METHODS

In this study, a set of 69 benzothiazinone derivatives, together with computational tools such as molecular descriptor analysis in chemometrics, metabolic prediction, and molecular coupling to 4 proteins: DprE1, InhA, PS, and DHFR important for the bacillus were investigated.

RESULTS

The chemometric model computed in the Volsurf+ program presented good predictive values for both amphiphilicity and molecular volume. These are essential for biological activity. Metabolites from the cytochrome isoforms CYP3A4 and 2D6 interactions revealed coupling divergences which, noting that the metabolites did not present changes to the QSAR proposed pharmacophore structures, may be due to the reaction medium and existing differences in the benzothiazinone structures. Similarly, molecular docking with the four TB enzymes presented good interactions for the more active compounds. The fragments found using QSAR (being essential for biological activity) also presented as being essential for ligand-protein site interactions.

CONCLUSION

From the benzothiazinone derivative series evaluated, compound 11026134 presented the best profile in all study analyses, noting that the trifluoromethyl, nitro group, and piperazine fragment with aliphatic hydrocarbon groups are likely pharmacophores for the benzothiazinones studied.

TPT sulfonate, a single, oral dose schistosomicidal prodrug: In vivo efficacy, disposition and metabolic profiling

Treatment of schistosomiasis relies precariously on just one drug, praziquantel (PZQ). In the search for alternatives, 15 S-[2-(alkylamino)alkane] thiosulfuric acids were obtained from a previous research program and profiled in mice for efficacy against both mature (>42-day-old) and juvenile (21-day-old) Schistosoma mansoni using a screening dose of 100 mg/kg PO QDx4. One compound, S-[2-(tert-butylamino)-1-phenylethane] thiosulfuric acid (TPT sulfonate), was the most effective by decreasing female and male worm burdens by ≥ 90% and ≥46% (mature), and ≥89% and ≥79% (juvenile), respectively. In contrast, PZQ decreased mature female and male worm burdens by 95% and 94%, respectively, but was ineffective against juvenile stages. Against 7-day-old lung-stage worms, TPT sulfonate was only effective at twice the dose decreasing female and male burdens by 95 and 80%, respectively. Single oral doses at 400 and/or 600 mg/kg across various developmental time-points (1-, 7-, 15-, 21- and/or 42 day-old) were consistent with the QD x4 data; efficacy was strongest once the parasites had completed lung migration, and female and male burdens were decreased by at least 90% and 80%, respectively. In vitro, TPT sulfonate is inactive against the parasite suggesting a pro-drug mechanism of action. In mice, TPT sulfonate is fully absorbed and subject to rapid, non-CYP-mediated, first-pass metabolism that is initiated by desulfation and yields a series of metabolites. The initially-formed free thiol-containing metabolite, termed TP thiol, was chemically synthesized; it dose-dependently decreased S. mansoni and Schistosoma haematobium motility in vitro. Also, when administered as a single 50 mg/kg IP dose, TP thiol decreased 33-day-old S. mansoni female and male burdens by 35% and 44%, with less severe organomegaly. Overall, TPT sulfonate's efficacy profile is competitive with that of PZQ. Also, the characterization of a parasiticidal metabolite facilitates an understanding and improvement of the chemistry, and identification of the mechanism of action and/or target.

•

TPT sulfonate provides single dose efficacy that is competitive with the current drug, praziquantel.

•

TPT sulfonate must be biotransformed to be active.

•

TPT sulfonate is fully absorbed and subject to rapid, non-CYP-mediated, first-pass metabolism.

•

One of the key metabolites, TP thiol, is anti-schistosomal in vitro and in vivo.

Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2

Vitamin K₁ (VK1) plays an important role in the modulation of bleeding disorders. It has been reported that ω-hydroxylation on the VK1 aliphatic chain is catalyzed by cytochrome P450 4F2 (CYP4F2), an enzyme responsible for the metabolism of eicosanoids. However, the mechanism of VK1 ω-hydroxylation by CYP4F2 has not been disclosed. In this study, we employed a combination of quantum mechanism (QM) calculations, homology modeling, molecular docking, molecular dynamics (MD) simulations, and combined quantum mechanism/molecular mechanism (QM/MM) calculations to investigate the metabolism profile of VK1 ω-hydroxylation. QM calculations based on the truncated VK1 model show that the energy barrier for ω-hydroxylation is about 6-25 kJ/mol higher than those at other potential sites of metabolism. However, results from the MD simulations indicate that hydroxylation at the ω-site is more favorable than at the other potential sites, which is in accordance with the experimental observation. The evaluation of MD simulations was further endorsed by the QM/MM calculation results. Our studies thus suggest that the active site residues of CYP4F2 play a determinant role in the ω-hydroxylation. Our results provide structural insights into the mechanism of VK1 ω-hydroxylation by CYP4F2 at the atomistic level and are helpful not only for characterizing the CYP4F2 functions but also for looking into the ω-hydroxylation mediated by other CYP4 enzymes.

The current limits in virtual screening and property prediction

Beyond finding inhibitors that show high binding affinity to the respective target, there is the challenge of optimizing their properties with respect to metabolic and toxicological issues, as well as further off-target effects. To reduce the experimental effort of synthesizing and testing actual substances in corresponding assays, virtual screening has become an indispensable toolbox in preclinical development. The scope of application covers the prediction of molecular properties including solubility, metabolic liability and binding to antitargets, such as the hERG channel. Furthermore, prediction of binding sites and drugable targets are emerging aspects of virtual screening. Issues involved with the currently applied computational models including machine learning algorithms are outlined, such as limitations to the accuracy of prediction and overfitting.

Biotransformation of 4-fluoro-N-(1-{2-[(propan-2-yl)phenoxy]ethyl}-8-azabicyclo[3.2.1]octan-3-yl)-benzenesulfonamide, a novel potent 5-HT7 receptor antagonist with antidepressant-like and anxiolytic properties: In vitro and in silico approach

The aim of the study was to investigate the metabolism of 4-fluoro-N-(1-{2-[(propan-2-yl)phenoxy]ethyl}-8-azabicyclo[3.2.1]octan-3-yl)-benzenesulfonamide (PZ-1150), a novel 5-HT₇ receptor antagonist with antidepressant-like and anxiolytic properties, by the following three ways: in vitro with microsomes; in vitro employing Cunninghamella echinulata, and in silico using MetaSite. Biotransformation of PZ-1150 with microsomes resulted in five metabolites, while transformation with C. echinulata afforded two metabolites. In both models, the predominant metabolite occurred due to hydroxylation of benzene ring. In silico data coincide with in vitro experiments, as three MetaSite metabolites matched compounds identified in microsomal samples. In human liver microsomes PZ-1150 exhibited in vitro half-life of 64 min, with microsomal intrinsic clearance of 54.1 μL/min/mg and intrinsic clearance of 48.7 mL/min/kg. Therefore, PZ-1150 is predicted to be a high-clearance agent. The study demonstrated the applicability of using microsomal model coupled with microbial model to elucidate the metabolic pathways of compounds and comparison with in silico metabolite predictions.

Plant organ cultures as masked mycotoxin biofactories: Deciphering the fate of zearalenone in micropropagated durum wheat roots and leaves

"Masked mycotoxins" senso strictu are conjugates of mycotoxins resulting from metabolic pathways activated by the interplay between pathogenic fungi and infected plants. Zearalenone, an estrogenic mycotoxin produced by Fusarium spp, was the first masked mycotoxin ever described in the literature, but its biotransformation has been studied to a lesser extent if compared to other compounds such as deoxynivalenol. We presented herein the first application of organ and tissue culture techniques to study the metabolic fate of zearalenone in durum wheat, using an untargeted HR-LCMS approach. A complete, quick absorption of zearalenone by uninfected plant organs was noticed, and its biotransformation into a large spectrum of phase I and phase II metabolites has been depicted. Therefore, wheat organ tissue cultures can be effectively used as a biocatalytic tool for the production of masked mycotoxins, as well as a replicable model for the investigation of the interplay between mycotoxins and wheat physiology.

Comparative Analysis of the Neurochemical Profile and MAO Inhibition Properties of N-(Furan-2-ylmethyl)-N-methylprop-2-yn-1-amine

The regulation of brain monoamine levels is paramount for cognitive functions, and the monoamine oxidase (MAO A and B) enzymes play a central role in these processes. The aim of this study was to evaluate whether the procognitive properties exerted by propargylamine N-(furan-2-ylmethyl)-N-methylprop-2-yn-1-amine (F2MPA) are related to changes in monoamine content via MAO inhibition. In vivo microdialysis and ex vivo amine metabolite measurement demonstrated region-specific alterations in monoamine metabolism that differ from both of the classic MAO A and MAO B inhibitors, clorgyline and l-deprenyl, respectively. Although all the inhibitors (1 and 4 mg/kg) increased cortical serotonin tissue content, only F2MPA increased the levels of cortical noradrenaline. In the striatum, clorgyline (1 mg/kg), but not F2MPA (1 mg/kg), reduced extracellular levels of dopamine metabolites at rest or stimulated by the intrastriatal application of the MAO substrate 3-methoxytyramine. In vitro, F2MPA exhibited a low affinity toward MAO B and MAO A. Nonetheless, it modified the B form of MAO, forming a flavin adduct structurally similar to that with deprenyl. F2MPA was rapidly metabolized in the presence of rat but not human microsomes, producing a hydroxylated derivative. In conclusion, the effect of F2MPA on cognition may arise from monoaminergic changes in the cortex, but the role of MAO in this process is likely to be negligible, consistent with the poor affinity of F2MPA for MAO.

Structure-metabolism relationships in human-AOX: Chemical insights from a large database of aza-aromatic and amide compounds

Aldehyde oxidase (AOX) is a metabolic enzyme catalyzing the oxidation of aldehyde and aza-aromatic compounds and the hydrolysis of amides, moieties frequently shared by the majority of drugs. Despite its key role in human metabolism, to date only fragmentary information about the chemical features responsible for AOX susceptibility are reported and only "very local" structure-metabolism relationships based on a small number of similar compounds have been developed. This study reports a more comprehensive coverage of the chemical space of structures with a high risk of AOX phase I metabolism in humans. More than 270 compounds were studied to identify the site of metabolism and the metabolite(s). Both electronic [supported by density functional theory (DFT) calculations] and exposure effects were considered when rationalizing the structure-metabolism relationship.

Deep Learning to Predict the Formation of Quinone Species in Drug Metabolism

Many adverse drug reactions are thought to be caused by electrophilically reactive drug metabolites that conjugate to nucleophilic sites within DNA and proteins, causing cancer or toxic immune responses. Quinone species, including quinone-imines, quinone-methides, and imine-methides, are electrophilic Michael acceptors that are often highly reactive and comprise over 40% of all known reactive metabolites. Quinone metabolites are created by cytochromes P450 and peroxidases. For example, cytochromes P450 oxidize acetaminophen to N-acetyl-p-benzoquinone imine, which is electrophilically reactive and covalently binds to nucleophilic sites within proteins. This reactive quinone metabolite elicits a toxic immune response when acetaminophen exceeds a safe dose. Using a deep learning approach, this study reports the first published method for predicting quinone formation: the formation of a quinone species by metabolic oxidation. We model both one- and two-step quinone formation, enabling accurate quinone formation predictions in nonobvious cases. We predict atom pairs that form quinones with an AUC accuracy of 97.6%, and we identify molecules that form quinones with 88.2% AUC. By modeling the formation of quinones, one of the most common types of reactive metabolites, our method provides a rapid screening tool for a key drug toxicity risk. The XenoSite quinone formation model is available at http://swami.wustl.edu/xenosite/p/quinone .

Modeling Reactivity to Biological Macromolecules with a Deep Multitask Network

Most small-molecule drug candidates fail before entering the market, frequently because of unexpected toxicity. Often, toxicity is detected only late in drug development, because many types of toxicities, especially idiosyncratic adverse drug reactions (IADRs), are particularly hard to predict and detect. Moreover, drug-induced liver injury (DILI) is the most frequent reason drugs are withdrawn from the market and causes 50% of acute liver failure cases in the United States. A common mechanism often underlies many types of drug toxicities, including both DILI and IADRs. Drugs are bioactivated by drug-metabolizing enzymes into reactive metabolites, which then conjugate to sites in proteins or DNA to form adducts. DNA adducts are often mutagenic and may alter the reading and copying of genes and their regulatory elements, causing gene dysregulation and even triggering cancer. Similarly, protein adducts can disrupt their normal biological functions and induce harmful immune responses. Unfortunately, reactive metabolites are not reliably detected by experiments, and it is also expensive to test drug candidates for potential to form DNA or protein adducts during the early stages of drug development. In contrast, computational methods have the potential to quickly screen for covalent binding potential, thereby flagging problematic molecules and reducing the total number of necessary experiments. Here, we train a deep convolution neural network-the XenoSite reactivity model-using literature data to accurately predict both sites and probability of reactivity for molecules with glutathione, cyanide, protein, and DNA. On the site level, cross-validated predictions had area under the curve (AUC) performances of 89.8% for DNA and 94.4% for protein. Furthermore, the model separated molecules electrophilically reactive with DNA and protein from nonreactive molecules with cross-validated AUC performances of 78.7% and 79.8%, respectively. On both the site- and molecule-level, the model's performances significantly outperformed reactivity indices derived from quantum simulations that are reported in the literature. Moreover, we developed and applied a selectivity score to assess preferential reactions with the macromolecules as opposed to the common screening traps. For the entire data set of 2803 molecules, this approach yielded totals of 257 (9.2%) and 227 (8.1%) molecules predicted to be reactive only with DNA and protein, respectively, and hence those that would be missed by standard reactivity screening experiments. Site of reactivity data is an underutilized resource that can be used to not only predict if molecules are reactive, but also show where they might be modified to reduce toxicity while retaining efficacy. The XenoSite reactivity model is available at http://swami.wustl.edu/xenosite/p/reactivity.

Discovery of in vitro antitubercular agents through in silico ligand-based approaches

The development of new anti-tubercular agents represents a constant challenge mostly due to the insurgency of resistance to the currently available drugs. In this study, a set of 60 molecules were selected by screening the Asinex and the ZINC collections and an in house library by means of in silico ligand-based approaches. Biological assays in Mycobacterium tuberculosis H37Ra ATCC 25177 strain highlighted (±)-1-(4-chlorophenyl)-2-(1H-imidazol-1-yl)ethyl-4-(3,4-dichlorophenyl)piperazine-1-carboxylate (5i) and 3-(4-chlorophenyl)-5-(2,4-dimethylpyrimidin-5-yl)-2-methylpyrazolo[1.5-a]pyrimidin-7(4H)-one (42) as the most potent compounds, having a Minimum Inhibitory Concentration (MIC) of 4 and 2 μg/mL respectively. These molecules represent a good starting point for further optimization of effective anti-TB agents.

Effects of different routes of application on ethylenediurea persistence in tobacco leaves

Ethylenediurea (EDU) is a common research tool for investigating ozone impacts on vegetation, although the role of different application routes (foliar spray vs soil drench) on EDU persistence in the leaves is unknown. We quantified EDU concentrations in leaves of the O3-sensitive Bel-W3 cultivar of tobacco treated with EDU as either foliar spray or soil drench. Foliar EDU concentrations were measured by Q-TOF LC/MS. When EDU was applied as foliar spray, 1 h was enough for reaching a measurable concentration within the leaf. EDU concentration increased over the 21-day period when the leaf was not washed after the application (treatment #1), while it decreased when the leaf was washed after the application (treatment #2). These results suggest that: a) dry deposition of EDU onto the leaf surface was gradually absorbed into the unwashed leaf, although the mechanisms of such uptake were unclear; b) concentration of EDU was decreased quickly (-35%) during the first 24 h from application and more slowly during the following three days (-20%) in the washed leaves. Degradation did not involve enzymatic reactions and was not affected by the presence of ROS. When EDU was applied as soil drench, foliar concentrations increased over time, likely due to adsorption onto soil organic matter and gradual re-solubilization by irrigation water. An analysis of EDU concentration in protoplast and intercellular washing fluid showed that EDU did not enter the cells, but was retained in the apoplast only. Possible implications of EDU in the apoplast and recommendations for EDU application are discussed.

Understanding Oxadiazolothiazinone Biological Properties: Negative Inotropic Activity versus Cytochrome P450-Mediated Metabolism

We present a series of oxadiazolothiazinones, selective inotropic agents on isolated cardiac tissues, devoid of chronotropy and vasorelaxant activity. Functional and binding data for the precursor of the series (compound 1) let us hypothesize LTCC blocking activity and the existence of a recognition site specific for this scaffold. We synthesized and tested 22 new derivatives: introducing a para-methoxyphenyl at C-8 led to compound 12 (EC50 = 0.022 μM), twice as potent as its para-bromo analogue (1). For 10 analogues, we extended the characterization of the biological properties by including the assessment of metabolic stability in human liver microsomes and cytochrome P450 inhibition potential. We observed that the methoxy group led to active compounds with low metabolic stability and high CYP inhibition, whereas the protective effect of bromine resulted in enhanced metabolic stability and reduced CYP inhibition. Thus, we identified two para-bromo benzothiazino-analogues as candidates for further studies.

Pyridofuropyrrolo[1,2-a]pyrimidines and pyridofuropyrimido[1,2-a]azepines: new chemical entities (NCE) with anticonvulsive and psychotropic properties

A series of pyridofuropyrrolo[1,2-a]pyrimidines4and of pyridofuropyrimido[1,2-a]azepines5having as precursors a some new condensed furo[2,3-b]pyridines3were synthesized and evaluated for their anticonvulsive and psychotropic properties.

Pyridofuropyrrolo[1,2-a]pyrimidines and pyridofuropyrimido[1,2-a]azepines: new chemical entities (NCE) with anticonvulsive and psychotropic properties

A series of pyridofuropyrrolo[1,2-a]pyrimidines4and pyridofuropyrimido[1,2-a]azepines5, with new condensed furo[2,3-b]pyridines3as precursors, were synthesized and evaluated for their anticonvulsive and psychotropic properties.

Modeling Epoxidation of Drug-like Molecules with a Deep Machine Learning Network

Drug toxicity is frequently caused by electrophilic reactive metabolites that covalently bind to proteins. Epoxides comprise a large class of three-membered cyclic ethers. These molecules are electrophilic and typically highly reactive due to ring tension and polarized carbon-oxygen bonds. Epoxides are metabolites often formed by cytochromes P450 acting on aromatic or double bonds. The specific location on a molecule that undergoes epoxidation is its site of epoxidation (SOE). Identifying a molecule's SOE can aid in interpreting adverse events related to reactive metabolites and direct modification to prevent epoxidation for safer drugs. This study utilized a database of 702 epoxidation reactions to build a model that accurately predicted sites of epoxidation. The foundation for this model was an algorithm originally designed to model sites of cytochromes P450 metabolism (called XenoSite) that was recently applied to model the intrinsic reactivity of diverse molecules with glutathione. This modeling algorithm systematically and quantitatively summarizes the knowledge from hundreds of epoxidation reactions with a deep convolution network. This network makes predictions at both an atom and molecule level. The final epoxidation model constructed with this approach identified SOEs with 94.9% area under the curve (AUC) performance and separated epoxidized and non-epoxidized molecules with 79.3% AUC. Moreover, within epoxidized molecules, the model separated aromatic or double bond SOEs from all other aromatic or double bonds with AUCs of 92.5% and 95.1%, respectively. Finally, the model separated SOEs from sites of sp(2) hydroxylation with 83.2% AUC. Our model is the first of its kind and may be useful for the development of safer drugs. The epoxidation model is available at http://swami.wustl.edu/xenosite.

Modeling of interactions between xenobiotics and cytochrome P450 (CYP) enzymes

The adverse effects to humans and environment of only few chemicals are well known. Absorption, distribution, metabolism, and excretion (ADME) are the steps of pharmaco/toxicokinetics that determine the internal dose of chemicals to which the organism is exposed. Of all the xenobiotic-metabolizing enzymes, the cytochrome P450 (CYP) enzymes are the most important due to their abundance and versatility. Reactions catalyzed by CYPs usually turn xenobiotics to harmless and excretable metabolites, but sometimes an innocuous xenobiotic is transformed into a toxic metabolite. Data on ADME and toxicity properties of compounds are increasingly generated using in vitro and modeling (in silico) tools. Both physics-based and empirical modeling approaches are used. Numerous ligand-based and target-based as well as combined modeling methods have been employed to evaluate determinants of CYP ligand binding as well as predicting sites of metabolism and inhibition characteristics of test molecules. In silico prediction of CYP–ligand interactions have made crucial contributions in understanding (1) determinants of CYP ligand binding recognition and affinity; (2) prediction of likely metabolites from substrates; (3) prediction of inhibitors and their inhibition potency. Truly predictive models of toxic outcomes cannot be created without incorporating metabolic characteristics; in silico methods help producing such information and filling gaps in experimentally derived data. Currently modeling methods are not mature enough to replace standard in vitro and in vivo approaches, but they are already used as an important component in risk assessment of drugs and other chemicals.

Synthesis and L-type calcium channel blocking activity of new chiral oxadiazolothiazinones

Oxadiazolo[3,4-c][1,4]thiazin-3-ones are cardiovascular agents that block L-type calcium channels. Previous data of cardiac and vasorelaxant activity on guinea-pig for several derivatives indicated the two positions ortho to the thiazine's sulphur as crucial for modulating the activity; but these positions are likely susceptible to metabolic biotransformations, as indicated by in silico predictions. We designed new derivatives, and obtained three negative inotropic agents with EC50 in the low nanomolar range, more potent than all the precursors published so far. In particular, benzocondensation at the thiazine ring led to 3a (EC50 = 0.013 μM) and 3b (EC50 = 0.006 μM). Besides negative inotropy, we also observed relaxant activity on nonvascular muscle in the micromolar range. We resolved the new derivatives by chiral chromatography, and determined their absolute configuration by comparing experimental and calculated chiroptical properties (VCD, ECD and ORD): they hold the same absolute configuration-optical rotation relationship, (S)-(+)/(R)-(-). Both cardiac and nonvascular activity are majorly or mostly retained in the R-form for all the compounds, but for the nonvascular activity we observed a strong stereoselectivity for 3a, with the R-form in the nanomolar range (IC50 = 0.020 μM) and 259-fold more potent than the S-one.

Biotransformation in vitro: An essential consideration in the quantitative in vitro-to-in vivo extrapolation (QIVIVE) of toxicity data

Early consideration of the multiplicity of factors that govern the biological fate of foreign compounds in living systems is a necessary prerequisite for the quantitative in vitro-in vivo extrapolation (QIVIVE) of toxicity data. Substantial technological advances in in vitro methodologies have facilitated the study of in vitro metabolism and the further use of such data for in vivo prediction. However, extrapolation to in vivo with a comfortable degree of confidence, requires continuous progress in the field to address challenges such as e.g., in vitro evaluation of chemical-chemical interactions, accounting for individual variability but also analytical challenges for ensuring sensitive measurement technologies. This paper discusses the current status of in vitro metabolism studies for QIVIVE extrapolation, serving today's hazard and risk assessment needs. A short overview of the methodologies for in vitro metabolism studies is given. Furthermore, recommendations for priority research and other activities are provided to ensure further widespread uptake of in vitro metabolism methods in 21st century toxicology. The need for more streamlined and explicitly described integrated approaches to reflect the physiology and the related dynamic and kinetic processes of the human body is highlighted i.e., using in vitro data in combination with in silico approaches.

Flavin monooxygenase metabolism: why medicinal chemists should matter

FMO enzymes (FMOs) play a key role in the processes of detoxification and/or bioactivation of specific pharmaceuticals and xenobiotics bearing nucleophilic centers. The N-oxide and S-oxide metabolites produced by FMOs are often active metabolites. The FMOs are more active than cytochromes in the brain and work in tandem with CYP3A4 in the liver. FMOs might reduce the risk of phospholipidosis of CAD-like drugs, although some FMOs metabolites seem to be neurotoxic and hepatotoxic. However, in silico methods for FMO metabolism prediction are not yet available. This paper reports, for the first time, a substrate-specificity and catalytic-activity model for FMO3, the most relevant isoform of the FMOs in humans. The application of this model to a series of compounds with unknown FMO metabolism is also reported. The model has also been very useful to design compounds with optimal clearance and in finding erroneous literature data, particularly cases in which substances have been reported to be FMO3 substrates when, in reality, the experimentally validated in silico model correctly predicts that they are not.

Model-based estimates of the effects of efavirenz on bedaquiline pharmacokinetics and suggested dose adjustments for patients coinfected with HIV and tuberculosis

Safe, effective concomitant treatment regimens for tuberculosis (TB) and HIV infection are urgently needed. Bedaquiline (BDQ) is a promising new anti-TB drug, and efavirenz (EFV) is a commonly used antiretroviral. Due to EFV's induction of cytochrome P450 3A4, the metabolic enzyme responsible for BDQ biotransformation, the drugs are expected to interact. Based on data from a phase I, single-dose pharmacokinetic study, a nonlinear mixed-effects model characterizing BDQ pharmacokinetics and interaction with multiple-dose EFV was developed. BDQ pharmacokinetics were best described by a 3-compartment disposition model with absorption through a dynamic transit compartment model. Metabolites M2 and M3 were described by 2-compartment models with clearance of BDQ and M2, respectively, as input. Impact of induction was described as an instantaneous change in clearance 1 week after initialization of EFV treatment and estimated for all compounds. The model predicts average steady-state concentrations of BDQ and M2 to be reduced by 52% (relative standard error [RSE], 3.7%) with chronic coadministration. A range of models with alternative structural assumptions regarding onset of induction effect and fraction metabolized resulted in similar estimates of the typical reduction and did not offer a markedly better fit to data. Simulations to investigate alternative regimens mitigating the estimated interaction effect were performed. The results suggest that simple adjustments of the standard regimen during EFV coadministration can prevent reduced exposure to BDQ without increasing exposures to M2. However, exposure to M3 would increase. Evaluation in clinical trials of adjusted regimens is necessary to ensure appropriate dosing for HIV-infected TB patients on an EFV-based regimen.