Application of Deep Neural Network Models in Drug Discovery Programs.

In silico driven optimization of compound properties related to pharmacokinetics, pharmacodynamics, and safety is a key requirement in modern drug discovery. Nowadays, large and harmonized datasets allow to implement deep neural networks (DNNs) as a framework for leveraging predictive models. Nevertheless, various available model architectures differ in their global applicability and performance in lead optimization projects, such as stability over time and interpretability of the results. Here, we describe and compare the value of established DNN-based methods for the prediction of key ADME property trends and biological activity in an industrial drug discovery environment, represented by microsomal lability, CYP3A4 inhibition and factor Xa inhibition. Three architectures are exemplified, our earlier described multilayer perceptron approach (MLP), graph convolutional network-based models (GCN) and a vector representation approach, Mol2Vec. From a statistical perspective, MLP and GCN were found to perform superior over Mol2Vec, when applied to external validation sets. Interestingly, GCN-based predictions are most stable over a longer period in a time series validation study. Apart from those statistical observations, DNN prove of value to guide local SAR. To illustrate this important aspect in pharmaceutical research projects, we discuss challenging applications in medicinal chemistry towards a more realistic picture of artificial intelligence in drug discovery.

Photosensitive Ru(II) Complexes as Inhibitors of the Major Human Drug Metabolizing Enzyme CYP3A4.

We report the synthesis and photochemical and biological characterization of the first selective and potent metal-based inhibitors of cytochrome P450 3A4 (CYP3A4), the major human drug metabolizing enzyme. Five Ru(II)-based derivatives were prepared from two analogs of the CYP3A4 inhibitor ritonavir, 4 and 6: [Ru(tpy)(L)(6)]Cl2 (tpy = 2,2':6',2″-terpyridine) with L = 6,6'-dimethyl-2,2'-bipyridine (Me2bpy; 8), dimethylbenzo[i]dipyrido[3,2-a:2',3'-c]phenazine (Me2dppn; 10) and 3,6-dimethyl-10,15-diphenylbenzo[i]dipyrido[3,2-a:2',3'-c]phenazine (Me2Ph2dppn; 11), [Ru(tpy)(Me2bpy)(4)]Cl2 (7) and [Ru(tpy)(Me2dppn)(4)]Cl2 (9). Photochemical release of 4 or 6 from 7-11 was demonstrated, and the spectrophotometric evaluation of 7 showed that it behaves similarly to free 4 (type II heme ligation) after irradiation with visible light but not in the dark. Unexpectedly, the intact Ru(II) complexes 7 and 8 were found to inhibit CYP3A4 potently and specifically through direct binding to the active site without heme ligation. Caged inhibitors 9-11 showed dual action properties by combining photoactivated dissociation of 4 or 6 with efficient 1O2 production. In prostate adenocarcinoma DU-145 cells, compound 9 had the best synergistic effect with vinblastine, the anticancer drug primarily metabolized by CYP3A4 in vivo. Thus, our study establishes a new paradigm in CYP inhibition using metalated complexes and suggests possible utilization of photoactive CYP3A4 inhibitory compounds in clinical applications, such as enhancement of therapeutic efficacy of anticancer drugs.

Rational Design of CYP3A4 Inhibitors: A One-Atom Linker Elongation in Ritonavir-Like Compounds Leads to a Marked Improvement in the Binding Strength.

