Molecular Simulation Studies on the Binding Selectivity of Type-I Inhibitors in the Complexes with ROS1 versus ALK

Current treatment and novel insights regarding ROS1-targeted therapy in malignant tumors

BACKGROUND

The proto-oncogene ROS1 encodes an intrinsic type I membrane protein of the tyrosine kinase/insulin receptor family. ROS1 facilitates the progression of various malignancies via self-mutations or rearrangements. Studies on ROS1-directed tyrosine kinase inhibitors have been conducted, and some have been approved by the FDA for clinical use. However, the adverse effects and mechanisms of resistance associated with ROS1 inhibitors remain unknown. In addition, next-generation ROS1 inhibitors, which have the advantage of treating central nervous system metastases and alleviating endogenous drug resistance, are still in the clinical trial stage.

METHOD

In this study, we searched relevant articles reporting the mechanism and clinical application of ROS1 in recent years; systematically reviewed the biological mechanisms, diagnostic methods, and research progress on ROS1 inhibitors; and provided perspectives for the future of ROS1-targeted therapy.

RESULTS

ROS1 is most expressed in malignant tumours. Only a few ROS1 kinase inhibitors are currently approved for use in NSCLC, the efficacy of other TKIs for NSCLC and other malignancies has not been ascertained. There is no effective standard treatment for adverse events or resistance to ROS1-targeted therapy. Next-generation TKIs appear capable of overcoming resistance and delaying central nervous system metastasis, but with a greater incidence of adverse effects.

CONCLUSIONS

Further research on next-generation TKIs regarding the localization of ROS1 and its fusion partners, binding sites for targeted drugs, and coadministration with other drugs is required. The correlation between TKIs and chemotherapy or immunotherapy in clinical practice requires further study.

Hot-Spot Residue-Based Virtual Screening of Novel Selective Estrogen-Receptor Degraders for Breast Cancer Treatment

The estrogen-receptor alfa (ERα) is considered pivotal for breast cancer treatment. Although selective estrogen-receptor degraders (SERDs) have been developed to induce ERα degradation and antagonism, their agonistic effect on the uterine tissue and poor pharmacokinetic properties limit further application of ERα; thus, discovering novel SERDs is necessary. The ligand preferentially interacts with several key residues of the protein (defined as hot-spot residues). Improving the interaction with hot-spot residues of ERα offers a promising avenue for obtaining novel SERDs. In this study, pharmacophore modeling, molecular mechanics/generalized Born surface area (MM/GBSA), and amino-acid mutation were combined to determine several hot-spot residues. Focusing on the interaction with these hot-spot residues, hit fragments A1-A3 and A9 were virtually screened from two fragment libraries. Finally, these hit fragments were linked to generate compounds B1-B3, and their biological activities were evaluated. Remarkably, compound B1 exhibited potent antitumor activity against MCF-7 cells (IC50 = 4.21 nM), favorable ERα binding affinity (Ki = 14.6 nM), and excellent ERα degradative ability (DC50 = 9.7 nM), which indicated its potential to evolve as a promising SERD for breast cancer treatment.

Theoretical Studies on Selectivity of HPK1/JAK1 Inhibitors by Molecular Dynamics Simulations and Free Energy Calculations

Hematopoietic progenitor kinase 1 (HPK1) is a negative regulator of T cell receptor, which has been regarded as a potential target for immunotherapy. Yu et al. observed the off-target effect of the high-throughput screening HPK1 kinase inhibitor hits on JAK1 kinase. The off-target effect is usually due to the lack of specificity of the drug, resulting in toxic side effects. Therefore, exploring the mechanisms to selectively inhibit HPK1 is critical for developing effective and safe inhibitors. In this study, two indazole compounds as HPK1 inhibitors with different selectivity towards JAK1 were used to investigate the selectivity mechanism using multiple computational methods, including conventional molecular dynamics simulations, binding free energy calculations and umbrella sampling simulations. The results indicate that the salt bridge between the inhibitor and residue Asp101 of HPK1 favors their selectivity towards HPK1 over JAK1. Information obtained from this study can be used to discover and design more potent and selective HPK1 inhibitors for immunotherapy.

Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation

Covalent modification of the oncogenic mutant epidermal growth factor receptor (EGFR) by small molecules is an efficient strategy for achieving an enhanced and sustained pharmacological effect in the treatment of non-small-cell lung cancer. NSP-037 (18), an irreversible inhibitor of the L858R/T790M double-mutant EGFR (EGFR^DM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFR^DM over the wild type (EGFR^WT), as compared to clinically approved osimertinib (7). Here, we employ multiple computational approaches to elucidate the mechanism underlining this improved selectivity, as well as the effect of CFA on the selectivity enhancement of inhibitor 18 over 7. We find that EGFR^DM undergoes significantly larger conformational changes than EGFR^WT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFR^DM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFR^WT and EGFR^DM and 18's higher selectivity than 7. The binding free energy of the 18-bound complexes is -6.38 kcal/mol greater than that of the 7-bound complexes, explaining the difference in selectivity of these inhibitors. Further, free energy decomposition analysis indicates that the electrostatic contribution of key residues plays an important role in the 18-bound complexes. QM/MM calculations show that the most favored mechanism for the Cys797 alkylation reaction is the direct displacement mechanism through a CFA-based inhibitor, producing a reaction with the lowest energy barrier and most stable product.

Molecular Mechanism Exploration of Potent Fluorinated PI3K Inhibitors with a Triazine Scaffold: Unveiling the Unusual Synergistic Effect of Pyridine-to-Pyrimidine Ring Interconversion and CF3 Defluorination

The phosphatidylinostitol-3-kinase (PI3K)/AKT/mammalian target of rapamycin signaling pathway is a vital regulator of cell proliferation, growth, and survival, which is frequently overactivated in many human cancers. To this effect, PI3K, which is an important mediator of this pathway, has been pinpointed as a crucial target in cancer therapy and hence the importance of PI3K inhibitors. It was recently reported that defluorination and pyridine-to-pyrimidine ring interconversion increase the potency of specific small-molecule inhibitors of PI3K. Compound 4, an inhibitor with the difluorinated pyrimidine motif, was found to be eight times more potent against PI3K than compound 1, an inhibitor with the trifluorinated pyridine motif. This observation presents the need to rationally resolve the differential inhibitory mechanisms exhibited by both compounds. In this present work, we employed multiple computational approaches to investigate and distinguish the binding modes of 1 and 4 in addition to the effects they mediate on the secondary structure of PI3K. Likewise, we evaluated two other derivatives, compounds 2 with the difluorinated pyridine motif and 3 with the trifluorinated pyrimidine motif, to investigate the cooperativity effect between the defluorination of CF₃ and pyridine-to-pyrimidine ring interconversion. Findings revealed that PI3K, upon interaction with 4, exhibited a series of structural changes that favored the binding of the inhibitor at the active-site region. Furthermore, a positive (synergistic) cooperativity effect was observed between CF₃ defluorination and pyridine-to-pyrimidine ring interconversion. Moreover, there was a good correlation between the binding free energy estimated and the biological activity reported experimentally. Energy decomposition analysis revealed that the major contributing force to binding affinity variations between 1 and 4 is the electrostatic energy. Per-residue energy-based hierarchical clustering analysis further identified four hot-spot residues ASP841, TYR867, ASP964, and LYS833 and four warm-spot residues ASP836, SER806, ASP837, and LYS808, which essentially mediate the optimal and higher-affinity binding of compound 4 to PI3K relative to 1. This study therefore provides rational insights into the mechanisms by which 4 exhibited superior PI3K-inhibitory activities over 1, which is vital for future structure-based drug discovery efforts in PI3K targeting.

Identification of Flavonoids as Putative ROS-1 Kinase Inhibitors Using Pharmacophore Modeling for NSCLC Therapeutics

Non-small cell lung cancer (NSCLC) is a lethal non-immunogenic malignancy and proto-oncogene ROS-1 tyrosine kinase is one of its clinically relevant oncogenic markers. The ROS-1 inhibitor, crizotinib, demonstrated resistance due to the Gly2032Arg mutation. To curtail this resistance, researchers developed lorlatinib against the mutated kinase. In the present study, a receptor-ligand pharmacophore model exploiting the key features of lorlatinib binding with ROS-1 was exploited to identify inhibitors against the wild-type (WT) and the mutant (MT) kinase domain. The developed model was utilized to virtually screen the TimTec flavonoids database and the retrieved drug-like hits were subjected for docking with the WT and MT ROS-1 kinase. A total of 10 flavonoids displayed higher docking scores than lorlatinib. Subsequent molecular dynamics simulations of the acquired flavonoids with WT and MT ROS-1 revealed no steric clashes with the Arg2032 (MT ROS-1). The binding free energy calculations computed via molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) demonstrated one flavonoid (Hit) with better energy than lorlatinib in binding with WT and MT ROS-1. The Hit compound was observed to bind in the ROS-1 selectivity pocket comprised of residues from the β-3 sheet and DFG-motif. The identified Hit from this investigation could act as a potent WT and MT ROS-1 inhibitor.

