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

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.

Markov State Models and Molecular Dynamics Simulations Provide Understanding of the Nucleotide-Dependent Dimerization-Based Activation of LRRK2 ROC Domain

Mutations in leucine-rich repeat kinase 2 (LRRK2) are recognized as the most frequent cause of Parkinson's disease (PD). As a multidomain ROCO protein, LRRK2 is characterized by the presence of both a Ras-of-complex (ROC) GTPase domain and a kinase domain connected through the C-terminal of an ROC domain (COR). The bienzymatic ROC-COR-kinase catalytic triad indicated the potential role of GTPase domain in regulating kinase activity. However, as a functional GTPase, the detailed intrinsic regulation of the ROC activation cycle remains poorly understood. Here, combining extensive molecular dynamics simulations and Markov state models, we disclosed the dynamic structural rearrangement of ROC's homodimer during nucleotide turnover. Our study revealed the coupling between dimerization extent and nucleotide-binding state, indicating a nucleotide-dependent dimerization-based activation scheme adopted by ROC GTPase. Furthermore, inspired by the well-known R1441C/G/H PD-relevant mutations within the ROC domain, we illuminated the potential allosteric molecular mechanism for its pathogenetic effects through enabling faster interconversion between inactive and active states, thus trapping ROC in a prolonged activated state, while the implicated allostery could provide further guidance for identification of regulatory allosteric pockets on the ROC complex. Our investigations illuminated the thermodynamics and kinetics of ROC homodimer during nucleotide-dependent activation for the first time and provided guidance for further exploiting ROC as therapeutic targets for controlling LRRK2 functionality in PD treatment.

Activation pathway of a G protein-coupled receptor uncovers conformational intermediates as targets for allosteric drug design

G protein-coupled receptors (GPCRs) are the most common proteins targeted by approved drugs. A complete mechanistic elucidation of large-scale conformational transitions underlying the activation mechanisms of GPCRs is of critical importance for therapeutic drug development. Here, we apply a combined computational and experimental framework integrating extensive molecular dynamics simulations, Markov state models, site-directed mutagenesis, and conformational biosensors to investigate the conformational landscape of the angiotensin II (AngII) type 1 receptor (AT1 receptor) - a prototypical class A GPCR-activation. Our findings suggest a synergistic transition mechanism for AT1 receptor activation. A key intermediate state is identified in the activation pathway, which possesses a cryptic binding site within the intracellular region of the receptor. Mutation of this cryptic site prevents activation of the downstream G protein signaling and β-arrestin-mediated pathways by the endogenous AngII octapeptide agonist, suggesting an allosteric regulatory mechanism. Together, these findings provide a deeper understanding of AT1 receptor activation at an atomic level and suggest avenues for the design of allosteric AT1 receptor modulators with a broad range of applications in GPCR biology, biophysics, and medicinal chemistry.

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.

Allosteric Mechanism of Human Mitochondrial Phenylalanyl-tRNA Synthetase: An Atomistic MD Simulation and a Mutual Information-Based Network Study

Aminoacyl-tRNA synthetases (aaRSs), a family of ubiquitous and essential enzymes, can bind target tRNAs and catalyze the aminoacylation reaction in genetic code translation. In this work, we explore the dynamic properties and allosteric communication of human mitochondrial phenylalanyl-tRNA synthetase (hmPheRS) in free and bound states to understand the mechanisms of its tRNA^Phe recognition and allostery using molecular dynamics simulations combined with the torsional mutual information-based network model. Our results reveal that hmPheRS's residue mobility and inter-residue motional coupling are significantly enhanced by tRNA^Phe binding, and there occurs a strong allosteric communication which is critical for the aminoacylation reaction, suggesting the vital role of tRNA^Phe binding in the enzyme's function. The identified signaling pathways mainly make the connections between the anticodon binding domain (ABD) and catalytic domain (CAD), as well as within the CAD composed of many functional fragments and active sites, revealing the co-regulation role of them to act coordinately and achieve hmPheRS's aminoacylation function. Besides, several key residues along the communication pathways are identified to be involved in mediating the coordinated coupling between anticodon recognition at the ABD and activation process at the CAD, showing their pivotal role in the allosteric network, which are well consistent with the experimental observation. This study sheds light on the allosteric communication mechanism in hmPheRS and can provide important information for the structure-based drug design targeting aaRSs.

Allosteric Type and Pathways Are Governed by the Forces of Protein-Ligand Binding

