Discovery of lirafugratinib (RLY-4008), a highly selective irreversible small-molecule inhibitor of FGFR2.

Fibroblast growth factor receptor (FGFR) kinase inhibitors have been shown to be effective in the treatment of intrahepatic cholangiocarcinoma and other advanced solid tumors harboring FGFR2 alterations, but the toxicity of these drugs frequently leads to dose reduction or interruption of treatment such that maximum efficacy cannot be achieved. The most common adverse effects are hyperphosphatemia caused by FGFR1 inhibition and diarrhea due to FGFR4 inhibition, as current therapies are not selective among the FGFRs. Designing selective inhibitors has proved difficult with conventional approaches because the orthosteric sites of FGFR family members are observed to be highly similar in X-ray structures. In this study, aided by analysis of protein dynamics, we designed a selective, covalent FGFR2 inhibitor. In a key initial step, analysis of long-timescale molecular dynamics simulations of the FGFR1 and FGFR2 kinase domains allowed us to identify differential motion in their P-loops, which are located adjacent to the orthosteric site. Using this insight, we were able to design orthosteric binders that selectively and covalently engage the P-loop of FGFR2. Our drug discovery efforts culminated in the development of lirafugratinib (RLY-4008), a covalent inhibitor of FGFR2 that shows substantial selectivity over FGFR1 (~250-fold) and FGFR4 (~5,000-fold) in vitro, causes tumor regression in multiple FGFR2-altered human xenograft models, and was recently demonstrated to be efficacious in the clinic at doses that do not induce clinically significant hyperphosphatemia or diarrhea.

RLY-4008, the First Highly Selective FGFR2 Inhibitor with Activity across FGFR2 Alterations and Resistance Mutations.

Oncogenic activation of fibroblast growth factor receptor 2 (FGFR2) drives multiple cancers and represents a broad therapeutic opportunity, yet selective targeting of FGFR2 has not been achieved. Although the clinical efficacy of pan-FGFR inhibitors (pan-FGFRi) validates FGFR2 driver status in FGFR2 fusion-positive intrahepatic cholangiocarcinoma, their benefit is limited by incomplete target coverage due to FGFR1- and FGFR4-mediated toxicities (hyperphosphatemia and diarrhea, respectively) and the emergence of FGFR2 resistance mutations. RLY-4008 is a highly selective, irreversible FGFR2 inhibitor designed to overcome these limitations. In vitro, RLY-4008 demonstrates >250- and >5,000-fold selectivity over FGFR1 and FGFR4, respectively, and targets primary alterations and resistance mutations. In vivo, RLY-4008 induces regression in multiple xenograft models-including models with FGFR2 resistance mutations that drive clinical progression on current pan-FGFRi-while sparing FGFR1 and FGFR4. In early clinical testing, RLY-4008 induced responses without clinically significant off-isoform FGFR toxicities, confirming the broad therapeutic potential of selective FGFR2 targeting. SIGNIFICANCE: Patients with FGFR2-driven cancers derive limited benefit from pan-FGFRi due to multiple FGFR1-4-mediated toxicities and acquired FGFR2 resistance mutations. RLY-4008 is a highly selective FGFR2 inhibitor that targets primary alterations and resistance mutations and induces tumor regression while sparing other FGFRs, suggesting it may have broad therapeutic potential. See related commentary by Tripathi et al., p. 1964. This article is featured in Selected Articles from This Issue, p. 1949.

Structural mechanism of a drug-binding process involving a large conformational change of the protein target.

Proteins often undergo large conformational changes when binding small molecules, but atomic-level descriptions of such events have been elusive. Here, we report unguided molecular dynamics simulations of Abl kinase binding to the cancer drug imatinib. In the simulations, imatinib first selectively engages Abl kinase in its autoinhibitory conformation. Consistent with inferences drawn from previous experimental studies, imatinib then induces a large conformational change of the protein to reach a bound complex that closely resembles published crystal structures. Moreover, the simulations reveal a surprising local structural instability in the C-terminal lobe of Abl kinase during binding. The unstable region includes a number of residues that, when mutated, confer imatinib resistance by an unknown mechanism. Based on the simulations, NMR spectra, hydrogen-deuterium exchange measurements, and thermostability measurements and estimates, we suggest that these mutations confer imatinib resistance by exacerbating structural instability in the C-terminal lobe, rendering the imatinib-bound state energetically unfavorable.