Inhibition of the major human drug-metabolizing cytochrome P450 3A4 (CYP3A4) by pharmaceuticals and other xenobiotics could lead to toxicity, drug-drug interactions and other adverse effects, as well as pharmacoenhancement. Despite serious clinical implications, the structural basis and attributes required for the potent inhibition of CYP3A4 remain to be established. We utilized a rational inhibitor design to investigate the structure-activity relationships in the analogues of ritonavir, the most potent CYP3A4 inhibitor in clinical use. This study elucidated the optimal length of the head-group spacer using eleven (series V) analogues with the R1/R2 side-groups as phenyls or R1-phenyl/R2-indole/naphthalene in various stereo configurations. Spectral, functional and structural characterization of the inhibitory complexes showed that a one-atom head-group linker elongation, from pyridyl-ethyl to pyridyl-propyl, was beneficial and markedly improved Ks, IC50 and thermostability of CYP3A4. In contrast, a two-atom linker extension led to a multi-fold decrease in the binding and inhibitory strength, possibly due to spatial and/or conformational constraints. The lead compound, 3h, was among the best inhibitors designed so far and overall, the strongest binder (Ks and IC50 of 0.007 and 0.090 µM, respectively). 3h was the fourth structurally simpler inhibitor superior to ritonavir, which further demonstrates the power of our approach.

The relationship between stereochemical and both, pharmacological and ADME-Tox, properties of the potent hydantoin 5-HT7R antagonist MF-8.

This study concerns synthesis and evaluation of pharmacodynamic and pharmacokinetic profile for all four stereoisomers of MF-8 (5-(4-fluorophenyl)-3-(2-hydroxy-3-(4-(2-methoxyphenyl)piperazin-1-yl)propyl)-5-methylimidazolidine-2,4-dione), the previously described, highly potent 5-HT7R ligand with antidepressant activity on mice. The combination of DFT calculations of 1H NMR chemical shifts with docking and dynamic simulations, in comparison to experimental screening results, provided prediction of the configuration for one of two present stereogenic centers. The experimental data for stereoisomers (MF-8A-MF-8D) confirmed the significant impact of stereochemistry on both, 5-HT7R affinity and antagonistic action, with Ki and Kb values in the range of 3-366 nM and 0.024-99 µM, respectively. We also indicated the stereochemistry-dependent influence of the tested compounds on P-glycoprotein efflux, absorption in Caco-2 model, metabolic pathway as well as CYP3A4 and CYP2C9 activities.

Nitrogen-Walk Approach to Explore Bioisosteric Replacements in a Series of Potent A2B Adenosine Receptor Antagonists.

A systematic exploration of bioisosteric replacements for furan and thiophene cores in a series of potent A2BAR antagonists has been carried out using the nitrogen-walk approach. A collection of 42 novel alkyl 4-substituted-2-methyl-1,4-dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-3-carboxylates, which contain 18 different pentagonal heterocyclic frameworks at position 4, was synthesized and evaluated. This study enabled the identification of new ligands that combine remarkable affinity (Ki < 30 nM) and exquisite selectivity. The structure-activity relationship (SAR) trends identified were substantiated by a molecular modeling study, based on a receptor-driven docking model and including a systematic free energy perturbation (FEP) study. Preliminary evaluation of the CYP3A4 and CYP2D6 inhibitory activity in optimized ligands evidenced weak and negligible activity, respectively. The stereospecific interaction between hA2BAR and the eutomer of the most attractive novel antagonist (S)-18g (Ki = 3.66 nM) was validated.

An increase in side-group hydrophobicity largely improves the potency of ritonavir-like inhibitors of CYP3A4.

Identification of structural determinants required for potent inhibition of drug-metabolizing cytochrome P450 3A4 (CYP3A4) could help develop safer drugs and more effective pharmacoenhancers. We utilize a rational inhibitor design to decipher structure-activity relationships in analogues of ritonavir, a highly potent CYP3A4 inhibitor marketed as pharmacoenhancer. Analysis of compounds with the R1 side-group as phenyl or naphthalene and R2 as indole or naphthalene in different stereo configuration showed that (i) analogues with the R2-naphthalene tend to bind tighter and inhibit CYP3A4 more potently than the R2-phenyl/indole containing counterparts; (ii) stereochemistry becomes a more important contributing factor, as the bulky side-groups limit the ability to optimize protein-ligand interactions; (iii) the relationship between the R1/R2 configuration and preferential binding to CYP3A4 is complex and depends on the side-group functionality/interplay and backbone spacing; and (iv) three inhibitors, 5a-b and 7d, were superior to ritonavir (IC50 of 0.055-0.085 µM vs. 0.130 µM, respectively).

