Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks

Simulation of the ligand-leaving process of the human heat shock protein

The human heat shock protein plays a critical role in various diseases and is an important target for pharmacological modulation. Simulation of conformational changes and free energy profiles of the human heat shock protein derived by the ligand-leaving process is a challenging issue. In this work, steered molecular dynamics simulation was adopted to simulate the ligand-leaving process. Two composite systems of heat shock protein NHSP90 and small molecules 6FJ and 6G7 are selected as research objects. The free energy during the leaving of ligand small molecules is calculated using conventional molecular dynamics simulation, steered molecular dynamics simulation (SMD), and the umbrella sampling method. We found that the a slower pulling velocity (0.001 nm ns-1) will result in 2.19 kcal mol-1, and the umbrella sampling method gives a value of 3.26 kcal mol-1 for the free energy difference for the two systems, which reasonably agrees with experimental results. A faster-pulling velocity (0.01 nm ns-1) leads to a large overestimation of free energy. At the same time, the conformational analysis indicated that the faster pulling velocity may lead to the conformational change of NHSP90, which was proved to be false by the slower pulling velocity and the umbrella sampling method.

Machine learning and protein allostery

The fundamental biological importance and complexity of allosterically regulated proteins stem from their central role in signal transduction and cellular processes. Recently, machine-learning approaches have been developed and actively deployed to facilitate theoretical and experimental studies of protein dynamics and allosteric mechanisms. In this review, we survey recent developments in applications of machine-learning methods for studies of allosteric mechanisms, prediction of allosteric effects and allostery-related physicochemical properties, and allosteric protein engineering. We also review the applications of machine-learning strategies for characterization of allosteric mechanisms and drug design targeting SARS-CoV-2. Continuous development and task-specific adaptation of machine-learning methods for protein allosteric mechanisms will have an increasingly important role in bridging a wide spectrum of data-intensive experimental and theoretical technologies.

Integrated metabolomics and phosphoproteomics reveal the protective role of exosomes from human umbilical cord mesenchymal stem cells in naturally aging mouse livers

BACKGROUND

Aging is characterized by a general decline in cellular function, which ultimately affects whole body homeostasis. This study aimed to investigate the effects and underlying mechanisms of exosomes derived from human umbilical cord mesenchymal stem cells (hUCMSC-exos) on the livers of naturally aging mice.

METHOD

Twenty-two-month-old C57BL6 mice were used as a natural aging animal model, divided into a saline-treated wild-type aged control group (WT-AC) and a hUCMSC-exo-treated group (WT-AEX), and then detected by morphology, metabolomics and phosphoproteomics.

RESULTS

Morphological analysis showed that hUCMSC-exos ameliorated structural disorder and decreased markers of senescence and genome instability in aging livers. Metabolomics showed that hUCMSC-exos decreased the contents of saturated glycerophospholipids, palmitoyl-glycerols and eicosanoid derivatives associated with lipotoxicity and inflammation, consistent with the decreased phosphorylation of metabolic enzymes, such as propionate-CoA ligase (Acss2), at S267 detected by phosphoproteomics. Moreover, phosphoproteomics indicated that hUCMSC-exos reduced the phosphorylation of proteins participating in nuclear transport and cancer signaling, such as heat shock protein HSP90-beta (Hsp90ab1) at S226 and nucleoprotein TPR (Tpr) at S453 and S379, while increasing those involved in intracellular communication, such as calnexin (Canx) at S563 and PDZ domain-containing protein 8 (Pdzd8). Finally, phosphorylated HSP90β and Tpr were verified predominantly in hepatocytes.

CONCLUSION

HUCMSC-exos improved metabolic reprogramming and genome stability mainly associated with phosphorylated HSP90β in hepatocytes in natural aging livers. This work provides a comprehensive resource of biological data by omics to support future investigations of hUCMSC-exos in aging.

Probing conformational landscapes and mechanisms of allosteric communication in the functional states of the ABL kinase domain using multiscale simulations and network-based mutational profiling of allosteric residue potentials

In the current study, multiscale simulation approaches and dynamic network methods are employed to examine the dynamic and energetic details of conformational landscapes and allosteric interactions in the ABL kinase domain that determine the kinase functions. Using a plethora of synergistic computational approaches, we elucidate how conformational transitions between the active and inactive ABL states can employ allosteric regulatory switches to modulate intramolecular communication networks between the ATP site, the substrate binding region, and the allosteric binding pocket. A perturbation-based network approach that implements mutational profiling of allosteric residue propensities and communications in the ABL states is proposed. Consistent with biophysical experiments, the results reveal functionally significant shifts of the allosteric interaction networks in which preferential communication paths between the ATP binding site and substrate regions in the active ABL state become suppressed in the closed inactive ABL form, which in turn features favorable allosteric coupling between the ATP site and the allosteric binding pocket. By integrating the results of atomistic simulations with dimensionality reduction methods and Markov state models, we analyze the mechanistic role of macrostates and characterize kinetic transitions between the ABL conformational states. Using network-based mutational scanning of allosteric residue propensities, this study provides a comprehensive computational analysis of long-range communications in the ABL kinase domain and identifies conserved regulatory hotspots that modulate kinase activity and allosteric crosstalk between the allosteric pocket, ATP binding site, and substrate binding regions.

Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery

Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.

Heat Shock Protein 90 (HSP90) Inhibitors as Anticancer Medicines: A Review on the Computer-Aided Drug Discovery Approaches over the Past Five Years

