Untangling Dual-Targeting Therapeutic Mechanism of Epidermal Growth Factor Receptor (EGFR) Based on Reversed Allosteric Communication

Mechanistic insights into the allosteric inactivation mechanism of ZAP-70 induced by the hot spot W165C mutation

Zeta chain-associated protein kinase 70 (ZAP-70) is a non-receptor tyrosine kinase that interacts with the activated T-cell receptor to transduce downstream signals, and thus plays an important role in the adaptive immune system. The biphosphorylated immunotyrosine-based activation motifs (ITAM-Y2P) binds to the N-SH2 and C-SH2 domains of ZAP-70 to promote the activation of ZAP-70. The present study explores molecular mechanisms of allosteric inactivation of ZAP-70 induced by the hot spot W165C mutation through atomically detailed molecular dynamics simulation approaches. We report microsecond-length simulations of two states of the tandem SH2 domains of ZAP-70 in complex with the ITAM-Y2P motif, including the wild-type and W165C mutant. Extensive analysis of local flexibility and dynamical correlated motions show that W165C mutation changes coupled motions of protein domains and community networks. The binding affinities of the ITAM-Y2P motif to the wild-type and W165C mutant of ZAP-70 are predicted using binding free energy calculations. The results suggest that the driving force to decrease the binding affinity in the W165C mutant derives from the difference in the protein-protein electrostatic interactions. Moreover, the per-residue free energy decomposition unravels that the contributions from residues in the phosphorylated Tyr315 (pY315) binding site, in particular pY315 of ITAM-Y2P, and Arg43, Tyr240 of ZAP-70, are the key determinants for the loss of binding affinity. This study may insights into our understanding of the pathological mechanism of ZAP-70.Communicated by Ramaswamy H. Sarma.

RGD-decorated nanoliposomes for combined delivery of arsenic trioxide and curcumin to prostate cancer cells

Nanotechnology and drug co-delivery offer a novel avenue in drug delivery research liposome-based co-delivery of anticancer drugs targeting the apoptosis pathway as a promising new approach to treat cancer. In this study, a co-delivery system of liposomes (arsenic trioxide/curcumin) modified with RGD peptide was designed to aim for enhancing the treatment of prostate cancer cells (PC3 cell line). Liposomal co-loaded curcumin and arsenic trioxide modified by RGD peptide (NLPs-RGD-Cur-ATO) were prepared by thin-layer lipid hydration techniques for the treatment of prostate cancer. The stability of the NLPs-RGD-Cur-ATO was evaluated by particle size analysis through dynamic light scattering (DLS) analysis and transmission electron microscopy (TEM). The percentage of cytotoxicity and apoptotic effect in PC3 cells treated with NLPs-RGD-Cur-ATO were detected by MTT and Annexin V-FITC (fluorescein isothiocyanate)/PI affinity assay, respectively. The particle size of NLPs-RGD-Cur-ATO was approximately 100 nm, with an encapsulation efficiency of about 99.52% and 70.61%, for ATO and Cur, respectively. Besides, NLPs-RGD-Cur-ATO displayed an enhanced anti-proliferative effect, increased the percentage of apoptotic cells 98 ± 1.85% (p < 0.0001), and significantly reduced EGFR gene expression level (p < 0.001) in the cell line tested. These results indicated that our NLPs-RGD-Cur-ATO co-delivery system was a promising strategy for prostate cancer therapy.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

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.

Decoupling the dynamic mechanism revealed by FGFR2 mutation-induced population shift

