Insight into the mechanism of allosteric activation of PI3Kα by oncoprotein K-Ras4B

Structural and energetic insights into the selective inhibition of PKMYT1 against WEE1

Protein kinase, membrane-associated tyrosine/threonine 1 (PKMYT1), a member of the WEE family and responsible for the regulation of CDK1 phosphorylation, has been considered a promising therapeutic target for cancer therapy. However, the highly structural conservation of the ATP-binding sites of the WEE family poses a challenge to the design of selective inhibitors for PKMYT1. Here, molecular docking, multiple microsecond-length molecular dynamics (MD) simulations and end-point free energy calculations were performed to uncover the molecular mechanism of the binding selectivity of RP-6306 toward PKMYT1 over its highly homologous kinase WEE1. The binding specificity of RP-6306 reported in previous experimental bioassays was clarified by MD simulations and binding free energy calculations. Further, the binding free energy prediction indicated that the binding selectivity of RP-6306 largely derived from the difference in the protein-ligand electrostatic interactions. The per-residue free energy decomposition suggested that the non-conserved gatekeeper residue in the hinge domain of PKMYT1/WEE1, Thr187/Asn376, is the critical factor responsible for the binding selectivity of RP-6306 toward PKMYT1. In addition, a water-mediated hydrogen bond was formed between RP-6306 and Gly191 at the hinge domain in the PKMYT1/RP-6306 complex, which was absent in the WEE1/RP-6306 complex. This study is expected to offer useful information for the design of more potent and selective PKMYT1 inhibitors.Communicated by Ramaswamy H. Sarma.

Selective inhibitors targeting Fis1/Mid51 protein-protein interactions protect against hypoxia-induced damage in cardiomyocytes

Cardiovascular diseases (CVDs) are the most common non-communicable diseases globally. An estimated 17.9 million people died from CVDs in 2019, representing 32% of all global deaths. Mitochondria play critical roles in cellular metabolic homeostasis, cell survival, and cell death, as well as producing most of the cell's energy. Protein-protein interactions (PPIs) have a significant role in physiological and pathological processes, and aberrant PPIs are associated with various diseases, therefore they are potential drug targets for a broad range of therapeutic areas. Due to their ability to mimic natural interaction motifs and cover relatively larger interaction region, peptides are very promising as PPI inhibitors. To expedite drug discovery, computational approaches are widely used for screening potential lead compounds. Here, we developed peptides that inhibit mitochondrial fission 1 (Fis1)/mitochondrial dynamics 51 kDa (Mid51) PPI to reduce the cellular damage that can lead to various human pathologies, such as CVDs. Based on a rational design approach we developed peptide inhibitors of the Fis1/Mid51 PPI. In silico and in vitro studies were done to evaluate the biological activity and molecular interactions of the peptides. Two peptides, CVP-241 and CVP-242 were identified based on low binding energy and molecular dynamics simulations. These peptides inhibit Fis1/Mid51 PPI (-1324.9 kcal mol-1) in docking calculations (CVP-241, -741.3 kcal mol-1, and CVP-242, -747.4 kcal mol-1), as well as in vitro experimental studies Fis1/Mid51 PPI (KD 0.054 µM) Fis1/Mid51 PPI + CVP-241 (KD 3.43 µM), and Fis1/Mid51 PPI + CVP-242 (KD 44.58 µM). Finally, these peptides have no toxicity to H9c2 cells, and they increase cell viability in cardiomyocytes (H9c2 cells). Consequently, the identified inhibitor peptides could serve as potent molecules in basic research and as leads for therapeutic development.

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 insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase

Calcium (Ca²⁺) signaling plays an important role in the regulation of many cellular functions. Ca²⁺-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca²⁺/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca²⁺ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca²⁺-bound and-unbound states to reveal the importance of Ca²⁺. MD simulations revealed that removal of Ca²⁺ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca²⁺-unbound state. Structural analysis further revealed that Ca²⁺ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca²⁺ in physiological processes.

In Silico Study of the Acquired Resistance Caused by the Secondary Mutations of KRAS G12C Protein Using Long Time Molecular Dynamics Simulation and Markov State Model Analysis