Targetable HER3 functions driving tumorigenic signaling in HER2-amplified cancers.

Effective inactivation of the HER2-HER3 tumor driver has remained elusive because of the challenging attributes of the pseudokinase HER3. We report a structure-function study of constitutive HER2-HER3 signaling to identify opportunities for targeting. The allosteric activation of the HER2 kinase domain (KD) by the HER3 KD is required for tumorigenic signaling and can potentially be targeted by allosteric inhibitors. ATP binding within the catalytically inactive HER3 KD provides structural rigidity that is important for signaling, but this is mimicked, not opposed, by small molecule ATP analogs, reported here in a bosutinib-bound crystal structure. Mutational disruption of ATP binding and molecular dynamics simulation of the apo KD of HER3 identify a conformational coupling of the ATP pocket with a hydrophobic AP-2 pocket, analogous to EGFR, that is critical for tumorigenic signaling and feasible for targeting. The value of these potential target sites is confirmed in tumor growth assays using gene replacement techniques.

Structural basis for ALK2/BMPR2 receptor complex signaling through kinase domain oligomerization.

Upon ligand binding, bone morphogenetic protein (BMP) receptors form active tetrameric complexes, comprised of two type I and two type II receptors, which then transmit signals to SMAD proteins. The link between receptor tetramerization and the mechanism of kinase activation, however, has not been elucidated. Here, using hydrogen deuterium exchange mass spectrometry (HDX-MS), small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations, combined with analysis of SMAD signaling, we show that the kinase domain of the type I receptor ALK2 and type II receptor BMPR2 form a heterodimeric complex via their C-terminal lobes. Formation of this dimer is essential for ligand-induced receptor signaling and is targeted by mutations in BMPR2 in patients with pulmonary arterial hypertension (PAH). We further show that the type I/type II kinase domain heterodimer serves as the scaffold for assembly of the active tetrameric receptor complexes to enable phosphorylation of the GS domain and activation of SMADs.

Structural Basis of AZD9291 Selectivity for EGFR T790M.

AZD9291 (Osimertinib) is highly effective in treating EGFR-mutated non-small-cell lung cancers (NSCLCs) with T790M-mediated drug resistance. Despite the remarkable success of AZD9291, its binding pose with EGFR T790M remains unclear. Here, we report unbiased, atomic-level molecular dynamics (MD) simulations in which spontaneous association of AZD9291 with EGFR kinases having WT and T790M mutant gatekeepers was observed. Simulation-generated structural models suggest that the binding pose of AZD9291 with T790M differs from its binding pose with the WT, and that AZD9291 interacts extensively with the gatekeeper residue (Met 790) in T790M but not with Thr 790 in the WT, which explains why AZD9291 binds T790M with higher affinity. The MD simulation-generated models were confirmed by experimentally determined EGFR/T790M complex crystal structures. This work may facilitate the rational design of drugs that can overcome resistance mutations to AZD9291, and more generally it suggests the potential of using unbiased MD simulation to elucidate small-molecule binding poses.

Conformational Adaption May Explain the Slow Dissociation Kinetics of Roniciclib (BAY 1000394), a Type I CDK Inhibitor with Kinetic Selectivity for CDK2 and CDK9.

Roniciclib (BAY 1000394) is a type I pan-CDK (cyclin-dependent kinase) inhibitor which has revealed potent efficacy in xenograft cancer models. Here, we show that roniciclib displays prolonged residence times on CDK2 and CDK9, whereas residence times on other CDKs are transient, thus giving rise to a kinetic selectivity of roniciclib. Surprisingly, variation of the substituent at the 5-position of the pyrimidine scaffold results in changes of up to 3 orders of magnitude of the drug-target residence time. CDK2 X-ray cocrystal structures have revealed a DFG-loop adaption for the 5-(trifluoromethyl) substituent, while for hydrogen and bromo substituents the DFG loop remains in its characteristic type I inhibitor position. In tumor cells, the prolonged residence times of roniciclib on CDK2 and CDK9 are reflected in a sustained inhibitory effect on retinoblastoma protein (RB) phosphorylation, indicating that the target residence time on CDK2 may contribute to sustained target engagement and antitumor efficacy.

