Chemical Perturbation of Oncogenic Protein Folding: from the Prediction of Locally Unstable Structures to the Design of Disruptors of Hsp90-Client Interactions

N-Glycosylation-Induced Pathologic Protein Conformations as a Tool to Guide the Selection of Biologically Active Small Molecules

Post-translational modifications such as protein N-glycosylation, significantly influence cellular processes. Dysregulated N-glycosylation, exemplified in Grp94, a member of the Hsp90 family, leads to structural changes and the formation of epichaperomes, contributing to pathologies. Targeting N-glycosylation-induced conformations offers opportunities for developing selective chemical tools and drugs for these pathologic forms of chaperones. We here demonstrate how a specific Grp94 conformation induced by N-glycosylation, identified previously via molecular dynamics simulations, rationalizes the distinct behavior of similar ligands. Integrating dynamic ligand unbinding information with SAR development, we differentiate ligands productively engaging the pathologic Grp94 conformers from those that are not. Additionally, analyzing binding site stereoelectronic properties and QSAR models using cytotoxicity data unveils relationships between chemical, conformational properties, and biological activities. These findings facilitate the design of ligands targeting specific Grp94 conformations induced by abnormal glycosylation, selectively disrupting pathogenic protein networks while sparing normal mechanisms.

Cracking the chaperone code through the computational microscope

The heat shock protein 90 kDa (Hsp90) chaperone machinery plays a crucial role in maintaining cellular homeostasis. Beyond its traditional role in protein folding, Hsp90 is integral to key pathways influencing cellular function in health and disease. Hsp90 operates through the modular assembly of large multiprotein complexes, with their composition, stability, and localization adapting to the cell's needs. Its functional dynamics are finely tuned by ligand binding and post-translational modifications (PTMs). Here, we discuss how to disentangle the intricacies of the complex code that governs the crosstalk between dynamics, binding, PTMs, and the functions of the Hsp90 machinery using computer-based approaches. Specifically, we outline the contributions of computational and theoretical methods to the understanding of Hsp90 functions, ranging from providing atomic-level insights into its dynamics to clarifying the mechanisms of interactions with protein clients, cochaperones, and ligands. The knowledge generated in this framework can be actionable for the design and development of chemical tools and drugs targeting Hsp90 in specific disease-associated cellular contexts. Finally, we provide our perspective on how computation can be integrated into the study of the fine-tuning of functions in the highly complex Hsp90 landscape, complementing experimental methods for a comprehensive understanding of this important chaperone system.

An update on the status of HSP90 inhibitors in cancer clinical trials

The evolutionary conserved molecular chaperone heat shock protein 90 (HSP90) plays an indispensable role in tumorigenesis by stabilizing client oncoproteins. Although the functionality of HSP90 is tightly regulated, cancer cells exhibit a unique dependence on this chaperone, leading to its overexpression, which has been associated with poor prognosis in certain malignancies. While various strategies targeting heat shock proteins (HSPs) involved in carcinogenesis have been explored, only inhibition of HSP90 has consistently and effectively resulted in proteasomal degradation of its client proteins. To date, a total of 22 HSP90 inhibitors (HSP90i) have been tested in 186 cancer clinical trials, as reported by clinicaltrials.gov. Among these trials, 60 % have been completed, 10 % are currently active, and 30 % have been suspended, terminated, or withdrawn. HSP90 inhibitors (HSP90i) have been used as single agents or in combination with other drugs for the treatment of various cancer types in clinical trials. Notably, improved clinical outcomes have been observed when HSP90i are used in combination therapies, as they exhibit a synergistic antitumor effect. However, as single agents, HSP90i have shown limited clinical activity due to drug-related toxicity or therapy resistance. Recently, active trials conducted in Japan evaluating TAS-116 (pimitespib) have demonstrated promising results with low toxicity as monotherapy and in combination with the immune checkpoint inhibitor nivolumab. Exploratory biomarker analyses performed in various trials have demonstrated target engagement that suggests the potential for identifying patient populations that may respond favorably to the therapy. In this review, we discuss the advances made in the past 5 years regarding HSP90i and their implications in anticancer therapeutics. Our focus lies in evaluating drug efficacy, prognosis forecast, pharmacodynamic biomarkers, and clinical outcomes reported in published trials. Through this comprehensive review, we aim to shed light on the progress and potential of HSP90i as promising therapeutic agents in cancer treatment.

Structure-based discovery of small molecule inhibitors of FKBP51-Hsp90 protein-protein interaction

The heat shock protein 90 kDa (Hsp90) molecular chaperone machinery is responsible for the folding and activation of hundreds of important clients such as kinases, steroid hormone receptors, transcription factors, etc. This process is dynamically regulated in an ATP-dependent manner by Hsp90 co-chaperones including a group of tetratricopeptide (TPR) motif proteins that bind to the C-terminus of Hsp90. Among these TPR containing co-chaperones, FK506-binding protein 51 kDa (FKBP51) is reported to play an important role in stress-related pathologies, psychiatric disorders, Alzheimer's disease, and cancer, making FKBP51-Hsp90 interaction a potential therapeutic target. In this study, we report identification of potent and selective inhibitors of FKBP51-Hsp90 protein-protein interaction using a structure-based virtual screening approach. Upon in vitro evaluation, the identified hits show a considerable degree of selectivity towards FKBP51 over other TPR proteins, particularly for highly homologous FKBP52. Tyr355 of FKBP51 emerged as an important contributor to inhibitor's specificity. Additionally, we demonstrate the impact of these inhibitors on cellular energy metabolism, and neurite outgrowth, which are subjects of FKBP51 regulation. Overall, the results from this study highlight a novel pharmacological approach towards regulation of FKBP51 function and more generally, Hsp90 function via its interaction with TPR co-chaperones.

