Finding the sweet spot: the role of nature and nurture in medicinal chemistry

Impact of preformulation on drug development

INTRODUCTION

Preformulation assists scientists in screening lead candidates based on their physicochemical and biopharmaceutical properties. This data is useful for selection of new chemical entities (NCEs) for preclinical efficacy/toxicity studies which is a major section under investigational new drug application. A strong collaboration between discovery and formulation group is essential for selecting right NCEs in order to reduce attrition rate in the late stage development.

AREAS COVERED

This article describes the significance of preformulation research in drug discovery and development. Various crucial preformulation parameters with case studies have been discussed.

EXPERT OPINION

Physicochemical and biopharmaceutical characterization of NCEs is a decisive parameter during product development. Early prediction of these properties helps in selecting suitable physical form (salt, polymorph, etc.) of the candidate. Based on pharmacokinetic and efficacy/toxicity studies, suitable formulation for Phase I clinical studies can be developed. Overall these activities contribute in streamlining efficacy/toxicology evaluation, allowing pharmacologically effective and developable molecules to reach the clinic and eventually to the market. In this review, the magnitude of understanding preformulation properties of NCEs and their utility in product development has been elaborated with case studies.

Discovery of O-Alkylamino Tethered Niclosamide Derivatives as Potent and Orally Bioavailable Anticancer Agents

Niclosamide has been identified to potently inhibit the activation, nuclear translocation, and transactivation of STAT3. Nevertheless, the poor aqueous solubility and bioavailability of niclosamide has hindered its further clinical development for cancer therapy. To discover new molecules with enhanced drug-like properties, a series of novel O-alkylamino tethered derivatives of niclosamide have been designed, synthesized, and biologically evaluated. Among them, compound 11 (HJC0152) has been demonstrated to significantly suppress MDA-MB-231 xenograft tumor growth in vivo (i.p. & p.o.), indicating its great potential as efficacious and orally bioavailable therapeutics for human cancer.

A combined bioinformatics and chemoinformatics approach for developing asymmetric bivalent AMPA receptor positive allosteric modulators as neuroprotective agents

PAMs new in town! An effective, combined bioinformatics and chemoinformatics approach was applied to the design of novel asymmetric bivalent α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor positive allosteric modulators (PAMs) with marked potency in vitro and efficacy in vivo for preventing neuroapoptosis. The novel chemotype could provide pharmacological probes and potential therapeutic agents for glutamatergic hypofunction and its related neurological and psychiatric disorders.

Investigation of the influence of molecular topology on ligand binding

Molecular topology class has previously been put forward as a new concept of describing compound quality and it has been shown that compared to general bioactive compounds, drugs is more similar to natural products and human metabolites in terms of molecular topology class distribution, in which they are enriched with compounds having only one ring system. To further understand how the molecular topology is influencing the drug discovery process, we have investigated the compound potency of different molecular topologies in published chemical patents. Our study shows that the potency for compounds having one ring system is higher compared to compounds that have more than one ring system. Compounds with one ring system are significantly less lipophilic and smaller compared to compounds with several ring systems. Further the influence of the molecular topology on ligand efficiency (LE), ligand lipophilic efficiency (LLE) and ligand-efficiency-dependent lipophilicity (LELP) was also analyzed and it was found that on average compounds with fewer ring systems and in particular compounds with only one ring system show consistently better LE, LLE and LELP. The results suggest that compounds with fewer ring systems and in particular compounds with only one ring system have good properties and that they might be useful starting point for drug discovery projects.

Contributions of molecular properties to drug promiscuity

In contrast to designed polypharmacology that can result in efficient drugs for complex disorders, unintended drug promiscuity has detrimental contribution to side effects and toxicology. Characterization of promiscuous compounds enhances the understanding of complex interaction patterns and aids the design of compounds with broader selectivity against off-targets that has a major impact on medicinal chemistry outcome. In this Miniperspective we provide insights to the effect of physicochemical parameters on promiscuity. Information collected from recent, large-scale in vitro studies enabled us to discuss the relationships between physicochemical properties and promiscuity in detail. In light of these data, lipophilicity and basic character have the highest influence. On the basis of the accumulated knowledge, we propose the extensive use of pre- and postsynthesis metrics, as well as strict control of physicochemical properties during medicinal chemistry optimizations.