Mechanistic Study of Potent Fluorinated EGFR Kinase Inhibitors with a Quinazoline Scaffold against L858R/T790M/C797S Resistance Mutation: Unveiling the Fluorine Substituent Cooperativity Effect on the Inhibitory Activity

Fluorination has considerable potential with regard to the design of kinase inhibitors for anticarcinoma therapy. It was recently reported that fluorination increases the potency of inhibitors of the epidermal growth factor receptor (EGFR), mutations of which have been linked specifically to nonsmall-cell lung cancer. For the L858R/T790M/C797S triplet mutant (EGFR^TM), a difluorinated inhibitor, 25g, was found to have 4.23 times greater potency against the EGFR^TM than an unfluorinated inhibitor, 25a. This discovery necessitates a rational explanation for the underlying inhibitory mechanisms. Here, we apply multiple computational approaches to explore, validate, and differentiate the binding modes of 25a and 25g in the EGFR^TM and investigate the cooperativity effect of fluorine substituents on the inhibitory activity. Our results showed that the EGFR^TM in the presence of 25g undergoes a series of conformational changes that favor inhibitor binding to both the active and allosteric sites. Further, the cooperativity effect of fluorine substituents is positive: the complex stability is increased by each additional fluorine substituent. Estimated binding free energies show good correlation with the experimental biological activity. Subsequently, the decomposition energy analysis revealed that the van der Waals interaction is the principal force contributing to variations in the binding affinities of 25a and 25g to the EGFR^TM. Per-residue energy-based hierarchical clustering analysis suggests that three hot-spot residues, L718, K745, and D855, are the key in achieving optimal binding modes for 25g with higher affinity in the EGFR^TM compared to 25a. This study provides a rationale for the superior EGFR^TM-inhibitory potency exhibited by 25g over 25a, which is expected to be useful for the future rational structure-based design of novel EGFR^TM inhibitors with improved potency and selectivity.

A Mechanistic Study of a Potent and Selective Epidermal Growth Factor Receptor Inhibitor against the L858R/T790M Resistance Mutation

