Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Gaining insight into crizotinib resistance mechanisms caused by L2026M and G2032R mutations in ROS1 via molecular dynamics simulations and free-energy calculations

ROS1 fusion kinase—highly expressed in a variety of human cancers—has emerged as an important and attractive target for anticancer drug design. Crizotinib, a well-known drug approved by the FDA as an ALK inhibitor to treat advanced NSCLC, also shows potent inhibitoy activity against ROS1. However, the development of serious resistance due to secondary mutations has been observed in clinical studies. To provide insight into the mechanisms of this drug resistance, molecular dynamics simulations and free-energy calculations were carried out for complexes of crizotinib with wild-type (WT) ROS1 as well as the mutated L2026M and G2032R forms. MD simulations indicated that the L2026M and G2032R systems are slightly less flexible than the WT system. Binding free energy calculations showed that the L2026M and G2032R mutations significantly reduce the binding affinity of crizotinib for ROS1, and that the resistance to crizotinib caused by the L2026M and G2032R mutations arises mostly from increases in entropic terms. Furthermore, calculations of per-residue binding free energies highlighted increased and decreased contributions of some residues in the L2026M and G2032R systems relative to those in the WT system. The present study therefore yielded detailed insight into the mechanisms of resistance to crizotinib caused by the L2026M and G2032R mutations, which should provide the basis for rational drug design to combat crizotinib resistance.

Open image in new window

Graphical Abstract Superposition of the average structures obtained from the last 10 ns of the molecular dynamics simulation trajectoriy for WT (green) and mutated ROS1 (cyan)

     

Exploring the crizotinib resistance mechanism of NSCLC with the L1196M mutation using molecular dynamics simulation

Crizotinib is an anticancer tyrosine kinase inhibitor that is approved for use as a first-line treatment for some non-small-cell lung cancers. L1196M is the most frequently observed mutation in NSCLC patients. This mutation, known as the gatekeeper mutation in the ALK kinase domain, confers resistance to crizotinib by sterically blocking the binding of the drug. However, the molecular mechanism of crizotinib resistance caused by the L1196M mutation is still unclear. Molecular dynamics simulation was therefore utilized in this study to investigate the mechanism by which the L1196M mutation may affect crizotinib resistance. Our results suggest that larger fluctuations in some important regions of the mutant complex compared to the wild-type complex may contribute to the resistance of the mutant complex to crizotinib. Also, mutation-induced alterations to the secondary structure of the complex as well as unstable hydrogen-bonding patterns in the A-loop and P-loop regions decrease the total binding energy of the complex. This study therefore provides a molecular explanation for the resistance to crizotinib caused by the L1196M mutation, which could aid the design of more efficient and selective drugs.     

Hydrogen bonding analysis of phosphoric acid–N,N-dimethylformamide mixtures

Crizotinib is the most effective and the only drug that has been approved for the treatment of anaplastic lymphoma kinase (ALK)-positive lung cancer. Reports suggest that there is a development of an acquired resistance against crizotinib action due to the emergence of several mutations in the ALK gene and F1174L is one such mutation. In this study, we used molecular docking and molecular dynamics (MD) approach to decipher the effect of F1174L mutation in drug–target binding. Docking results suggest that crizotinib was found to adopt the most promising conformations to the native-type ALK by identifying the M1199 residue as a prospective partner for making a hydrogen bond as compared to the mutant-type ALK. MD results showed that the average atom, especially atoms of the native-type ALK-crizotinib complex, movements were less, displayed less fluctuation, fast convergence of energy, and changes in geometry. This shows the stable binding of crizotinib with the native-type ALK in comparison to the mutant-type ALK. We believe that this study could be useful for the logical design of stronger, more selective, and more consistent ALK inhibitor against drug-resistant F1174L mutation.     

Insight into drug resistance mechanisms and discovery of potential inhibitors against wild-type and L1196M mutant ALK from FDA-approved drugs

