Discovery of BI-9508, a Brain-Penetrant GPR88-Receptor-Agonist Tool Compound for In Vivo Mouse Studies.

Decreased activity and expression of the G-protein coupled receptor GPR88 is linked to many behavior-linked neurological disorders. Published preclinical GPR88 allosteric agonists all have in vivo pharmacokinetic properties that preclude their progression to the clinic, including high lipophilicity and poor brain penetration. Here, we describe our attempts to improve GPR88 agonists' drug-like properties and our analysis of the trade-offs required to successfully target GPR88's allosteric pocket. We discovered two new GPR88 agonists: One that reduced morphine-induced locomotor activity in a murine proof-of-concept study, and the atropoisomeric BI-9508, which is a brain penetrant and has improved pharmacokinetic properties and dosing that recommend it for future in vivo studies in rodents. BI-9508 still suffers from high lipophilicity, and research on this series was halted. Because of its utility as a tool compound, we now offer researchers access to BI-9508 and a negative control free of charge via Boehringer Ingelheim's open innovation portal opnMe.com.

Action of Dipeptidyl Peptidase-4 Inhibitors on SARS-CoV-2 Main Protease.

In a recent publication, Eleftheriou et al. proposed that inhibitors of dipeptidyl peptidase-4 (DPP-4) are functional inhibitors of the main protease (Mpro ) of SARS-CoV-2. Their predictions prompted the authors to suggest linagliptin, a DPP-4 inhibitor and approved anti-diabetes drug, as a repurposed drug candidate against the ongoing COVID-19 pandemic. We used an enzymatic assay measuring the inhibition of Mpro catalytic activity in the presence of four different commercially available gliptins (linagliptin, sitagliptin, alogliptin and saxagliptin) and several structural analogues of linagliptin to study the binding of DPP-4 inhibitors to Mpro and their functional activity. We show here that DPP-4 inhibitors like linagliptin, other gliptins and structural analogues are inactive against Mpro .

Multiple Binding Sites Contribute to the Mechanism of Mixed Agonistic and Positive Allosteric Modulators of the Cannabinoid CB1 Receptor.

The cannabinoid CB1 receptor (CB1R) is an abundant metabotropic G-protein-coupled receptor that has been difficult to address therapeutically because of CNS side effects exerted by orthosteric drug candidates. Recent efforts have focused on developing allosteric modulators that target CB1R. Compounds from the recently discovered class of mixed agonistic and positive allosteric modulators (Ago-PAMs) based on 2-phenylindoles have shown promising functional and binding properties as CB1R ligands. Here, we identify binding modes of both the CP 55,940 agonist and GAT228, a 2-phenylindole allosteric modulator, by using our metadynamics simulation protocol, and quantify their affinity and cooperativity by atomistic simulations. We demonstrate the involvement of multiple adjunct binding sites in the Ago-PAM characteristics of the 2-phenylindole modulators and explain their ability to compete with orthosteric agonists at higher concentrations. We validate these results experimentally by showing the contribution of multiple sites on the allosteric binding of ZCZ011, another homologous member of the class, together with the orthosteric agonist.

Enantiomeric atropisomers inhibit HCV polymerase and/or HIV matrix: characterizing hindered bond rotations and target selectivity.

An anthranilic acid series of allosteric thumb pocket 2 HCV NS5B polymerase inhibitors exhibited hindered rotation along a covalent bond axis, and the existence of atropisomer chirality was confirmed by NMR, HPLC analysis on chiral supports, and computational studies. A thorough understanding of the concerted rotational properties and the influence exerted by substituents involved in this steric phenomenon was attained through biophysical studies on a series of truncated analogues. The racemization half-life of a compound within this series was determined to be 69 min, which was consistent with a class 2 atropisomer (intermediate conformational exchange). It was further found by X-ray crystallography that one enantiomer of a compound bound to the intended HCV NS5B polymerase target whereas the mirror image atropisomer was able to bind to an unrelated HIV matrix target. Analogues were then identified that selectively inhibited the former. These studies highlight that atropisomer chirality can lead to distinct entities with specific properties, and the phenomenon of atropisomerism in drug discovery should be evaluated and appropriately managed.

Molecular dynamics simulations and structure-based rational design lead to allosteric HCV NS5B polymerase thumb pocket 2 inhibitor with picomolar cellular replicon potency.

