Atropisomerism in Drug Discovery: A Medicinal Chemistry Perspective Inspired by Atropisomeric Class I PI3K Inhibitors.

Atropisomerism is a type of axial chirality resulting from hindered rotation about a σ bond that gives rise to nonsuperimposable stereoisomers (termed "atropisomers"). The inversion of chirality of an atropisomeric axis is a time- and temperature-dependent dynamic process occurring by simple bond rotation. For this reason, the rotational energy barrier (ΔErot) and the interconversion rate between an atropisomeric pair of biologically active molecules are important parameters to consider in drug discovery.Many compounds with atropisomeric axes advance into development every year. The vast majority of them have low rotational energy barriers (ΔErot lower than 20 kcal/mol), meaning they are rapidly equilibrating conformers and considered achiral (class 1 atropisomers). Compounds in class 2 (ΔErot = 20 to 30 kcal/mol) can be challenging to develop given that the stereochemical integrity of the atropisomeric axes can be compromised over time. It has been recommended that small molecule drug candidates containing one or more atropisomeric axes with rotational energy barriers greater than 30 kcal/mol (class 3 atropisomers) should be developed as single atropisomers.In medicinal chemistry, a σ bond with restricted rotation is engineered into a bioactive molecule primarily to limit its number of accessible conformations, thereby minimizing entropic and/or enthalpic energy penalties associated with biological target binding. In addition to enhanced pharmacology, potential positive outcomes of introducing atropisomerism include improved physicochemical properties and superior pharmacokinetics/ADME profiles. The application of atropisomerism in medicinal chemistry has become increasingly enabled due to recent advances in synthesis, purification, and analysis, as described in this special issue and recent review articles.Herein, we discuss two case studies from our own work in which restricting rotation about axes of atropisomerism led to significant improvements in pharmacological, physicochemical, and ADME properties for different series of PI3K inhibitors. In the first instance, a restricted axis of rotation was designed to mitigate an acid-mediated hydrolytic degradation pathway observed in a series of PI3Kδ inhibitors. The conformational constraint disrupts conjugation between a quinazolinone and a pyridine, leading to improved chemical stability under acidic conditions. In the second case study, introduction of a restricted axis of rotation between two heteroaromatic systems in a series of PI3Kß inhibitors generated pairs of atropisomeric compounds with significantly different biological activities. Advanced profiling also demonstrated clear substrate stereospecificity in regard to metabolism by aldehyde oxidase. Gratifyingly, the eutomer (more active atropisomer) shows significantly less susceptibility for oxidative metabolism relative to the distomer (less active atropisomer). The improvements in potency, selectivity, chemical stability, and metabolic stability discussed in this manuscript are all directly related to the concept of atropisomerism.

Discovery of an Atropisomeric PI3Kß Selective Inhibitor through Optimization of the Hinge Binding Motif.

A series of PI3Kß selective inhibitors derived from a novel 4-(1H-benzo[d]imidazol-1-yl)quinoline chemotype has been rationally designed. Crucial to achieving the desired selectivity over the other class I PI3K isoforms, including the challenging δ-isoform, was the identification of a subset of substituted pyridine hinge binders. This work led to the discovery of (P)-14, a highly selective and orally bioavailable PI3Kß inhibitor displaying an excellent pharmacokinetic profile in addition to great cellular potency in various PTEN-deficient tumor cell lines. Results from a dog toxicology study revealing structure-related, off-target ocular toxicity are also briefly discussed.

Discovery of Lanraplenib (GS-9876): A Once-Daily Spleen Tyrosine Kinase Inhibitor for Autoimmune Diseases.

Spleen tyrosine kinase (SYK) is a critical regulator of signaling in a variety of immune cell types such as B-cells, monocytes, and macrophages. Accordingly, there have been numerous efforts to identify compounds that selectively inhibit SYK as a means to treat autoimmune and inflammatory diseases. We previously disclosed GS-9973 (entospletinib) as a selective SYK inhibitor that is under clinical evaluation in hematological malignancies. However, a BID dosing regimen and drug interaction with proton pump inhibitors (PPI) prevented development of entospletinib in inflammatory diseases. Herein, we report the discovery of a second-generation SYK inhibitor, GS-9876 (lanraplenib), which has human pharmacokinetic properties suitable for once-daily administration and is devoid of any interactions with PPI. Lanraplenib is currently under clinical evaluation in multiple autoimmune indications.

Atropisomerism by Design: Discovery of a Selective and Stable Phosphoinositide 3-Kinase (PI3K) ß Inhibitor.

