The SOS1 Inhibitor MRTX0902 Blocks KRAS Activation and Demonstrates Antitumor Activity in Cancers Dependent on KRAS Nucleotide Loading.

KRAS is the most frequently mutated oncogene in human cancer and facilitates uncontrolled growth through hyperactivation of the receptor tyrosine kinase (RTK)/mitogen-activated protein kinase (MAPK) pathway. The Son of Sevenless homolog 1 (SOS1) protein functions as a guanine nucleotide exchange factor (GEF) for the RAS subfamily of small GTPases and represents a druggable target in the pathway. Using a structure-based drug discovery approach, MRTX0902 was identified as a selective and potent SOS1 inhibitor that disrupts the KRAS:SOS1 protein-protein interaction to prevent SOS1-mediated nucleotide exchange on KRAS and translates into an anti-proliferative effect in cancer cell lines with genetic alterations of the KRAS-MAPK pathway. MRTX0902 augmented the antitumor activity of the KRAS G12C inhibitor adagrasib when dosed in combination in eight out of 12 KRAS G12C-mutant human non-small cell lung cancer and colorectal cancer xenograft models. Pharmacogenomic profiling in preclinical models identified cell cycle genes and the SOS2 homolog as genetic co-dependencies and implicated tumor suppressor genes (NF1 and PTEN) in resistance following combination treatment. Lastly, combined vertical inhibition of RTK/MAPK pathway signaling by MRTX0902 with inhibitors of EGFR or RAF/MEK led to greater downregulation of pathway signaling and improved antitumor responses in KRAS-MAPK pathway-mutant models. These studies demonstrate the potential clinical application of dual inhibition of SOS1 and KRAS G12C and additional SOS1 combination strategies that will aide in the understanding of SOS1 and RTK/MAPK biology in targeted cancer therapy.

Discovery of Pyridopyrimidinones that Selectively Inhibit the H1047R PI3Kα Mutant Protein.

The H1047R mutation of PIK3CA is highly prevalent in breast cancers and other solid tumors. Selectively targeting PI3KαH1047R over PI3KαWT is crucial due to the role that PI3KαWT plays in normal cellular processes, including glucose homeostasis. Currently, only one PI3KαH1047R-selective inhibitor has progressed into clinical trials, while three pan mutant (H1047R, H1047L, H1047Y, E542K, and E545K) selective PI3Kα inhibitors have also reached the clinical stage. Herein, we report the design and discovery of a series of pyridopyrimidinones that inhibit PI3KαH1047R with high selectivity over PI3KαWT, resulting in the discovery of compound 17. When dosed in the HCC1954 tumor model in mice, 17 provided tumor regressions and a clear pharmacodynamic response. X-ray cocrystal structures from several PI3Kα inhibitors were obtained, revealing three distinct binding modes within PI3KαH1047R including a previously reported cryptic pocket in the C-terminus of the kinase domain wherein we observe a ligand-induced interaction with Arg1047.

Discovery of Five SOS2 Fragment Hits with Binding Modes Determined by SOS2 X-Ray Cocrystallography.

SOS1 and SOS2 are guanine nucleotide exchange factors that mediate RTK-stimulated RAS activation. Selective SOS1:KRAS PPI inhibitors are currently under clinical investigation, whereas there are no reports to date of SOS2:KRAS PPI inhibitors. SOS2 activity is implicated in MAPK rebound when divergent SOS1 mutant cell lines are treated with the SOS1 inhibitor BI-3406; therefore, SOS2:KRAS inhibitors are of therapeutic interest. In this report, we detail a fragment-based screening strategy to identify X-ray cocrystal structures of five diverse fragment hits bound to SOS2.

Fragment optimization and elaboration strategies - the discovery of two lead series of PRMT5/MTA inhibitors from five fragment hits.