Kirsten rat sarcoma viral oncogene homolog (KRAS) is a small GTPase protein which plays an important role in the treatment of KRAS mutant cancers. The FDA-approved AMG510 and MRTX849 (phase III clinical trials) are two potent KRASG12C-selective inhibitors that target KRAS G12C. However, the drug resistance caused by the second-site mutation in KRAS has emerged, and the mechanisms of drug resistance at atom level are still unclear. To clarify the mechanisms of drug resistance, we conducted long time molecular dynamics simulations (75 μs in total) to study the structural and energetic features of KRAS G12C and its four drug resistant variants to inhibitors. The combined binding free energy calculation and protein-ligand interaction fingerprint revealed that these second-site mutations indeed caused KRAS to produce different degrees of resistance to AMG510 and MRTX849. Furthermore, Markov State Models and 2D-free energy landscapes analysis revealed the difference in conformational changes of mutated KRAS bound with and without inhibitors. Furthermore, the comparative analysis of these systems showed that there were differences in their allosteric signal pathways. These findings provide the molecular mechanism of drug resistance, which helps to guide novel KRAS G12C inhibitor design to overcome drug resistance.

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.

Computer-Aided Drug Design Boosts RAS Inhibitor Discovery

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.

Mechanism of zinc ejection by disulfiram in nonstructural protein 5A

Hepatitis C virus (HCV) is a notorious member of the Flaviviridae family of enveloped, positive-strand RNA viruses. Non-structural protein 5A (NS5A) plays a key role in HCV replication and assembly. NS5A is a multi-domain protein which includes an N-terminal amphipathic membrane anchoring alpha helix, a highly structured domain-1, and two intrinsically disordered domains 2-3. The highly structured domain-1 contains a zinc finger (Zf)-site, and binding of zinc stabilizes the overall structure, while ejection of this zinc from the Zf-site destabilizes the overall structure. Therefore, NS5A is an attractive target for anti-HCV therapy by disulfiram, through ejection of zinc from the Zf-site. However, the zinc ejection mechanism is poorly understood. To disclose this mechanism based on three different states, A-state (NS5A protein), B-state (NS5A + Zn), and C-state (NS5A + Zn + disulfiram), we have performed molecular dynamics (MD) simulation in tandem with DFT calculations in the current study. The MD results indicate that disulfiram triggers Zn ejection from the Zf-site predominantly through altering the overall conformation ensemble. On the other hand, the DFT assessment demonstrates that the Zn adopts a tetrahedral configuration at the Zf-site with four Cys residues, which indicates a stable protein structure morphology. Disulfiram binding induces major conformational changes at the Zf-site, introduces new interactions of Cys39 with disulfiram, and further weakens the interaction of this residue with Zn, causing ejection of zinc from the Zf-site. The proposed mechanism elucidates the therapeutic potential of disulfiram and offers theoretical guidance for the advancement of drug candidates.

The Importance of Being PI3K in the RAS Signaling Network

Ras proteins are essential mediators of a multitude of cellular processes, and its deregulation is frequently associated with cancer appearance, progression, and metastasis. Ras-driven cancers are usually aggressive and difficult to treat. Although the recent Food and Drug Administration (FDA) approval of the first Ras G12C inhibitor is an important milestone, only a small percentage of patients will benefit from it. A better understanding of the context in which Ras operates in different tumor types and the outcomes mediated by each effector pathway may help to identify additional strategies and targets to treat Ras-driven tumors. Evidence emerging in recent years suggests that both oncogenic Ras signaling in tumor cells and non-oncogenic Ras signaling in stromal cells play an essential role in cancer. PI3K is one of the main Ras effectors, regulating important cellular processes such as cell viability or resistance to therapy or angiogenesis upon oncogenic Ras activation. In this review, we will summarize recent advances in the understanding of Ras-dependent activation of PI3K both in physiological conditions and cancer, with a focus on how this signaling pathway contributes to the formation of a tumor stroma that promotes tumor cell proliferation, migration, and spread.

Effects of miR-202-5p silencing PIK3CA gene expression on proliferation, invasion, and epithelial-mesenchymal transition of cervical cancer SiHa cells through inhibiting PI3K/Akt/mTOR signaling pathway activation