Molecules inducing specific cyclin-dependent kinase degradation and their possible use in cancer therapy

Cyclin-dependent kinases (CDKs) play an important role in the regulation of cell proliferation, and many CDK inhibitors were developed. However, pan-CDK inhibitors failed to be approved due to intolerant toxicity or low efficacy and the use of selective CDK4/6 inhibitors is limited by resistance. Protein degraders have the potential to increase selectivity, efficacy and overcome resistance, which provides a novel strategy for regulating CDKs. In this review, we summarized the function of CDKs in regulating the cell cycle and transcription, and introduced the representative CDK inhibitors. Then we made a detailed introduction about four types of CDKs degraders, including their action mechanisms, research status and application prospects, which could help the development of novel CDKs degraders.

Molecular mechanisms of chaperone-directed protein folding: Insights from atomistic simulations

Molecular chaperones, a family of proteins of which Hsp90 and Hsp70 are integral members, form an essential machinery to maintain healthy proteomes by controlling the folding and activation of a plethora of substrate client proteins. This is achieved through cycles in which Hsp90 and Hsp70, regulated by task-specific co-chaperones, process ATP and become part of a complex network that undergoes extensive compositional and conformational variations. Despite impressive advances in structural knowledge, the mechanisms that regulate the dynamics of functional assemblies, their response to nucleotides, and their relevance for client remodeling are still elusive. Here, we focus on the glucocorticoid receptor (GR):Hsp90:Hsp70:co-chaperone Hop client-loading and the GR:Hsp90:co-chaperone p23 client-maturation complexes, key assemblies in the folding cycle of glucocorticoid receptor (GR), a client strictly dependent upon Hsp90/Hsp70 for activity. Using a combination of molecular dynamics simulation approaches, we unveil with unprecedented detail the mechanisms that underpin function in these chaperone machineries. Specifically, we dissect the processes by which the nucleotide-encoded message is relayed to the client and how the distinct partners of the assemblies cooperate to (pre)organize partially folded GR during Loading and Maturation. We show how different ligand states determine distinct dynamic profiles for the functional interfaces defining the interactions in the complexes and modulate their overall flexibility to facilitate progress along the chaperone cycle. Finally, we also show that the GR regions engaged by the chaperone machinery display peculiar energetic signatures in the folded state, which enhance the probability of partial unfolding fluctuations. From these results, we propose a model where a dynamic cross-talk emerges between the chaperone dynamics states and remodeling of client-interacting regions. This factor, coupled to the highly dynamic nature of the assemblies and the conformational heterogeneity of their interactions, provides the basis for regulating the functions of distinct assemblies during the chaperoning cycle.

Pazopanib ameliorates rotenone-induced Parkinsonism in rats by suppressing multiple regulated cell death mechanisms

Parkinson's disease (PD) is characterized by motor impairments and progressive dopaminergic neuronal death in the substantia nigra (SN). Recently, the involvement of other regulated cell death (RCD) machineries has been highlighted in PD. Necroptosis is controlled by p-RIPK1, p-RIPK3, and p-MLKL and negatively regulated by caspase-8. Ferroptosis is characterized by iron overload and accumulation of reactive oxygen species. Interestingly, the molecular chaperone complex HSP90/CDC37 has been reported to directly regulate necroptosis, ferroptosis, and some PD-associated proteins. We investigated the potential anti-necroptotic and anti-ferroptotic effects of the anti-cancer drug pazopanib, uncovering the HSP90/CDC37 complex as a master RCD modulator in rotenone-induced Parkinsonism in rats. Oral administration of 15 mg/kg pazopanib to rotenone-intoxicated rats for three weeks improved motor deficits, debilitated histopathological changes, and increased striatal dopaminergic levels. Pazopanib suppressed LRRK2 and c-Abl. Pazopanib displayed an anti-necroptotic effect through inhibition of the p-RIPK1/p-RIPK3/p-MLKL pathway and activation of caspase-8. Moreover, pazopanib inhibited the ferroptotic p-VEGFR2-PKCβII-PLC-γ-ACSL-4 pathway, iron, 4-HNE, and PTGS2 while increasing GPX-4 and GSH levels. Taken together, the current research sheds light on the repositioning of pazopanib targeting HSP90/CDC37 and its multiple RCD mechanisms, which would offer a new perspective for therapeutic strategies in PD.

Lapatinib ditosylate rescues motor deficits in rotenone-intoxicated rats: Potential repurposing of anti-cancer drug as a disease-modifying agent in Parkinson's disease