Synthesis of Spin-Labelled Bergamottin: A Potent CYP3A4 Inhibitor with Antiproliferative Activity.

Bergamottin (BM, 1), a component of grapefruit juice, acts as an inhibitor of some isoforms of the cytochrome P450 (CYP) enzyme, particularly CYP3A4. Herein, a new bergamottin containing a nitroxide moiety (SL-bergamottin, SL-BM, 10) was synthesized; chemically characterized, evaluated as a potential inhibitor of the CYP2C19, CYP3A4, and CYP2C9 enzymes; and compared to BM and known inhibitors such as ketoconazole (KET) (3A4), warfarin (WAR) (2C9), and ticlopidine (TIC) (2C19). The antitumor activity of the new SL-bergamottin was also investigated. Among the compounds studied, BM showed the strongest inhibition of the CYP2C9 and 2C19 enzymes. SL-BM is a more potent inhibitor of CYP3A4 than the parent compound; this finding was also supported by docking studies, suggesting that the binding positions of BM and SL-BM to the active site of CYP3A4 are very similar, but that SL-BM had a better ∆Gbind value than that of BM. The nitroxide moiety markedly increased the antitumor activity of BM toward HeLa cells and marginally increased its toxicity toward a normal cell line. In conclusion, modification of the geranyl sidechain of BM can result in new CYP3A4 enzyme inhibitors with strong antitumor effects.

Tofacitinib Is a Mechanism-Based Inactivator of Cytochrome P450 3A4.

Tofacitinib (TFT) is an oral JAK inhibitor which has been approved for the treatment of moderately and severely active rheumatoid arthritis. TFT was found to show concentration-, time-, and NADPH-dependent inhibition of CYP3A4, and irreversibility of the inactivation was also observed. Incubation (40 min, 37 °C) of recombinant CYP3A4 with TFT at 200 µM resulted in >70% loss of CYP3A4 activity. Estimated kinact and KI were 0.037 min-1 and 93.2 µM, respectively. GSH and superoxide dismutase/catalase revealed minor or little protection against the CYP3A4 inactivation. Furthermore, ketoconazole attenuated TFT-mediated CYP3A4 inactivation. Epoxide and α-keto-aldehyde intermediates of TFT were trapped and characterized in microsomal incubations, respectively. The aldehyde intermediate is believed to be the key for the enzyme inactivation. Multiple P450 enzymes, including CYPs2C19, 3A4, 2D6, and 1A2, participated in the metabolism of TFT to the epoxide, while the formation of the aldehyde was mainly catalyzed by CYP3A4. In conclusion, TFT was proven to be a mechanism-based inactivator of CYP3A4.

Discovery of 2,4-diamino-5-cyanopyrimidine derivatives as protein kinase C theta inhibitors with mitigated time-dependent drug-drug interactions.

Protein kinase C theta (PKCθ) plays a critical role in T cell signaling and has therapeutic potential for T cell-mediated diseases such as transplant rejection and rheumatoid arthritis. PKCθ inhibitors have emerged as effective immunomodulative agents for the prevention of transplant rejection. We previously reported that the 2,4-diamino-5-cyanopyrimidine derivative 2 was a potent PKCθ inhibitor; however, it exhibited CYP3A4 time-dependent inhibition (TDI). Here, we report the structural modification of compound 2 into 34 focusing on mitigating CYP3A4 TDI. Compound 34 exhibited potent in vitro activity with mitigated CYP3A4 TDI and efficacy in vivo transplant model.

Design, Synthesis and Biological Evaluation of Spiro Cyclohexane-1,2- Quinazoline Derivatives as Potent Dipeptidyl Peptidase IV Inhibitors.