Covalent targeting is a promising strategy for increasing the potency and selectivity of potential drug candidates. This therapeutic approach was recently reported for the epidermal growth factor receptor (EGFR), wherein a covalent binder, 20g [N-(3-{7-[2-methoxy-4-(4-methylpiperazin-1-yl)phenylamino]-3,4-dihydro-3-isopropyl-2,4-dioxopyrimido[4,5-d]pyrimidin-1(2H)-yl}phenyl)acrylamide], demonstrated significant selectivity and inhibitory activity toward the EGFR L858R/T790M double mutant (EGFR^DM) relative to the EGFR wild-type form (EGFR^WT). The enhanced therapeutic potency of 20g against EGFR^DM is 263 times greater than that against EGFR^WT, which necessitates a rational explanation for the underlying selective and inhibitory mechanisms. In this work, we investigate the differential binding modes of 20g with EGFR^WT and EGFR^DM using molecular dynamics simulations coupled with free energy calculations and further identify key residues involved in the selective targeting, binding, and inhibitory mechanisms mediated by 20g. We find that systematic orientational and conformational changes in the α-loop, p-loop, active loop, and αC-helix are responsible for the disparate binding mechanisms and inhibitory prowess of 20g with respect to EGFR^WT and EGFR^DM. The calculated binding free energies show good correlation with the experimental biological activity. The total binding free energy difference between EGFR^WT-20g and EGFR^DM-20g is -11.47 kcal/mol, implying that 20g binds more strongly to EGFR^DM. This enhanced binding affinity of 20g for EGFR^DM is a result of a large increase in the van der Waals and electrostatic interactions with three critical residues (Met790, Gln791, and Met793) that are chiefly responsible for the high-affinity interactions mediated by 20g with EGFR^DM relative to EGFR^WT.

Importance of Incorporating Protein Flexibility in Molecule Modeling: A Theoretical Study on Type I1/2 NIK Inhibitors

NF-κB inducing kinase (NIK), which is considered as the central component of the non-canonical NF-κB pathway, has been proved to be an important target for the regulation of the immune system. In the past few years, NIK inhibitors with various scaffolds have been successively reported, among which type I^1/2 inhibitors that can not only bind in the ATP-binding pocket at the DFG-in state but also extend into an additional back pocket, make up the largest proportion of the NIK inhibitors, and are worthy of more attention. In this study, an integration protocol that combines molecule docking, MD simulations, ensemble docking, MM/GB(PB)SA binding free energy calculations, and decomposition was employed to understand the binding mechanism of 21 tricyclic type I^1/2 NIK inhibitors. It is found that the docking accuracy is largely dependent on the selection of docking protocols as well as the crystal structures. The predictions given by the ensemble docking based on multiple receptor conformations (MRCs) and the MM/GB(PB)SA calculations based on MD simulations showed higher linear correlations with the experimental data than those given by conventional rigid receptor docking (RRD) methods (Glide, GOLD, and Autodock Vina), highlighting the importance of incorporating protein flexibility in predicting protein–ligand interactions. Further analysis based on MM/GBSA demonstrates that the hydrophobic interactions play the most essential role in the ligand binding to NIK, and the polar interactions also make an important contribution to the NIK-ligand recognition. A deeper comparison of several pairs of representative derivatives reveals that the hydrophobic interactions are vitally important in the structural optimization of analogs as well. Besides, the H-bond interactions with some key residues and the large desolvation effect in the back pocket devote to the affinity distinction. It is expected that our study could provide valuable insights into the design of novel and potent type I^1/2 NIK inhibitors.

Covalent vs. Non-Covalent Inhibition: Tackling Drug Resistance in EGFR - A Thorough Dynamic Perspective

A persistent challenge in the treatment of non-small cell lung cancer (NSCLC) with EGFR is the emergence of drug-resistant caused by somatic mutations. The EGFR L858R/T790 M double mutant (EGFR^DM ) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFR^DM effectively. Recently, 11h and 45a, have been found to be potent inhibitors against EGFR^DM through two distinctive mechanisms, non-covalent and covalent binding, respectively. However, the structural and dynamic implications of the two modes of inhibitions remain unexplored. Herein, two molecular dynamics simulation protocols, coupled with free-energy calculations, were applied to gain insight into the atomistic nature of each binding mode. The comparative analysis confirmed that there is a significant difference in the binding free energy between 11h and 45a (ΔΔG_bind =-21.17 kcal/mol). The main binding force that governs the binding of both inhibitors is vdW, with a higher contribution for 45a. Two residues ARG841 and THR854 were found to have curtailed role in the binding of 45a to EGFR^DM by stabilizing its flexible alcohol chain. The 45a binding to EGFR^DM induces structural rearrangement in the active site to allow easier accessibility of 45a to target residue CYS797. The findings of this work can substantially shed light on new strategies for developing novel classes of covalent and non-covalent inhibitors with increased specificity and potency.

