Vinylpyridine as a Tunable Covalent Warhead Targeting C797 in EGFR.

To further facilitate the discovery of cysteine reactive covalent inhibitors, there is a need to develop new reactive groups beyond the traditional acrylamide-type warheads. Herein we describe the design and synthesis of covalent EGFR inhibitors that use vinylpyridine as the reactive group. Structure-based design identified the quinazoline-containing vinylpyridine 6 as a starting point. Further modifications focused on reducing reactivity resulted in substituted vinyl compound 12, which shows high EGFR potency and good kinase selectivity, as well as significantly reduced reactivity compared to the starting compound 6, confirming that vinylpyridines can be applied as an alternative cysteine reactive warhead with tunable reactivity.

Augmenting DMTA using predictive AI modelling at AstraZeneca.

Design-Make-Test-Analyse (DMTA) is the discovery cycle through which molecules are designed, synthesised, and assayed to produce data that in turn are analysed to inform the next iteration. The process is repeated until viable drug candidates are identified, often requiring many cycles before reaching a sweet spot. The advent of artificial intelligence (AI) and cloud computing presents an opportunity to innovate drug discovery to reduce the number of cycles needed to yield a candidate. Here, we present the Predictive Insight Platform (PIP), a cloud-native modelling platform developed at AstraZeneca. The impact of PIP in each step of DMTA, as well as its architecture, integration, and usage, are discussed and used to provide insights into the future of drug discovery.

Identification and analyses of inhibitors targeting apolipoprotein(a) kringle domains KIV-7, KIV-10, and KV provide insight into kringle domain function.

Increased plasma concentrations of lipoprotein(a) (Lp(a)) are associated with an increased risk for cardiovascular disease. Lp(a) is composed of apolipoprotein(a) (apo(a)) covalently bound to apolipoprotein B of low-density lipoprotein (LDL). Many of apo(a)'s potential pathological properties, such as inhibition of plasmin generation, have been attributed to its main structural domains, the kringles, and have been proposed to be mediated by their lysine-binding sites. However, available small-molecule inhibitors, such as lysine analogs, bind unselectively to kringle domains and are therefore unsuitable for functional characterization of specific kringle domains. Here, we discovered small molecules that specifically bind to the apo(a) kringle domains KIV-7, KIV-10, and KV. Chemical synthesis yielded compound AZ-05, which bound to KIV-10 with a Kd of 0.8 µm and exhibited more than 100-fold selectivity for KIV-10, compared with the other kringle domains tested, including plasminogen kringle 1. To better understand and further improve ligand selectivity, we determined the crystal structures of KIV-7, KIV-10, and KV in complex with small-molecule ligands at 1.6-2.1 Å resolutions. Furthermore, we used these small molecules as chemical probes to characterize the roles of the different apo(a) kringle domains in in vitro assays. These assays revealed the assembly of Lp(a) from apo(a) and LDL, as well as potential pathophysiological mechanisms of Lp(a), including (i) binding to fibrin, (ii) stimulation of smooth-muscle cell proliferation, and (iii) stimulation of LDL uptake into differentiated monocytes. Our results indicate that a small-molecule inhibitor targeting the lysine-binding site of KIV-10 can combat the pathophysiological effects of Lp(a).

Structural analysis identifies an escape route from the adverse lipogenic effects of liver X receptor ligands.

Liver X receptors (LXRs) are attractive drug targets for cardiovascular disease treatment due to their role in regulating cholesterol homeostasis and immunity. The anti-atherogenic properties of LXRs have prompted development of synthetic ligands, but these cause major adverse effects-such as increased lipogenesis-which are challenging to dissect from their beneficial activities. Here we show that LXR compounds displaying diverse functional responses in animal models induce distinct receptor conformations. Combination of hydrogen/deuterium exchange mass spectrometry and multivariate analysis allowed identification of LXR regions differentially correlating with anti-atherogenic and lipogenic activities of ligands. We show that lipogenic compounds stabilize active states of LXRα and LXRß while the anti-atherogenic expression of the cholesterol transporter ABCA1 is associated with the ligand-induced stabilization of LXRα helix 3. Our data indicates that avoiding ligand interaction with the activation helix 12 while engaging helix 3 may provide directions for development of ligands with improved therapeutic profiles.

An Angle on MK2 Inhibition-Optimization and Evaluation of Prevention of Activation Inhibitors.