The fibroblast growth factor receptor 2 (FGFR2) is a key component in cellular signaling networks, and its dysfunctional activation has been implicated in various diseases including cancer and developmental disorders. Mutations at the activation loop (A-loop) have been suggested to trigger an increased basal kinase activity. However, the molecular mechanism underlying this highly dynamic process has not been fully understood due to the limitation of static structural information. Here, we conducted multiple, large-scale Gaussian accelerated molecular dynamics simulations of five (K659E, K659N, K659M, K659Q, and K659T) FGFR2 mutants at the A-loop, and comprehensively analyzed the dynamic molecular basis of FGFR2 activation. The results quantified the population shift of each system, revealing that all mutants had a higher proportion of active-like states. Using Markov state models, we extracted the representative structure of different conformational states and identified key residues related to the increased kinase activity. Furthermore, community network analysis showed enhanced information connections in the mutants, highlighting the long-range allosteric communication between the A-loop and the hinge region. Our findings may provide insights into the dynamic mechanism for FGFR2 dysfunctional activation and allosteric drug discovery.Communicated by Ramaswamy H. Sarma.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches

BACKGROUND

Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear.

RESULTS

This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin.

CONCLUSION

Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

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.

Decoding molecular mechanism of species-selective targeting of fungal versus human HSP90 using multiple replica molecular dynamics simulations and binding free energy calculations

As a highly evolutionarily conserved molecular chaperone, heat shock protein (HSP90), plays an important role in virulence traits, representing a therapeutic target for the treatment of fungal infections. The close evolutionary relationship between fungi and their human hosts poses a key challenge for the development of selective antifungal agents. In this work, molecular docking, multiple replica microsecond-based molecular dynamics (MD) simulations, and binding free energy calculations were performed to decode molecular mechanism of species-selective targeting of fungal versus human HSP90 triggered by the compound A11. MD simulations reveal that binding of compound A11 to human HSP90 nucleotide-binding domain (NBD) leads to obvious conformational changes relative to fungal HSP90 NBD. Binding free energy calculations show that the binding of compound A11 to fungal HSP90 NBD is stronger than that to human HSP90 NBD. Per residue-based free energy decomposition analysis was used to evaluate the inhibitor - residue interaction profile. The results efficiently identify the hot spot residues that play vital roles in favorable binding of compound A11 to fungal HSP90 NBD. This study is expected to provide a useful guidance for the development of selective inhibitors toward fungal HSP90.Communicated by Ramaswamy H. Sarma.

Mechanistic Elucidation of Activation/Deactivation Signal Transduction within Neurotensin Receptor 1 Triggered by 'Driver Chemical Groups' of Modulators: A Comparative Molecular Dynamics Simulation

Small-molecule modulators of neurotensin receptor 1 (NTSR1), a class A G-protein-coupled receptor (GPCR), has emerged as promising therapeutic agent for psychiatric disorders and cancer. Interestingly, a chemical group substitution in NTSR1 modulators can launch different types of downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted a synergistic application of a Gaussian accelerated molecular dynamics simulation, a conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underlying mechanism of 'driver chemical groups' of modulators triggering inverse signaling. The results indicated that the flexibility of the leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of the adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R322^6.54, Y319^6.51, F353^7.42, R148^3.32, S356^7.45, and S357^7.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs but also lay a foundation for the development of targeted drugs harboring different regulatory functions of NTSR1.

Deciphering the Shared and Specific Drug Resistance Mechanisms of Anaplastic Lymphoma Kinase via Binding Free Energy Computation

Anaplastic lymphoma kinase (ALK), a tyrosine receptor kinase, has been proven to be associated with the occurrence of numerous malignancies. Although there have been already at least 3 generations of ALK inhibitors approved by FDA or in clinical trials, the occurrence of various mutations seriously attenuates the effectiveness of the drugs. Unfortunately, most of the drug resistance mechanisms still remain obscure. Therefore, it is necessary to reveal the bottom reasons of the drug resistance mechanisms caused by the mutations. In this work, on the basis of verifying the accuracy of 2 main kinds of binding free energy calculation methodologies [end-point method of Molecular Mechanics with Poisson-Boltzmann/Generalized Born and Surface Area (MM/PB(GB)SA) and alchemical method of Thermodynamic Integration (TI)], we performed a systematic analysis on the ALK systems to explore the underlying shared and specific drug resistance mechanisms, covering the one-drug-multiple-mutation and multiple-drug-one-mutation cases. Through conventional molecular dynamics (cMD) simulation in conjunction with MM/PB(GB)SA and umbrella sampling (US) in conjunction with contact network analysis (CNA), the resistance mechanisms of the in-pocket, out-pocket, and multiple-site mutations were revealed. Especially for the out-pocket mutation, a possible transfer chain of the mutation effect was revealed, and the reason why different drugs exhibited various sensitivities to the same mutation was also uncovered. The proposed mechanisms may be prevalent in various drug resistance cases.

