Structure and energy based quantitative missense variant effect analysis provides insights into drug resistance mechanisms of anaplastic lymphoma kinase mutations

A computational examination of the therapeutic advantages of fourth-generation ALK inhibitors TPX-0131 and repotrectinib over third-generation lorlatinib for NSCLC with ALK F1174C/L/V mutations

Background: In non-small-cell lung cancer (NSCLC), a pivotal factor in promoting cancer development is the rearrangement in the anaplastic lymphoma kinase ALK gene, resulting in elevated ALK protein expression. F1174C/L/V is the acquired secondary resistant mutation in ALK. Significant survival improvements have been seen while tyrosine kinase inhibitors specifically target ALK. Nevertheless, the emergence of drug resistance hinders the clinical effectiveness of these drugs. Objective: This research sought to find the binding affinity/inhibitory effects of the existing drug lorlatinib (LOR) and upcoming TPX-0131 (zotizalkib/TPX) and repotrectinib (TPX-0005/REP) inhibitors against ALK F1174C/L/V mutations using computational approaches to identify potential strategies over resistance. Methods: We conducted molecular docking, molecular dynamics simulation, and MMPBSA calculations to investigate how compact macrocyclic inhibitors, such as TPX-0131 and repotrectinib, fit within the ATP-binding boundary and differ from LOR. Results: Our results demonstrated that TPX-0131 and repotrectinib contributed to higher binding energy in F1174C and F1174L mutations than LOR. Repotrectinib showed greater binding energy in the F1174V mutation, whereas LOR and TPX-0131 exhibited similar binding energy. However, all three inhibitors showed significant binding energy toward F1174C/L/V mutations found in NSCLC. Conclusion: This comparative study of the potential binding effects of fourth-generation inhibitors TPX-0131 and repotrectinib and third-generation inhibitor LOR for ALK F1174C/L/V mutations revealed the atomistic insights of the binding mechanism. These computational findings enable us to carry out further research for the clinical implementation of fourth-generation ALK inhibitors on ALK-positive NSCLC.

Polyphyllin VII induces autophagy-dependent ferroptosis in human gastric cancer through targeting T-lymphokine-activated killer cell-originated protein kinase

T-lymphokine-activated killer cell-originated protein kinase (TOPK) is a serine-threonine kinase that is overexpressed in gastric cancer (GC) and promotes tumor progression. Polyphyllin VII (PPVII), a pennogenin isolated from the rhizomes of Paris polyphylla, shows anticancer effects. Here, we explored the antitumor activity and mechanism of PPVII in GC. Ferroptosis was detected by transmission electron microscope, malondialdehyde, and iron determination assays. Autophagy and its upstream signaling pathway were detected by Western blot, and gene alterations. The binding of PPVII and TOPK was examined through microscale thermophoresis and drug affinity responsive target stability assays. An in vivo mouse model was performed to evaluate the therapeutic of PPVII. PPVII inhibits GC by inducing autophagy-mediated ferroptosis. PPVII promotes the degradation of ferritin heavy chain 1, which is responsible for autophagy-mediated ferroptosis. PPVII activates the Unc-51-like autophagy-activating kinase 1 (ULK1) upstream of autophagy. PPVII inhibits the activity of TOPK, thereby weakening the inhibition of downstream ULK1. PPVII stabilizes the dimer of the inactive form of TOPK by direct binding. PPVII inhibits tumor growth without causing obvious toxicity in vivo. Collectively, this study suggests that PPVII is a potential agent for the treatment of GC by targeting TOPK to activate autophagy-mediated ferroptosis.

Unraveling Extremely Damaging IRAK4 Variants and Their Potential Implications for IRAK4 Inhibitor Efficacy

