Mechanistic Insights Into Co-Administration of Allosteric and Orthosteric Drugs to Overcome Drug-Resistance in T315I BCR-ABL1

Exploring the allosteric effect of SHP2 Tyr62 phosphorylation on the emergence of acquired resistance to allosteric inhibitor SHP099

The Src homology-2 (SH2)-containing phosphatase 2 (SHP2), a non-receptor protein tyrosine phosphatase, is a key regulator modulating various signaling pathways. Recent studies have revealed that phosphorylation of Tyr62 (pY62) on the N-SH2 domain of SHP2 causes the emergence of acquired resistance to the allosteric inhibitor of SHP2 (SHP099) that occupies the PTP catalytic domain. However, the allosteric mechanism underlying the insensitivity of the allosteric inhibitor SHP099 to the phosphorylated SHP2 (pSHP2) remains unexplored. In this study, multiple replica molecular dynamics (MD) simulations and the post-trajectory analyses (principal component analysis, dynamics cross-correlation matrix analysis, allosteric community analysis, and binding free energy calculations) were performed for the SHP2, pSHP2, SHP2-SHP099, and pSHP2-SHP099 complexes. MD results showed that SHP099 binding contributed to stabilize SHP2, but pY62 had a detrimental role in the stability of the pSHP2-SHP099 complex. Domain correlation analysis showed that pY62 increased the anti-correlated motions between the C-SH2 and N-SH2/PTP domains. Binding free energy calculations revealed that the protein-ligand interactions in the SHP2 - SHP099 complex were stronger than that of the pSHP2 - SHP099 complex. Further, Thr108, Phe113, and Glu250 might be the critical residues responsible for the loss of the binding affinity in the pSHP2 - SHP099 complex through a per-residue decomposition analysis and H-bond occupancy time analysis. Overall, this study may provide a mechanistic insight into the mechanism how the allosteric effect of pY62 of SHP2 on SHP099 binding.

Structure and Dynamics of the ABL1 Tyrosine Kinase and Its Important Role in Chronic Myeloid Leukemia

Abelson (ABL1) tyrosine kinase is an essential component of non-receptor tyrosine kinases and is associated with numerous cellular processes, including differentiation and proliferation. The structural features of ABL1 include a distinct N-terminal cap region, a C-terminal tail, a bilobed kinase, SH2, and SH3 domains. These domains enable its engagement in several signaling cascades and dynamic control. The pathophysiology of chronic myeloid leukemia (CML) is mainly driven by the BCR-ABL1 oncoprotein, arising from dysregulation of ABL1 kinase, namely through its fusion to the breakpoint cluster region (BCR) gene. ABL1 is a crucial target in the treatment of CML as the BCR-ABL1 fusion causes uncontrolled cellular proliferation and resistance to apoptosis. Tyrosine kinase inhibitors (TKIs) targeting the ABL1 tyrosine kinase are playing a critical role in the treatment of CML through the inhibition of persistently activated signaling pathways mediated by the BCR-ABL1 fusion protein. The article examines the structural characteristics of ABL1, how they relate to CML, and the interactions between ABL1 and the current FDA-approved TKIs, emphasizing the kinase's critical function in carcinogenesis and its possible target status for tyrosine kinase inhibitors.

Nanohydrogel of Curcumin/Berberine Co-Crystals Induces Apoptosis via Dual Covalent/Noncovalent Inhibition of Caspases in Endometrial Cancer Cell Lines: The Synergy Between Pharmacokinetics and Pharmacodynamics