Toward the next generation EGFR inhibitors: an overview of osimertinib resistance mediated by EGFR mutations in non-small cell lung cancer

Epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI) is currently the standard first-line therapy for EGFR-mutated advanced non-small cell lung cancer (NSCLC). The life quality and survival of this subgroup of patients were constantly improving owing to the continuous iteration and optimization of EGFR-TKI. Osimertinib, an oral, third-generation, irreversible EGFR-TKI, was initially approved for the treatment of NSCLC patients carrying EGFR T790M mutations, and has currently become the dominant first-line targeted therapy for most EGFR mutant lung cancer. Unfortunately, resistance to osimertinib inevitably develops during the treatment and therefore limits its long-term effectiveness. For both fundamental and clinical researchers, it stands for a major challenge to reveal the mechanism, and a dire need to develop novel therapeutics to overcome the resistance. In this article, we focus on the acquired resistance to osimertinib caused by EGFR mutations which account for approximately 1/3 of all reported resistance mechanisms. We also review the proposed therapeutic strategies for each type of mutation conferring resistance to osimertinib and give an outlook to the development of the next generation EGFR inhibitors. Video Abstract.

The allosteric effect of the upper half of SENP1 contributes to its substrate selectivity for SUMO1 over SUMO2

SUMOylation regulates various cellular process and SENP1 (SUMO-specific protease 1) serves as a SUMO (small ubiquitin-related modifier) specific protease that participates in the SUMO cycle. Given its extensive influences on metabolic activities, SENP1 has gained more and more attentions in clinical treatments. However, there remains a question on why does the SENP1 prefer to process SUMO1 rather than SUMO2. Here, we performed molecular dynamics simulations of SENP1-SUMO1, SENP1-SUMO2, and apo SENP1 systems and observed distinct conformational dynamics in the upper half of the clamp and the three loops in the catalytic center of the SENP1. Principal component analysis revealed that the most prominent canonical variable represented the spatial distribution of the upper half of the clamp, while the openness of the cleft was closely related to the catalytic ability of SENP1. Further analysis of the SENP1-SUMO interactions revealed that the extensive and strong interactions between the SENP1 and SUMO1 were both in the interface of the upper half region and the catalytic center. Dynamic cross-correlation matrix analysis demonstrated that the inter-residue correlations in the SUMO1 system was much stronger, especially in the two essential regions belonging to the upper and lower half of cleft. Based on these observations, we proposed an allosteric propagation model and further testified it using the community analysis. These results revealed the propagation pathway of allosteric communication that contributed to the substrate discrimination of SENP1 upon SUMO1 and SUMO2.Communicated by Ramaswamy H. Sarma.

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 Understanding of the Palmitoylation of Go Protein in the Allosteric Regulation of Adhesion Receptor GPR97

Adhesion G-protein-coupled receptors (aGPCRs)—a major family of GPCRs—play critical roles in the regulation of tissue development and cancer progression. The orphan receptor GPR97, activated by glucocorticoid stress hormones, is a prototypical aGPCR. Although it has been established that the palmitoylation of the C-terminal G_o protein is essential for G_o’s efficient engagement with the active GPR97, the detailed allosteric mechanism remains to be clarified. Hence, we performed extensive large-scale molecular dynamics (MD) simulations of the GPR97−G_o complex in the presence or absence of G_o palmitoylation. The conformational landscapes analyzed by Markov state models revealed that the overall conformation of GPR97 is preferred to be fully active when interacting with palmitoylated G_o protein. Structural and energetic analyses indicated that the palmitoylation of G_o can allosterically stabilize the critical residues in the ligand-binding pocket of GPR97 and increase the affinity of the ligand for GPR97. Furthermore, the community network analysis suggests that the palmitoylation of G_o not only allosterically strengthens the internal interactions between G_αo and G_βγ, but also enhances the coupling between G_o and GPR97. Our study provides mechanistic insights into the regulation of aGPCRs via post-translational modifications of the G_o protein, and offers guidance for future drug design of aGPCRs.