Interleukin-1-receptor-associated kinase 4 (IRAK4) possesses a crucial function in the toll-like receptor (TLR) signaling pathway, and the dysfunction of this molecule could lead to various infectious and immune-related diseases in addition to cancers. IRAK4 genetic variants have been linked to various types of diseases. Therefore, we conducted a comprehensive analysis to recognize the missense variants with the most damaging impacts on IRAK4 with the employment of diverse bioinformatics tools to study single-nucleotide polymorphisms' effects on function, stability, secondary structures, and 3D structure. The residues' location on the protein domain and their conservation status were investigated as well. Moreover, docking tools along with structural biology were engaged in analyzing the SNPs' effects on one of the developed IRAK4 inhibitors. By analyzing IRAK4 gene SNPs, the analysis distinguished ten variants as the most detrimental missense variants. All variants were situated in highly conserved positions on an important protein domain. L318S and L318F mutations were linked to changes in IRAK4 secondary structures. Eight SNPs were revealed to have a decreasing effect on the stability of IRAK4 via both I-Mutant 2.0 and Mu-Pro tools, while Mu-Pro tool identified a decreasing effect for the G198E SNP. In addition, detrimental effects on the 3D structure of IRAK4 were also discovered for the selected variants. Molecular modeling studies highlighted the detrimental impact of these identified SNP mutant residues on the druggability of the IRAK4 ATP-binding site towards the known target inhibitor, HG-12-6, as compared to the native protein. The loss of important ligand residue-wise contacts, altered protein global flexibility, increased steric clashes, and even electronic penalties at the ligand-binding site interfaces were all suggested to be associated with SNP models for hampering the HG-12-6 affinity towards IRAK4 target protein. This given model lays the foundation for the better prediction of various disorders relevant to IRAK4 malfunction and sheds light on the impact of deleterious IRAK4 variants on IRAK4 inhibitor efficacy.

Targeting ALK Rearrangements in NSCLC: Current State of the Art

Anaplastic lymphoma kinase (ALK) alterations in non-small cell lung cancer (NSCLC) can be effectively treated with a variety of ALK-targeted drugs. After the approval of the first-generation ALK inhibitor crizotinib which achieved better results in prolonging the progression-free survival (PFS) compared with chemotherapy, a number of next-generation ALK inhibitors have been developed including ceritinib, alectinib, brigatinib, and ensartinib. Recently, a potent, third-generation ALK inhibitor, lorlatinib, has been approved by the Food and Drug Administration (FDA) for the first-line treatment of ALK-positive (ALK+) NSCLC. These drugs have manageable toxicity profiles. Responses to ALK inhibitors are however often not durable, and acquired resistance can occur as on-target or off-target alterations. Studies are underway to explore the mechanisms of resistance and optimal treatment options beyond progression. Efforts have also been undertaken to develop further generations of ALK inhibitors. This review will summarize the current situation of targeting the ALK signaling pathway.

Deciphering the Mechanism of Gilteritinib Overcoming Lorlatinib Resistance to the Double Mutant I1171N/F1174I in Anaplastic Lymphoma Kinase

Anaplastic lymphoma kinase (ALK) is validated as a therapeutic molecular target in multiple malignancies, such as non-small cell lung cancer (NSCLC). However, the feasibility of targeted therapies exerted by ALK inhibitors is inevitably hindered owing to drug resistance. The emergence of clinically acquired drug mutations has become a major challenge to targeted therapies and personalized medicines. Thus, elucidating the mechanism of resistance to ALK inhibitors is helpful for providing new therapeutic strategies for the design of next-generation drug. Here, we used molecular docking and multiple molecular dynamics simulations combined with correlated and energetical analyses to explore the mechanism of how gilteritinib overcomes lorlatinib resistance to the double mutant ALK I1171N/F1174I. We found that the conformational dynamics of the ALK kinase domain was reduced by the double mutations I1171N/F1174I. Moreover, energetical and structural analyses implied that the double mutations largely disturbed the conserved hydrogen bonding interactions from the hinge residues Glu1197 and Met1199 in the lorlatinib-bound state, whereas they had no discernible adverse impact on the binding affinity and stability of gilteritinib-bound state. These discrepancies created the capacity of the double mutant ALK I1171N/F1174I to confer drug resistance to lorlatinib. Our result anticipates to provide a mechanistic insight into the mechanism of drug resistance induced by ALK I1171N/F1174I that are resistant to lorlatinib treatment in NSCLC.

NGS-based liquid biopsy profiling identifies mechanisms of resistance to ALK inhibitors: a step toward personalized NSCLC treatment

Despite impressive and durable responses, nonsmall cell lung cancer (NSCLC) patients treated with anaplastic lymphoma kinase (ALK) inhibitors (ALK-Is) ultimately progress due to development of resistance. Here, we have evaluated the clinical utility of circulating tumor DNA (ctDNA) profiling by next-generation sequencing (NGS) upon disease progression. We collected 26 plasma and two cerebrospinal fluid samples from 24 advanced ALK-positive NSCLC patients at disease progression to an ALK-I. These samples were analyzed by NGS and digital PCR. A tool to retrieve variants at the ALK locus was developed (VALK tool). We identified at least one resistance mutation in the ALK locus in ten (38.5%) plasma samples; the G1269A and G1202R mutations were the most prevalent among patients progressing to first- and second-generation ALK-Is, respectively. Overall, 61 somatic mutations were detected in 14 genes: TP53, ALK, PIK3CA, SMAD4, MAP2K1 (MEK1), FGFR2, FGFR3, BRAF, EGFR, IDH2, MYC, MET, CCND3, and CCND1. Specifically, a deletion in exon 19 in EGFR, a non-V600 BRAF mutation (G466V), and the F129L mutation in MAP2K1 were identified in four patients who showed no objective survival benefit from ALK-Is. Potential ALK-I-resistance mutations were also found in PIK3CA and IDH2. Finally, a c-MYC gain, along with a loss of CCND1 and FGFR3, was detected in a patient progressing on a first-line treatment with crizotinib. We conclude that NGS analysis of liquid biopsies upon disease progression identified different putative ALK-I-resistance mutations in most cases and could be a valuable approach for therapy decision making.