OBJECTIVE: Inhibition of dipeptidyl peptidase IV (DPP-4) is currently one of the most valuable and potential chemotherapeutic regimes for the medication of Type 2 Diabetes Mellitus (T2DM). METHOD: Based on linagliptin, this study discusses the design, synthesis and biological evaluation of spiro cyclohexane-1,2'-quinazoline scaffold hybridized with various heterocyclic ring systems through different atomic spacers as a highly potent DPP-4 inhibitors. DPP-4 enzyme assay represented that most of the target compounds are 102-103 folds more active than the reference drug linagliptin (IC50: 0.0005-0.0089 nM vs 0.77 nM; respectively). Moreover, in vivo oral hypoglycemic activity assay revealed that most of the tested candidates were more potent than the reference drug, sitagliptin, producing rapid onset with long duration of activity that extends to 24 h. Interestingly, the derivatives 11, 16, 18a and 23 showed evidence of mild cytochrome P450 3A4 (CYP3A4) inhibition (IC50; > 210 µM) and their acute toxicity (LD50) was more than 1.9 gm/kg. Molecular simulation study of the new quinazoline derivatives explained the obtained biological results. CONCLUSION: Finally, we conclude that our target compounds could be highly beneficial for diabetic patients in the clinic.

Amelioration of mechanism-based inactivation of CYP3A4 by a H-PGDS inhibitor.

This work describes the rational amelioration of mechanism-based inactivation (MBI) of Cytochrome P450 (CYP) 3A4 in a human hematopoietic prostaglandin D synthase (hH-PGDS) inhibitor (cpd 1). We utilized metabolism reports in order to check if patterns in the metabolism of 1 and similar compounds by CYP3A4 could be deciphered. Then we used structure based design, first modifying the CYP3A4 crystal structure (pdb code: 4NY4) by adding an oxyferryl moiety to the heme, followed by validating the modified structure to obtain the 1' and 4 position oxidation products of midazolam and then recapitulating the metabolism patterns deciphered previously for 1 and analogs. We checked if the pattern deciphered could lead to a putative reactive moiety. Finally we used the docking pose of 1 into this model of the modified CYP3A4 crystal structure to guide transformation of 1 into MBI-free H-PGDS inhibitors.

Novel hydrazide-hydrazone and amide substituted coumarin derivatives: Synthesis, cytotoxicity screening, microarray, radiolabeling and in vivo pharmacokinetic studies.

The current work presents the synthesis and biological evaluation of new series of coumarin hydrazide-hydrazone derivatives that showed in vitro broad spectrum antitumor activities against resistant pancreatic carcinoma (Panc-1), hepatocellular carcinoma (HepG2) and leukemia (CCRF) cell lines using doxorubicin as reference standard. Bromocoumarin hydrazide-hydrazone derivative (BCHHD) 11b showed excellent anticancer activity against all tested cancer cell lines. Enzyme assays showed that BCHHD 11b induced apoptosis due to activation of caspases 3/7. Moreover, 11b inhibited GST and CYP3A4 in a dose dependent manner and the induced cell death could be attributed to metabolic inhibition. Moreover, 11b microarray analysis showed significant up- and down-regulation of many genes in the treated cells related to apoptosis, cell cycle, tumor growth and suppressor genes. All of the above presents BCHHD 11b as a potent anticancer agent able to overcome drug resistance. In addition, compound 11b was able to serve as a chemical carrier for 99mTc and the in vivo biodistribution study of 99mTc-11b complex revealed a remarkable targeting ability of 99mTc into solid tumor showing that 99mTc-11b might be used as a promising radiopharmaceutical imaging agent for cancer.

2-Aminopyridine-Based Mitogen-Activated Protein Kinase Kinase Kinase Kinase 4 (MAP4K4) Inhibitors: Assessment of Mechanism-Based Safety.

Studies have linked the serine-threonine kinase MAP4K4 to the regulation of a number of biological processes and/or diseases, including diabetes, cancer, inflammation, and angiogenesis. With a majority of the members of our lead series (e.g., 1) suffering from time-dependent inhibition (TDI) of CYP3A4, we sought design avenues that would eliminate this risk. One such approach arose from the observation that carboxylic acid-based intermediates employed in our discovery efforts retained high MAP4K4 inhibitory potency and were devoid of the TDI risk. The medicinal chemistry effort that led to the discovery of this central nervous system-impaired inhibitor together with its preclinical safety profile is described.