Understanding gilteritinib resistance to FLT3-F691L mutation through an integrated computational strategy

FMS-like tyrosine kinase 3 (FLT3) serves as an important drug target for acute myeloid leukemia (AML), and gene mutations of FLT3 have been closely associated with AML patients with an incidence rate of ~ 30%. However, the mechanism of the clinically relevant F691L gatekeeper mutation conferred resistance to the drug gilteritinib remained poorly understood. In this study, multiple microsecond molecular dynamics (MD) simulations, end-point free energy calculations, and dynamic correlated and network analyses were performed to investigate the molecular basis of gilteritinib resistance to the FLT3-F691L mutation. The simulations revealed that the resistant mutation largely induced the conformational changes of the activation loop (A-loop), the phosphate-binding loop, and the helix αC of the FLT3 protein. The binding abilities of the gilteritinib to the wild-type and the F691L mutant were different through the binding free energy prediction. The simulation results further indicated that the driving force to determine the binding affinity of gilteritinib was derived from the differences in the energy terms of electrostatic and van der Waals interactions. Moreover, the per-residue free energy decomposition suggested that the four residues (Phe803, Gly831, Leu832, and Ala833) located at the A-loop of FLT3 had a significant impact on the binding affinity of gilteritinib to the F691L mutant. This study may provide useful information for the design of novel FLT3 inhibitors specially targeting the F691L gatekeeper mutant.

Mechanistic Insights into the Mechanism of Inhibitor Selectivity toward the Dark Kinase STK17B against Its High Homology STK17A

As a member of the death-associated protein kinase (DAPK) family, STK17B plays an important role in the regulation of cellular apoptosis and has been considered as a promising drug target for hepatocellular carcinoma. However, the highly conserved ATP-binding site of protein kinases represents a challenge to design selective inhibitors for a specific DAPK isoform. In this study, molecular docking, multiple large-scale molecular dynamics (MD) simulations, and binding free energy calculations were performed to decipher the molecular mechanism of the binding selectivity of PKIS43 toward STK17B against its high homology STK17A. MD simulations revealed that STK17A underwent a significant conformational arrangement of the activation loop compared to STK17B. The binding free energy predictions suggested that the driving force to control the binding selectivity of PKIS43 was derived from the difference in the protein–ligand electrostatic interactions. Furthermore, the per-residue free energy decomposition unveiled that the energy contribution from Arg41 at the phosphate-binding loop of STK17B was the determinant factor responsible for the binding specificity of PKIS43. This study may provide useful information for the rational design of novel and potent selective inhibitors toward STK17B.

Emerging strategies to overcome resistance to third-generation EGFR inhibitors

Epidermal growth factor receptor (EGFR), the receptor for members of the epidermal growth factor family, regulates cell proliferation and signal transduction; moreover, EGFR is related to the inhibition of tumor cell proliferation, angiogenesis, invasion, metastasis, and apoptosis. Therefore, EGFR has become an important target for the treatment of cancer, including non-small cell lung cancer, head and neck cancer, breast cancer, glioma, cervical cancer, and bladder cancer. First- to third-generation EGFR inhibitors have shown considerable efficacy and have significantly improved disease prognosis. However, most patients develop drug resistance after treatment. The challenge of overcoming intrinsic and acquired resistance in primary and recurrent cancer mediated by EGFR mutations is thus driving the search for alternative strategies in the design of new therapeutic agents. In view of resistance to third-generation inhibitors, understanding the intricate mechanisms of resistance will offer insight for the development of more advanced targeted therapies. In this review, we discuss the molecular mechanisms of resistance to third-generation EGFR inhibitors and review recent strategies for overcoming resistance, new challenges, and future development directions.