To explore the mechanism of miR-202-5p targeting the expression of PIK3CA and mediating the activation of PI3K/Akt/mTOR signaling pathway on the proliferation, invasion, and epithelial-mesenchymal transition (EMT) of cervical cancer. The objects of study were 105 cases of cervical cancer and their corresponding normal tissues. qRT-PCR was used to detect the expression of miR-202-5p and PIK3CA in adjacent normal tissue and cervical cancer tissue. Dual luciferase reporter assay was used to verify the targeting relationship between miR-202-5p and PIK3CA gene. Human cervical cancer cell lines HPV-16E6, SiHa, HeLa, and CaSki were purchased for our cell experiments. The expression levels of PIK3CA in the cells were detected by qRT-PCR. The cell line with higher expression levels was selected to complete the follow-up experiment. The cultured cells were transfected and divided into the miR-202-5p mimic NC group, miR-202-5p mimic group, miR-202-5p inhibitor NC group, miR-202-5p inhibitor group, siRNA-PIK3CA NC group, siRNA-PIK3CA group, miR-202-5p inhibitor NC + siRNA-PIK3CA NC group, miR-202-5p inhibitor + siRNA-PIK3CA NC group, and miR-202-5p inhibitor + siRNA-PIK3CA group. QRT-PCR was used to detect the expression of miR-202-5p. Western blot and qRT-PCR were applied to detect the mRNA and protein expression levels of related pathway proteins (PIK3CA, PI3K, PTEN, p-Akt1, and p-mTOR) and epithelial-mesenchymal transition-related factors (N-cadherin, E-cadherin, and vimentin). Cell proliferation was detected by plate colony formation assay. Transwell assay was used to detect the invasion ability of each group. When compared with the adjacent tissues, PIK3CA mRNA expression level was significantly increased and miR-202-5p expression level was significantly decreased in cervical cancer tissues (all P < 0.05). PIK3CA was a target gene of miR-202-5p. The mRNA expression level of PIK3CA in SiHa cervical cancer cells was significantly higher than that in CaSki, HeLa, and HPV-16E6 cells (all P < 0.05), and SiHa cervical cancer cells were selected to complete the follow-up experiments. When compared with the corresponding NC group, the expression of miR-202-5p in miR-202-5p mimic group was increased. In addition, the mRNA and protein expression levels of E-cadherin and PTEN in miR-202-5p mimic and siRNA-PIK3CA groups were increased, and the protein expression of p-Akt1 and p-mTOR was decreased, and also, the mRNA and protein expression levels of PIK3CA, PI3K, N-cadherin, and vimentin were decreased (all P < 0.05); in miR-202-5p inhibitor group, the expression levels of miR-202-5p, E-cadherin, and PTEN decreased, the protein expression of p-Akt1 and p-mTOR increased, and the mRNA and protein expression of PIK3CA, PI3K, N-cadherin, and vimentin increased in miR-202-5p inhibitor group (all P < 0.05); in miR-202-5p inhibitor + siRNA-PIK3CA group, the expression of miR-202-5p decreased (P < 0.05), but the mRNA and protein expression of PIK3CA, PI3K, p-Akt1, p-mTOR, N-cadherin, E-cadherin, and vimentin had no significant changes (all P > 0.05). When compared with the corresponding NC group, the number of cell clones in miR-202-5p mimic group and siRNA-PIK3CA group was decreased, and the invasion ability of miR-202-5p inhibitor group was increased, and the invasion ability was enhanced (all P < 0.05); miR-202-5p inhibitor + siRNA-PIK3CA group showed no significant change in the number of cell clones and the rate of invasion (P > 0.05). In conclusion, the overexpression of miR-202-5p can suppress PIK3CA gene expression and the activation of PI3K/Akt/mTOR signaling pathway to suppress the proliferation, invasion, and EMT of cervical cancer.

Rational Design of Peptide-Based Inhibitors Disrupting Protein-Protein Interactions

Protein-protein interactions (PPIs) are well-established as a class of promising drug targets for their implications in a wide range of biological processes. However, drug development toward PPIs is inevitably hampered by their flat and wide interfaces, which generally lack suitable pockets for ligand binding, rendering most PPI systems "undruggable." Here, we summarized drug design strategies for developing peptide-based PPI inhibitors. Importantly, several quintessential examples toward well-established PPI targets such as Bcl-2 family members, p53-MDM2, as well as APC-Asef are presented to illustrate the detailed schemes for peptide-based PPI inhibitor development and optimizations. This review supplies a comprehensive overview of recent progresses in drug discovery targeting PPIs through peptides or peptidomimetics, and will shed light on future therapeutic agent development toward the historically "intractable" PPI systems.

Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators

The recent discovery of activator compounds binding to an allosteric site on the NAD⁺-dependent protein lysine deacetylase, sirtuin 6 (SIRT6) has attracted interest and presents a pharmaceutical target for aging-related and cancer diseases. However, the mechanism underlying allosteric activation of SIRT6 by the activator MDL-801 remains largely elusive because no major conformational changes are observed upon activator binding. By combining molecular dynamics simulations with biochemical and kinetic analyses of wild-type SIRT6 and its variant M136A, we show that conformational rotation of 2-methyl-4-fluoro-5-bromo substituent on the right phenyl ring (R-ring) of MDL-801, which uncovers previously unseen hydrophobic interactions, contributes to increased activating deacetylation activity of SIRT6. This hypothesis is further supported by the two newly synthesized MDL-801 derivatives through the removal of the 5-Br atom on the R-ring (MDL-801-D1) or the restraint of the rotation of the R-ring (MDL-801-D2). We further propose that the 5-Br atom serves as an allosteric driver that controls the ligand allosteric efficacy. Our study highlights the effect of allosteric enzyme catalytic activity by activator binding and provides a rational approach for enhancing deacetylation activity.

Deciphering the selective binding mechanisms of anaplastic lymphoma kinase-derived neuroblastoma tumor neoepitopes to human leukocyte antigen

Neuroblastoma (NB), as a metastatic form of solid tumor, has a high fatality rate found in early childhood. The two anaplastic lymphoma kinase (ALK) neoepitopes nonamer and decamer used in cancer immunotherapy against NB cancer can selectively bind to the human leukocyte antigen (HLA-B*15:01) groove with high affinities, whereas the native self-peptide is unable to interact with the HLA-B*15:01. Here, we performed molecular dynamics (MD) simulations and subsequent molecular mechanics-generalized born surface area (MM-GBSA) binding free energy calculations to explore the selective binding mechanisms of nonamer and decamer to the HLA-B*15:01 against the self-peptide. MD simulations revealed the significant conformational dynamics of the self-peptide in the HLA-B*15:01 groove against the nonamer and decamer. Binding free energy calculations showed that the binding affinities of HLA-B*15:01-neoepitope complexes were followed in the order decamer > nonamer > self-peptide. Detailed analysis of HLA-B*15:01-neoepitope structural complexes showed that compared to the nonamer, the self-peptide tended to move outward to the solvent, whereas the decamer moved deep to the HLA-B*15:01 groove. These different dynamic observations of the ALK neoepitopes can explain the distinct binding affinities of self-peptide, nonamer, and decamer to the HLA-B*15:01. The results may be useful for the design of more selective ALK neoepitopes.

Conformational Selection Mechanism Provides Structural Insights into the Optimization of APC-Asef Inhibitors

Metastasis is the major cause of death in colorectal cancer and it has been proven that inhibiting an interaction between adenomatous polyposis coli (APC) and Rho guanine nucleotide exchange factor 4 (Asef) efficaciously restrain metastasis. However, current inhibitors cannot achieve a satisfying effect in vivo and need to be optimized. In the present study, we applied molecular dynamics (MD) simulations and extensive analyses to apo and holo APC systems in order to reveal the inhibitor mechanism in detail and provide insights into optimization. MD simulations suggested that apo APC takes on a broad array of conformations and inhibitors stabilize conformation selectively. Representative structures in trajectories show specific APC-ligand interactions, explaining the different binding process. The stability and dynamic properties of systems elucidate the inherent factors of the conformation selection mechanism. Binding free energy analysis quantitatively confirms key interface residues and guide optimization. This study elucidates the conformation selection mechanism in APC-Asef inhibition and provides insights into peptide-based drug design.

Mechanistic insights into the effect of phosphorylation on Ras conformational dynamics and its interactions with cell signaling proteins

Ras undergoes interconversion between the active GTP-bound state and the inactive GDP-bound state. This GTPase cycle, which controls the activities of Ras, is accelerated by Ras GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (SOS). Oncogenic Ras mutations could affect the GTPase cycle and impair Ras functions. Additionally, Src-induced K-Ras Y32/64 dual phosphorylation has been reported to disrupt GTPase cycle and hinder Ras downstream signaling. However, the underlying mechanisms remain unclear. To address this, we performed molecular dynamics simulations (~30 μs in total) on unphosphorylated and phosphorylated K-Ras4B in GTP- and GDP-bound states, and on their complexes with GTPase cycle regulators (GAP and SOS) and the effector protein Raf. We found that K-Ras4B dual phosphorylation mainly alters the conformation at the nucleotide binding site and creates disorder at the catalytic site, resulting in the enlargement of GDP binding pocket and the retard of Ras-GTP intrinsic hydrolysis. We observed phosphorylation-induced shift in the distribution of Ras-GTP inactive-active sub-states and recognized potential druggable pockets in the phosphorylated Ras-GTP. Moreover, decreased catalytic competence or signal delivery abilities due to reduced binding affinities and/or distorted catalytic conformations of GAP, SOS and Raf were observed. In addition, the allosteric pathway from Ras/Raf interface to the distal Raf L4 loop was compromised by Ras phosphorylation. These results reveal the mechanisms by which phosphorylation influences the intrinsic or GAP/SOS catalyzed transformations between GTP- and GDP-bound states of Ras and its signal transduction to Raf. Our findings project Ras phosphorylation as a target for cancer drug discovery.