Anaplastic lymphoma kinase (ALK) plays a crucial role in multiple malignant cancers. It is known as a well-established target for the treatment of ALK-dependent cancers. Even though substantial efforts have been made to develop ALK inhibitors, only crizotinib, ceritinib, and alectinib had been approved by the U.S. Food and Drug Administration for patients with ALK-positive non-small cell lung cancer (NSCLC). The secondary mutations with drug-resistance bring up difficulties to develop effective drugs for ALK-positive cancers. To give a comprehensive understanding of molecular mechanism underlying inhibitor response to ALK tyrosine kinase mutations, we established an accurate assessment for the extensive profile of drug against ALK mutations by means of computational approaches. The molecular mechanics-generalized Born surface area (MM-GBSA) method based on molecular dynamics (MD) simulation was carried out to calculate relative binding free energies for receptor-drug systems. In addition, the structure-based virtual screening was utilized to screen effective inhibitors targeting wild-type ALK and the gatekeeper mutation L1196M from 3180 approved drugs. Finally, the mechanism of drug resistance was discussed, several novel potential wild-type and L1196M mutant ALK inhibitors were successfully identified.     

Resistance to ROS1 Inhibition Mediated by EGFR Pathway Activation in Non-Small Cell Lung Cancer

The targeting of oncogenic ‘driver’ kinases with small molecule inhibitors has proven to be a highly effective therapeutic strategy in selected non-small cell lung cancer (NSCLC) patients. However, acquired resistance to targeted therapies invariably arises and is a major limitation to patient care. ROS1 fusion proteins are a recently described class of oncogenic driver, and NSCLC patients that express these fusions generally respond well to ROS1-targeted therapy. In this study, we sought to determine mechanisms of acquired resistance to ROS1 inhibition. To accomplish this, we analyzed tumor samples from a patient who initially responded to the ROS1 inhibitor crizotinib but eventually developed acquired resistance. In addition, we generated a ROS1 inhibition-resistant derivative of the initially sensitive NSCLC cell line HCC78. Previously described mechanisms of acquired resistance to tyrosine kinase inhibitors including target kinase-domain mutation, target copy number gain, epithelial-mesenchymal transition, and conversion to small cell lung cancer histology were found to not underlie resistance in the patient sample or resistant cell line. However, we did observe a switch in the control of growth and survival signaling pathways from ROS1 to EGFR in the resistant cell line. As a result of this switch, ROS1 inhibition-resistant HCC78 cells became sensitive to EGFR inhibition, an effect that was enhanced by co-treatment with a ROS1 inhibitor. Our results suggest that co-inhibition of ROS1 and EGFR may be an effective strategy to combat resistance to targeted therapy in some ROS1 fusion-positive NSCLC patients.     

A molecular dynamics investigation on the crizotinib resistance mechanism of C1156Y mutation in ALK

Crizotinib is an anaplastic lymphoma kinase (ALK) inhibitor that has recently been approved in the US for the treatment of non-small cell lung carcinoma (NSCLC). Despite its outstanding safety and efficacy, several resistant mutations against crizotinib have been detected in the treatment of NSCLC. However, in contrast to the widely accepted mechanism of steric hindrance by mutations at the active site, the mechanism by which the C1156Y non-active site mutation confers resistance against crizotinib remains unclear. In the present study, the resistance mechanism of C1156Y in ALK was investigated using molecular dynamics simulations. The results suggest that despite the non-active site mutation, C1156Y causes the dislocation of crizotinib as well as the indirect conformational changes in the binding cavity, which results in a marked decrease in the van der Waals and electrostatic interactions between crizotinib and ALK. The obtained results provide a detailed explanation of the resistance caused by C1156Y and may give a vital clue for the design of drugs to combat crizotinib resistance.     

Crizotinib-Induced Abnormal Signal Processing in the Retina

Molecular target therapy for cancer is characterized by unique adverse effects that are not usually observed with cytotoxic chemotherapy. For example, the anaplastic lymphoma kinase (ALK)-tyrosine kinase inhibitor crizotinib causes characteristic visual disturbances, whereas such effects are rare when another ALK-tyrosine kinase inhibitor, alectinib, is used. To elucidate the mechanism responsible for these visual disturbances, the responses to light exhibited by retinal ganglion cells treated with these agents were evaluated using a C57BL6 mouse ex vivo model. Both crizotinib and alectinib changed the firing rate of ON and OFF type retinal ganglion cells. However, the ratio of alectinib-affected cells (15.7%) was significantly lower than that of crizotinib-affected cells (38.6%). Furthermore, these drugs changed the response properties to light stimuli of retinal ganglion cells in some of the affected cells, i.e., OFF cells responded to both ON and OFF stimuli, etc. Finally, the expressions of ALK (a target receptor of both crizotinib and alectinib) and of MET and ROS1 (additional target receptors of crizotinib) were observed at the mRNA level in the retina. Our findings suggest that these drugs might target retinal ganglion cells and that the potency of the drug actions on the light responses of retinal ganglion cells might be responsible for the difference in the frequencies of visual disturbances observed between patients treated with crizotinib and those treated with alectinib. The present experimental system might be useful for screening new molecular target agents prior to their use in clinical trials.     