Development of (6 R)-2-Nitro-6-[4-(trifluoromethoxy)phenoxy]-6,7-dihydro-5 H-imidazo[2,1- b][1,3]oxazine (DNDI-8219): A New Lead for Visceral Leishmaniasis.

Discovery of the potent antileishmanial effects of antitubercular 6-nitro-2,3-dihydroimidazo[2,1- b][1,3]oxazoles and 7-substituted 2-nitro-5,6-dihydroimidazo[2,1- b][1,3]oxazines stimulated the examination of further scaffolds (e.g., 2-nitro-5,6,7,8-tetrahydroimidazo[2,1- b][1,3]oxazepines), but the results for these seemed less attractive. Following the screening of a 900-compound pretomanid analogue library, several hits with more suitable potency, solubility, and microsomal stability were identified, and the superior efficacy of newly synthesized 6 R enantiomers with phenylpyridine-based side chains was established through head-to-head assessments in a Leishmania donovani mouse model. Two such leads ( R-84 and R-89) displayed promising activity in the more stringent Leishmania infantum hamster model but were unexpectedly found to be potent inhibitors of hERG. An extensive structure-activity relationship investigation pinpointed two compounds ( R-6 and pyridine R-136) with better solubility and pharmacokinetic properties that also provided excellent oral efficacy in the same hamster model (>97% parasite clearance at 25 mg/kg, twice daily) and exhibited minimal hERG inhibition. Additional profiling earmarked R-6 as the favored backup development candidate.

Synthesis of Furanocoumarin, Benzofuran and Coumarin Derivatives Possessing an Inhibitory Effect on Human CYP, and Elucidation of the Inhibitory Mechanism.

Grapefruit juice (GFJ) consumption has been shown to increase the bioavailability of certain orally administered drugs. The furanocoumarin derivatives Paradisin A and bergamottin, which are present in GFJ, are potent mechanism-based inhibitors of CYP3A4. The primary aim of this work was to synthesize a series of furanocoumarin derivatives with a view to determining the relationship between the structure of the inhibitors and their inhibitory CYP3A4 activity. Furanocoumarin derivatives that were more stable and accessible than the furanocoumarin derivatives in GFJ were prepared, and their ability to inhibit CYP3A4 was examined. Synthesized furanocoumarin monomers showed strong mechanism-based inhibition of CYP3A4. The furanocoumarin dimers are also mechanism-based inhibitors of CYP3A4. These monomers and dimers are more potent inhibitors of CYP3A4 than bergamottin and Paradisin A, respectively.

Design, Synthesis, and Evaluation of the Highly Selective and Potent G-Protein-Coupled Receptor Kinase 2 (GRK2) Inhibitor for the Potential Treatment of Heart Failure.

A novel class of therapeutic drug candidates for heart failure, highly potent and selective GRK2 inhibitors, exhibit potentiation of ß-adrenergic signaling in vitro studies. Hydrazone derivative 5 and 1,2,4-triazole derivative 24a were identified as hit compounds by HTS. New scaffold generation and SAR studies of all parts resulted in a 4-methyl-1,2,4-triazole derivative with an N-benzylcarboxamide moiety with highly potent activity toward GRK2 and selectivity over other kinases. In terms of subtype selectivity, these compounds showed enough selectivity against GRK1, 5, 6, and 7 with almost equipotent inhibition to GRK3. Our medicinal chemistry efforts led to the discovery of 115h (GRK2 IC50 = 18 nM), which was obtained the cocrystal structure with human GRK2 and an inhibitor of GRK2 that potentiates ß-adrenergic receptor (ßAR)-mediated cAMP accumulation and prevents internalization of ßARs in ß2AR-expressing HEK293 cells treated with isoproterenol. Therefore, 115h appears to be a novel class of therapeutic for heart failure treatment.