Atropisomerism is a type of axial chirality in which enantiomers or diastereoisomers arise due to hindered rotation around a bond axis. In this manuscript, we report a case in which torsional scan studies guided the thoughtful creation of a restricted axis of rotation between two heteroaromatic systems of a phosphoinositide 3-kinase (PI3K) ß inhibitor, generating a pair of atropisomeric compounds with significantly different pharmacological and pharmacokinetic profiles. Emblematic of these differences, the metabolism of inactive ( M)-28 is primarily due to the cytosolic enzyme aldehyde oxidase, while active ( P)-28 has lower affinity for aldehyde oxidase, resulting in substantially better metabolic stability. Additionally, we report torsional scan and experimental studies used to determine the barriers of rotation of this novel PI3Kß inhibitor.

Discovery of a Phosphoinositide 3-Kinase (PI3K) ß/δ Inhibitor for the Treatment of Phosphatase and Tensin Homolog (PTEN) Deficient Tumors: Building PI3Kß Potency in a PI3Kδ-Selective Template by Targeting Nonconserved Asp856.

Phosphoinositide 3-kinase (PI3K) ß signaling is required to sustain cancer cell growth in which the tumor suppressor phosphatase and tensin homolog (PTEN) has been deactivated. This manuscript describes the discovery, optimization, and in vivo evaluation of a novel series of PI3Kß/δ inhibitors in which PI3Kß potency was built in a PI3Kδ-selective template. This work led to the discovery of a highly selective PI3Kß/δ inhibitor displaying excellent pharmacokinetic profile and efficacy in a human PTEN-deficient LNCaP prostate carcinoma xenograft tumor model.

Interaction of 2,4-Diaminopyrimidine-Containing Drugs Including Fedratinib and Trimethoprim with Thiamine Transporters.

Inhibition of thiamine transporters has been proposed as a putative mechanism for the observation of Wernicke's encephalopathy and subsequent termination of clinical development of fedratinib, a Janus kinase inhibitor (JAKi). This study aimed to determine the potential for other JAKi to inhibit thiamine transport using human epithelial colorectal adenocarcinoma (Caco-2) and thiamine transporter (THTR) overexpressing cells and to better elucidate the structural basis for interacting with THTR. Only JAKi containing a 2,4-diaminopyrimidine were observed to inhibit thiamine transporters. Fedratinib inhibited thiamine uptake into Caco-2 cells (IC50 = 0.940 µM) and THTR-2 (IC50 = 1.36 µM) and, to a lesser extent, THTR-1 (IC50 = 7.10 µM) overexpressing cells. Two other JAKi containing this moiety, AZD1480 and cerdulatinib, were weaker inhibitors of the thiamine transporters. Other JAKi-including monoaminopyrimidines, such as momelotinib, and nonaminopyrimidines, such as filgotinib-did not have any inhibitory effects on thiamine transport. A pharmacophore model derived from the minimized structure of thiamine suggests that 2,4-diaminopyrimidine-containing compounds can adopt a conformation matching several key features of thiamine. Further studies with drugs containing a 2,4-diaminopyrimidine resulted in the discovery that the antibiotic trimethoprim also potently inhibits thiamine uptake mediated by THTR-1 (IC50 = 6.84 µM) and THTR-2 (IC50 = 5.56 µM). Fedratinib and trimethoprim were also found to be substrates for THTR, a finding with important implications for their disposition in the body. In summary, our results show that not all JAKi have the potential to inhibit thiamine transport and further establish the interaction of these transporters with xenobiotics.

Selectivity switch between FAK and Pyk2: Macrocyclization of FAK inhibitors improves Pyk2 potency.

Herein, we describe the synthesis of Pyk2 inhibitors via macrocyclization of FAK and dual Pyk2-FAK inhibitors. We identified macrocycle 25a as a highly potent Pyk2 inhibitor (IC50=0.7nM), with ∼175-fold improvement in Pyk2 potency as compared to its acyclic counterpart. In many cases, macrocyclization improved Pyk2 potency while weakening FAK potency, thereby improving the Pyk2/FAK selectivity ratio for this structural class of inhibitors. Various macrocyclic linkers were studied in an attempt to optimize Pyk2 selectivity. We observed macrocyclic atropisomerism during the synthesis of 19-membered macrocycles 10a-d, and successfully obtained crystallographic evidence of one atropisomer (10a-AtropB) preferentially bound to Pyk2.

Discovery of Orally Efficacious Phosphoinositide 3-Kinase δ Inhibitors with Improved Metabolic Stability.

Aberrant signaling of phosphoinositide 3-kinase δ (PI3Kδ) has been implicated in numerous pathologies including hematological malignancies and rheumatoid arthritis. Described in this manuscript are the discovery, optimization, and in vivo evaluation of a novel series of pyridine-containing PI3Kδ inhibitors. This work led to the discovery of 35, a highly selective inhibitor of PI3Kδ which displays an excellent pharmacokinetic profile and is efficacious in a rodent model of rheumatoid arthritis.