Here we describe the early stages of a fragment-based lead discovery (FBLD) project for a recently elucidated synthetic lethal target, the PRMT5/MTA complex, for the treatment of MTAP-deleted cancers. Starting with five fragment/PRMT5/MTA X-ray co-crystal structures, we employed a two-phase fragment elaboration process encompassing optimization of fragment hits and subsequent fragment growth to increase potency, assess synthetic tractability, and enable structure-based drug design. Two lead series were identified, one of which led to the discovery of the clinical candidate MRTX1719.

Design and evaluation of achiral, non-atropisomeric 4-(aminomethyl)phthalazin-1(2H)-one derivatives as novel PRMT5/MTA inhibitors.

MRTX1719 is an inhibitor of the PRMT5/MTA complex and recently entered clinical trials for the treatment of MTAP-deleted cancers. MRTX1719 is a class 3 atropisomeric compound that requires a chiral synthesis or a chiral separation step in its preparation. Here, we report the SAR and medicinal chemistry design strategy, supported by structural insights from X-ray crystallography, to discover a class 1 atropisomeric compound from the same series that does not require a chiral synthesis or a chiral separation step in its preparation.

Design and Discovery of MRTX0902, a Potent, Selective, Brain-Penetrant, and Orally Bioavailable Inhibitor of the SOS1:KRAS Protein-Protein Interaction.

SOS1 is one of the major guanine nucleotide exchange factors that regulates the ability of KRAS to cycle through its "on" and "off" states. Disrupting the SOS1:KRASG12C protein-protein interaction (PPI) can increase the proportion of GDP-loaded KRASG12C, providing a strong mechanistic rationale for combining inhibitors of the SOS1:KRAS complex with inhibitors like MRTX849 that target GDP-loaded KRASG12C. In this report, we detail the design and discovery of MRTX0902─a potent, selective, brain-penetrant, and orally bioavailable SOS1 binder that disrupts the SOS1:KRASG12C PPI. Oral administration of MRTX0902 in combination with MRTX849 results in a significant increase in antitumor activity relative to that of either single agent, including tumor regressions in a subset of animals in the MIA PaCa-2 tumor mouse xenograft model.

Fragment-Based Discovery of MRTX1719, a Synthetic Lethal Inhibitor of the PRMT5•MTA Complex for the Treatment of MTAP-Deleted Cancers.

The PRMT5•MTA complex has recently emerged as a new synthetically lethal drug target for the treatment of MTAP-deleted cancers. Here, we report the discovery of development candidate MRTX1719. MRTX1719 is a potent and selective binder to the PRMT5•MTA complex and selectively inhibits PRMT5 activity in MTAP-deleted cells compared to MTAP-wild-type cells. Daily oral administration of MRTX1719 to tumor xenograft-bearing mice demonstrated dose-dependent inhibition of PRMT5-dependent symmetric dimethylarginine protein modification in MTAP-deleted tumors that correlated with antitumor activity. A 4-(aminomethyl)phthalazin-1(2H)-one hit was identified through a fragment-based screen, followed by X-ray crystallography, to confirm binding to the PRMT5•MTA complex. Fragment growth supported by structural insights from X-ray crystallography coupled with optimization of pharmacokinetic properties aided the discovery of development candidate MRTX1719.

Structure-guided optimization of a novel class of ASK1 inhibitors with increased sp3 character and an exquisite selectivity profile.

Apoptosis Signal-Regulating Kinase-1 (ASK1) is a known member of the Mitogen-Activated Protein Kinase Kinase Kinase (MAP3K) family and upon stimulation will activate the p38- and JNK-pathways leading to cardiac apoptosis, fibrosis, and hypertrophy. Using Structure-Based Drug Design (SBDD) in parallel with deconstruction of a published compound, a novel series of ASK1 inhibitors was optimized, which incorporated a saturated heterocycle proximal to the hinge-binding motif. This yielded a unique chemical series with excellent selectivity across the broader kinome, and desirable drug-like properties. The lead compound (10) is highly soluble and permeable, and exhibits a cellular EC50 = 24 nM and Kd < 1 nM. Of the 350 kinases tested, 10 has an IC50 ≤ 500 nM for only eight of them. This paper will describe the design hypotheses behind this series, key data points during the optimization phase, as well as a possible structural rationale for the kinome selectivity. Based on crystallographic data, the presence of an aliphatic cycle adjacent to the hinge-binder in the active site of the protein kinase showed up in <1% of the >5000 structures in the Protein Data Bank, potentially conferring the selectivity seen in this series.