Parkinson's disease (PD) is a neurodegenerative disorder characterized by motor deficits induced by dopaminergic neuronal death in the substantia nigra (SN). Finding a successful neuroprotective therapy is still challenging despite improved knowledge of the etiology of PD and a variety of medications intended to reduce symptoms. Lapatinib (LAP), an FDA-approved anti-cancer medication, has been stated to exert its effect through the modulation of oxidative stress. Furthermore, recent studies display the neuroprotective effects of LAP in epilepsy, encephalomyelitis, and Alzheimer's disease in rodent models through the modulation of oxidative stress and ferroptosis. Nevertheless, it is questionable whether LAP exerts neuroprotective effects in PD. In the current study, administration of 100 mg/kg LAP in rotenone-treated rats for 21 days ameliorates motor impairment, debilitated histopathological alterations, and revived dopaminergic neurons by increasing tyrosine hydroxylase (TH) expression in SN, along with increased dopamine level. LAP remarkably restored the antioxidant defense mechanism system, GPX4/GSH/NRF2 axis, inhibiting oxidative markers, including iron, TfR1, PTGS2, and 4-HNE, along with suppression of p-EGFR/c-SRC/PKCβII/PLC-γ/ACSL-4 pathway. Moreover, LAP modulates HSP90/CDC37 chaperone complex, regulating many key pathological markers of PD, including LRRK2, c-ABL, and α-syn. It is concluded that LAP has neuroprotective effects in PD via modulation of many key parameters implicated in PD pathogenesis. Taken together, the current study offers insights into the potential repositioning of LAP as a disease-modifying drug in PD.

Activation of autophagy depends on Atg1/Ulk1-mediated phosphorylation and inhibition of the Hsp90 chaperone machinery

Cellular homeostasis relies on both the chaperoning of proteins and the intracellular degradation system that delivers cytoplasmic constituents to the lysosome, a process known as autophagy. The crosstalk between these processes and their underlying regulatory mechanisms is poorly understood. Here, we show that the molecular chaperone heat shock protein 90 (Hsp90) forms a complex with the autophagy-initiating kinase Atg1 (yeast)/Ulk1 (mammalian), which suppresses its kinase activity. Conversely, environmental cues lead to Atg1/Ulk1-mediated phosphorylation of a conserved serine in the amino domain of Hsp90, inhibiting its ATPase activity and altering the chaperone dynamics. These events impact a conformotypic peptide adjacent to the activation and catalytic loop of Atg1/Ulk1. Finally, Atg1/Ulk1-mediated phosphorylation of Hsp90 leads to dissociation of the Hsp90:Atg1/Ulk1 complex and activation of Atg1/Ulk1, which is essential for initiation of autophagy. Our work indicates a reciprocal regulatory mechanism between the chaperone Hsp90 and the autophagy kinase Atg1/Ulk1 and consequent maintenance of cellular proteostasis.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

Emerging approaches to CDK inhibitor development, a structural perspective

Aberrant activity of the cyclin-dependent kinase family is frequently noted in a number of diseases identifying them as potential targets for drug development. However, current CDK inhibitors lack specificity owing to the high sequence and structural conservation of the ATP binding cleft across family members, highlighting the necessity of finding novel modes of CDK inhibition. The wealth of structural information regarding CDK assemblies and inhibitor complexes derived from X-ray crystallographic studies has been recently complemented through the use of cryo-electron microscopy. These recent advances have provided insights into the functional roles and regulatory mechanisms of CDKs and their interaction partners. This review explores the conformational malleability of the CDK subunit, the importance of SLiM recognition sites in CDK complexes, the progress made in chemically induced CDK degradation and how these studies can contribute to CDK inhibitor design. Additionally, fragment-based drug discovery can be utilised to identify small molecules that bind to allosteric sites on the CDK surface employing interactions which mimic those of native protein-protein interactions. These recent structural advances in CDK inhibitor mechanisms and in chemical probes which do not occupy the orthosteric ATP binding site can provide important insights for targeted CDK therapies.

Computational Epitope Prediction and Design for Antibody Development and Detection

The design of optimized protein antigens is a fundamental step in the development of new vaccine candidates and in the detection of therapeutic antibodies. A fundamental prerequisite is the identification of antigenic regions that are most prone to interact with antibodies, namely, B-cell epitopes. Here, we describe an efficient structure-based computational method for epitope prediction, called MLCE. In this approach, all that is required is the 3D structure of the antigen of interest. MLCE can be applied to glycosylated proteins, facilitating the identification of immunoreactive versus immune-shielding carbohydrates.

Second Virtual International Symposium on Cellular and Organismal Stress Responses, September 8-9, 2022

The Second International Symposium on Cellular and Organismal Stress Responses took place virtually on September 8-9, 2022. This meeting was supported by the Cell Stress Society International (CSSI) and organized by Patricija Van Oosten-Hawle and Andrew Truman (University of North Carolina at Charlotte, USA) and Mehdi Mollapour (SUNY Upstate Medical University, USA). The goal of this symposium was to continue the theme from the initial meeting in 2020 by providing a platform for established researchers, new investigators, postdoctoral fellows, and students to present and exchange ideas on various topics on cellular stress and chaperones. We will summarize the highlights of the meeting here and recognize those that received recognition from the CSSI.

Impact of Co-chaperones and Posttranslational Modifications Toward Hsp90 Drug Sensitivity