The mitogen-activated protein kinase p38α pathway has been an attractive target for the treatment of inflammatory conditions such as rheumatoid arthritis. While a number of p38α inhibitors have been taken to the clinic, they have been limited by their efficacy and toxicological profile. A lead identification program was initiated to selectively target prevention of activation (PoA) of mitogen-activated protein kinase-activated protein kinase 2 (MK2) rather than mitogen- and stress-activated protein kinase 1 (MSK1), both immediate downstream substrates of p38α, to improve the efficacy/safety profile over direct p38α inhibition. Starting with a series of pyrazole amide PoA MK2 inhibitor leads, and guided by structural chemistry and rational design, a highly selective imidazole 9 (2-(3'-(2-amino-2-oxoethyl)-[1,1'-biphenyl]-3-yl)-N-(5-(N,N-dimethylsulfamoyl)-2-methylphenyl)-1-propyl-1H-imidazole-5-carboxamide) and the orally bioavailable imidazole 18 (3-methyl-N-(2-methyl-5-sulfamoylphenyl)-2-(o-tolyl)imidazole-4-carboxamide) were discovered. The PoA concept was further evaluated by protein immunoblotting, which showed that the optimized PoA MK2 compounds, despite their biochemical selectivity against MSK1 phosphorylation, behaved similarly to p38 inhibitors in cellular signaling. This study highlights the importance of selective tool compounds in untangling complex signaling pathways, and although 9 and 18 were not differentiated from p38α inhibitors in a cellular context, they are still useful tools for further research directed to understand the role of MK2 in the p38α signaling pathway.

A hierarchical screening approach to enantiomeric separation.

The screening of a number of chiral stationary phases (CSPs) with different modifiers in supercritical fluid chromatography to find a chromatographic method for separation of enantiomers can be time-consuming. Computational methods for data analysis were utilized to establish a hierarchical screening strategy, using a dataset of 110 drug-like chiral compounds with diverse structures tested on 15 CSPs with two different modifiers. This dataset was analyzed using a combinatorial algorithm, principal component analysis (PCA), and a correlation matrix. The primary goal was to find a set of eight columns resolving a large number of compounds, but also having complementary enantioselective properties. In addition to the hereby defined hierarchical experimental strategy, quantitative structure enantioselective models (QSERs) were evaluated. The diverse chemical space and relatively limited size of the training set reduced the accuracy of the QSERs. However, including separation factors from other CSPs increased the accuracies of the QSERs substantially. Hence, such combined models can support the experimental strategy in prioritizing the CSPs of the second screening phase, when a compound is not separated by the primary set of columns.

Matched Molecular Pair Analysis in Short: Algorithms, Applications and Limitations.

Molecular matched pair (MMP) analysis has been used for more than 40 years within molecular design and is still an important tool to analyse potency data and other compound properties. The methods used to find matched pairs range from manual inspection, through supervised methods to unsupervised methods, which are able to find previously unknown molecular pairs. Recent publications demonstrate the value of automatic MMP analysis of publicly available bioactivity databases. The MMP concept has its limitations, but because of its easy to use and intuitive nature, it will remain one of the most important tools in the toolbox of many drug designers.

Predicting the relative binding affinity of mineralocorticoid receptor antagonists by density functional methods.

In drug discovery, prediction of binding affinity ahead of synthesis to aid compound prioritization is still hampered by the low throughput of the more accurate methods and the lack of general pertinence of one method that fits all systems. Here we show the applicability of a method based on density functional theory using core fragments and a protein model with only the first shell residues surrounding the core, to predict relative binding affinity of a matched series of mineralocorticoid receptor (MR) antagonists. Antagonists of MR are used for treatment of chronic heart failure and hypertension. Marketed MR antagonists, spironolactone and eplerenone, are also believed to be highly efficacious in treatment of chronic kidney disease in diabetes patients, but is contra-indicated due to the increased risk for hyperkalemia. These findings and a significant unmet medical need among patients with chronic kidney disease continues to stimulate efforts in the discovery of new MR antagonist with maintained efficacy but low or no risk for hyperkalemia. Applied on a matched series of MR antagonists the quantum mechanical based method gave an R(2) = 0.76 for the experimental lipophilic ligand efficiency versus relative predicted binding affinity calculated with the M06-2X functional in gas phase and an R(2) = 0.64 for experimental binding affinity versus relative predicted binding affinity calculated with the M06-2X functional including an implicit solvation model. The quantum mechanical approach using core fragments was compared to free energy perturbation calculations using the full sized compound structures.

Experimental and Quantum Chemical Evaluations of Pyridine Oxidation Under Drug Development Stress Test Conditions.

The oxidation reaction of pyridine by hydrogen peroxides in water media was investigated by combining quantum chemical calculations and laboratory experiments. Pyridine was selected as a model system for aromatic amines that frequently occurs in drug molecules. Several different reaction conditions, commonly used in stress testing of drug molecules during drug development, were investigated to increase mechanistic insight to this class of oxidation reactions. Of special interest is to note that small amounts of acetonitrile, a regularly used cosolvent to keep poorly soluble drug molecules in water solution, could catalyze the oxidation reaction in the presence of hydrogen peroxide. Consequently, attention needs to be taken when comparing data from different stress test studies of amine oxidation by hydrogen peroxides at different pH, and with and without acetonitrile. In particular, they need to be controlled when identifying the proper intrinsic stability of the drug molecule.

A qualitative method for prediction of amine oxidation in methanol and water.

