A combined computational and experimental strategy identifies mutations conferring resistance to drugs targeting the BCR-ABL fusion protein

Predicting Mutational Effects on Ca2+-Activated Chloride Conduction of TMEM16A Based on a Simulation Study

The dysfunction and defects of ion channels are associated with many human diseases, especially for loss-of-function mutations in ion channels such as cystic fibrosis transmembrane conductance regulator mutations in cystic fibrosis. Understanding ion channels is of great current importance for both medical and fundamental purposes. Such an understanding should include the ability to predict mutational effects and describe functional and mechanistic effects. In this work, we introduce an approach to predict mutational effects based on kinetic information (including reaction barriers and transition state locations) obtained by studying the working mechanism of target proteins. Specifically, we take the Ca2+-activated chloride channel TMEM16A as an example and utilize the computational biology model to predict the mutational effects of key residues. Encouragingly, we verified our predictions through electrophysiological experiments, demonstrating a 94% prediction accuracy regarding mutational directions. The mutational strength assessed by Pearson's correlation coefficient is -0.80 between our calculations and the experimental results. These findings suggest that the proposed methodology is reliable and can provide valuable guidance for revealing functional mechanisms and identifying key residues of the TMEM16A channel. The proposed approach can be extended to a broad scope of biophysical systems.

Nanodrug regulates ROS homeostasis via enhancing fatty acid oxidation and inhibiting autophagy to overcome tumor drug resistance

The treatment of drug-resistant tumors poses a significant challenge in the field of tumor therapy. Disrupting the homeostasis of reactive oxygen species (ROS) within tumor cells may represent a pivotal strategy for overcoming the prevalent issue of drug resistance. However, the restricted sustainability of ROS generation and the increased autophagy capacity exhibited by tumor cells hinder the application of ROS-based therapies. In this study, we developed liposome nanoparticles (Ato/CQ@L) for co-encapsulation of atorvastatin (Ato), an activator of AMP-activated protein kinase (AMPK), and chloroquine (CQ), an autophagy inhibitor. Upon internalization by tumor cells, Ato upregulated carnitine palmitoyltransferase 1(CPT1) concentration and promoted fatty acid oxidation (FAO) within the tumor cells. The process of FAO coupled with an abundance of fatty acid substrates, facilitates a sustained generation of ROS production. Concurrently, a positive feedback loop is established between escalated concentration of ROS and AMPK protein levels, resulting in a persistent elevation in ROS levels. In addition, CQ disrupted lysosomes, leading to an increased lysosomal pH and reducing autophagy in tumor cells. In both in vivo and in vitro experiments, the Ato/CQ@L treatment group exhibited a considerable enhancement in tumor cell apoptosis, validating the efficacy of this combined therapy. In summary, the combined therapy involving Ato and CQ addresses the inherent limitations of conventional ROS therapy, which include insufficient ROS production and increased autophagy. This approach holds significant potential as a treatment strategy for drug-resistant triple-negative breast cancer.

In silico resources help combat cancer drug resistance mediated by target mutations

Drug resistance causes catastrophic cancer treatment failures. Mutations in target proteins with altered drug binding indicate a main mechanism of cancer drug resistance (CDR). Global research has generated considerable CDR-related data and well-established knowledge bases and predictive tools. Unfortunately, these resources are fragmented and underutilized. Here, we examine computational resources for exploring CDR caused by target mutations, analyzing these tools based on their functional characteristics, data capacity, data sources, methodologies and performance. We also discuss their disadvantages and provide examples of how potential inhibitors of CDR have been discovered using these resources. This toolkit is designed to help specialists explore resistance occurrence effectively and to explain resistance prediction to non-specialists easily.

Prediction of Amino Acid Substitutions in ABL1 Protein Leading to Tumor Drug Resistance Based on "Structure-Property" Relationship Classification Models

Drug resistance to anticancer drugs is a serious complication in patients with cancer. Typically, drug resistance occurs due to amino acid substitutions (AAS) in drug target proteins. The study aimed at developing and validating a new approach to the creation of structure-property relationships (SPR) classification models to predict AASs leading to drug resistance to inhibitors of tyrosine-protein kinase ABL1. The approach was based on the representation of AASs as peptides described in terms of structural formulas. The data on drug-resistant and non-resistant variants of AAS for two isoforms of ABL1 were extracted from the COSMIC database. The given training sets (approximately 700 missense variants) were used for the creation of SPR models in MultiPASS software based on substructural atom-centric multiple neighborhoods of atom (MNA) descriptors for the description of the structural formula of protein fragments and a Bayesian-like algorithm for revealing structure-property relationships. It was found that MNA descriptors of the 6th level and peptides from 11 amino acid residues were the best combination for ABL1 isoform 1 with the prediction accuracy (AUC) of resistance to imatinib (0.897) and dasatinib (0.996). For ABL1 isoform 2 (resistance to imatinib), the best combination was MNA descriptors of the 6th level, peptides form 15 amino acids (AUC value was 0.909). The prediction of possible drug-resistant AASs was made for dbSNP and gnomAD data. The six selected most probable imatinib-resistant AASs were additionally validated by molecular modeling and docking, which confirmed the possibility of resistance for the E334V and T392I variants.