How Parkinson's disease-related mutations disrupt the dimerization of WD40 domain in LRRK2: a comparative molecular dynamics simulation study

The multidomain kinase enzyme leucine-rich-repeat kinase 2 (LRRK2), activated through a homodimerization manner, has been identified as an important pathogenic factor in Parkinson's disease (PD), the second most common neurodegenerative disease wordwide. The Trp-Asp-40 (WD40) domain, located in the C-terminal LRRK2, harbours one of the most frequent PD-related variants, G2385R. However, the detailed dynamics of WD40 during LRRK2 dimerization and the underlying mechanism through which the pathogenic mutations disrupt the formation of the WD40 dimer have remained elusive. Here, microsecond-scale molecular dynamics simulations were employed to provide a mechanistic view underlying the WD40 dimerization and unveil the structural basis by which the interface-based mutations G2385R, H2391D and R2394E compromise the corresponding process. The simulation results identified important residues, D2351, R2394, E2395, R2413, and R2443, involved in establishing the complex binding network along the dimerization interface, which was significantly weakened in the presence of interfacial mutations. A "sag-bulge" model was proposed to explain the unfavorable dimer formation in the mutant systems. In addition, mutations altered the community configuration in the wild-type system in which inter-monomeric interplay is prominent, leading to the destabilization of the WD40 dimer under mutation.

PI3K Driver Mutations: A Biophysical Membrane-Centric Perspective

Ras activates its effectors at the membrane. Active PI3Kα and its associated kinases/phosphatases assemble at membrane regions enriched in signaling lipids. In contrast, the Raf kinase domain extends into the cytoplasm and its assembly is away from the crowded membrane surface. Our structural membrane-centric outlook underscores the spatiotemporal principles of membrane and signaling lipids, which helps clarify PI3Kα activation. Here we focus on mechanisms of activation driven by PI3Kα driver mutations, spotlighting the PI3Kα double (multiple) activating mutations. Single mutations can be potent, but double mutations are stronger: their combination is specific, a single strong driver cannot fully activate PI3K, and two weak drivers may or may not do so. In contrast, two strong drivers may successfully activate PI3K, where one, for example, H1047R, modulates membrane interactions facilitating substrate binding at the active site (km) and the other, for example, E542K and E545K, reduces the transition state barrier (ka), releasing autoinhibition by nSH2. Although mostly unidentified, weak drivers are expected to be common, so we ask here how common double mutations are likely to be and why PI3Kα with double mutations responds effectively to inhibitors. We provide a structural view of hotspot and weak driver mutations in PI3Kα activation, explain their mechanisms, compare these with mechanisms of Raf activation, and point to targeting cell-specific, chromatin-accessible, and parallel (or redundant) pathways to thwart the expected emergence of drug resistance. Collectively, our biophysical outlook delineates activation and highlights the challenges of drug resistance.

Computational methods-guided design of modulators targeting protein-protein interactions (PPIs)

Protein-protein interactions (PPIs) play a pivotal role in extensive biological processes and are thus crucial to human health and the development of disease states. Due to their critical implications, PPIs have been spotlighted as promising drug targets of broad-spectrum therapeutic interests. However, owing to the general properties of PPIs, such as flat surfaces, featureless conformations, difficult topologies, and shallow pockets, previous attempts were faced with serious obstacles when targeting PPIs and almost portrayed them as "intractable" for decades. To date, rapid progress in computational chemistry and structural biology methods has promoted the exploitation of PPIs in drug discovery. These techniques boost their cost-effective and high-throughput traits, and enable the study of dynamic PPI interfaces. Thus, computational methods represent an alternative strategy to target "undruggable" PPI interfaces and have attracted intense pharmaceutical interest in recent years, as exemplified by the accumulating number of successful cases. In this review, we first introduce a diverse set of computational methods used to design PPI modulators. Herein, we focus on the recent progress in computational strategies and provide a comprehensive overview covering various methodologies. Importantly, a list of recently-reported successful examples is highlighted to verify the feasibility of these computational approaches. Finally, we conclude the general role of computational methods in targeting PPIs, and also discuss future perspectives on the development of such aids.