Development and Comprehensive Benchmark of a High-Quality AMBER-Consistent Small Molecule Force Field with Broad Chemical Space Coverage for Molecular Modeling and Free Energy Calculation.

Biomolecular simulations have become an essential tool in contemporary drug discovery, and molecular mechanics force fields (FFs) constitute its cornerstone. Developing a high quality and broad coverage general FF is a significant undertaking that requires substantial expert knowledge and computing resources, which is beyond the scope of general practitioners. Existing FFs originate from only a limited number of groups and organizations, and they either suffer from limited numbers of training sets, lower than desired quality because of oversimplified representations, or are costly for the molecular modeling community to access. To address these issues, in this work, we developed an AMBER-consistent small molecule FF with extensive chemical space coverage, and we provide Open Access parameters for the entire modeling community. To validate our FF, we carried out benchmarks of quantum mechanics (QM)/molecular mechanics conformer comparison and free energy perturbation calculations on several benchmark data sets. Our FF achieves a higher level of performance at reproducing QM energies and geometries than two popular open-source FFs, OpenFF2 and GAFF2. In relative binding free energy calculations for 31 protein-ligand data sets, comprising 1079 pairs of ligands, the new FF achieves an overall root-mean-square error of 1.19 kcal/mol for ΔΔG and 0.92 kcal/mol for ΔG on a subset of 463 ligands without bespoke fitting to the data sets. The results are on par with those of the leading commercial series of OPLS FFs.

1,2-Difunctionalized bicyclo[1.1.1]pentanes: Long-sought-after mimetics for ortho/meta-substituted arenes.

The development of a versatile platform for the synthesis of 1,2-difunctionalized bicyclo[1.1.1]pentanes to potentially mimic ortho/meta-substituted arenes is described. The syntheses of useful building blocks bearing alcohol, amine, and carboxylic acid functional handles have been achieved from a simple common intermediate. Several ortho- and meta-substituted benzene analogs, as well as simple molecular matched pairs, have also been prepared using this platform. The results of in-depth ADME (absorption, distribution, metabolism, and excretion) investigations of these systems are presented, as well as computational studies which validate the ortho- or meta-character of these bioisosteres.

Characterization of Specific N-α-Acetyltransferase 50 (Naa50) Inhibitors Identified Using a DNA Encoded Library.

Two novel compounds were identified as Naa50 binders/inhibitors using DNA-encoded technology screening. Biophysical and biochemical data as well as cocrystal structures were obtained for both compounds (3a and 4a) to understand their mechanism of action. These data were also used to rationalize the binding affinity differences observed between the two compounds and a MLGP peptide-containing substrate. Cellular target engagement experiments further confirm the Naa50 binding of 4a and demonstrate its selectivity toward related enzymes (Naa10 and Naa60). Additional analogs of inhibitor 4a were also evaluated to study the binding mode observed in the cocrystal structures.

Sequential ALK Inhibitors Can Select for Lorlatinib-Resistant Compound ALK Mutations in ALK-Positive Lung Cancer.

The cornerstone of treatment for advanced ALK-positive lung cancer is sequential therapy with increasingly potent and selective ALK inhibitors. The third-generation ALK inhibitor lorlatinib has demonstrated clinical activity in patients who failed previous ALK inhibitors. To define the spectrum of ALK mutations that confer lorlatinib resistance, we performed accelerated mutagenesis screening of Ba/F3 cells expressing EML4-ALK. Under comparable conditions, N-ethyl-N-nitrosourea (ENU) mutagenesis generated numerous crizotinib-resistant but no lorlatinib-resistant clones harboring single ALK mutations. In similar screens with EML4-ALK containing single ALK resistance mutations, numerous lorlatinib-resistant clones emerged harboring compound ALK mutations. To determine the clinical relevance of these mutations, we analyzed repeat biopsies from lorlatinib-resistant patients. Seven of 20 samples (35%) harbored compound ALK mutations, including two identified in the ENU screen. Whole-exome sequencing in three cases confirmed the stepwise accumulation of ALK mutations during sequential treatment. These results suggest that sequential ALK inhibitors can foster the emergence of compound ALK mutations, identification of which is critical to informing drug design and developing effective therapeutic strategies.Significance: Treatment with sequential first-, second-, and third-generation ALK inhibitors can select for compound ALK mutations that confer high-level resistance to ALK-targeted therapies. A more efficacious long-term strategy may be up-front treatment with a third-generation ALK inhibitor to prevent the emergence of on-target resistance. Cancer Discov; 8(6); 714-29. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 663.

Resensitization to Crizotinib by the Lorlatinib ALK Resistance Mutation L1198F.

In a patient who had metastatic anaplastic lymphoma kinase (ALK)-rearranged lung cancer, resistance to crizotinib developed because of a mutation in the ALK kinase domain. This mutation is predicted to result in a substitution of cysteine by tyrosine at amino acid residue 1156 (C1156Y). Her tumor did not respond to a second-generation ALK inhibitor, but it did respond to lorlatinib (PF-06463922), a third-generation inhibitor. When her tumor relapsed, sequencing of the resistant tumor revealed an ALK L1198F mutation in addition to the C1156Y mutation. The L1198F substitution confers resistance to lorlatinib through steric interference with drug binding. However, L1198F paradoxically enhances binding to crizotinib, negating the effect of C1156Y and resensitizing resistant cancers to crizotinib. The patient received crizotinib again, and her cancer-related symptoms and liver failure resolved. (Funded by Pfizer and others; ClinicalTrials.gov number, NCT01970865.).

