Novel Substrate Prediction for the TAM Family of RTKs Using Phosphoproteomics and Structure-Based Modeling.

The TAM family of receptor tyrosine kinases is implicated in multiple distinct oncogenic signaling pathways. However, to date, there are no FDA-approved small molecule inhibitors for the TAM kinases. Inhibitor design and screening rely on tools to study the kinase activity. Our goal was to address this gap by designing a set of synthetic peptide substrates for each of the TAM family members: Tyro3, Axl, and Mer. We used an in vitro phosphoproteomics workflow to determine the substrate profile of each TAM kinase and input the identified substrates into our data processing pipeline, KINATEST-ID, producing a position-specific scoring matrix for each target kinase and generating a list of candidate synthetic peptide substrates. We synthesized and characterized a set of those substrate candidates, systematically measuring their initial phosphorylation rate with each TAM kinase by LC-MS. We also used the multimer modeling function of AlphaFold2 (AF2) to predict peptide-kinase interactions at the active site for each of the novel candidate peptide sequences against each of the TAM family kinases and observed that, remarkably, every sequence for which it predicted a putative catalytically competent interaction was also demonstrated biochemically to be a substrate for one or more of the TAM kinases. This work shows that kinase substrate design can be achieved using a combination of preference motifs and structural modeling, and it provides the first demonstration of peptide-protein interaction modeling with AF2 for predicting the likelihood of constructive catalytic interactions.

Source, co-occurrence, and prognostic value of PTEN mutations or loss in colorectal cancer.

Somatic PTEN mutations are common and have driver function in some cancer types. However, in colorectal cancers (CRCs), somatic PTEN-inactivating mutations occur at a low frequency (~8-9%), and whether these mutations are actively selected and promote tumor aggressiveness has been controversial. Analysis of genomic data from ~53,000 CRCs indicates that hotspot mutation patterns in PTEN partially reflect DNA-dependent selection pressures, but also suggests a strong selection pressure based on protein function. In microsatellite stable (MSS) tumors, PTEN alterations co-occur with mutations activating BRAF or PI3K, or with TP53 deletions, but not in CRC with microsatellite instability (MSI). Unexpectedly, PTEN deletions are associated with poor survival in MSS CRC, whereas PTEN mutations are associated with improved survival in MSI CRC. These and other data suggest use of PTEN as a prognostic marker is valid in CRC, but such use must consider driver mutation landscape, tumor subtype, and category of PTEN alteration.

Efficient Hit-to-Lead Searching of Kinase Inhibitor Chemical Space via Computational Fragment Merging.

In early-stage drug discovery, the hit-to-lead optimization (or "hit expansion") stage entails starting from a newly identified active compound and improving its potency or other properties. Traditionally, this process relies on synthesizing and evaluating a series of analogues to build up structure-activity relationships. Here, we describe a computational strategy focused on kinase inhibitors, intended to expedite the process of identifying analogues with improved potency. Our protocol begins from an inhibitor of the target kinase and generalizes the synthetic route used to access it. By searching for commercially available replacements for the individual building blocks used to make the parent inhibitor, we compile an enumerated library of compounds that can be accessed using the same chemical transformations; these huge libraries can exceed many millions─or billions─of compounds. Because the resulting libraries are much too large for explicit virtual screening, we instead consider alternate approaches to identify the top-scoring compounds. We find that contributions from individual substituents are well described by a pairwise additivity approximation, provided that the corresponding fragments position their shared core in precisely the same way relative to the binding site. This key insight allows us to determine which fragments are suitable for merging into single new compounds and which are not. Further, the use of pairwise approximation allows interaction energies to be assigned to each compound in the library without the need for any further structure-based modeling: interaction energies instead can be reliably estimated from the energies of the component fragments, and the reduced computational requirements allow for flexible energy minimizations that allow the kinase to respond to each substitution. We demonstrate this protocol using libraries built from six representative kinase inhibitors drawn from the literature, which target five different kinases: CDK9, CHK1, CDK2, EGFRT790M, and ACK1. In each example, the enumerated library includes additional analogues reported by the original study to have activity, and these analogues are successfully prioritized within the library. We envision that the insights from this work can facilitate the rapid assembly and screening of increasingly large libraries for focused hit-to-lead optimization. To enable adoption of these methods and to encourage further analyses, we disseminate the computational tools needed to deploy this protocol.

Rationalizing PROTAC-Mediated Ternary Complex Formation Using Rosetta.

Proteolysis-targeting chimaeras (PROTACs) are molecules that combine a target-binding warhead with an E3 ligase-recruiting moiety; by drawing the target protein into a ternary complex with the E3 ligase, PROTACs induce target protein degradation. While PROTACs hold exciting potential as chemical probes and as therapeutic agents, development of a PROTAC typically requires synthesis of numerous analogs to thoroughly explore variations on the chemical linker; without extensive trial and error, it is unclear how to link the two protein-recruiting moieties to promote formation of a productive ternary complex. Here, we describe a structure-based computational method for evaluating the suitability of a given linker for ternary complex formation. Our method uses Rosetta to dock the protein components and then builds the PROTAC from its component fragments into each binding mode; complete models of the ternary complex are then refined. We apply this approach to retrospectively evaluate multiple PROTACs from the literature, spanning diverse target proteins. We find that modeling ternary complex formation is sufficient to explain both activity and selectivity reported for these PROTACs, implying that other cellular factors are not key determinants of activity in these cases. We further find that interpreting PROTAC activity is best approached using an ensemble of structures of the ternary complex rather than a single static conformation and that members of a structurally conserved protein family can be recruited by the same PROTAC through vastly different binding modes. To encourage adoption of these methods and promote further analyses, we disseminate both the computational methods and the models of ternary complexes.