Oxazolidinone resistance mutations in 23S rRNA of Escherichia coli reveal the central region of domain V as the primary site of drug action

Oxazolidinone antibiotics inhibit bacterial protein synthesis by interacting with the large ribosomal subunit. The structure and exact location of the oxazolidinone binding site remain obscure, as does the manner in which these drugs inhibit translation. To investigate the drug-ribosome interaction, we selected Escherichia coli oxazolidinone-resistant mutants, which contained a randomly mutagenized plasmid-borne rRNA operon. The same mutation, G2032 to A, was identified in the 23S rRNA genes of several independent resistant isolates. Engineering of this mutation by site-directed mutagenesis in the wild-type rRNA operon produced an oxazolidinone resistance phenotype, establishing that the G2032A substitution was the determinant of resistance. Engineered U and C substitutions at G2032, as well as a G2447-to-U mutation, also conferred resistance to oxazolidinone. All the characterized resistance mutations were clustered in the vicinity of the central loop of domain V of 23S rRNA, suggesting that this rRNA region plays a major role in the interaction of the drug with the ribosome. Although the central loop of domain V is an essential integral component of the ribosomal peptidyl transferase, oxazolidinones do not inhibit peptide bond formation, and thus these drugs presumably interfere with another activity associated with the peptidyl transferase center.     

Molecular dynamics simulations provide insights into the origin of gleevec’s selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec’s selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec’s selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases.

Communicated by Ramaswamy H. Sarma     

Insight on mutation-induced resistance to anaplastic lymphoma kinase inhibitor ceritinib from molecular dynamics simulations

Ceritinib, an advanced anaplastic lymphoma kinase (ALK) next-generation inhibitor, has been proved excellent antitumor activity in the treatment of ALK-associated cancers. However, the accumulation of acquired resistance mutations compromise the therapeutic efficacy of ceritinib. Despite abundant mutagenesis data, the structural determinants for reduced ceritinib binding in mutants remains elusive. Focusing on the G1123S and F1174C mutations, we applied molecular dynamics (MD) simulations to study possible reasons for drug resistance caused by these mutations. The MD simulations predict that the studied mutations allosterically impact the configurations of the ATP-binding pocket. An important hydrophobic cluster is identified that connects P-loop and the αC-helix, which has effects on stabilizing the conformation of ATP-binding pocket. It is suggested, in this study, that the G1123S and F1174C mutations can induce the conformational change of P-loop thereby causing the reduced ceritinib affinity and causing drug resistance.     

Influence of Chirality of Crizotinib on Its MTH1 Protein Inhibitory Activity: Insight from Molecular Dynamics Simulations and Binding Free Energy Calculations

As a promising target for the treatment of lung cancer, the MutT Homolog 1 (MTH1) protein can be inhibited by crizotinib. A recent work shows that the inhibitory potency of (S)-crizotinib against MTH1 is about 20 times over that of (R)-crizotinib. But the detailed molecular mechanism remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of chirality of crizotinib on the inhibitory activity against MTH1. The binding free energy of (S)-crizotinib predicted by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Adaptive biasing force (ABF) methodologies is much lower than that of (R)-crizotinib, which is consistent with the experimental data. The analysis of the individual energy terms suggests that the van der Waals interactions are important for distinguishing the binding of (S)-crizotinib and (R)-crizotinib. The binding free energy decomposition analysis illustrated that residues Tyr7, Phe27, Phe72 and Trp117 were important for the selective binding of (S)-crizotinib to MTH1. The adaptive biasing force (ABF) method was further employed to elucidate the unbinding process of (S)-crizotinib and (R)-crizotinib from the binding pocket of MTH1. ABF simulation results suggest that the reaction coordinates of the (S)-crizotinib from the binding pocket is different from (R)-crizotinib. The results from our study can reveal the details about the effect of chirality on the inhibition activity of crizotinib to MTH1 and provide valuable information for the design of more potent inhibitors.     

Structural Mechanisms Determining Inhibition of the Collagen Receptor DDR1 by Selective and Multi-Targeted Type II Kinase Inhibitors