The Role of Conformational Dynamics and Allostery in the Control of Distinct Efficacies of Agonists to the Glucocorticoid Receptor

Glucocorticoid receptor (GR) regulates various cellular functions. Given its broad influence on metabolic activities, it has been the target of drug discovery for decades. However, how drugs induce conformational changes in GR has remained elusive. Herein, we used five GR agonists (dex, AZ938, pred, cor, and dibC) with different efficacies to investigate which aspect of the ligand induced the differences in efficacy. We performed molecular dynamics simulations on the five systems (dex-, AZ938-, pred-, cor-, and dibC-bound systems) and observed a distinct discrepancy in the conformation of the cofactor TIF2. Moreover, we discovered ligand-induced differences regarding the level of conformational changes posed by the binding of cofactor TIF2 and identified a pair of essential residues D590 and T39. We further found a positive correlation between the efficacies of ligands and the interaction of the two binding pockets' domains, where D590 and T739 were involved, implying their significance in the participation of allosteric communication. Using community network analysis, two essential communities containing D590 and T739 were identified with their connectivity correlating to the efficacy of ligands. The potential communication pathways between these two residues were revealed. These results revealed the underlying mechanism of allosteric communication between the ligand-binding and cofactor-binding pockets and identified a pair of important residues in the allosteric communication pathway, which can serve as a guide for future drug discovery.

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.

Deciphering the conformational transitions of LIMK2 active and inactive states to ponder specific druggable states through microsecond scale molecular dynamics simulation

LIMK2 inhibitors are one of the potential therapeutic modalities for treating various diseases. In the current scenario, there is a paucity of effective LIMK inhibitors that are highly specific with minimal off-target effects. To date, the conformational transitions of LIMK2 from DFG_inαC_in (CIDI) (active) to DFG_outαC_out (CODO) (inactive) states are yet to be probed and are essential for capturing the unique, druggable conformations. Therefore, this study was intended to capture the diverse conformational states of LIMK2 for accelerating the rational identification of conformation specific inhibitors through high-end structural bioinformatics protocols. Hence, in this study, molecular modelling followed by an extensive microsecond timescale of molecular dynamics simulation was performed encompassing perturbation response scanning, metapath, and community analysis towards the conformational sampling of LIMK2. Overall this study precisely identifies the conformational ensemble of LIMK2 the intermediate inactive states namely, CIDO, CinterDinter, CIDinter, CinterDI, CinterDO, CODI, CODinter apart from CIDI and CODO. This also facilitated observing that β8 preceding XDFG, αC (F373, L374), and αD (L413) as the major effectors that may facilitate the regulation of varying conformational transitions among the states. Additionally, the conserved β sheets and the loops namely, C.l, b.l, and G/P.loop were observed to be involved in the metapath for allosteric communication among the intermediates with CIDI and CODO state. Moreover, only the CODO state was observed to have closed type A.l, while the CIDI and other intermediate states except for CIDO were observed to have open-DFG out type A.l, thereby enabling the binding of substrate. Apart from these, the druggable site analysis inferred that the CIDI and CODO states harbor prominent druggable sites spanning the conserved N-lobe, while the intermediates were observed to have unraveled allosteric druggable sites distal from the ATP binding site, majorly spanning the C-lobe of LIMK2. Thus, this study provides potential insights into the intermediate conformational druggable states of LIMK2 and also the druggable conformations, especially the inactive states of LIMK2, as a specific therapeutic targeting mode. Thus, providing a widened avenue to ponder the allosteric sites or the isoform selectivity conformations for targeting LIMK2 in various disease conditions.

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.

Allosteric pluripotency: challenges and opportunities