Allostery is central to many cellular processes, by up- or down-regulating target function. However, what determines the allosteric type remains elusive and currently it is impossible to predict whether the allosteric compounds would activate or inhibit target function before experimental studies. We demonstrated that the allosteric type and allosteric pathways are governed by the forces imposed by ligand binding to target protein using the anisotropic network model and developed an allosteric type prediction method (AlloType). AlloType correctly predicted 13 of the 16 allosteric systems in the data set with experimentally determined protein and complex structures as well as verified allosteric types, which was also used to identify allosteric pathways. When applied to glutathione peroxidase 4, a protein with no complex structure information, AlloType could still be able to predict the allosteric type of the recently reported allosteric activators, demonstrating its potential application in designing specific allosteric drugs and uncovering allosteric mechanisms.

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.

Emerging roles of SIRT6 in human diseases and its modulators

The biological functions of sirtuin 6 (SIRT6; e.g., deacetylation, defatty-acylation, and mono-ADP-ribosylation) play a pivotal role in regulating lifespan and several fundamental processes controlling aging such as DNA repair, gene expression, and telomeric maintenance. Over the past decades, the aberration of SIRT6 has been extensively observed in diverse life-threatening human diseases. In this comprehensive review, we summarize the critical roles of SIRT6 in the onset and progression of human diseases including cancer, inflammation, diabetes, steatohepatitis, arthritis, cardiovascular diseases, neurodegenerative diseases, viral infections, renal and corneal injuries, as well as the elucidation of the related signaling pathways. Moreover, we discuss the advances in the development of small molecule SIRT6 modulators including activators and inhibitors as well as their pharmacological profiles toward potential therapeutics for SIRT6-mediated diseases.

Predicting Multi-Interfacial Binding Mechanisms of NLRP3 and ASC Pyrin Domains in Inflammasome Activation

NLRP3-PYD inflammasome activates an inflammatory pathway in response to a wide variety of cell damage or infections. Dysregulated NLRP3 inflammatory signaling has many chronic inflammatory and autoimmune disorders. NLRP3 and ASC have a PYD, a superfamily member of the Death Domain, which plays a key role in inflammatory assembly. The ASC interacts with NLRP3 through a homotypic PYD and recruits the procaspase-1 through a homotypic caspase recruitment domain interaction. Here, we used several computational approaches to reveal the interactions of the NLRP3 and ASC PYD domains that lead to the activation of the inflammasome complex. We have characterized ASC and NLRP3-PYD intermolecular interactions by protein-protein docking, and further molecular dynamics (MD) simulations were conducted to evaluate the stability of NLRP3/ASC-PYD complex. Subsequently, we have identified several residues that stabilize the NLRP3/ASC-PYD complex in different faces (i.e., Face-1 to Face-4). The research framework offers new insights into the molecular mechanisms of inflammasome and apoptosis signaling as well as the ease of the drug discovery process.

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.

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

Allostery, which is one of the most direct and efficient methods to fine-tune protein functions, has gained increasing recognition in drug discovery. However, there are several challenges associated with the identification of allosteric sites, which is the fundamental cornerstone of drug design. Previous studies on allosteric site predictions have focused on communication signals propagating from the allosteric sites to the orthosteric sites. However, recent biochemical studies have revealed that allosteric coupling is bidirectional and that orthosteric perturbations can modulate allosteric sites through reversed allosteric communication. Here, we proposed a new framework for the prediction of allosteric sites based on reversed allosteric communication using a combination of computational and experimental strategies (molecular dynamics simulations, Markov state models, and site-directed mutagenesis). The desirable performance of our approach was demonstrated by predicting the known allosteric site of the small molecule MDL-801 in nicotinamide dinucleotide (NAD⁺)-dependent protein lysine deacetylase sirtuin 6 (Sirt6). A potential novel cryptic allosteric site located around the L116, R119, and S120 residues within the dynamic ensemble of Sirt6 was identified. The allosteric effect of the predicted site was further quantified and validated using both computational and experimental approaches. This study proposed a state-of-the-art computational pipeline for detecting allosteric sites based on reversed allosteric communication. This method enabled the identification of a previously uncharacterized potential cryptic allosteric site on Sirt6, which provides a starting point for allosteric drug design that can aid the identification of candidate pockets in other therapeutic targets.

Using reversed allosteric communication, we performed MD simulations, MSMs, and mutagenesis experiments, to discover allosteric sites. It reproduced the known allosteric site for MDL-801 on Sirt6 and uncovered a novel cryptic allosteric Pocket X.