2,4,6-Triaminopyrimidine as a Novel Hinge Binder in a Series of PI3Kδ Selective Inhibitors.

Inhibition of phosphoinositide 3-kinase δ (PI3Kδ) is an appealing target for several hematological malignancies and inflammatory diseases. Herein, we describe the discovery and optimization of a series of propeller shaped PI3Kδ inhibitors comprising a novel triaminopyrimidine hinge binder. Combinations of electronic and structural strategies were employed to mitigate aldehyde oxidase mediated metabolism. This medicinal chemistry effort culminated in the identification of 52, a potent and highly selective inhibitor of PI3Kδ that demonstrates efficacy in a rat model of arthritis.

1-Benzylbenzimidazoles: the discovery of a novel series of bradykinin B(1) receptor antagonists.

The design, synthesis, and structure-activity studies of a novel series of BK B(1) receptor antagonists based on a 1-benzylbenzimidazole chemotype are described. A number of compounds, for example, 38g, with excellent affinity for the cynomolgus macaque and rat bradykinin B(1) receptor were discovered.

2-Arylpyrimidines: novel CRF-1 receptor antagonists.

The design, synthesis and structure-activity relationship studies of a novel series of CRF-1 receptor antagonists, the 2-arylpyrimidines, are described. The effects of substitution on the aromatic ring and the pyrimidine core on CRF-1 receptor binding were investigated. A number of compounds with K(i) values below 10 nM and lipophilicity in a minimally acceptable range for a CNS drug (cLogP<5) were discovered.

Aminoquinazolines as TRPV1 antagonists: modulation of drug-like properties through the exploration of 2-position substitution.

A focused SAR exploration of the lead 4-aminoquinazoline TRPV1 antagonist 2 led to the discovery of compound 18. In rats, compound 18 is readily absorbed following oral dosing and demonstrates excellent in vivo potency and efficacy in an acute inflammatory pain model.

Strategies and tactics for optimizing the Hit-to-Lead process and beyond--a computational chemistry perspective.

The Hit-to-Lead-to-Candidate process continues to evolve rapidly, and while technological advances offer much potential, the reality often pales to the promise. Conversely, strategies and tactics implementing existing technologies may result in more benefit in the end. This article focuses on some of the thinking and approaches that may improve the efficiency and effectiveness of the beginnings of the drug discovery path. From the perspective of computational chemists, different types of strategy and philosophy of approach will be treated including: considerations of early lead choices, strategies for improving poor leads, multivariate optimization, opportunities for informatics, and engineering good decisions.

The design, synthesis and structure-activity relationships of 1-aryl-4-aminoalkylisoquinolines: a novel series of CRF-1 receptor antagonists.

The design, synthesis and structure-activity relationships of a novel series of CRF-1 receptor antagonist, the 1-aryl-4-alkylaminoisoquinolines, is described. The effects of substitution on the aromatic ring, the amino group and the isoquinoline core on CRF-1 receptor binding were investigated.

From arylureas to biarylamides to aminoquinazolines: discovery of a novel, potent TRPV1 antagonist.

Bioisosteric replacement of piperazine with an aryl ring in lead VR1 antagonist 1 led to the biarylamide series. The development of B-ring SAR led to the conformationally constrained analog 70. The resulting aminoquinazoline 70 represents a novel VR1 antagonist with improved in vitro potency and oral bioavailability vs the analogous compounds from the lead series.

Contributions of conformational compression and preferential transition state stabilization to the rate enhancement by chorismate mutase.

The rate enhancement provided by the chorismate mutase (CM) enzyme for the Claisen rearrangement of chorismate to prephenate has been investigated by application of the concept of near attack conformations (NACs). Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82% and 100% of chorismate conformers were found to be NAC structures in water and in the CM active site, respectively. Consequently, the conversion of non-NACs to NACs does not contribute to the free energy of activation from preorganization of the substrate into NACs. The FEP calculations yielded differences in free energies of activation that well reproduce the experimental data. Additional calculations indicate that the rate enhancement by CM over the aqueous phase results primarily from conformational compression of NACs by the enzyme and that this process is enthalpically controlled. This suggests that preferential stabilization of the transition state in the enzyme environment relative to water plays a secondary role in the catalysis by CM.

Investigation of solvent effects for the Claisen rearrangement of chorismate to prephenate: mechanistic interpretation via near attack conformations.