Evolution strategy of ROS1 kinase inhibitors for use in cancer therapy

The abnormal expression of c-ros oncogene1 receptor tyrosine kinase (ROS1) has been identified as clinically actionable oncogenic driver in non-small-cell lung cancer. Since crizotinib was approved by the US FDA for the treatment of advanced ROS1-positive non-small-cell lung cancer, ROS1 kinase has become a promising therapeutic target. Under the guidance of some advanced computer-assisted technologies, such as structure-based drug design, homology modeling and lipophilic efficiency parameters, several potent and selective inhibitors against wild-type and mutant ROS1 were designed and synthesized. In this article, we will review a series of scaffolds targeting ROS1 kinase from the hit-to-drug evolution strategies of their representative compounds and it is hoped that these design strategies would facilitate medicinal chemists to optimize the process of drug design.

Sulfonated Compounds Bind with Prostatic Acid Phosphatase (PAP248-286) to Inhibit the Formation of Amyloid Fibrils

The peptide segment of prostatic acid phosphatase (PAP248-286) aggregates to form SEVI (semen-derived enhancer of virus infection) amyloid fibrils. These are characteristic seminal amyloids that have the ability to promote the effect of HIV infection. In this paper, we explore the binding of sulfonated compounds with PAP248-286 through an in silico study. Three derivatives of suramin, NF110, NF279, and NF340, are selected. All of these sulfonated molecules bind to PAP248-286 and alter the conformation of the peptide, even though they have various structures, sizes, and configurations. The compounds bind with PAP248-286 through multiple interactions, such as hydrogen-bonding interactions, hydrophobic interactions, π-π stacking interactions, and electrostatic interactions. However, NF110, which has an X-shaped configuration, has the highest binding affinity of the three derivatives investigated. We also perform surface plasmon resonance and a Congo red assay to validate the results. The interactions between PAP248-286 and the sulfonated compounds are proposed to depend on the orientations of the sulfonate groups and the specific configurations of the compounds instead of the number of sulfonate groups. NF110 molecules occupy the exposed binding sites of PAP248-286, blocking interactions between the peptides. Therefore, these compounds are important in inhibiting the aggregation of PAP248-286. Herein, we provide useful information to develop new efficient microbicides to antagonize seminal amyloid fibrils and to block HIV transmission.

Design, synthesis, biological evaluation and molecular modeling of novel 2-amino-4-(1-phenylethoxy) pyridine derivatives as potential ROS1 inhibitors

With the aim of discovering potential and selective inhibitors targeting ROS1 kinase, we rationally designed, synthesized and evaluated two series of novel 2-amino-pyridine derivatives with 1-phenylethoxy at C-3 and C-4 position. The enzymic assays results indicated that six of the new compounds 13b-13d and 14a-14c showed remarkably higher inhibitory activities against ROS1 kinase. The most promising compounds, 13d and 14c displayed the most desired ROS1 inhibitory activity with IC₅₀ values of 440 nM and 370 nM respectively. Furthermore, 13d and 14c displayed ROS1 inhibitory selectivity of about 7-fold and 12-fold, relative to that of ALK sharing about 49% amino acid sequence homology in the kinase domains. They also showed good anti-proliferative effects against ROS1-addicted HCC78 cell lines with the IC₅₀ values of 8.1 μM and 65.3 μM, respectively. Moreover, molecular docking and molecular dynamics simulation studies disclosed that compound 14c and 13d shared similar binding poses with Crizotinib except the selective binding site of ROS1. It also gave a probable molecular explanation for their activity and selectivity, which the methoxyl group in benzene ring was the crucial to the selectivity to ROS1 versus ALK.