Allosteric pluripotency arises when the functional response of an allosteric receptor to an allosteric stimulus depends on additional allosteric modulators. Here, we discuss allosteric pluripotency as observed in the prototypical Protein Kinase A (PKA) as well as in other signaling systems, from typical multidomain signaling proteins to bacterial enzymes. We identify key drivers of pluripotent allostery and illustrate how hypothesizing allosteric pluripotency may solve apparent discrepancies currently present in the literature regarding the dual nature of known allosteric modulators. We also outline the implications of allosteric pluripotency for cellular signaling and allosteric drug design, and analyze the challenges and opportunities opened by the pluripotent nature of allostery.

Autopromotion of K-Ras4B Feedback Activation Through an SOS-Mediated Long-Range Allosteric Effect

The Ras-specific guanine nucleotide exchange factors Son of Sevenless (SOS) regulates Ras activation by converting inactive GDP-bound to active GTP-bound states. The catalytic activity of Ras is further allosterically regulated by GTP-Ras bound to a distal site through a positive feedback loop. To address the mechanism underlying the long-range allosteric activation of the catalytic K-Ras4B by an additional allosteric GTP-Ras through SOS, we employed molecular dynamics simulation of the K-Ras4BG13D•SOScat complex with and without an allosteric GTP-bound K-Ras4BG13D. We found that the binding of an allosteric GTP-K-Ras4BG13D enhanced the affinity between the catalytic K-Ras4BG13D and SOScat, forming a more stable conformational state. The peeling away of the switch I from the nucleotide binding site facilitated the dissociation of GDP, thereby contributing to the increased nucleotide exchange rate. The community networks further showed stronger edge connection upon allosteric GTP-K-Ras4BG13D binding, which represented an increased interaction between catalytic K-Ras4BG13D and SOScat. Moreover, GTP-K-Ras4BG13D binding transmitted allosteric signaling pathways though the Cdc25 domain of SOS that enhanced the allosteric regulatory from the K-Ras4BG13D allosteric site to the catalytic site. This study may provide an in-depth mechanism for abnormal activation and allosteric regulation of K-Ras4BG13D.

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

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm, driven by the BCR-ABL1 fusion oncoprotein. The discovery of orthosteric BCR-ABL1 tyrosine kinase inhibitors (TKIs) targeting its active ATP-binding pocket, such as first-generation Imatinib and second-generation Nilotinib (NIL), has profoundly revolutionized the therapeutic landscape of CML. However, currently targeted therapeutics still face considerable challenges with the inevitable emergence of drug-resistant mutations within BCR-ABL1. One of the most common resistant mutations in BCR-ABL1 is the T315I gatekeeper mutation, which confers resistance to most current TKIs in use. To resolve such conundrum, co-administration of orthosteric TKIs and allosteric drugs offers a novel paradigm to tackle drug resistance. Remarkably, previous studies have confirmed that the dual targeting BCR-ABL1 utilizing orthosteric TKI NIL and allosteric inhibitor ABL001 resulted in eradication of the CML xenograft tumors, exhibiting promising therapeutic potential. Previous studies have demonstrated the cooperated mechanism of two drugs. However, the conformational landscapes of synergistic effects remain unclear, hampering future efforts in optimizations and improvements. Hence, extensive large-scale molecular dynamics (MD) simulations of wide type (WT), WT-NIL, T315I, T315I-NIL, T315I-ABL001 and T315I-ABL001-NIL systems were carried out in an attempt to address such question. Simulation data revealed that the dynamic landscape of NIL-bound BCR-ABL1 was significantly reshaped upon ABL001 binding, as it shifted from an active conformation towards an inactive conformation. The community network of allosteric signaling was analyzed to elucidate the atomistic overview of allosteric regulation within BCR-ABL1. Moreover, binding free energy analysis unveiled that the affinity of NIL to BCR-ABL1 increased by the induction of ABL001, which led to its favorable binding and the release of drug resistance. The findings uncovered the in-depth structural mechanisms underpinning dual-targeting towards T315I BCR-ABL1 to overcome its drug resistance and will offer guidance for the rational design of next generations of BCR-ABL1 modulators and future combinatory therapeutic regimens.

Elucidation of the conformational dynamics and assembly of Argonaute-RNA complexes by distinct yet coordinated actions of the supplementary microRNA