Posttranslational modifications (PTMs) regulate myriad cellular processes by modulating protein function and protein-protein interaction. Heat shock protein 90 (Hsp90) is an ATP-dependent molecular chaperone whose activity is responsible for the stabilization and maturation of more than 300 client proteins. Hsp90 is a substrate for numerous PTMs, which have diverse effects on Hsp90 function. Interestingly, many Hsp90 clients are enzymes that catalyze PTM, demonstrating one of the several modes of regulation of Hsp90 activity. Approximately 25 co-chaperone regulatory proteins of Hsp90 impact structural rearrangements, ATP hydrolysis, and client interaction, representing a second layer of influence on Hsp90 activity. A growing body of literature has also established that PTM of these co-chaperones fine-tune their activity toward Hsp90; however, many of the identified PTMs remain uncharacterized. Given the critical role of Hsp90 in supporting signaling in cancer, clinical evaluation of Hsp90 inhibitors is an area of great interest. Interestingly, differential PTM and co-chaperone interaction have been shown to impact Hsp90 binding to its inhibitors. Therefore, understanding these layers of Hsp90 regulation will provide a more complete understanding of the chaperone code, facilitating the development of new biomarkers and combination therapies.

Synthesis and Evaluation of Small Molecule Disruptors of the Aha1/Hsp90 Complex for the Reduction of Tau Aggregation

KU-177 was recently shown to disrupt interactions between Hsp90 and Aha1 in vitro. Subsequent studies in recombinant thioflavin T (ThT) assays demonstrated that KU-177 ablates Aha1-driven enhancement of Hsp90-dependent tau aggregation, which was confirmed by TEM. Using KU-177 as a lead compound, derivatives of KU-177 were synthesized and evaluated for their ability to disrupt Aha1/Hsp90 interactions and inhibit P301L tau aggregation. Preliminary structure-activity relationships were revealed, which led to the identification of a new lead compound that contains a cis-like amide bond. The new lead compounds retain the ability to disrupt Aha1/Hsp90 interactions in SH-SY5Y and SK-BR-3 cells without direct inhibition of Hsp90, providing a new scaffold for subsequent drug discovery efforts.

Therapeutic potential of CDK4/6 inhibitors in renal cell carcinoma

The treatment of advanced and metastatic kidney cancer has entered a golden era with the addition of more therapeutic options, improved survival and new targeted therapies. Tyrosine kinase inhibitors, mammalian target of rapamycin (mTOR) inhibitors and immune checkpoint blockade have all been shown to be promising strategies in the treatment of renal cell carcinoma (RCC). However, little is known about the best therapeutic approach for individual patients with RCC and how to combat therapeutic resistance. Cancers, including RCC, rely on sustained replicative potential. The cyclin-dependent kinases CDK4 and CDK6 are involved in cell-cycle regulation with additional roles in metabolism, immunogenicity and antitumour immune response. Inhibitors of CDK4 and CDK6 are now commonly used as approved and investigative treatments in breast cancer, as well as several other tumours. Furthermore, CDK4/6 inhibitors have been shown to work synergistically with other kinase inhibitors, including mTOR inhibitors, as well as with immune checkpoint inhibitors in preclinical cancer models. The effect of CDK4/6 inhibitors in kidney cancer is relatively understudied compared with other cancers, but the preclinical studies available are promising. Collectively, growing evidence suggests that targeting CDK4 and CDK6 in kidney cancer, alone and in combination with current therapeutics including mTOR and immune checkpoint inhibitors, might have therapeutic benefit and should be further explored.

Discovery of potent Camkk1 kinase inhibitors through e-pharmacophore and molecular screening approaches

Diabetes is recognized as a major health problem and according to WHO estimates global prevalence of diabetes is expected to increase from 171 million in 2000 to 366 million in 2030, among which 21.7% will be Indians. The chronic nature of diabetes leads to several metabolic complications like kidney failure, cardiac problems and hypertension, etc.Camkk family members are attractive and emerging targets for the development of anti-diabetic drugs. However, the selectivity of inhibitors is a crucial property as a lack of selectivity could lead to serious adverse effects. STO-609 recently reported the role of Camkks a as potent inhibitor. In this study, Combined Molecular Docking and Pharmacophore Mapping were employed to identify potent lead molecules. E-Pharmacophore based virtual screening was performed against commercially available databases to identify the best lead molecule which was docked with the targets and analyzed for the binding pattern. Also, ADME and density function theory (DFT) studies of the compound were performed and the hits that showed good binding to the active sites and that matched with the pharmacophore models were considered as possible functional molecules against Camkk1. The results from e-pharmacophore based virtual screening and MD simulations evidenced that the top three compounds namely (Lifechemicals_1, Zinc_0910993 and Binding_10131) will be a promising inhibitor for Camkk1 family.Communicated by Ramaswamy H. Sarma.

Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q²) of 0.68, 0.641 and linear regression values (r²) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets.

The tumor suppressor folliculin inhibits lactate dehydrogenase A and regulates the Warburg effect

Aerobic glycolysis in cancer cells, also known as the 'Warburg effect', is driven by hyperactivity of lactate dehydrogenase A (LDHA). LDHA is thought to be a substrate-regulated enzyme, but it is unclear whether a dedicated intracellular protein also regulates its activity. Here, we identify the human tumor suppressor folliculin (FLCN) as a binding partner and uncompetitive inhibitor of LDHA. A flexible loop within the amino terminus of FLCN controls movement of the LDHA active-site loop, tightly regulating its enzyme activity and, consequently, metabolic homeostasis in normal cells. Cancer cells that experience the Warburg effect show FLCN dissociation from LDHA. Treatment of these cells with a decapeptide derived from the FLCN loop region causes cell death. Our data suggest that the glycolytic shift of cancer cells is the result of FLCN inactivation or dissociation from LDHA. Together, FLCN-mediated inhibition of LDHA provides a new paradigm for the regulation of glycolysis.