The design and preliminary SAR of a new series of 1H-quinazolin-4-one (QAZ) allosteric HCV NS5B thumb pocket 2 (TP-2) inhibitors was recently reported. To support optimization efforts, a molecular dynamics (MD) based modeling workflow was implemented, providing information on QAZ binding interactions with NS5B. This approach predicted a small but critical ligand-binding induced movement of a protein backbone region which increases the pocket size and improves access to the backbone carbonyl groups of Val 494 and Pro 495. This localized backbone shift was consistent with key SAR results and was subsequently confirmed by X-ray crystallography. The MD protocol guided the design of inhibitors, exploiting novel H-bond interactions with the two backbone carbonyl groups, leading to the first thumb pocket 2 NS5B inhibitor with picomolar antiviral potency in genotype (gt) 1a and 1b replicons (EC50 = 120 and 110 pM, respectively) and with EC50 ≤ 80 nM against gt 2-6.

Structure-based design of novel HCV NS5B thumb pocket 2 allosteric inhibitors with submicromolar gt1 replicon potency: discovery of a quinazolinone chemotype.

We describe the structure-based design of a novel lead chemotype that binds to thumb pocket 2 of HCV NS5B polymerase and inhibits cell-based gt1 subgenomic reporter replicons at sub-micromolar concentrations (EC50<200nM). This new class of potent thumb pocket 2 inhibitors features a 1H-quinazolin-4-one scaffold derived from hybridization of a previously reported, low affinity thiazolone chemotype with our recently described anthranilic acid series. Guided by X-ray structural information, a key NS5B-ligand interaction involving the carboxylate group of anthranilic acid based inhibitors was replaced by a neutral two-point hydrogen bonding interaction between the quinazolinone scaffold and the protein backbone. The in vitro ADME and in vivo rat PK profile of representative analogs are also presented and provide areas for future optimization of this new class of HCV polymerase inhibitors.

Design, synthesis and biological evaluation of novel aminothiazoles as antiviral compounds acting against human rhinovirus.

We describe here the design, synthesis and biological evaluation of antiviral compounds acting against human rhinovirus (HRV). A series of aminothiazoles demonstrated pan-activity against the HRV genotypes screened and productive structure-activity relationships. A comprehensive investigational library was designed and performed allowing the identification of potent compounds with lower molecular weight and improved ADME profile. 31d-1, 31d-2, 31f showed good exposures in CD-1 mice. The mechanism of action was discovered to be a host target: the lipid kinase phosphatidylinositol 4-kinase III beta (PI4KIIIß). The identification of the pan-HRV active compound 31f combined with a structurally distinct literature compound T-00127-HEV1 allowed the assessment of target related tolerability of inhibiting this kinase for a short period of time in order to prevent HRV replication.

Phosphatidylinositol 4-kinase III beta is essential for replication of human rhinovirus and its inhibition causes a lethal phenotype in vivo.

Human rhinovirus (HRV) is the predominant cause of the common cold, but more importantly, infection may have serious repercussions in asthmatics and chronic obstructive pulmonary disorder (COPD) patients. A cell-based antiviral screen against HRV was performed with a subset of our proprietary compound collection, and an aminothiazole series with pan-HRV species and enteroviral activity was identified. The series was found to act at the level of replication in the HRV infectious cycle. In vitro selection and sequencing of aminothiazole series-resistant HRV variants revealed a single-nucleotide mutation leading to the amino acid change I42V in the essential HRV 3A protein. This same mutation has been previously implicated in resistance to enviroxime, a former clinical-stage antipicornavirus agent. Enviroxime-like compounds have recently been shown to target the lipid kinase phosphatidylinositol 4-kinase III beta (PI4KIIIß). A good correlation between PI4KIIIß activity and HRV antiviral potency was found when analyzing the data over 80 compounds of the aminothiazole series, covering a 750-fold potency range. The mechanism of action through PI4KIIIß inhibition was further demonstrated by small interfering RNA (siRNA) knockdown of PI4KB, which reduced HRV replication and also increased the potency of the PI4KIIIß inhibitors. Inhibitors from two different structural classes with promising pharmacokinetic profiles and with very good selectivity for PI4KIIIß were used to dissociate compound-related toxicity from target-related toxicity. Mortality was seen in all dosing groups of mice treated with either compound, therefore suggesting that short-term inhibition of PI4KIIIß is deleterious.