A universal homogeneous assay for high-throughput determination of binding kinetics.

There is an increasing demand for assay technologies that enable accurate, cost-effective, and high-throughput measurements of drug-target association and dissociation rates. Here we introduce a universal homogeneous kinetic probe competition assay (kPCA) that meets these requirements. The time-resolved fluorescence energy transfer (TR-FRET) procedure combines the versatility of radioligand binding assays with the advantages of homogeneous nonradioactive techniques while approaching the time resolution of surface plasmon resonance (SPR) and related biosensors. We show application of kPCA for three important target classes: enzymes, protein-protein interactions, and G protein-coupled receptors (GPCRs). This method is capable of supporting early stages of drug discovery with large amounts of kinetic information.

A tethered delivery mechanism explains the catalytic action of a microtubule polymerase.

Stu2p/XMAP215 proteins are essential microtubule polymerases that use multiple αß-tubulin-interacting TOG domains to bind microtubule plus ends and catalyze fast microtubule growth. We report here the structure of the TOG2 domain from Stu2p bound to yeast αß-tubulin. Like TOG1, TOG2 binds selectively to a fully 'curved' conformation of αß-tubulin, incompatible with a microtubule lattice. We also show that TOG1-TOG2 binds non-cooperatively to two αß-tubulins. Preferential interactions between TOGs and fully curved αß-tubulin that cannot exist elsewhere in the microtubule explain how these polymerases localize to the extreme microtubule end. We propose that these polymerases promote elongation because their linked TOG domains concentrate unpolymerized αß-tubulin near curved subunits already bound at the microtubule end. This tethering model can explain catalyst-like behavior and also predicts that the polymerase action changes the configuration of the microtubule end.

A TOG:αß-tubulin complex structure reveals conformation-based mechanisms for a microtubule polymerase.

Stu2p/XMAP215/Dis1 family proteins are evolutionarily conserved regulatory factors that use αß-tubulin-interacting tumor overexpressed gene (TOG) domains to catalyze fast microtubule growth. Catalysis requires that these polymerases discriminate between unpolymerized and polymerized forms of αß-tubulin, but the mechanism by which they do so has remained unclear. Here, we report the structure of the TOG1 domain from Stu2p bound to yeast αß-tubulin. TOG1 binds αß-tubulin in a way that excludes equivalent binding of a second TOG domain. Furthermore, TOG1 preferentially binds a curved conformation of αß-tubulin that cannot be incorporated into microtubules, contacting α- and ß-tubulin surfaces that do not participate in microtubule assembly. Conformation-selective interactions with αß-tubulin explain how TOG-containing polymerases discriminate between unpolymerized and polymerized forms of αß-tubulin and how they selectively recognize the growing end of the microtubule.

Design, overexpression, and purification of polymerization-blocked yeast αß-tubulin mutants.

Microtubule dynamics play essential roles in intracellular organization and cell division. They result from structural and biochemical properties of αß-tubulin heterodimers and how these polymerizing subunits interact with themselves and with regulatory proteins. A broad understanding of the underlying mechanisms has been established, but fundamental questions remain unresolved. The lack of routine access to recombinant αß-tubulin represents an obstacle to deeper insight into αß-tubulin structure, biochemistry, and recognition. Indeed, the widespread reliance on animal brain αß-tubulin means that very few in vitro studies have taken advantage of powerful and ordinarily routine techniques like site-directed mutagenesis. Here we report new methods for purifying wild-type or mutant yeast αß-tubulin from inducibly overexpressing strains of Saccharomyces cerevisiae. Inducible overexpression is an improvement over existing approaches that rely on constitutive expression: it provides higher yields while also allowing otherwise lethal mutants to be purified. We also designed and purified polymerization-blocked αß-tubulin mutants. These "blocked" forms of αß-tubulin give a dominant lethal phenotype when expressed in cells; they cannot form microtubules in vitro and when present in mixtures inhibit the polymerization of wild-type αß-tubulin. The effects of blocking mutations are very specific, because purified mutants exhibit normal hydrodynamic properties, bind GTP, and interact with a tubulin-binding domain. The ability to overexpress and purify wild-type αß-tubulin, or mutants like the ones we report here, creates new opportunities for structural studies of αß-tubulin and its complexes with regulatory proteins, and for biochemical and functional studies of microtubule dynamics and its regulation.