Head to head evaluation of second generation ALK inhibitors brigatinib and alectinib as first-line treatment for ALK+ NSCLC using an in silico systems biology-based approach

Around 3-7% of patients with non-small cell lung cancer (NSCLC), which represent 85% of diagnosed lung cancers, have a rearrangement in the ALK gene that produces an abnormal activity of the ALK protein cell signaling pathway. The developed ALK tyrosine kinase inhibitors (TKIs), such as crizotinib, ceritinib, alectinib, brigatinib and lorlatinb present good performance treating ALK+ NSCLC, although all patients invariably develop resistance due to ALK secondary mutations or bypass mechanisms. In the present study, we compare the potential differences between brigatinib and alectinib's mechanisms of action as first-line treatment for ALK+ NSCLC in a systems biology-based in silico setting. Therapeutic performance mapping system (TPMS) technology was used to characterize the mechanisms of action of brigatinib and alectinib and the impact of potential resistances and drug interferences with concomitant treatments. The analyses indicate that brigatinib and alectinib affect cell growth, apoptosis and immune evasion through ALK inhibition. However, brigatinib seems to achieve a more diverse downstream effect due to a broader cancer-related kinase target spectrum. Brigatinib also shows a robust effect over invasiveness and central nervous system metastasis-related mechanisms, whereas alectinib seems to have a greater impact on the immune evasion mechanism. Based on this in silico head to head study, we conclude that brigatinib shows a predicted efficacy similar to alectinib and could be a good candidate in a first-line setting against ALK+ NSCLC. Future investigation involving clinical studies will be needed to confirm these findings. These in silico systems biology-based models could be applied for exploring other unanswered questions.

Next-generation sequencing to dynamically detect mechanisms of resistance to ALK inhibitors in ALK-positive NSCLC patients: a case report

Tyrosine kinase inhibitors (TKIs) of the anaplastic lymphoma kinase gene (ALK) have significantly improved the quality of life and survival of non-small cell lung cancer (NSCLC) patients whose tumors harbor an ALK translocation. However, most of these patients relapse within 2 to 3 years as the tumor acquires resistance mutations. Unlike beaming and digital PCR (dPCR), which only allow a few mutations to be analyzed, next-generation sequencing (NGS) approaches enable the simultaneous screening of multiple genetic alterations even when the frequencies of the variants are very low. We present the case of a 52-year-old man who was diagnosed with an ALK-positive NSCLC and was treated with crizotinib and, subsequently, ceritinib. The analysis of serial liquid biopsies by NGS detected two asynchronous mutations arising in the ALK locus during disease progression, namely p.Gly1269Ala (c.3806G>C) and p.Gly1202Arg (c.3604G>A), that conferred resistance to crizotinib and ceritinib, respectively. The resistance mutations were detected independently at different times, and could be imputed to different metastatic lesions, thereby highlighting the importance of heterogeneity in advance disease. Plasma levels of ALK resistance mutations correlated well with tumor responses assessed by CT scans and bone scintigraphy, demonstrating that non-invasive tumor molecular profiling by NGS allows the efficient dynamic monitoring of ALK-positive NSCLC patients, and outperforms dPCR and beaming because more somatic mutations can be tracked over the course of the treatment. In conclusion, this case report illustrates the usefulness NGS to guide therapeutic decisions in ALK-positive NSCLC patients based tumor molecular profile upon disease progression.

ZX-29, a novel ALK inhibitor, induces apoptosis via ER stress in ALK rearrangement NSCLC cells and overcomes cell resistance caused by an ALK mutation