Optimization of a 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of HIV capsid assembly inhibitors 2: structure-activity relationships (SAR) of the C3-phenyl moiety.

Detailed structure-activity relationships of the C3-phenyl moiety that allow for the optimization of antiviral potency of a series of 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione inhibitors of HIV capsid (CA) assembly are described. Combination of favorable substitutions gave additive SAR and allowed for the identification of the most potent compound in the series, analog 27. Productive SAR also transferred to the benzotriazepine and spirobenzodiazepine scaffolds, providing a solution to the labile stereocenter at the C3 position. The molecular basis of how compound 27 inhibits mature CA assembly is rationalized using high-resolution structural information. Our understanding of how compound 27 may inhibit immature Gag assembly is also discussed.

Optimization of a 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of HIV capsid assembly inhibitors 1: addressing configurational instability through scaffold modification.

The optimization of a 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of inhibitors of HIV-1 capsid assembly that possess a labile stereocenter at C3 is described. Quaternization of the C3 position of compound 1 in order to prevent racemization gave compound 2, which was inactive in our capsid disassembly assay. A likely explanation for this finding was revealed by in silico analysis predicting a dramatic increase in energy of the bioactive conformation upon quaternization of the C3 position. Replacement of the C3 of the diazepine ring with a nitrogen atom to give the 1,5-dihydro-benzo[f][1,3,5]triazepine-2,4-dione analog 4 was well tolerated. Introduction of a rigid spirocyclic system at the C3 position gave configurationally stable 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione analog 5, which was able to access the bioactive conformation without a severe energetic penalty and inhibit capsid assembly. Preliminary structure-activity relationships (SAR) and X-ray crystallographic data show that knowledge from the 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of inhibitors of HIV-1 capsid assembly can be transferred to these new scaffolds.

Discovery and structural characterization of a new inhibitor series of HIV-1 nucleocapsid function: NMR solution structure determination of a ternary complex involving a 2:1 inhibitor/NC stoichiometry.

The nucleocapsid (NC) protein is an essential factor with multiple functions within the human immunodeficiency virus type 1 (HIV-1) replication cycle. In this study, we describe the discovery of a novel series of inhibitors that targets HIV-1 NC protein by blocking its interaction with nucleic acids. This series was identified using a previously described capsid (CA) assembly assay, employing a recombinant HIV-1 CA-NC protein and immobilized TG-rich deoxyoligonucleotides. Using visible absorption spectroscopy, we were able to demonstrate that this new inhibitor series binds specifically and reversibly to the NC with a peculiar 2:1 stoichiometry. A fluorescence-polarization-based binding assay was also developed in order to monitor the inhibitory activities of this series of inhibitors. To better characterize the structural aspect of inhibitor binding onto NC, we performed NMR studies using unlabeled and (13)C,(15)N-double-labeled NC(1-55) protein constructs. This allowed the determination of the solution structure of a ternary complex characterized by two inhibitor molecules binding to the two zinc knuckles of the NC protein. To the best of our knowledge, this represents the first report of a high-resolution structure of a small-molecule inhibitor bound to NC, demonstrating sub-micromolar potency and moderate antiviral potency with one analogue of the series. This structure was compared with available NC/oligonucleotide complex structures and further underlined the high flexibility of the NC protein, allowing it to adopt many conformations in order to bind its different oligonucleotide/nucleomimetic targets. In addition, analysis of the interaction details between the inhibitor molecules and NC demonstrated how this novel inhibitor series is mimicking the guanosine nucleobases found in many reported complex structures.

Human intestinal transporter database: QSAR modeling and virtual profiling of drug uptake, efflux and interactions.

PURPOSE: Membrane transporters mediate many biological effects of chemicals and play a major role in pharmacokinetics and drug resistance. The selection of viable drug candidates among biologically active compounds requires the assessment of their transporter interaction profiles. METHODS: Using public sources, we have assembled and curated the largest, to our knowledge, human intestinal transporter database (>5,000 interaction entries for >3,700 molecules). This data was used to develop thoroughly validated classification Quantitative Structure-Activity Relationship (QSAR) models of transport and/or inhibition of several major transporters including MDR1, BCRP, MRP1-4, PEPT1, ASBT, OATP2B1, OCT1, and MCT1. RESULTS: QSAR models have been developed with advanced machine learning techniques such as Support Vector Machines, Random Forest, and k Nearest Neighbors using Dragon and MOE chemical descriptors. These models afforded high external prediction accuracies of 71-100% estimated by 5-fold external validation, and showed hit retrieval rates with up to 20-fold enrichment in the virtual screening of DrugBank compounds. CONCLUSIONS: The compendium of predictive QSAR models developed in this study can be used for virtual profiling of drug candidates and/or environmental agents with the optimal transporter profiles.