Argonaute (AGO) proteins, the core of RNA-induced silencing complex, are guided by microRNAs (miRNAs) to recognize target RNA for repression. The miRNA-target RNA recognition forms initially through pairing at the seed region while the additional supplementary pairing can enhance target recognition and compensate for seed mismatch. The extension of miRNA lengths can strengthen the target affinity when pairing both in the seed and supplementary regions. However, the mechanism underlying the effect of the supplementary pairing on the conformational dynamics and the assembly of AGO-RNA complex remains poorly understood. To address this, we performed large-scale molecular dynamics simulations of AGO-RNA complexes with different pairing patterns and miRNA lengths. The results reveal that the additional supplementary pairing can not only strengthen the interaction between miRNA and target RNA, but also induce the increased plasticity of the PAZ domain and enhance the domain connectivity among the PAZ, PIWI, N domains of the AGO protein. The strong community network between these domains tightens the mouth of the supplementary chamber of AGO protein, which prevents the escape of target RNA from the complex and shields it from solvent water attack. Importantly, the inner stronger matching pairs between the miRNA and target RNA can compensate for weaker mismatches at the edge of supplementary region. These findings provide guidance for the design of miRNA mimics and anti-miRNAs for both clinical and experimental use and open the way for further engineering of AGO proteins as a new tool in the field of gene regulation.

Outliers in SAR and QSAR: 4. effects of allosteric protein-ligand interactions on the classical quantitative structure-activity relationships

Effects of allosteric interactions on the classical structure-activity relationship (SAR) and quantitative SAR (QSAR) have been investigated. Apprehending the outliers in SAR and QSAR studies can improve the quality, predictability, and use of QSAR in designing unknown compounds in drug discovery research. We explored allosteric protein-ligand interactions as a possible source of outliers in SAR/QSAR. We used glycogen phosphorylase as an example of a protein that has an allosteric site. Examination of the ligand-bound x-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site. The case provided an example of constructive use of QSAR identifying outliers with alternative binding modes. Other allosteric QSARs that captured our attention were the inverted parabola/bilinear QSARs. The x-ray crystal structures and the QSAR analyses indicated that the inverted parabola QSARs could be associated with the conformational changes in the allosteric interactions. Our results showed that the normal parabola, as well as the inverted parabola QSARs, can describe the allosteric interactions. Examination of the ligand-bound X-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site.

Targeting a cryptic allosteric site of SIRT6 with small-molecule inhibitors that inhibit the migration of pancreatic cancer cells

SIRT6 belongs to the conserved NAD⁺-dependent deacetylase superfamily and mediates multiple biological and pathological processes. Targeting SIRT6 by allosteric modulators represents a novel direction for therapeutics, which can overcome the selectivity problem caused by the structural similarity of orthosteric sites among deacetylases. Here, developing a reversed allosteric strategy AlloReverse, we identified a cryptic allosteric site, Pocket Z, which was only induced by the bi-directional allosteric signal triggered upon orthosteric binding of NAD⁺. Based on Pocket Z, we discovered an SIRT6 allosteric inhibitor named JYQ-42. JYQ-42 selectively targets SIRT6 among other histone deacetylases and effectively inhibits SIRT6 deacetylation, with an IC₅₀ of 2.33 μmol/L. JYQ-42 significantly suppresses SIRT6-mediated cancer cell migration and pro-inflammatory cytokine production. JYQ-42, to our knowledge, is the most potent and selective allosteric SIRT6 inhibitor. This study provides a novel strategy for allosteric drug design and will help in the challenging development of therapeutic agents that can selectively bind SIRT6.

Delineating the activation mechanism and conformational landscape of a class B G protein-coupled receptor glucagon receptor