Correction to "Structure-Based Design of ASK1 Inhibitors as Potential Agents for Heart Failure".

[This corrects the article DOI: 10.1021/acsmedchemlett.6b00481.].

1,2,4-Triazolo-[1,5-a]pyridine HIF Prolylhydroxylase Domain-1 (PHD-1) Inhibitors With a Novel Monodentate Binding Interaction.

Herein we describe the identification of 4-{[1,2,4]triazolo[1,5-a]pyridin-5-yl}benzonitrile-based inhibitors of the hypoxia-inducible factor prolylhydroxylase domain-1 (PHD-1) enzyme. These inhibitors were shown to possess a novel binding mode by X-ray crystallography, in which the triazolo N1 atom coordinates in a hitherto unreported monodentate interaction with the active site Fe2+ ion, while the benzonitrile group accepts a hydrogen-bonding interaction from the side chain residue of Asn315. Further optimization led to potent PHD-1 inhibitors with good physicochemical and pharmacokinetic properties.

Structure-Based Design of ASK1 Inhibitors as Potential Agents for Heart Failure.

Apoptosis signal-regulating kinase 1 (ASK1/MAP3K) is a mitogen-activated protein kinase family member shown to contribute to acute ischemia/reperfusion injury. Using structure-based drug design, deconstruction, and reoptimization of a known ASK1 inhibitor, a lead compound was identified. This compound displayed robust MAP3K pathway inhibition and reduction of infarct size in an isolated perfused heart model of cardiac injury.

High-resolution structure of the human GPR40 receptor bound to allosteric agonist TAK-875.

Human GPR40 receptor (hGPR40), also known as free fatty-acid receptor 1 (FFAR1), is a G-protein-coupled receptor that binds long-chain free fatty acids to enhance glucose-dependent insulin secretion. Novel treatments for type-2 diabetes mellitus are therefore possible by targeting hGPR40 with partial or full agonists. TAK-875, or fasiglifam, is an orally available, potent and selective partial agonist of hGPR40 receptor, which reached phase III clinical trials for the potential treatment of type-2 diabetes mellitus. Data from clinical studies indicate that TAK-875, which is an ago-allosteric modulator of hGPR40 (ref. 3), demonstrates improved glycaemic control and low hypoglycaemic risk in diabetic patients. Here we report the crystal structure of hGPR40 receptor bound to TAK-875 at 2.3 Å resolution. The co-complex structure reveals a unique binding mode of TAK-875 and suggests that entry to the non-canonical binding pocket most probably occurs via the lipid bilayer. The atomic details of the extensive charge network in the ligand binding pocket reveal additional interactions not identified in previous studies and contribute to a clear understanding of TAK-875 binding to the receptor. The hGPR40-TAK-875 structure also provides insights into the plausible binding of multiple ligands to the receptor, which has been observed in radioligand binding and Ca(2+) influx assay studies. Comparison of the transmembrane helix architecture with other G-protein-coupled receptors suggests that the crystallized TAK-875-bound hGPR40 complex is in an inactive-like state.