Cancer is a disease caused by the uncontrolled, abnormal growth of cells in different anatomic sites. In 2018, it was predicted that the worldwide cancer burden would rise to 18.1 million new cases and 9.6 million deaths. Anticancer compounds, often known as chemotherapeutic medicines, have gained much interest in recent cancer research. These medicines work through various biological processes in targeting cells at various stages of the cell's life cycle. One of the most significant roadblocks to developing anticancer drugs is that traditional chemotherapy affects normal cells and cancer cells, resulting in substantial side effects. Recently, advancements in new drug development methodologies and the prediction of the targeted interatomic and intermolecular ligand interaction sites have been beneficial. This has prompted further research into developing and discovering novel chemical species as preferred therapeutic compounds against specific cancer types. Identifying new drug molecules with high selectivity and specificity for cancer is a prerequisite in the treatment and management of the disease. The overexpression of HSP90 occurs in patients with cancer, and the HSP90 triggers unstable harmful kinase functions, which enhance carcinogenesis. Therefore, the development of potent HSP90 inhibitors with high selectivity and specificity becomes very imperative. The activities of HSP90 as chaperones and cochaperones are complex due to the conformational dynamism, and this could be one of the reasons why no HSP90 drugs have made it beyond the clinical trials. Nevertheless, HSP90 modulations appear to be preferred due to the competitive inhibition of the targeted N-terminal adenosine triphosphate pocket. This study, therefore, presents an overview of the various computational models implored in the development of HSP90 inhibitors as anticancer medicines. We hereby suggest an extensive investigation of advanced computational modelling of the three different domains of HSP90 for potent, effective inhibitor design with minimal off-target effects.

Dissecting mutational allosteric effects in alkaline phosphatases associated with different Hypophosphatasia phenotypes: An integrative computational investigation

Hypophosphatasia (HPP) is a rare inherited disorder characterized by defective bone mineralization and is highly variable in its clinical phenotype. The disease occurs due to various loss-of-function mutations in ALPL, the gene encoding tissue-nonspecific alkaline phosphatase (TNSALP). In this work, a data-driven and biophysics-based approach is proposed for the large-scale analysis of ALPL mutations-from nonpathogenic to severe HPPs. By using a pipeline of synergistic approaches including sequence-structure analysis, network modeling, elastic network models and atomistic simulations, we characterized allosteric signatures and effects of the ALPL mutations on protein dynamics and function. Statistical analysis of molecular features computed for the ALPL mutations showed a significant difference between the control, mild and severe HPP phenotypes. Molecular dynamics simulations coupled with protein structure network analysis were employed to analyze the effect of single-residue variation on conformational dynamics of TNSALP dimers, and the developed machine learning model suggested that the topological network parameters could serve as a robust indicator of severe mutations. The results indicated that the severity of disease-associated mutations is often linked with mutation-induced modulation of allosteric communications in the protein. This study suggested that ALPL mutations associated with mild and more severe HPPs can exert markedly distinct effects on the protein stability and long-range network communications. By linking the disease phenotypes with dynamic and allosteric molecular signatures, the proposed integrative computational approach enabled to characterize and quantify the allosteric effects of ALPL mutations and role of allostery in the pathogenesis of HPPs.

Anti-Cancer Properties of Ginkgolic Acids in Human Nasopharyngeal Carcinoma CNE-2Z Cells via Inhibition of Heat Shock Protein 90

Ginkgo biloba L. has been used in traditional Chinese medicine (TCM) for thousands of years. However, the anti-cancer properties of ginkgolic acids (GAS) isolated from G. biloba have not been investigated in human nasopharyngeal carcinoma cells. In this study, GAS exhibited an inhibitory effect on the ATPase activity of heat shock protein 90 (Hsp90) and anti-proliferative activities against four human cancer cell lines, with IC₅₀ values ranging from 14.91 to 23.81 μg·mL⁻¹. In vivo experiments confirmed that GAS inhibited tumor growth in CNE-2Z cell-xenografted nude mice with low hepatotoxicity. We further demonstrated that GAS suppressed migration and invasion and induced the apoptosis of CNE-2Z cells by inducing the degradation of Hsp90 client proteins (MMP-2, MMP-9, Her-2, c-Raf, Akt, and Bcl-2). Together, GAS are new Hsp90 inhibitors by binding to Hsp90 (hydrogen bond and hydrophobic interaction). Thus, GAS from G. biloba might represent promising Hsp90 inhibitors for the development of anti-nasopharyngeal carcinoma agents.

Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins

The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.

Post-translational modifications of Hsp90 and translating the chaperone code

Cells have a remarkable ability to synthesize large amounts of protein in a very short period of time. Under these conditions, many hydrophobic surfaces on proteins may be transiently exposed, and the likelihood of deleterious interactions is quite high. To counter this threat to cell viability, molecular chaperones have evolved to help nascent polypeptides fold correctly and multimeric protein complexes assemble productively, while minimizing the danger of protein aggregation. Heat shock protein 90 (Hsp90) is an evolutionarily conserved molecular chaperone that is involved in the stability and activation of at least 300 proteins, also known as clients, under normal cellular conditions. The Hsp90 clients participate in the full breadth of cellular processes, including cell growth and cell cycle control, signal transduction, DNA repair, transcription, and many others. Hsp90 chaperone function is coupled to its ability to bind and hydrolyze ATP, which is tightly regulated both by co-chaperone proteins and post-translational modifications (PTMs). Many reported PTMs of Hsp90 alter chaperone function and consequently affect myriad cellular processes. Here, we review the contributions of PTMs, such as phosphorylation, acetylation, SUMOylation, methylation, O-GlcNAcylation, ubiquitination, and others, toward regulation of Hsp90 function. We also discuss how the Hsp90 modification state affects cellular sensitivity to Hsp90-targeted therapeutics that specifically bind and inhibit its chaperone activity. The ultimate challenge is to decipher the comprehensive and combinatorial array of PTMs that modulate Hsp90 chaperone function, a phenomenon termed the "chaperone code."