Chemical tools for epichaperome-mediated interactome dysfunctions of the central nervous system

Diseases are a manifestation of how thousands of proteins interact. In several diseases, such as cancer and Alzheimer’s disease, proteome-wide disturbances in protein-protein interactions are caused by alterations to chaperome scaffolds termed epichaperomes. Epichaperome-directed chemical probes may be useful for detecting and reversing defective chaperomes. Here we provide structural, biochemical, and functional insights into the discovery of epichaperome probes, with a focus on their use in central nervous system diseases. We demonstrate on-target activity and kinetic selectivity of a radiolabeled epichaperome probe in both cells and mice, together with a proof-of-principle in human patients in an exploratory single group assignment diagnostic study (ClinicalTrials.gov Identifier: NCT03371420). The clinical study is designed to determine the pharmacokinetic parameters and the incidence of adverse events in patients receiving a single microdose of the radiolabeled probe administered by intravenous injection. In sum, we introduce a discovery platform for brain-directed chemical probes that specifically modulate epichaperomes and provide proof-of-principle applications in their use in the detection, quantification, and modulation of the target in complex biological systems.

Here, the authors show structural, biochemical, and functional insights into the discovery of epichaperome‐ directed chemical probes for use in central nervous system diseases. Probes emerging from this work have translated to human clinical studies in Alzheimer’s disease and cancer.

Design and Synthesis of TRAP1 Selective Inhibitors: H-Bonding with Asn171 Residue in TRAP1 Increases Paralog Selectivity

Tumor necrosis factor receptor-associated protein 1 (TRAP1) is overexpressed in the mitochondria of various cancer cells, reprograms cellular metabolism to enable cancer cells to adapt to harsh tumor environments. As inactivation of TRAP1 induces massive apoptosis in cancer cells in vitro and in vivo, the development of TRAP1-selective inhibitors has become an attractive approach. A series of purine-8-one and pyrrolo[2,3-d]pyrimidine derivatives was developed based on TRAP1 structure and identified to be highly selective in vitro for TRAP1 over the paralogous enzymes, Hsp90α and Grp94. The TRAP1-selective inhibition strategy via utilization of the Asn171 residue of the ATP-lid was investigated using X-ray crystallography and molecular dynamics simulation studies. Among various synthesized potent TRAP1 inhibitors, 5f possessed a 65-fold selectivity over Hsp90α and a 13-fold selectivity over Grp94. Additionally, 6f had a half-maximal inhibitory concentration (IC₅₀) of 63.5 nM for TRAP1, with a 78-fold and 30-fold selectivity over Hsp90α and Grp94, respectively.

Unraveling the Allosteric Cross-Talk between the Coactivator Peptide and the Ligand-Binding Site in the Glucocorticoid Receptor

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known allosteric cross-talk between the ligand and coregulator binding sites within the GR ligand-binding domain that is crucial for the control of the functional response. However, the molecular mechanisms underlying such an allosteric control remain elusive. Here, molecular dynamics (MD) simulations, bioinformatic analysis, and biophysical measurements are integrated to capture the structural and dynamic features of the allosteric cross-talk within the GR. We identified a network of evolutionarily conserved residues that enables the allosteric signal transduction, in agreement with experimental data. MD simulations clarify how such a network is dynamically interconnected and offer a mechanistic explanation of how different peptides affect the intensity of the allosteric signal. This study provides useful insights to elucidate the GR allosteric regulation, ultimately providing a foundation for designing novel drugs.

Deciphering the Molecular Terms of Arp2/3 Allosteric Regulation from All-Atom Simulations and Dynamical Network Theory

The Arp2/3 molecular machine stimulates the generation of branched actin networks at the cytosolic surface of cellular membranes. Arp2/3 is thus pivotal for cell motility and migration, and its aberrant function is implicated in cancer invasion and metastasis. Here, all-atom multi μs-long molecular dynamics simulations and dynamical NetWork Analysis (NWA) unprecedentedly disclose the molecular terms of Arp2/3 regulation (activation/inhibition) by positive/negative allosteric modulators. After identifying the crucial structural elements underlying Arp2/3's conformational transition toward its active actin-polymerization-competent state, we decrypt the activating signaling paths heading from the allosteric effector (ATP) binding sites to these pivotal regions, also elucidating how small-molecule inhibitors scramble this signal-exchange. As a result, while ATP-induced signaling triggers a harmonious conformational transition toward active Arp2/3, the inhibitors disturb these information channels, desynchronizing Arp2/3 functional movements, thus hindering its activation. Our outcomes supply a conceptual basis for devising small-molecule inhibitors to block infiltrative cancer migration.

Discovery of pseudolaric acid A as a new Hsp90 inhibitor uncovers its potential anticancer mechanism

Pseudolaric acid A (PAA), one of the main bioactive ingredients in traditional medicine Pseudolarix cortex, exhibits remarkable anticancer activities. Yet its mechanism of action and molecular target have not been investigated and remain unclear. In this work, mechanistic study showed that PAA induced cell cycle arrest at G2/M phase and promoted cell death through caspase-8/caspase-3 pathway, demonstrating potent antiproliferation and anticancer activities. PAA was discovered to be a new Hsp90 inhibitor and multiple biophysical experiments confirmed that PAA directly bind to Hsp90. Active PAA-probe was designed, synthesized and biological evaluated. It was subsequently employed to verify the cellular interaction with Hsp90 in HeLa cells through photoaffinity labeling approach. Furthermore, NMR experiments showed that N-terminal domain of Hsp90 and essential groups in PAA are important for the protein-inhibitor recognition. Structure-activity relationship studies revealed the correlation between its Hsp90 inhibitory activity with anticancer activity. This work proposed a potential mechanism involved with the anticancer activity of PAA and will improve the appreciation of PAA as a potential cancer therapy candidate.