Discovery of novel inhibitors of HIV-1 reverse transcriptase through virtual screening of experimental and theoretical ensembles.

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are potent anti-HIV chemotherapeutics. Although there are FDA-approved NNRTIs, challenges such as the development of resistance have limited their utility. Here, we describe the identification of novel NNRTIs through a combination of computational and experimental approaches. Based on the known plasticity of the NNRTI binding pocket (NNIBP), we adopted an ensemble-based virtual screening strategy: coupling receptor conformations from 10 X-ray crystal structures with 120 snapshots from a total of 480 ns of molecular dynamics (MD) trajectories. A screening library of 2864 National Cancer Institute (NCI) compounds was built and docked against the ensembles in a hierarchical fashion. Sixteen diverse compounds were tested for their ability to block HIV infection in human tissue cultures using a luciferase-based reporter assay. Three promising compounds were further characterized, using a HIV-1 RT-based polymerase assay, to determine the specific mechanism of inhibition. We found that 2 of the three compounds inhibited the polymerase activity of RT (with potency similar to the positive control, the FDA-approved drug nevirapine). Through a computational approach, we were able to discover two compounds which inhibit HIV replication and block the activity of RT, thus offering the potential for optimization into mature inhibitors.

Simulations of biased agonists in the ß(2) adrenergic receptor with accelerated molecular dynamics.

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signaling pathway a GPCR promotes intracellular signals though ß-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signaling through the G protein and ß-arrestin. Here we report on the dynamics of the ß2 adrenergic receptor bound to the ß-arrestin and G protein-biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring of the transition within the nanosecond time scale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the ß-arrestin-biased agonist N-cyclopentylbutanepherine, we observe a different pattern of motions in helix 7 when compared to simulations with the G protein-biased agonist salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs.

Locating a plausible binding site for an open-channel blocker, GlyH-101, in the pore of the cystic fibrosis transmembrane conductance regulator.

High-throughput screening has led to the identification of small-molecule blockers of the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel, but the structural basis of blocker binding remains to be defined. We developed molecular models of the CFTR channel on the basis of homology to the bacterial transporter Sav1866, which could permit blocker binding to be analyzed in silico. The models accurately predicted the existence of a narrow region in the pore that is a likely candidate for the binding site of an open-channel pore blocker such as N-(2-naphthalenyl)-[(3,5-dibromo-2,4-dihydroxyphenyl)methylene]glycine hydrazide (GlyH-101), which is thought to act by entering the channel from the extracellular side. As a more-stringent test of predictions of the CFTR pore model, we applied induced-fit, virtual, ligand-docking techniques to identify potential binding sites for GlyH-101 within the CFTR pore. The highest-scoring docked position was near two pore-lining residues, Phe337 and Thr338, and the rates of reactions of anionic, thiol-directed reagents with cysteines substituted at these positions were slowed in the presence of the blocker, consistent with the predicted repulsive effect of the net negative charge on GlyH-101. When a bulky phenylalanine that forms part of the predicted binding pocket (Phe342) was replaced with alanine, the apparent affinity of the blocker was increased ∼200-fold. A molecular mechanics-generalized Born/surface area analysis of GlyH-101 binding predicted that substitution of Phe342 with alanine would substantially increase blocker affinity, primarily because of decreased intramolecular strain within the blocker-protein complex. This study suggests that GlyH-101 blocks the CFTR channel by binding within the pore bottleneck.

A combinatorial TIR1/AFB-Aux/IAA co-receptor system for differential sensing of auxin.

The plant hormone auxin regulates virtually every aspect of plant growth and development. Auxin acts by binding the F-box protein transport inhibitor response 1 (TIR1) and promotes the degradation of the AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) transcriptional repressors. Here we show that efficient auxin binding requires assembly of an auxin co-receptor complex consisting of TIR1 and an Aux/IAA protein. Heterologous experiments in yeast and quantitative IAA binding assays using purified proteins showed that different combinations of TIR1 and Aux/IAA proteins form co-receptor complexes with a wide range of auxin-binding affinities. Auxin affinity seems to be largely determined by the Aux/IAA. As there are 6 TIR1/AUXIN SIGNALING F-BOX proteins (AFBs) and 29 Aux/IAA proteins in Arabidopsis thaliana, combinatorial interactions may result in many co-receptors with distinct auxin-sensing properties. We also demonstrate that the AFB5-Aux/IAA co-receptor selectively binds the auxinic herbicide picloram. This co-receptor system broadens the effective concentration range of the hormone and may contribute to the complexity of auxin response.