Class B G protein-coupled receptors (GPCRs) are important targets in the treatment of metabolic syndrome and diabetes. Although multiple structures of class B GPCRs-G protein complexes have been elucidated, the detailed activation mechanism of the receptors remains unclear. Here, we combine Gaussian accelerated molecular dynamics simulations and Markov state models (MSM) to investigate the activation mechanism of a canonical class B GPCR, human glucagon receptor-GCGR, including the negative allosteric modulator-bound inactive state, the agonist glucagon-bound active state, and both glucagon- and Gs-bound fully active state. The free-energy landscapes of GCGR show the conformational ensemble consisting of three activation-associated states: inactive, active, and fully active. The structural analysis indicates the high dynamics of GCGR upon glucagon binding with both active and inactive conformations in the ensemble. Significantly, the H8 and TM6 exhibits distinct features from the inactive to the active states. The additional simulations demonstrate the role of H8 in the recruitment of Gs. Gs binding presents a crucial function of stabilizing the glucagon binding site and MSM highlights the absolute requirement of Gs to help the GCGR reach the fully active state. Together, our results reveal the detailed activation mechanism of GCGR from the view of conformational dynamics.

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.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Hidden allosteric sites and De-Novo drug design

INTRODUCTION

Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge.

AREAS COVERED

This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator.

EXPERT OPINION

In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.

Atomic-scale insights into allosteric inhibition and evolutional rescue mechanism of Streptococcus thermophilus Cas9 by the anti-CRISPR protein AcrIIA6

CRISPR-Cas systems are prokaryotic adaptive immunity against invading phages and plasmids. Phages have evolved diverse protein inhibitors of CRISPR-Cas systems, called anti-CRISPR (Acr) proteins, to neutralize this CRISPR machinery. In response, bacteria have co-evolved Cas variants to escape phage's anti-CRISPR strategies, called anti-anti-CRISPR systems. Here we explore the anti-CRISPR allosteric inhibition and anti-anti-CRISPR rescue mechanisms between Streptococcus thermophilus Cas9 (St1Cas9) and the anti-CRISPR protein AcrIIA6 at the atomic level, by generating mutants of key residues in St1Cas9. Extensive unbiased molecular dynamics simulations show that the functional motions of St1Cas9 in the presence of AcrIIA6 differ substantially from those of St1Cas9 alone. AcrIIA6 binding triggers a shift of St1Cas9 conformational ensemble towards a less catalytically competent state; this state significantly compromises protospacer adjacent motif (PAM) recognition and nuclease activity by altering interdependently conformational dynamics and allosteric signals among nuclease domains, PAM-interacting (PI) regions, and AcrIIA6 binding motifs. Via in vitro DNA cleavage assays, we further elucidate the rescue mechanism of efficiently escaping AcrIIA6 inhibition harboring St1Cas9 triple mutations (G993K/K1008M/K1010E) in the PI domain and identify the evolutionary landscape of such mutational escape within species. Our results provide mechanistic insights into Acr proteins as natural brakes for the CRISPR-Cas systems and a promising potential for the design of allosteric Acr peptidomimetics.

How does nintedanib overcome cancer drug-resistant mutation of RET protein-tyrosine kinase: insights from molecular dynamics simulations

Targeted drug therapies represent a therapeutic breakthrough in the treatment of human cancer. However, the emergence of acquired resistance inevitably compromises therapeutic drugs. Rearranged during transfection (RET) proto-oncogene, which encodes a receptor tyrosine kinase, is a target for several kinds of human cancer such as thyroid, breast, and colorectal carcinoma. A single mutation L881V at the RET kinase domain was found in familial medullary thyroid carcinoma. Nintedanib can effectively inhibit the RET L881V mutant, whereas its analog compound 1 is unable to combat this mutant. However, the underlying mechanism was still unexplored. Here, molecular dynamics (MD) simulations, binding free energy calculations, and structural analysis were performed to uncover the mechanism of overcoming the resistance of RET L881V mutant to nintedanib. Energetic analysis revealed that the L881V mutant remained sensitive to the treatment of nintedanib, whereas it was insensitive to the compound 1. Structural analysis further showed that the distribution of K758, D892, and N879 network had a detrimental effect on the binding of compound 1 to the L881V mutant. The obtained results may provide insight into the mechanism of overcoming resistance in the RET kinase.