γ-Secretase Modulators: Can We Combine Potency with Safety?

γ-Secretase modulation has been proposed as a potential disease modifying anti-Alzheimer's approach. γ-Secretase modulators (GSMs) cause a product shift from the longer amyloid-beta (Aβ) peptide isoforms to shorter, more soluble, and less amyloidogenic isoforms, without inhibiting APP or Notch proteolytic processing. As such, modulating γ-secretase may avoid some of the adverse effects observed with γ-secretase inhibitors. Since the termination of the GSM tarenfurbil in 2008 due to negative phase III trial results, a considerable progress has been made towards more potent and better brain penetrable compounds. However, an analysis of their lipophilic efficiency indices indicates that their increased potency can be largely attributed to their increased lipophilicity. The need for early and chronic dosing with GSMs will require high-safety margins. This will be a challenge to achieve with the current, highly lipophilic GSMs. We will demonstrate that by focusing on the drug-like properties of GSMs, a combination of high in vitro potency and reduced lipophilicity can be achieved and does result in better tolerated compounds. The next hurdle will be to translate this knowledge into GSMs which are highly efficacious and safe in vivo.

Förster resonance energy transfer competitive displacement assay for human soluble epoxide hydrolase

The soluble epoxide hydrolase (sEH), responsible for the hydrolysis of various fatty acid epoxides to their corresponding 1,2-diols, is becoming an attractive pharmaceutical target. These fatty acid epoxides, particularly epoxyeicosatrienoic acids (EETs), play an important role in human homeostatic and inflammation processes. Therefore, inhibition of human sEH, which stabilizes EETs in vivo, brings several beneficial effects to human health. Although there are several catalytic assays available to determine the potency of sEH inhibitors, measuring the in vitro inhibition constant (K(i)) for these inhibitors using catalytic assay is laborious. In addition, k(off), which has been recently suggested to correlate better with the in vivo potency of inhibitors, has never been measured for sEH inhibitors. To better measure the potency of sEH inhibitors, a reporting ligand, 1-(adamantan-1-yl)-3-(1-(2-(7-hydroxy-2-oxo-2H-chromen-4-yl)acetyl) piperidin-4-yl)urea (ACPU), was designed and synthesized. With ACPU, we have developed a Förster resonance energy transfer (FRET)-based competitive displacement assay using intrinsic tryptophan fluorescence from sEH. In addition, the resulting assay allows us to measure the K(i) values of very potent compounds to the picomolar level and to obtain relative k(off) values of the inhibitors. This assay provides additional data to evaluate the potency of sEH inhibitors.

3D-QSAR-assisted drug design: identification of a potent quinazoline-based Aurora kinase inhibitor

We describe the 3D-QSAR-assisted design of an Aurora kinase A inhibitor with improved physicochemical properties, in vitro activity, and in vivo pharmacokinetic profiles over those of the initial lead. Three different 3D-QSAR models were built and validated by using a set of 66 pyrazole (Model I) and furanopyrimidine (Model II) compounds with IC(50) values toward Aurora kinase A ranging from 33 nM to 10.5 μM. The best 3D-QSAR model, Model III, constructed with 24 training set compounds from both series, showed robustness (r(2) (CV) =0.54 and 0.52 for CoMFA and CoMSIA, respectively) and superior predictive capacity for 42 test set compounds (R(2) (pred) =0.52 and 0.67, CoMFA and CoMSIA). Superimposition of CoMFA and CoMSIA Model III over the crystal structure of Aurora kinase A suggests the potential to improve the activity of the ligands by decreasing the steric clash with Val147 and Leu139 and by increasing hydrophobic contact with Leu139 and Gly216 residues in the solvent-exposed region of the enzyme. Based on these suggestions, the rational redesign of furanopyrimidine 24 (clog P=7.41; Aurora A IC(50) =43 nM; HCT-116 IC(50) =400 nM) led to the identification of quinazoline 67 (clog P=5.28; Aurora A IC(50) =25 nM; HCT-116 IC(50) =23 nM). Rat in vivo pharmacokinetic studies showed that 67 has better systemic exposure after i.v. administration than 24, and holds potential for further development.

Lysosomotropic properties of weakly basic anticancer agents promote cancer cell selectivity in vitro