Endometrial cancer remains a significant therapeutic challenge due to drug resistance and heterogeneity. This study leverages the synergistic potential of curcumin (CUR) and berberine (BBR) co-crystals encapsulated in a nanohydrogel to address these challenges through a pharmacokinetically and pharmacodynamically targeted therapeutic strategy. The nanohydrogel formulation significantly improves the solubility, stability, and bioavailability of CUR/BBR co-crystals, optimizing their therapeutic delivery and sustained release under physiological and tumor microenvironment conditions. On the other hand, the dual inhibitory mechanism of CUR and BBR, with CUR covalently binding to the active site of caspase-3 and BBR non-covalently targeting the allosteric site, achieves enhanced apoptotic activity by disrupting both the catalytic and conformational functions of caspase-3. In vitro cytotoxicity assays demonstrate remarkable efficacy of the CUR/BBR nanohydrogel, achieving an IC50 of 12.36 μg/mL against HEC-59 endometrial cancer cells, significantly outperforming the individual components and the standard drug Camptothecin (IC50: 17.27 μg/mL). Caspase-3/7 assays confirm enhanced apoptosis induction for the nanohydrogel formulation compared to co-crystals alone and Camptothecin. Molecular dynamics simulations and binding free energy analyses further validate the synergistic interaction of CUR and BBR in their dual binding mode. This study introduces a novel therapeutic approach by enhancing drug delivery and dual targeting mechanisms, demonstrating the potential of CUR-BBR nanohydrogel as a robust therapy for EC. This strategy offers a promising platform for addressing drug resistance and improving outcomes in endometrial cancer therapy.

Clinical Pharmacology of Asciminib: A Review

Asciminib is a first-in-class allosteric inhibitor of the kinase activity of BCR::ABL1, specifically targeting the ABL myristoyl pocket (STAMP). This review focuses on the pharmacokinetic (PK) and pharmacodynamic data of asciminib, which is approved at a total daily dose of 80 mg for the treatment of adult patients with chronic myeloid leukemia in chronic phase who are either resistant or intolerant to ≥ 2 tyrosine kinase inhibitors or those harboring the T315I mutation (at a dose of 200 mg twice daily). Asciminib is predicted to be almost completely absorbed from the gut, with an absolute bioavailability (F) of approximately 73%. It should be administered in a fasted state, as food (particularly high-fat meals) reduces exposure. Asciminib displays a slightly greater than dose-proportional increase in exposure, with no time-dependent changes in PK observed following repeated dosing. This drug shows low clearance (6.31 L/h), with a moderate volume of distribution (111 L) and high human plasma protein binding (97.3%). The apparent terminal elimination half-life (t1/2) across studies was estimated to be between 7 and 15 h. The PK of asciminib is not substantially affected by body weight, age, gender, race, or renal or hepatic impairment. Asciminib is primarily metabolized via CYP3A4-mediated oxidation (36.0%) and UGT2B7- and UGT2B17-mediated glucuronidation (13.3% and 7.8%, respectively); biliary secretion via breast cancer resistance protein contributes to about 31.1% to total systemic clearance, which is mainly through hepatic metabolism and biliary secretion through the fecal pathway, with renal excretion playing a minor role. The potential for PK drug interaction for asciminib both as a victim and a perpetrator has been summarized here based on clinical and predicted drug-drug interaction studies. Robust exposure-response models characterized asciminib exposure-efficacy and exposure-safety relationships. In patients without the T315I mutation, the exposure-efficacy analysis of the time course of BCR::ABL1IS percentages highlighted the existence of a slightly positive, albeit not clinically significant, relationship. Higher exposure was required for efficacy in patients harboring the T315I mutation compared with those who did not. The exposure-safety relationship analysis showed no apparent association between exposure and adverse events of interest over the broad range of exposure or dose levels investigated. Asciminib has also been shown to have no clinically relevant effect on cardiac repolarization. Here, we review the clinical pharmacology data available to date for asciminib that supported its clinical development program and regulatory applications.

Predictive identification and design of potent inhibitors targeting resistance-inducing candidate genes from E. coli whole-genome sequences

Introduction: This work utilizes predictive modeling in drug discovery to unravel potential candidate genes from Escherichia coli that are implicated in antimicrobial resistance; we subsequently target the gidB, MacB, and KatG genes with some compounds from plants with reported antibacterial potentials. Method: The resistance genes and plasmids were identified from 10 whole-genome sequence datasets of E. coli; forty two plant compounds were selected, and their 3D structures were retrieved and optimized for docking. The 3D crystal structures of KatG, MacB, and gidB were retrieved and prepared for molecular docking, molecular dynamics simulations, and ADMET profiling. Result: Hesperidin showed the least binding energy (kcal/mol) against KatG (-9.3), MacB (-10.7), and gidB (-6.7); additionally, good pharmacokinetic profiles and structure-dynamics integrity with their respective protein complexes were observed. Conclusion: Although these findings suggest hesperidin as a potential inhibitor against MacB, gidB, and KatG in E. coli, further validations through in vitro and in vivo experiments are needed. This research is expected to provide an alternative avenue for addressing existing antimicrobial resistances associated with E. coli's MacB, gidB, and KatG.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations

CONTEXT

ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation.

METHODS

In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.

Allosteric regulation and inhibition of protein kinases

The human genome encodes more than 500 different protein kinases: signaling enzymes with tightly regulated activity. Enzymatic activity within the conserved kinase domain is influenced by numerous regulatory inputs including the binding of regulatory domains, substrates, and the effect of post-translational modifications such as autophosphorylation. Integration of these diverse inputs occurs via allosteric sites that relate signals via networks of amino acid residues to the active site and ensures controlled phosphorylation of kinase substrates. Here, we review mechanisms of allosteric regulation of protein kinases and recent advances in the field.

Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations

OBJECTIVE

RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood.

METHODS

Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant.

RESULTS

The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant.

CONCLUSIONS

The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.

Mechanistic Insights into the Protection Effect of Argonaute-RNA Complex on the HCV Genome

While host miRNA usually plays an antiviral role, the relentless tides of viral evolution have carved out a mechanism to recruit host miRNA as a viral protector. By complementing miR-122 at the 5' end of the genome, the hepatitis C virus (HCV) gene can form a complex with Argonaute 2 (Ago2) protein to protect the 5' end of HCV RNA from exonucleolytic attacks. Experiments showed that the disruption of the stem-loop 1(SL1) structure and the 9th nucleotide (T9) of HCV site 1 RNA could enhance the affinity of the Ago2 protein to the HCV site 1 RNA (target RNA). However, the underlying mechanism of how the conformation and dynamics of the Ago2: miRNA: target RNA complex is affected by the SL1 and T9 remains unclear. To address this, we performed large-scale molecular dynamics simulations on the AGO2-miRNA complex binding with the WT target, T9-abasic target and SL1-disruption target, respectively. The results revealed that the T9 and SL1 structures could induce the departing motion of the PAZ, PIWI and N domains, propping up the mouth of the central groove which accommodates the target RNA, causing the instability of the target RNA and disrupting the Ago2 binding. The coordinated motion among the PAZ, PIWI and N domains were also weakened by the T9 and SL1 structures. Moreover, we proposed a new model wherein the Ago2 protein could adopt a more constraint conformation with the proximity and more correlated motions of the PAZ, N and PIWI domains to protect the target RNA from dissociation. These findings reveal the mechanism of the Ago2-miRNA complex's protective effect on the HCV genome at the atomic level, which will offer guidance for the design of drugs to confront the protection effect and engineering of Ago2 as a gene-regulation tool.

Mechanistic insights into the clinical Y96D mutation with acquired resistance to AMG510 in the KRASG12C

Special oncogenic mutations in the RAS proteins lead to the aberrant activation of RAS and its downstream signaling pathways. AMG510, the first approval drug for KRAS, covalently binds to the mutated cysteine 12 of KRASG12C protein and has shown promising antitumor activity in clinical trials. Recent studies have reported that the clinically acquired Y96D mutation could severely affect the effectiveness of AMG510. However, the underlying mechanism of the drug-resistance remains unclear. To address this, we performed multiple microsecond molecular dynamics simulations on the KRASG12C−AMG510 and KRASG12C/Y96D−AMG510 complexes at the atomic level. The direct interaction between the residue 96 and AMG510 was impaired owing to the Y96D mutation. Moreover, the mutation yielded higher flexibility and more coupled motion of the switch II and α3-helix, which led to the departing motion of the switch II and α3-helix. The resulting departing motion impaired the interaction between the switch II and α3-helix and subsequently induced the opening and loosening of the AMG510 binding pocket, which further disrupted the interaction between the key residues in the pocket and AMG510 and induced an increased solvent exposure of AMG510. These findings reveal the resistance mechanism of AMG510 to KRASG12C/Y96D, which will help to offer guidance for the development of KRAS targeted drugs to overcome acquired resistance.

Mechanistic Insights into the Long-range Allosteric Regulation of KRAS Via Neurofibromatosis Type 1 (NF1) Scaffold Upon SPRED1 Loading

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1-SUMO1 complex.

Understanding the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.