The discoidin domain receptors (DDRs), DDR1 and DDR2, form a unique subfamily of receptor tyrosine kinases that are activated by the binding of triple-helical collagen. Excessive signaling by DDR1 and DDR2 has been linked to the progression of various human diseases, including fibrosis, atherosclerosis and cancer. We report the inhibition of these unusual receptor tyrosine kinases by the multi-targeted cancer drugs imatinib and ponatinib, as well as the selective type II inhibitor DDR1-IN-1. Ponatinib is identified as the more potent molecule, which inhibits DDR1 and DDR2 with an IC₅₀ of 9 nM. Co-crystal structures of human DDR1 reveal a DFG-out conformation (DFG, Asp-Phe-Gly) of the kinase domain that is stabilized by an unusual salt bridge between the activation loop and αD helix. Differences to Abelson kinase (ABL) are observed in the DDR1 P-loop, where a β-hairpin replaces the cage-like structure of ABL. P-loop residues in DDR1 that confer drug resistance in ABL are therefore accommodated outside the ATP pocket. Whereas imatinib and ponatinib bind potently to both the DDR and ABL kinases, the hydrophobic interactions of the ABL P-loop appear poorly satisfied by DDR1-IN-1 suggesting a structural basis for its DDR1 selectivity. Such inhibitors may have applications in clinical indications of DDR1 and DDR2 overexpression or mutation, including lung cancer.     

Structural basis for drug resistance mechanisms against anaplastic lymphoma kinase

Drug resistance to anaplastic lymphoma kinase (ALK) inhibitors (crizotinib and ceritinib) is caused by mutation in the region encoding kinase domain of ALK. Compounds with potential ability to inhibit all strains of ALK are a solution to tackle the problem of drug resistance. In this study, we delineated positions of residues possessing the ability to make ALK drug resistant upon mutation by assessing them using five parameters (conservation index, binding-site root-mean-square deviation, protein structure stability, change in ATP, and drug-binding affinity). Four residual positions (Leu 1122, Thr 1151, Phe 1245, and Gly 1269) were ascertained. This study will be beneficial for designing drugs with better proficiency against ALK and the issues of drug resistance. This study can be taken as a pipeline for investigating drug-resistant mutations in other diseases as well.     

Insights into the drug resistance induced by the BaDHPS mutations: molecular dynamic simulations and MM/GBSA studies

Dihydropteroate synthase (DHPS) is essential for the folic acid biosynthetic pathway in prokaryotes; the mutation forms for DHPS are found to be relative to the urgent drug resistance problems. In our study, the Bacillus anthracis DHPS (BaDHPS) was selected for molecular dynamics and binding free energy studies to investigate the biochemistry behaviors of the wild-type and mutation form BaDHPS proteins (D184N and K220Q). It is found that the conformational change of the ligand dihydropteroate sulfathiazole binding site in mutation D184N and K220Q systems is mainly attributed from the Loop 1, Loop 2, and Loop 7 regions, and the binding free energy of these mutation systems is lower than that of the wild-type system. Additionally, some important hydrogen bonds of the mutation systems are disrupted during the simulations. But the shortening of the distance between residue Thr67 and the ligand would cause significant change of the binding pose in the K220Q system. These studies of DHPS family will be helpful for further drug resistance investigations.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:24     

序列变异对miR-378生物发生以及靶标关系的影响

本研究旨在探索序列变异对miR-378的结构、表达水平以及靶标关系的影响.利用PCR测序比对不同猪种miR-378的序列突变,预测突变对miR-378二级结构和自由能的改变;构建miR-378表达载体检测突变对其加工过程中各级产物表达水平的作用;运用TargetScan和TargetRank对突变引起的靶标关系变化进行...     

Isolation of a new melanoma antigen, MART-2, containing a mutated epitope recognized by autologous tumor-infiltrating T lymphocytes