Inhibition of HIV-1 capsid assembly: optimization of the antiviral potency by site selective modifications at N1, C2 and C16 of a 5-(5-furan-2-yl-pyrazol-1-yl)-1H-benzimidazole scaffold.

A uHTS campaign led to the discovery of a 5-(5-furan-2-ylpyrazol-1-yl)-1H-benzimidazole series that inhibits assembly of HIV-1 capsid. Synthetic manipulations at N1, C2 and C16 positions improved the antiviral potency by a . The X-ray structure of 33 complexed with the capsid N-terminal domain allowed identification of major interactions between the inhibitor and the protein.

Revealing atropisomer axial chirality in drug discovery.

An often overlooked source of chirality is atropisomerism, which results from slow rotation along a bond axis due to steric hindrance and/or electronic factors. If undetected or not managed properly, this time-dependent chirality has the potential to lead to serious consequences, because atropisomers can be present as distinct enantiomers or diastereoisomers with their attendant different properties. Herein we introduce a strategy to reveal and classify compounds that have atropisomeric chirality. Energy barriers to axial rotation were calculated using quantum mechanics, from which predicted high barriers could be experimentally validated. A calculated rotational energy barrier of 20 kcal mol(-1) was established as a suitable threshold to distinguish between atropisomers and non-atropisomers with a prediction accuracy of 86%. This methodology was applied to subsets of drug databases in the course of which atropisomeric drugs were identified. In addition, some drugs were exposed that were not yet known to have this chiral attribute. The most valuable utility of this tool will be to predict atropisomerism along the drug discovery pathway. When used in concert with our compound classification scheme, decisions can be made during early discovery stages such as "hit-to-lead" and "lead optimization," to foresee and validate the presence of atropisomers and to exercise options of removing, further stabilizing, or rendering the chiral axis of interest more freely rotatable via SAR design, thereby decreasing this potential liability within a compound series. The strategy can also improve drug development plans, such as determining whether a drug or series should be developed as a racemic mixture or as an isolated single compound. Moreover, the work described herein can be extended to other chemical fields that require the assessment of potential chiral axes.

Combined X-ray, NMR, and kinetic analyses reveal uncommon binding characteristics of the hepatitis C virus NS3-NS4A protease inhibitor BI 201335.

Hepatitis C virus infection, a major cause of liver disease worldwide, is curable, but currently approved therapies have suboptimal efficacy. Supplementing these therapies with direct-acting antiviral agents has the potential to considerably improve treatment prospects for hepatitis C virus-infected patients. The critical role played by the viral NS3 protease makes it an attractive target, and despite its shallow, solvent-exposed active site, several potent NS3 protease inhibitors are currently in the clinic. BI 201335, which is progressing through Phase IIb trials, contains a unique C-terminal carboxylic acid that binds noncovalently to the active site and a bromo-quinoline substitution on its proline residue that provides significant potency. In this work we have used stopped flow kinetics, x-ray crystallography, and NMR to characterize these distinctive features. Key findings include: slow association and dissociation rates within a single-step binding mechanism; the critical involvement of water molecules in acid binding; and protein side chain rearrangements, a bromine-oxygen halogen bond, and profound pK(a) changes within the catalytic triad associated with binding of the bromo-quinoline moiety.

Discovery of a 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of inhibitors of HIV-1 capsid assembly.

The discovery of a 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione series of inhibitors of HIV-1 capsid assembly is described. Synthesis of analogs of the 1,5-dihydrobenzo[b][1,4]diazepine-2,4-dione hit established structure-activity relationships. Replacement of the enamine functionality of the hit series with either an imidazole or a pyrazole ring led to compounds that inhibited both capsid assembly and reverse transcriptase. Optimization of the bicyclic benzodiazepine scaffold to include a 3-phenyl substituent led to lead compound 48, a pure capsid assembly inhibitor with improved antiviral activity.