Solvent effects on the rate of the Claisen rearrangement of chorismate to prephenate have been examined in water and methanol. The preequilibrium free-energy differences between diaxial and diequatorial conformers of chorismate, which had previously been implicated as the sole basis for the observed 100-fold rate increase in water over methanol, have been reframed using the near attack conformation (NAC) concept of Bruice and co-workers. Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82%, 57%, and 1% of chorismate conformers were found to be NAC structures (NACs) in water, methanol, and the gas phase, respectively. As a consequence, the conversion of non-NACs to NACs provides no free-energy contributions to the overall relative reaction rates in water versus methanol. Free-energy perturbation calculations yielded differences in free energies of activation for the two polar protic solvents and the gas phase. The rate enhancement in water over the gas phase arises from preferential hydration of the transition state (TS) relative to the reactants via increased hydrogen bonding and long-range electrostatic interactions, which accompany bringing the two negatively charged carboxylates into closer proximity. More specifically, there is an increase of 1.3 and 0.6 hydrogen bonds to the carboxylate groups and the ether oxygen, respectively, in going from the reactant to the TS in water. In methanol, the corresponding changes in hydrogen bonding with first shell solvent molecules are small; the rate enhancement arises primarily from the enhanced long-range interactions with solvent molecules. Thus, the reaction occurs faster in water than in methanol due to greater stabilization of the TS in water by specific interactions with first shell solvent molecules.

Design, synthesis, and discovery of 5-piperazinyl-1,2,6,7-tetrahydro-5H-azepino[3,2,1-hi]indol-4-one derivatives: a novel series of mixed dopamine D2/D4 receptor antagonist.

5-piperazinyl-1,2,6,7-tetrahydro-5H-azepino[3,2,1-hi]indol-4-one derivatives were designed, synthesized, and identified as a new series of mixed dopamine D(2)/D(4) receptor antagonists. This series featured a rigid tricyclic ring system as an important pharmacophore core structure for high binding affinity. Molecular modeling studies are also described.

PDDG/PM3 and PDDG/MNDO: improved semiempirical methods.

Two new semiempirical methods employing a Pairwise Distance Directed Gaussian modification have been developed: PDDG/PM3 and PDDG/MNDO; they are easily implemented in existing software, and yield heats of formation for compounds containing C, H, N, and O atoms with significantly improved accuracy over the standard NDDO schemes, PM5, PM3, AM1, and MNDO. The PDDG/PM3 results for heats of formation also show substantial improvement over density functional theory with large basis sets. The PDDG modifications consist of a single function, which is added to the existing pairwise core repulsion functions within PM3 and MNDO, a reparameterized semiempirical parameter set, and modified computation of the energy of formation of a gaseous atom. The PDDG addition introduces functional group information via pairwise atomic interactions using only atom-based parameters. For 622 diverse molecules containing C, H, N, and O atoms, mean absolute errors in calculated heats of formation are reduced from 4.4 to 3.2 kcal/mol and from 8.4 to 5.2 kcal/mol using the PDDG modified versions of PM3 and MNDO over the standard versions, respectively. Several specific problems are overcome, including the relative stability of hydrocarbon isomers, and energetics of small rings and molecules containing multiple heteroatoms. The internal consistency of PDDG energies is also significantly improved, enabling more reliable analysis of isomerization energies and trends across series of molecules; PDDG isomerization energies show significant improvement over B3LYP/6-31G* results. Comparison of heats of formation, ionization potentials, dipole moments, isomer, and conformer energetics, intermolecular interaction energies, activation energies, and molecular geometries from the PDDG techniques is made to experimental data and values from other semiempirical and ab initio methods.

Improved semiempirical heats of formation through the use of bond and group equivalents.

Deficiencies in energetics obtained using the common semiempirical methods, AM1, PM3, and MNDO, may partly be traced to the use of pseudoatomic equivalents for conversion of molecular energies to heats of formation at 298 K. We present an alternative scheme based on the use of bond and group equivalents. Values for the 61 bond and group equivalents necessary for treatment of molecules containing the common organic elements, hydrogen, carbon, nitrogen, and oxygen have been derived. For a set of 583 neutral, closed-shell molecules mean absolute errors in AM1, PM3, and MNDO heats of formation are reduced from 6.6, 4.2, and 8.2 kcal/mol to 2.3, 2.2, and 3.0 kcal/mol, respectively. Several systematic problems are overcome in the present scheme including relative stabilities of branched hydrocarbons, energetics of conjugated systems, heats of formation of long chain hydrocarbons, and enthalpies of molecules containing multiple heteroatoms. Although the approach is restricted to molecules with well-defined functional groups, the equivalents are easy to incorporate and are chemically relevant. This revised procedure allows semiempirical methods to be used for far more reliable evaluations of heats of reactions. Estimates are made of the errors inherent in these semiempirical formalisms, arising from integral approximations and the neglect of explicit treatment of electron correlation effects, while excluding those from inadequate parameterization.