Using cDNA expression cloning, a cDNA encoding a novel human melanoma Ag, MART-2 (melanoma Ag recognized by T cells-2), recognized by HLA-A1-restricted CD8(+) T cells from tumor-infiltrating lymphocytes (TIL1362) was isolated from an autologous melanoma cell line, 1362 mel. Homologous sequences to the cDNA had been registered in the EST database. This gene encoded an uncharacterized protein expressed ubiquitously in most normal and cancer cells. A mutation (A to G transition) was found in the cDNA obtained from the1362 mel melanoma cell line in the sequences encoding the phosphate binding loop (P-loop) that resulted in loss of the ability to bind GTP. Transfection of NIH-3T3 with the mutated MART-2 did not result in the development of significant foci. By screening 36 various cancer cell lines using single-strand conformation polymorphism, a possible mutation in the P-loop of MART-2 was found in one squamous cell lung cancer cell line, EBC1. The T cell epitope for TIL1362, FLEGNEVGKTY, was identified to be encoded by the mutated sequence of the MART-2 Ag. The mutation substituted glycine in the normal peptide with glutamic acid at the third amino acid of the epitope, which is an important primary anchor amino acid for HLA-A1 peptide binding. The normal peptide, FLGGNEVGKTY, was not recognized by TIL1362, suggesting that this T cell response was specific for the autologous tumor. Although transforming activity was not detected in the NIH-3T3 assay, MART-2 with the mutation in the P-loop may be involved in the generation of melanoma through a loss of GTP binding activity.     

How does the novel T315L mutation of breakpoint cluster region-abelson (BCR-ABL) kinase confer resistance to ponatinib: a comparative molecular dynamics simulation study

Abstract

Acute lymphocytic leukemia (ALL) is one of the most dangerous types of leukemia, and about 40% of them is Philadelphia chromosome-positive acute lymphocytic leukemia (Ph?+?ALL). Ph?+?ALL is caused by the fusion of the breakpoint cluster region (BCR) and the Ableson (ABL) genes, named the BCR-ABL fused gene that codes for an autonomously active tyrosine kinase. Tyrosine kinase inhibitors (TKIs) are among the first-line therapeutic agents for the treatment of Ph?+?ALL. Drug resistance are the major obstacle, limiting their clinical utility. The latest third-generation TKIs, ponatinib, can tackle most abnormal BCR-ABL kinases, including the T315I mutant that is resistant to first- and second-generations TKIs such as imatinib. However, drug resistance still emerges with the novel T315L mutation and the underlying mechanisms remain elusive. Here, using molecular dynamics (MD) simulations, we explored into the detailed interactions between ponatinib and BCR-ABL in the wild-type (WT), T315I, and T315L systems. The simulations revealed the significant conformational changes of ponatinib in its binding site due to the T315L mutation and the underlying structural mechanisms. Binding free energy analysis unveiled that the affinity of ponatinib to BCR-ABL decreased upon T315L mutation, which resulted in its unfavorable binding and drug resistance. Key residues responsible for the unfavored unbinding were also identified. This study elucidates the detailed mechanisms for the resistance of ponatinib in Ph?+?ALL triggered by the T315L mutation and will provide insights for future drug development and optimization.     

Identification of direct target engagement biomarkers for kinase-targeted therapeutics

Pharmacodynamic (PD) biomarkers are an increasingly valuable tool for decision-making and prioritization of lead compounds during preclinical and clinical studies as they link drug-target inhibition in cells with biological activity. They are of particular importance for novel, first-in-class mechanisms, where the ability of a targeted therapeutic to impact disease outcome is often unknown. By definition, proximal PD biomarkers aim to measure the interaction of a drug with its biological target. For kinase drug discovery, protein substrate phosphorylation sites represent candidate PD biomarkers. However, substrate phosphorylation is often controlled by input from multiple converging pathways complicating assessment of how potently a small molecule drug hits its target based on substrate phoshorylation measurements alone. Here, we report the use of quantitative, differential mass-spectrometry to identify and monitor novel drug-regulated phosphorylation sites on target kinases. Autophosphorylation sites constitute clinically validated biomarkers for select protein tyrosine kinase inhibitors. The present study extends this principle to phosphorylation sites in serine/threonine kinases looking beyond the T-loop autophosphorylation site. Specifically, for the 3'-phosphoinositide-dependent protein kinase 1 (PDK1), two phospho-residues p-PDK1(Ser410) and p-PDK1(Thr513) are modulated by small-molecule PDK1 inhibitors, and their degree of dephosphorylation correlates with inhibitor potency. We note that classical, ATP-competitive PDK1 inhibitors do not modulate PDK1 T-loop phosphorylation (p-PDK1(Ser241)), highlighting the value of an unbiased approach to identify drug target-regulated phosphorylation sites as these are complementary to pathway PD biomarkers. Finally, we extend our analysis to another protein Ser/Thr kinase, highlighting a broader utility of our approach for identification of kinase drug-target engagement biomarkers.