Evolutionary algorithm-based generation of optimum peptide sequences with dengue virus inhibitory activity

Background: There is currently no effective dengue virus (DENV) therapeutic. We aim to develop a genetic algorithm-based framework for the design of peptides with possible DENV inhibitory activity. Methods & results: A Python-based tool (denominated AutoPepGEN) based on a DENV support vector machine classifier as the objective function was implemented. AutoPepGEN was applied to the design of three- to seven-amino acid sequences and ten peptides were selected. Peptide-protease (DENV) docking and Molecular Mechanics-Generalized Born Surface Area calculations were performed for the selected sequences and favorable binding energies were observed. Conclusion: It is hoped that AutoPepGEN will serve as an in silico alternative to the experimental design of positional scanning combinatorial libraries, known to be prone to a combinatorial explosion. AutoPepGEN is available at: https://github.com/sjbarigye/AutoPepGEN.

Machine Learning Prediction of Allosteric Drug Activity from Molecular Dynamics

Allosteric drugs have been attracting increasing interest over the past few years. In this context, it is common practice to use high-throughput screening for the discovery of non-natural allosteric drugs. While the discovery stage is supported by a growing amount of biological information and increasing computing power, major challenges still remain in selecting allosteric ligands and predicting their effect on the target protein's function. Indeed, allosteric compounds can act both as inhibitors and activators of biological responses. Computational approaches to the problem have focused on variations on the theme of molecular docking coupled to molecular dynamics with the aim of recovering information on the (long-range) modulation typical of allosteric proteins.

Molecular Mechanisms of the Blockage of Glioblastoma Motility

Glioblastoma (GBM) is the most common and lethal brain tumor. GBM has a remarkable degree of motility and is able to infiltrate the healthy brain. In order to perform a rationale-based drug-repositioning study, we have used known inhibitors of two small Rho GTPases, Rac1 and Cdc42, which are upregulated in GBM and are involved in the signaling processes underlying the orchestration of the cytoskeleton and cellular motility. The selected inhibitors (R-ketorolac and ML141 for Cdc42 and R-ketorolac and EHT 1864 for Rac1) have been successfully employed to reduce the infiltration propensity of GBM in live cell imaging studies. Complementarily, all-atom simulations have elucidated the molecular basis of their inhibition mechanism, identifying the binding sites targeted by the inhibitors and dissecting their impact on the small Rho GTPases' function. Our results demonstrate the potential of targeting the Rac1 and Cdc42 proteins with small molecules to contrast GBM infiltration growth and supply precious information for future drug discovery studies aiming to fight GBM and other infiltrative cancer types.

Targeting the HSP90-CDC37-kinase chaperone cycle: A promising therapeutic strategy for cancer

Heat shock protein 90 (HSP90) is an indispensable molecular chaperone that facilitates the maturation of numerous oncoproteins in cancer cells, including protein kinases, ribonucleoproteins, steroid hormone receptors, and transcription factors. Although over 30 HSP90 inhibitors have steadily entered clinical trials, further clinical advancement has been restricted by their limited efficacy, inevitable heat shock response, and multiple side-effects, likely induced via an ATP inhibition mechanism. Since both ATP and various co-chaperones play essential roles in the HSP90 chaperone cycle to achieve integrated function, optimal therapeutics require an understanding of the dynamic interactions among HSP90, ATP, and cochaperones. To date, continuous research has promoted the exploration of the cochaperone cell division cycle 37 (CDC37) as a kinase-specific recognizer and has shown that the HSP90-CDC37-kinase complex is particularly relevant in cancers. Indeed, disrupting the HSP90-CDC37-kinase complex, rather than totally blocking the ATP function of HSP90, is emerging as an alternative way to avoid the limitations of current inhibitors. In this review, we first briefly introduce the HSP90-CDC37-kinase cycle and present the currently available approaches for inhibitor development targeting this cycle and provide insights into selective regulation of the kinase clients of HSP90 by more directional ways.

Structural Model for Recruitment of RIT1 to the LZTR1 E3 Ligase: Evidences from an Integrated Computational Approach

Leucine-zipper transcription regulator 1 (LZTR1) is a highly mutated tumor suppressor gene, involved in the pathogenesis of several cancer types and developmental disorders. In proteasomal degradation, it acts as an adaptor protein responsible for the recognition and recruitment of substrates to be ubiquitinated in Cullin3-RING ligase E3 (CRL3) machinery. LZTR1 belongs to the BTB-Kelch family, a multi-domain protein where the Kelch propeller plays as the substrate recognition region and for which no experimental structure has been solved. Recently, large effort mutational analyses pointed to the role of disease-associated LZTR1 mutations in the RAS/MAPK signaling pathway and RIT1, a small Ras-related GTPase protein, has been identified by mass spectroscopy to interact with LZTR1. Hence, a better understanding of native structure, molecular mechanism, and substrate specificity would help clarifying the role of LZTR1 in pathological diseases, thus promoting advancement in the development of novel therapeutic strategies. Here, we address the interaction model between adaptor LZTR1 and substrate RIT1 by applying an integrated computational approach, including molecular modeling and docking techniques. We observe that the interaction model LZTR1-RIT1 is stabilized by an electrostatic bond network established between the two protein surfaces, which is reminiscent of homologous ubiquitin ligases complexes. Then, running MD simulations, we characterize differential conformational dynamics of the multi-domain LZTR1, offering interesting implications on the mechanistic role of specific point mutations. We identify G248R and R283Q as damaging mutations involved in the recognition process of the substrate RIT1 and R412C as a possible allosteric mutation from the Kelch to the C-term BTB-domain. Our findings provide important structural insights on targeting CRL3s for drug discovery.