We have developed a predictive method, based on quantum chemical calculations, that qualitatively predicts N-oxidation by hydrogen peroxides in drug structures. The method uses linear correlations of two complementary approaches to estimate the activation barrier without calculating it explicitly. This method can therefore be automated as it avoids demanding transition state calculations. As such, it may be used by chemists without experience in molecular modeling and provide additional understanding to experimental findings. The predictive method gives relative rates for N,N-dimethylbenzylamine and N-methylmorpholine in good agreement with experiments. In water, the experimental rate constants show that N,N-dimethylbenzylamine is oxidized three times faster than N-methylmorpholine and in methanol it is two times faster. The method suggests it to be two and five times faster, respectively. The method was also used to correlate experimental with predicted activation barriers, linear free-energy relationships, for a test set of tertiary amines. A correlation coefficient R(2) = 0.74 was obtained, where internal diagnostics in the method itself allowed identification of outliers. The method was applied to four drugs: caffeine, azelastine, buspirone, and clomipramine, all possessing several nitrogens. Both overall susceptibility and selectivity of oxidation were predicted, and verified by experiments.

Discovery and characterization of MAPK-activated protein kinase-2 prevention of activation inhibitors.

Two structurally distinct series of novel, MAPK-activated kinase-2 prevention of activation inhibitors have been discovered by high throughput screening. Preliminary structure-activity relationship (SAR) studies revealed substructural features that influence the selective inhibition of the activation by p38α of the downstream kinase MK2 in preference to an alternative substrate, MSK1. Enzyme kinetics, surface plasmon resonance (SPR), 2D protein NMR, and X-ray crystallography were used to determine the binding mode and the molecular mechanism of action. The compounds bind competitively to the ATP binding site of p38α but unexpectedly with higher affinity in the p38α-MK2 complex compared with p38α alone. This observation is hypothesized to be the origin of the substrate selectivity. The two lead series identified are suitable for further investigation for their potential to treat chronic inflammatory diseases with improved tolerability over previously studied p38α inhibitors.

Discovery of the Fibrinolysis Inhibitor AZD6564, Acting via Interference of a Protein-Protein Interaction.

A class of novel oral fibrinolysis inhibitors has been discovered, which are lysine mimetics containing an isoxazolone as a carboxylic acid isostere. As evidenced by X-ray crystallography the inhibitors bind to the lysine binding site in plasmin thus preventing plasmin from binding to fibrin, hence blocking the protein-protein interaction. Optimization of the series, focusing on potency in human buffer and plasma clotlysis assays, permeability, and GABAa selectivity, led to the discovery of AZD6564 (19) displaying an in vitro human plasma clot lysis IC50 of 0.44 µM, no detectable activity against GABAa, and with DMPK properties leading to a predicted dose of 340 mg twice a day oral dosing in humans.

Prediction of drug candidates' sensitivity toward autoxidation: computational estimation of C-H dissociation energies of carbon-centered radicals.

A method to predict a compound's sensitivity toward autoxidation using bond dissociation energies for hydrogen abstraction is described. The methodology is based on quantum mechanics and has been validated with a small molecule test set. Through this work, it has been observed that stabilization of an incipient radical by more than a single functional group is normally required to trigger autoxidation. The method has also been used to understand salt effects, wherein protonation of a basic amine stabilizes proximal C-H bonds to autoxidation. It can be used to support understanding of autoxidation processes and can form a predictive role for propensity to form potentially genotoxic and other degradation products. An automated protocol has been developed that allows the nonspecialist to perform quantum chemical calculations. The protocol is robust to enable general usage such that drug-like molecules can be handled by the tool and give an answer in hours (up to some days) depending on the size of the molecule. The efficiency of the tool makes it possible to perform risk assessment for autoxidation of small lists of molecules and could typically be used for shortlisted candidates before drug nomination, during drug formulation development, and during due diligence for in-licensing compounds.

Synthesis and functionalization of cyclic sulfonimidamides: a novel chiral heterocyclic carboxylic Acid bioisostere.

An efficient synthesis of aryl substituted cyclic sulfonimidamides designed as chiral nonplanar heterocyclic carboxylic acid bioisosteres is described. The cyclic sulfonimidamide ring system could be prepared in two steps from a trifluoroacetyl protected sulfinamide and methyl ester protected amino acids. By varying the amino acid, a range of different C-3 substituted sulfonimidamides could be prepared. The compounds could be further derivatized in the aryl ring using standard cross-coupling reactions to yield highly substituted cyclic sulfonimidamides in excellent yields. The physicochemical properties of the final compounds were examined and compared to those of the corresponding carboxylic acid and tetrazole derivatives. The unique nonplanar shape in combination with the relatively strong acidity (pK a 5-6) and the ease of modifying the chemical structure to fine-tune the physicochemical properties suggest that this heterocycle can be a valuable addition to the range of available carboxylic acid isosteres.

Discovery of cyclopentane- and cyclohexane-trans-1,3-diamines as potent melanin-concentrating hormone receptor 1 antagonists.

We herein report the optimization of cyclopentane- and cyclohexane-1,3-diamine derivatives as novel and potent MCH-R1 antagonists. Structural modifications of the 2-amino-quinoline and thiophene moieties found in the initial lead compound served to improve its metabolic stability profile and MCH-R1 affinity, and revealed unprecedented SAR when compared to other 2-amino-quinoline-containing MCH-R1 antagonists.