Effect of water solvation on the lipophilicity of isomeric pyrimidine-carboxamides.

Incorporation of nitrogen is a common medicinal chemistry tactic to reduce logD values. Neighboring group participation influences logD, so the results are isomer dependent. The logD and logP differences observed between isomeric pyrimidines 1, 2 and 3 presumably result when the carbonyl or ether lone pairs are in close proximity to a heterocyclic nitrogen lone pair, recruiting water to bridge between the electron rich atoms. Various lipophilicity calculators did not discriminate between 1 (logD=2.6) and 3 (logD=1.0), but solvation energies using Poisson-Boltzmann and 3D-RISM methods rationalize the observed differences in lipophilicity among pyrimidine carboxamide isomers.

Discovery of (10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(metheno)pyrazolo[4,3-h][2,5,11]-benzoxadiazacyclotetradecine-3-carbonitrile (PF-06463922), a macrocyclic inhibitor of anaplastic lymphoma kinase (ALK) and c-ros oncogene 1 (ROS1) with preclinical brain exposure and broad-spectrum potency against ALK-resistant mutations.

Although crizotinib demonstrates robust efficacy in anaplastic lymphoma kinase (ALK)-positive non-small-cell lung carcinoma patients, progression during treatment eventually develops. Resistant patient samples revealed a variety of point mutations in the kinase domain of ALK, including the L1196M gatekeeper mutation. In addition, some patients progress due to cancer metastasis in the brain. Using structure-based drug design, lipophilic efficiency, and physical-property-based optimization, highly potent macrocyclic ALK inhibitors were prepared with good absorption, distribution, metabolism, and excretion (ADME), low propensity for p-glycoprotein 1-mediated efflux, and good passive permeability. These structurally unusual macrocyclic inhibitors were potent against wild-type ALK and clinically reported ALK kinase domain mutations. Significant synthetic challenges were overcome, utilizing novel transformations to enable the use of these macrocycles in drug discovery paradigms. This work led to the discovery of 8k (PF-06463922), combining broad-spectrum potency, central nervous system ADME, and a high degree of kinase selectivity.

Design of potent and selective inhibitors to overcome clinical anaplastic lymphoma kinase mutations resistant to crizotinib.

Crizotinib (1), an anaplastic lymphoma kinase (ALK) receptor tyrosine kinase inhibitor approved by the U.S. Food and Drug Administration in 2011, is efficacious in ALK and ROS positive patients. Under pressure of crizotinib treatment, point mutations arise in the kinase domain of ALK, resulting in resistance and progressive disease. The successful application of both structure-based and lipophilic-efficiency-focused drug design resulted in aminopyridine 8e, which was potent across a broad panel of engineered ALK mutant cell lines and showed suitable preclinical pharmacokinetics and robust tumor growth inhibition in a crizotinib-resistant cell line (H3122-L1196M).

Discovery of novel, potent, and selective inhibitors of 3-phosphoinositide-dependent kinase (PDK1).

Analogues substituted with various amines at the 6-position of the pyrazine ring on (4-amino-7-isopropyl-7H-pyrrolo[2,3-d]pyrimidin-5-yl)pyrazin-2-ylmethanone were discovered as potent and selective inhibitors of PDK1 with potential as anticancer agents. An early lead with 2-pyridine-3-ylethylamine as the pyrazine substituent showed moderate potency and selectivity. Structure-based drug design led to improved potency and selectivity against PI3Kα through a combination of cyclizing the ethylene spacer into a saturated, five-membered ring and substituting on the 4-position of the aryl ring with a fluorine. ADME properties were improved by lowering the lipophilicity with heteroatom replacements in the saturated, five-membered ring. The optimized analogues have a PDK1 Ki of 1 nM and >100-fold selectivity against PI3K/AKT-pathway kinases. The cellular potency of these analogues was assessed by the inhibition of AKT phosphorylation (T308) and by their antiproliferation activity against a number of tumor cell lines.

Discovery of aryloxy tetramethylcyclobutanes as novel androgen receptor antagonists.

An aryloxy tetramethylcyclobutane was identified as a novel template for androgen receptor (AR) antagonists via cell-based high-throughput screening. Follow-up to the initial "hit" established 5 as a viable lead. Further optimization to achieve full AR antagonism led to the discovery of 26 and 30, both of which demonstrated excellent in vivo tumor growth inhibition upon oral administration in a castration-resistant prostate cancer (CRPC) animal model.

Quantifying and predicting the promiscuity and isoform specificity of small-molecule cytochrome P450 inhibitors.

Drug promiscuity (i.e., inhibition of multiple enzymes by a single compound) is increasingly recognized as an important pharmacological consideration in the drug development process. However, systematic studies of functional or physicochemical characteristics that correlate with drug promiscuity are handicapped by the lack of a good way of quantifying promiscuity. In this article, we present a new entropy-based index of drug promiscuity. We apply this index to two high-throughput data sets describing inhibition of cytochrome P450 isoforms by small-molecule drugs and drug candidates, and we demonstrate how drug promiscuity or specificity can be quantified. For these drug-metabolizing enzymes, we find that there is essentially no correlation between a drug's potency and specificity. We also present an index to quantify the susceptibilities of different enzymes to inhibition by diverse substrates. Finally, we use partial least-squares regression to successfully predict isoform specificity and promiscuity of small molecules, using a set of fingerprint-based descriptors.