Hsp90 modulates human sperm capacitation via the Erk1/2 and p38 MAPK signaling pathways

BACKGROUND

Heat shock protein 90 (Hsp90) is a highly abundant eukaryotic molecular chaperone that plays important roles in client protein maturation, protein folding and degradation, and signal transduction. Previously, we found that both Hsp90 and its co-chaperone cell division cycle protein 37 (Cdc37) were expressed in human sperm. Hsp90 is known to be involved in human sperm capacitation via unknown underlying mechanism(s). As Cdc37 was a kinase-specific co-chaperone of Hsp90, Hsp90 may regulate human sperm capacitation via other kinases. It has been reported that two major mitogen-activated protein kinases (MAPKs), extracellular signal-regulated kinase 1/2 (Erk1/2) and p38, are expressed in human sperm in the same locations as Hsp90 and Cdc37. Phosphorylated Erk1/2 has been shown to promote sperm hyperactivated motility and acrosome reaction, while phosphorylated p38 inhibits sperm motility. Therefore, in this study we explored whether Hsp90 modulates human sperm capacitation via the Erk1/2 and p38 MAPK signaling pathways.

METHODS

Human sperm was treated with the Hsp90-specific inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG) during capacitation. Computer-assisted sperm analyzer (CASA) was used to detect sperm motility and hyperactivation. The sperm acrosome reaction was analyzed by using fluorescein isothiocyanate-conjugated Pisum sativum agglutinin (PSA-FITC) staining. The interactions between Hsp90, Cdc37, Erk1/2 and p38 were assessed using co-immunoprecipitation (Co-IP) experiments. Western blotting analysis was used to evaluate the levels of protein expression and phosphorylation.

RESULTS

Human sperm hyperactivation and acrosome reaction were inhibited by 17-AAG, suggesting that Hsp90 is involved in human sperm capacitation. In addition, Co-IP experiments revealed that 17-AAG reduced the interaction between Hsp90 and Cdc37, leading to the dissociation of Erk1/2 from the Hsp90-Cdc37 protein complex. Western blotting analysis revealed that levels of Erk1/2 and its phosphorylated form were subsequently decreased. Decreasing of Hsp90-Cdc37 complex also affected the interaction between Hsp90 and p38. Nevertheless, p38 dissociated from the Hsp90 protein complex and was activated by autophosphorylation.

CONCLUSIONS

Taken together, our findings indicate that Hsp90 is involved in human sperm hyperactivation and acrosome reaction. In particular, Hsp90 and its co-chaperone Cdc37 form a protein complex with Erk1/2 and p38 to regulate their kinase activity. These results suggest that Hsp90 regulates human sperm capacitation via the Erk1/2 and p38 MAPK signaling pathways.

Exploiting Folding and Degradation Machineries To Target Undruggable Proteins: What Can a Computational Approach Tell Us?

Advances in genomics and proteomics have unveiled an ever-growing number of key proteins and provided mechanistic insights into the genesis of pathologies. This wealth of data showed that changes in expression levels of specific proteins, mutations, and post-translational modifications can result in (often subtle) perturbations of functional protein-protein interaction networks, which ultimately determine disease phenotypes. Although many such validated pathogenic proteins have emerged as ideal drug targets, there are also several that escape traditional pharmacological regulation; these proteins have thus been labeled "undruggable". The challenges posed by undruggable targets call for new sorts of molecular intervention. One fascinating solution is to perturb a pathogenic protein's expression levels, rather than blocking its activities. In this Concept paper, we shall discuss chemical interventions aimed at recruiting undruggable proteins to the ubiquitin proteasome system, or aimed at disrupting protein-protein interactions in the chaperone-mediated cellular folding machinery: both kinds of intervention lead to a decrease in the amount of active pathogenic protein expressed. Specifically, we shall discuss the role of computational strategies in understanding the molecular determinants characterizing the function of synthetic molecules typically designed for either type of intervention. Finally, we shall provide our perspectives and views on the current limitations and possibilities to expand the scope of rational approaches to the design of chemical regulators of protein levels.

Markov state models and NMR uncover an overlooked allosteric loop in p53