Cystic fibrosis transmembrane conductance regulator: a molecular model defines the architecture of the anion conduction path and locates a "bottleneck" in the pore.

We developed molecular models for the cystic fibrosis transmembrane conductance regulator chloride channel based on the prokaryotic ABC transporter, Sav1866. Here we analyze predicted pore geometry and side-chain orientations for TM3, TM6, TM9, and TM12, with particular attention being paid to the location of the rate-limiting barrier for anion conduction. Side-chain orientations assayed by cysteine scanning were found to be from 77 to 90% in accord with model predictions. The predicted geometry of the anion conduction path was defined by a space-filling model of the pore and confirmed by visualizing the distribution of water molecules from a molecular dynamics simulation. The pore shape is that of an asymmetric hourglass, comprising a shallow outward-facing vestibule that tapers rapidly toward a narrow "bottleneck" linking the outer vestibule to a large inner cavity extending toward the cytoplasmic extent of the lipid bilayer. The junction between the outer vestibule and the bottleneck features an outward-facing rim marked by T338 in TM6 and I1131 in TM12, consistent with the observation that cysteines at both of these locations reacted with both channel-permeant and channel-impermeant, thiol-directed reagents. Conversely, cysteines substituted for S341 in TM6 or T1134 in TM12, predicted by the model to lie below the rim of the bottleneck, were found to react exclusively with channel-permeant reagents applied from the extracellular side. The predicted dimensions of the bottleneck are consistent with the demonstrated permeation of Cl(-), pseudohalide anions, water, and urea.

A molecular dynamics ensemble-based approach for the mapping of druggable binding sites.

An expanding repertoire of "allosteric" drugs is revealing that structure-based drug design (SBDD) is not restricted to the "active site" of the target protein. Such compounds have been shown to bind distant regions of the protein topography, potentially providing higher levels of target specificity, reduced toxicity and access to new regions of chemical space. Unfortunately, the location of such allosteric pockets is not obvious in the absence of a bound crystal structure and the ability to predict their presence would be useful in the discovery of novel therapies. Here, we describe a method for the prediction of "druggable" binding sites that takes protein flexibility into account through the use of molecular dynamics (MD) simulation. By using a dynamic representation of the target, we are able to sample multiple protein conformations that may expose new drug-binding surfaces. We perform a fragment-based mapping analysis of individual structures in the MD ensemble using the FTMAP algorithm and then rank the most prolific probe-binding protein residues to determine potential "hot-spots" for further examination. This approach has recently been applied to a pair of human G-protein-coupled receptors (GPCRs), resulting in the detection of five potential allosteric sites.

Molecular recognition in the case of flexible targets.

A protein's flexibility is well recognized to underlie its capacity to engage in critical functions, such as signal transduction, biomolecular transport and biochemical reactivity. Molecular recognition is also tightly linked to the dynamics of the binding partners, yet protein flexibility has largely been ignored by the growing field of structure-based drug design (SBDD). In combination with experimentally determined structures, a number of computational methods have been proposed to model protein movements, which may be important for small molecule binding. Such techniques have the ability to expose new binding site conformations, which may in turn recognize and lead to the discovery of more potent and selective drugs through molecular docking. In this article, we discuss various methods and focus on the Relaxed Complex Scheme (RCS), which uses Molecular Dynamics (MD) simulations to model full protein flexibility and enhance virtual screening programmes. We review practical applications of the RCS and use a recent study of the HIV-1 reverse transcriptase to illustrate the various phases of the scheme. We also discuss some encouraging developments, aimed at addressing current weaknesses of the RCS.

Predictive power of molecular dynamics receptor structures in virtual screening.

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range.