Development of specific "drug-like property" rules for carboxylate-containing oral drug candidates.

The carboxylate moiety is an important pharmacophore in the medicinal chemist's arsenal and is sometimes an irreplaceable functionality in drug-target interactions. Thus, practical guidance on its use in the most optimized manner would be a welcome addition to rational drug design. Key physicochemical and ADMET-PK properties from a dataset of drugs containing a carboxylate (COOH) moiety were assembled and compared with those of a broader, general drug dataset. Our main objective was to identify features specific to COOH-containing oral drugs that could be converted into simple rules delineating the boundaries within which prospective COOH-containing chemical series and COOH-containing drug candidates would be reasonably expected to possess properties suitable for oral administration. These specific "drug-like" property rules include molecular weight, the number of rotatable bonds, the number of hydrogen bond donors and acceptors, predictions of lipophilic character (calculated log P and log D values), topological polar surface area (TPSA), and the pK(a) value of the carboxylate moiety. Similar to the various sets of criteria that have emerged over the past decade and which have significantly reshaped the way medicinal chemists think about preferred drug chemical space, we propose these specific COOH "drug-like" property rules as a guide for the design of superior COOH-containing drug candidates and as a tool to better manage the liabilities generally associated with the presence of a COOH moiety.

The effect of the P1 side chain on the binding of optimized carboxylate and activated carbonyl inhibitors of the hepatitis C virus NS3 protease.

Peptidyl inhibitors of the hepatitis C virus NS3 protease hold much promise as direct-acting antiviral agents against hepatitis C infection. The optimization of N-terminal cleavage products, found to exhibit activity (product inhibition) against the enzyme, has led to potent tripeptide inhibitors that bear free C-terminal carboxylate groups. An analogous activated carbonyl compound (pentafluoroethyl ketone) bearing a P1 norvaline (Nva) was found to possess comparable activity against hepatitis C virus protease. However, an analogue bearing an aminocyclopropylcarboxylic acid (Acca) P1 residue exhibited very poor activity. (19)F-NMR studies indicate that the propensity of the Acca-derived activated carbonyl to form hemiketals is only slightly reduced compared with that of a P1 Nva equivalent. These results, as well as molecular modeling studies, argue against steric hindrance of the nucleophilic attack of Ser-139 accounting for the poor mechanism-based inhibition by the former. We hypothesize that the conformational properties of the respective C-termini in the context of an adaptable active site account for the divergent P1 structure-activity relationships.

Efficacy, pharmacokinetics, and metabolism of tetrahydroquinoline inhibitors of Plasmodium falciparum protein farnesyltransferase.

New antimalarials are urgently needed. We have shown that tetrahydroquinoline (THQ) protein farnesyltransferase (PFT) inhibitors (PFTIs) are effective against the Plasmodium falciparum PFT and are effective at killing P. falciparum in vitro. Previously described THQ PFTIs had limitations of poor oral bioavailability and rapid clearance from the circulation of rodents. In this paper, we validate both the Caco-2 cell permeability model for predicting THQ intestinal absorption and the in vitro liver microsome model for predicting THQ clearance in vivo. Incremental improvements in efficacy, oral absorption, and clearance rate were monitored by in vitro tests; and these tests were followed up with in vivo absorption, distribution, metabolism, and excretion studies. One compound, PB-93, achieved cure when it was given orally to P. berghei-infected rats every 8 h for a total of 72 h. However, PB-93 was rapidly cleared, and dosing every 12 h failed to cure the rats. Thus, the in vivo results corroborate the in vitro pharmacodynamics and demonstrate that 72 h of continuous high-level exposure to PFTIs is necessary to kill plasmodia. The metabolism of PB-93 was demonstrated by a novel technique that relied on double labeling with a radiolabel and heavy isotopes combined with radiometric liquid chromatography and mass spectrometry. The major liver microsome metabolite of PB-93 has the PFT Zn-binding N-methyl-imidazole removed; this metabolite is inactive in blocking PFT function. By solving the X-ray crystal structure of PB-93 bound to rat PFT, a model of PB-93 bound to malarial PFT was constructed. This model suggests areas of the THQ PFTIs that can be modified to retain efficacy and protect the Zn-binding N-methyl-imidazole from dealkylation.