A comprehensive protein design protocol to identify resistance mutations and signatures of adaptation in pathogens

Most pathogens mutate and evolve over time to escape immune and drug pressure. To achieve this, they alter specific hotspot residues in their intracellular proteins to render the targeted drug(s) ineffective and develop resistance. Such hotspot residues may be located as a cluster or uniformly as a signature of adaptation in a protein. Identifying the hotspots and signatures is extremely important to comprehensively understand the disease pathogenesis and rapidly develop next-generation therapeutics. As experimental methods are time-consuming and often cumbersome, there is a need to develop efficient computational protocols and adequately utilize them. To address this issue, we present a unique computational protein design protocol that identifies hotspot residues, resistance mutations and signatures of adaptation in a pathogen’s protein against a bound drug. Using the protocol, the binding affinity between the designed mutants and drug is computed quickly, which offers predictions for comparison with biophysical experiments. The applicability and accuracy of the protocol are shown using case studies of a few protein–drug complexes. As a validation, resistance mutations in severe acute respiratory syndrome coronavirus 2 main protease (Mpro) against narlaprevir (an inhibitor of hepatitis C NS3/4A serine protease) are identified. Notably, a detailed methodology and description of the working principles of the protocol are presented. In conclusion, our protocol will assist in providing a first-hand explanation of adaptation, hotspot-residue variations and surveillance of evolving resistance mutations in a pathogenic protein.

Structure-guided machine learning prediction of drug resistance mutations in Abelson 1 kinase

Kinases play crucial roles in cellular signalling and biological processes with their dysregulation associated with diseases, including cancers. Kinase inhibitors, most notably those targeting ABeLson 1 (ABL1) kinase in chronic myeloid leukemia, have had a significant impact on cancer survival, yet emergence of resistance mutations can reduce their effectiveness, leading to therapeutic failure. Limited effort, however, has been devoted to developing tools to accurately identify ABL1 resistance mutations, as well as providing insights into their molecular mechanisms. Here we investigated the structural basis of ABL1 mutations modulating binding affinity of eight FDA-approved drugs. We found mutations impair affinity of type I and type II inhibitors differently and used this insight to developed a novel web-based diagnostic tool, SUSPECT-ABL, to pre-emptively predict resistance profiles and binding free-energy changes (ΔΔG) of all possible ABL1 mutations against inhibitors with different binding modes. Resistance mutations in ABL1 were successfully identified, achieving a Matthew's Correlation Coefficient of up to 0.73 and the resulting change in ligand binding affinity with a Pearson's correlation of up to 0.77, with performances consistent across non-redundant blind tests. Through an in silico saturation mutagenesis, our tool has identified possibly emerging resistance mutations, which offers opportunities for in vivo experimental validation. We believe SUSPECT-ABL will be an important tool not just for improving precision medicine efforts, but for facilitating the development of next-generation inhibitors that are less prone to resistance. We have made our tool freely available at http://biosig.unimelb.edu.au/suspect_abl/.

Acquired rare recurrent EGFR mutations as mechanisms of resistance to Osimertinib in lung cancer and in silico structural modelling

A growing number of progression on Osimertinib among EGFR-mutated lung cancers represents a great challenge clinically. Our study aims to gain insights into novel mechanisms of acquired resistance to Osimertinib. We performed genomic studies on 2 large independent cohorts of lung cancer patients with progressed diseases on different tyrosine kinase inhibitors (TKIs). In silico modeling was used to study the structural mechanism of selected EGFR mutations. Compared with the 1^st-TKIs-resistant group, EGFR mutations C797S/G, L718Q/V, L792F/H were significantly more enriched in the Osimertinib-resistant cohort, whose sensitivities to Osimertinib were successfully predicted. Importantly, a total of 14 low-frequency EGFR mutations were exclusively or significantly observed in the Osimertinib-resistant group, 7 were predicted to dramatically reduce the binding affinity of EGFR to Osimertinib (G796S, V802F, T725M, Q791L/H, P794S/R). Analysis of pre-Osimertinib treatment samples of two patients supported that EGFR V802F and G796S were acquired during the treatment. In addition, EGFR G796S was predicted to be susceptible to gefitinib. This study represented the largest real-world data so far investigating Osimertinib resistance in EGFR-mutated lung cancer. We identified a collection of coexistent EGFR rare mutations and provided possible guidance for those patients who progressed on the first-line treatment of Osimertinib.