Discovery of Clinical Candidate 4-[2-(5-Amino-1H-pyrazol-4-yl)-4-chlorophenoxy]-5-chloro-2-fluoro-N-1,3-thiazol-4-ylbenzenesulfonamide (PF-05089771): Design and Optimization of Diaryl Ether Aryl Sulfonamides as Selective Inhibitors of NaV1.7.

A series of acidic diaryl ether heterocyclic sulfonamides that are potent and subtype selective NaV1.7 inhibitors is described. Optimization of early lead matter focused on removal of structural alerts, improving metabolic stability and reducing cytochrome P450 inhibition driven drug-drug interaction concerns to deliver the desired balance of preclinical in vitro properties. Concerns over nonmetabolic routes of clearance, variable clearance in preclinical species, and subsequent low confidence human pharmacokinetic predictions led to the decision to conduct a human microdose study to determine clinical pharmacokinetics. The design strategies and results from preclinical PK and clinical human microdose PK data are described leading to the discovery of the first subtype selective NaV1.7 inhibitor clinical candidate PF-05089771 (34) which binds to a site in the voltage sensing domain.

Scaffold Hopping and Optimization of Maleimide Based Porcupine Inhibitors.

Porcupine is an O-acyltransferase that regulates Wnt secretion. Inhibiting porcupine may block the Wnt pathway which is often dysregulated in various cancers. Consequently porcupine inhibitors are thought to be promising oncology therapeutics. A high throughput screen against porcupine revealed several potent hits that were confirmed to be Wnt pathway inhibitors in secondary assays. We developed a pharmacophore model and used the putative bioactive conformation of a xanthine inhibitor for scaffold hopping. The resulting maleimide scaffold was optimized to subnanomolar potency while retaining good physical druglike properties. A preclinical development candidate was selected for which extensive in vitro and in vivo profiling is reported.

Discovery of Highly Potent Dual Orexin Receptor Antagonists via a Scaffold-Hopping Approach.

Starting from suvorexant (trade name Belsomra), we successfully identified interesting templates leading to potent dual orexin receptor antagonists (DORAs) via a scaffold-hopping approach. Structure-activity relationship optimization allowed us not only to improve the antagonistic potency on both orexin 1 and orexin 2 receptors (Ox1 and Ox2, respectively), but also to increase metabolic stability in human liver microsomes (HLM), decrease time-dependent inhibition of cytochrome P450 (CYP) 3A4, and decrease P-glycoprotein (Pgp)-mediated efflux. Compound 80 c [{(1S,6R)-3-(6,7-difluoroquinoxalin-2-yl)-3,8-diazabicyclo[4.2.0]octan-8-yl}(4-methyl-[1,1'-biphenyl]-2-yl)methanone] is a potent and selective DORA that inhibits the stimulating effects of orexin peptides OXA and OXB at both Ox1 and Ox2. In calcium-release assays, 80 c was found to exhibit an insurmountable antagonistic profile at both Ox1 and Ox2, while displaying a sleep-promoting effect in rat and dog models, similar to that of the benchmark compound suvorexant.

Quinoxaline-based inhibitors of Ebola and Marburg VP40 egress.

We prepared a series of quinoxalin-2-mercapto-acetyl-urea analogs and evaluated them for their ability to inhibit viral egress in our Marburg and Ebola VP40 VLP budding assays in HEK293T cells. We also evaluated selected compounds in our bimolecular complementation assay (BiMC) to detect and visualize a Marburg mVP40-Nedd4 interaction in live mammalian cells. Antiviral activity was assessed for selected compounds using a live recombinant vesicular stomatitis virus (VSV) (M40 virus) that expresses the EBOV VP40 PPxY L-domain. Finally selected compounds were evaluated in several ADME assays to have an early assessment of their drug properties. Our compounds had low nM potency in these assays (e.g., compounds 21, 24, 26, 39), and had good human liver microsome stability, as well as little or no inhibition of P450 3A4.