Although anaplastic lymphoma kinase (ALK) inhibitors have good clinical efficacy, the inevitable development of drug resistance is the most common obstacle to their clinical application. There is an urgent need to develop more effective and selective ALK inhibitors to overcome the problem of drug resistance. Here, we screened a series of ALK inhibitors and found that ZX-29 displayed potent cytotoxic activity against ALK rearrangement non-small cell lung cancer (NSCLC) NCI-H2228 cells. Then, we investigated the antitumor effects of ZX-29. We demonstrated that ZX-29 time- and dose-dependently inhibited the viability of NCI-H2228 cells, induced cell cycle arrest in the G1 phase, and then they subsequently progressed into cell death. The type of cell death induced by ZX-29 was apoptosis through endoplasmic reticulum (ER) stress. Interestingly, ZX-29 induced protective autophagy, and inhibiting autophagy could enhance the antitumor effect of ZX-29. Furthermore, ZX-29 suppressed tumor growth in a mouse xenograft model. More importantly, ZX-29 could overcome the drug resistance caused by the ALK G1202R mutation. In conclusion, we demonstrated that ZX-29 showed excellent anti-ALK rearrangement NSCLC activity in vitro and in vivo and overcame the drug resistance caused by an ALK mutation. Therefore, ZX-29 is a promising antitumor drug targeting ALK rearrangement or ALK G1202R mutation NSCLC.

Insight into the structural stability of wild-type and histidine mutants in Pin1 by experimental and computational methods

Pin1, a polypeptide proline isomerase parvulin, plays a key role in Alzheimer's disease (AD), common tumors and cancers. Two conservative histidine residues, His59 and His157, are important for maintaining the stability of the PPIase domain. Hence multiple spectral and computational techniques were performed to investigate the potential mechanism of two histidine residues. Thermal denaturation indicated that both residues His59 and His157 are not sensitive to the lower temperatures, while residue His59 is more sensitive to the higher temperatures than residue His157. Acidic denaturation suggested that influences of both residues His59 and His157 to acidic stability were the difference from Pin1-WT. ANS and RLS spectra hinted that there was no significant effect on hydrophobic change and aggregation by histidine mutations. The GndHCl-induced denaturation implied that residues His59 and His157 contributed the most to the chemical stability. MD simulations revealed that residues His59 and His157 mutations resulted in that the hydrogen bond network of the dual histidine motif was destroyed wholly. In summary, these histidine residues play an important role in maintaining the structural stability of the PPIase domain.

SNP2SIM: a modular workflow for standardizing molecular simulation and functional analysis of protein variants

Background

Molecular simulations are used to provide insight into protein structure and dynamics, and have the potential to provide important context when predicting the impact of sequence variation on protein function. In addition to understanding molecular mechanisms and interactions on the atomic scale, translational applications of those approaches include drug screening, development of novel molecular therapies, and targeted treatment planning. Supporting the continued development of these applications, we have developed the SNP2SIM workflow that generates reproducible molecular dynamics and molecular docking simulations for downstream functional variant analysis. The Python workflow utilizes molecular dynamics software (NAMD (Phillips et al., J Comput Chem 26(16):1781-802, 2005), VMD (Humphrey et al., J Mol Graph 14(1):33-8, 27-8, 1996)) to generate variant specific scaffolds for simulated small molecule docking (AutoDock Vina (Trott and Olson, J Comput Chem 31(2):455-61, 2010)).

Results

SNP2SIM is composed of three independent modules that can be used sequentially to generate the variant scaffolds of missense protein variants from the wildtype protein structure. The workflow first generates the mutant structure and configuration files required to execute molecular dynamics simulations of solvated protein variant structures. The resulting trajectories are clustered based on the structural diversity of residues involved in ligand binding to produce one or more variant scaffolds of the protein structure. Finally, these unique structural conformations are bound to small molecule ligand libraries to predict variant induced changes to drug binding relative to the wildtype protein structure.

Conclusions

SNP2SIM provides a platform to apply molecular simulation based functional analysis of sequence variation in the protein targets of small molecule therapies. In addition to simplifying the simulation of variant specific drug interactions, the workflow enables large scale computational mutagenesis by controlling the parameterization of molecular simulations across multiple users or distributed computing infrastructures. This enables the parallelization of the computationally intensive molecular simulations to be aggregated for downstream functional analysis, and facilitates comparing various simulation options, such as the specific residues used to define structural variant clusters. The Python scripts that implement the SNP2SIM workflow are available (SNP2SIM Repository. https://github.com/mccoymd/SNP2SIM, Accessed 2019 February ), and individual SNP2SIM modules are available as apps on the Seven Bridges Cancer Genomics Cloud (Lau et al., Cancer Res 77(21):e3-e6, 2017; Cancer Genomics Cloud [www.cancergenomicscloud.org; Accessed 2018 November]).