The tumor suppressor p53 is the most frequently mutated gene in human cancer, and thus reactivation of mutated p53 is a promising avenue for cancer therapy. Analysis of wildtype p53 and the Y220C cancer mutant long-timescale molecular dynamics simulations with Markov state models and validation by NMR relaxation studies has uncovered the involvement of loop L6 in the slowest motions of the protein. Due to its distant location from the DNA-binding surface, the conformational dynamics of this loop has so far remained largely unexplored. We observe mutation-induced stabilization of alternate L6 conformations, distinct from all experimentally-determined structures, in which the loop is both extended and located further away from the DNA-interacting surface. Additionally, the effect of the L6-adjacent Y220C mutation on the conformational landscape of the functionally-important loop L1 suggests an allosteric role to this dynamic loop and the inactivation mechanism of the mutation. Finally, the simulations reveal a novel Y220C cryptic pocket that can be targeted for p53 rescue efforts. Our approach exemplifies the power of the MSM methodology for uncovering intrinsic dynamic and kinetic differences among distinct protein ensembles, such as for the investigation of mutation effects on protein function.

Wildtype and Y220C L1 and L6 loops conformational landscape, with MSM-identified L6 states highlighted on the right.

Allosteric inhibition explained through conformational ensembles sampling distinct "mixed" states

Allosteric modulation provides an effective avenue for selective and potent enzyme inhibition. Here, we summarize and critically discuss recent advances on the mechanisms of allosteric partial agonists for three representative signalling enzymes activated by cyclic nucleotides: the cAMP-dependent protein kinase (PKA), the cGMP-dependent protein kinase (PKG), and the exchange protein activated by cAMP (EPAC). The comparative analysis of partial agonism in PKA, PKG and EPAC reveals a common emerging theme, i.e. the sampling of distinct “mixed” conformational states, either within a single domain or between distinct domains. Here, we show how such “mixed” states play a crucial role in explaining the observed functional response, i.e. partial agonism and allosteric pluripotency, as well as in maximizing inhibition while minimizing potency losses. In addition, by combining Nuclear Magnetic Resonance (NMR), Molecular Dynamics (MD) simulations and Ensemble Allosteric Modeling (EAM), we also show how to map the free-energy landscape of conformational ensembles containing “mixed” states. By discussing selected case studies, we illustrate how MD simulations and EAM complement NMR to quantitatively relate protein dynamics to function. The resulting NMR- and MD-based EAMs are anticipated to inform not only the design of new generations of highly selective allosteric inhibitors, but also the choice of multidrug combinations.

Rapid in silico Design of Potential Cyclic Peptide Binders Targeting Protein-Protein Interfaces

Rational design of specific inhibitors of protein-protein interactions is desirable for drug design to control cellular signal transduction but also for studying protein-protein interaction networks. We have developed a rapid computational approach to rationally design cyclic peptides that potentially bind at desired regions of the interface of protein-protein complexes. The methodology is based on comparing the protein backbone structure of short peptide segments (epitopes) at the protein-protein interface with a collection of cyclic peptide backbone structures. A cyclic peptide that matches the backbone structure of the segment is used as a template for a binder by adapting the amino acid side chains to the side chains found in the target complex. For a small library of cyclic peptides with known high resolution structures we found for the majority (~82%) of 154 protein-protein complexes at least one very well fitting match for a cyclic peptide template to a protein-protein interface segment. The majority of the constructed protein-cyclic peptide complexes was very stable during Molecular Dynamics simulations and showed an interaction energy score that was typically more favorable compared to interaction scores of typical peptide-protein complexes. Our cPEPmatch approach could be a promising approach for rapid suggestion of cyclic peptide binders that could be tested experimentally and further improved by chemical modification.

The Answer Lies in the Energy: How Simple Atomistic Molecular Dynamics Simulations May Hold the Key to Epitope Prediction on the Fully Glycosylated SARS-CoV-2 Spike Protein

SARS-CoV-2 is a health threat with dire socioeconomical consequences. As the crucial mediator of infection, the viral glycosylated spike protein (S) has attracted the most attention and is at the center of efforts to develop therapeutics and diagnostics. Herein, we use an original decomposition approach to identify energetically uncoupled substructures as antibody binding sites on the fully glycosylated S. Crucially, all that is required are unbiased MD simulations; no prior knowledge of binding properties or ad hoc parameter combinations is needed. Our results are validated by experimentally confirmed structures of S in complex with anti- or nanobodies. We identify poorly coupled subdomains that are poised to host (several) epitopes and potentially involved in large functional conformational transitions. Moreover, we detect two distinct behaviors for glycans: those with stronger energetic coupling are structurally relevant and protect underlying peptidic epitopes, and those with weaker coupling could themselves be prone to antibody recognition.

ATP-Driven Nonequilibrium Activation of Kinase Clients by the Molecular Chaperone Hsp90

The molecular chaperone 90-kDa heat-shock protein (Hsp90) assists the late-stage folding and activation of diverse types of protein substrates (called clients), including many kinases. Previous studies have established that the Hsp90 homodimer undergoes an ATP-driven cycle through open and closed conformations. Here, I propose a model of client activation by Hsp90 that predicts that this cycle enables Hsp90 to use ATP energy to drive a client out of thermodynamic equilibrium toward its active conformation. My model assumes that an Hsp90-bound client can transition between a deactivating conformation and an activating conformation. It suggests that the cochaperone Cdc37 aids Hsp90 to activate kinase clients by differentiating between these two intermediate conformations. My model makes experimentally testable predictions, including how modulating the stepwise kinetics of the Hsp90 cycle-for example, by various cochaperones-affects the activation of different clients. My model may inform client-specific and cell-type-specific therapeutic intervention of Hsp90-mediated protein activation.

Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins

A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

mTOR–mLST8 interaction: hot spot identification through quantum biochemistry calculations

Quantum calculation of mTOR–mLST8 interaction.