Drug distribution in cells is a fundamentally important, yet often overlooked, variable in drug efficacy. Many weakly basic anticancer agents accumulate extensively in the acidic lysosomes of normal cells through ion trapping. Lysosomal trapping reduces the activity of anticancer drugs, since anticancer drug targets are often localized in the cell cytosol or nucleus. Some cancer cells have defective acidification of lysosomes, which causes a redistribution of trapped drugs from the lysosomes to the cytosol. We have previously established that such differences in drug localization between normal and cancer cells can contribute to the apparent selectivity of weakly basic drugs to cancer cells in vitro. In this work, we tested whether this intracellular distribution-based drug selectivity could be optimized based on the acid dissociation constant (pKa) of the drug, which is one of the determinants of lysosomal sequestration capacity. We synthesized seven weakly basic structural analogs of the Hsp90 inhibitor geldanamycin (GDA) with pKa values ranging from 5 to 12. The selectivity of each analog was expressed by taking ratios of anti-proliferative IC₅₀ values of the inhibitors in normal fibroblasts to the IC₅₀ values in human leukemic HL-60 cells. Similar selectivity assessments were performed in a pair of cancer cell lines that differed in lysosomal pH as a result of siRNA-mediated alteration of vacuolar proton ATPase subunit expression. Optimal selectivity was observed for analogs with pKa values near 8. Similar trends were observed with commercial anticancer agents with varying weakly basic pKa values. These evaluations advance our understanding of how weakly basic properties can be optimized to achieve maximum anticancer drug selectivity towards cancer cells with defective lysosomal acidification in vitro. Additional in vivo studies are needed to examine the utility of this approach for enhancing selectivity.

Structural sweet spot for A1 adenosine receptor activation by truncated (N)-methanocarba nucleosides: receptor docking and potent anticonvulsant activity

A(1) adenosine receptor (AR) agonists display antiischemic and antiepileptic neuroprotective activity, but peripheral cardiovascular side effects impeded their development. SAR study of N(6)-cycloalkylmethyl 4'-truncated (N)-methanocarba-adenosines identified 10 (MRS5474, N(6)-dicyclopropylmethyl, K(i) = 47.9 nM) as a moderately A(1)AR-selective full agonist. Two stereochemically defined N(6)-methynyl group substituents displayed narrow SAR; groups larger than cyclobutyl greatly reduced AR affinity, and those larger or smaller than cyclopropyl reduced A(1)AR selectivity. Nucleoside docking to A(1)AR homology model characterized distinct hydrophobic cyclopropyl subpockets, the larger "A" forming contacts with Thr270 (7.35), Tyr271 (7.36), Ile274 (7.39), and carbon chains of glutamates (EL2) and the smaller subpocket "B" forming contacts between TM6 and TM7. 10 suppressed minimal clonic seizures (6 Hz mouse model) without typical rotarod impairment of A(1)AR agonists. Truncated nucleosides, an appealing preclinical approach, have more druglike physicochemical properties than other A(1)AR agonists. Thus, we identified highly restricted regions for substitution around N(6) suitable for an A(1)AR agonist with anticonvulsant activity.

Anti-influenza drug discovery: identification of an orally bioavailable quinoline derivative through activity- and property-guided lead optimization

From a high-throughput screening (HTS) hit with inhibitory activity against virus-induced cytophathic in the low micromolar range, we have developed a potent anti-influenza lead through careful optimization without compromising the drug-like properties of the compound. An orally bioavailable compound was identified as a lead agent with nanomolar activity against influenza, representing a 140-fold improvement over the initial hit.

AtlasCBS: a web server to map and explore chemico-biological space

New approaches are needed that can help decrease the unsustainable failure in small-molecule drug discovery. Ligand Efficiency Indices (LEI) are making a great impact on early-stage compound selection and prioritization. Given a target-ligand database with chemical structures and associated biological affinities/activities for a target, the AtlasCBS server generates two-dimensional, dynamical representations of its contents in terms of LEI. These variables allow an effective decoupling of the chemical (angular) and biological (radial) components. BindingDB, PDBBind and ChEMBL databases are currently implemented. Proprietary datasets can also be uploaded and compared. The utility of this atlas-like representation in the future of drug design is highlighted with some examples. The web server can be accessed at http://ub.cbm.uam.es/atlascbs and https://www.ebi.ac.uk/chembl/atlascbs.