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

Advancing Molecular Simulations: Merging Physical Models, Experiments, and AI to Tackle Multiscale Complexity

Proteins and protein complexes form adaptable networks that regulate essential biochemical pathways and define cell phenotypes through dynamic mechanisms and interactions. Advances in structural biology and molecular simulations have revealed how protein systems respond to changes in their environments, such as ligand binding, stress conditions, or perturbations like mutations and post-translational modifications, influencing signal transduction and cellular phenotypes. Here, we discuss how computational approaches, ranging from molecular dynamics (MD) simulations to AI-driven methods, are instrumental in studying protein dynamics from isolated molecules to large assemblies. These techniques elucidate conformational landscapes, ligand-binding mechanisms, and protein-protein interactions and are starting to support the construction of multiscale realistic representations of highly complex systems, ranging up to whole cell models. With cryo-electron microscopy, cryo-electron tomography, and AlphaFold accelerating the structural characterization of protein networks, we suggest that integrating AI and Machine Learning with multiscale MD methods will enhance fundamental understating for systems of ever-increasing complexity, usher in exciting possibilities for predictive modeling of the behavior of cell compartments or even whole cells. These advances are indeed transforming biophysics and chemical biology, offering new opportunities to study biomolecular mechanisms at atomic resolution.

Allosteric modulation by the fatty acid site in the glycosylated SARS-CoV-2 spike

The spike protein is essential to the SARS-CoV-2 virus life cycle, facilitating virus entry and mediating viral-host membrane fusion. The spike contains a fatty acid (FA) binding site between every two neighbouring receptor-binding domains. This site is coupled to key regions in the protein, but the impact of glycans on these allosteric effects has not been investigated. Using dynamical nonequilibrium molecular dynamics (D-NEMD) simulations, we explore the allosteric effects of the FA site in the fully glycosylated spike of the SARS-CoV-2 ancestral variant. Our results identify the allosteric networks connecting the FA site to functionally important regions in the protein, including the receptor-binding motif, an antigenic supersite in the N-terminal domain, the fusion peptide region, and another allosteric site known to bind heme and biliverdin. The networks identified here highlight the complexity of the allosteric modulation in this protein and reveal a striking and unexpected link between different allosteric sites. Comparison of the FA site connections from D-NEMD in the glycosylated and non-glycosylated spike revealed that glycans do not qualitatively change the internal allosteric pathways but can facilitate the transmission of the structural changes within and between subunits.

Pathogenic mutation impairs functional dynamics of Hsp60 in mono- and oligomeric states

Mitochondrial chaperonin Heat Shock Protein 60 kDa (Hsp60) oversees the correct folding of client proteins in cooperation with Hsp10. Hsp60 monomers M first form 7-meric Single rings (S), which then pair into 14-meric Double rings (D) that accommodate clients in their lumen. Recruitment of 7 Hsp10 molecules per pole yields a sealed 28-meric Football-shaped complex (F). ATP hydrolysis in each Hsp60 unit drives client folding and F disassembly. The V72I mutation in hereditary spastic paraplegia form SPG13 impairs Hsp60 function despite being distant from the active site. We here investigate this impairment with atomistic molecular dynamics (MD) simulations of M, S, D, and F for both WT and mutant Hsp60, considering catalytic aspartates in D and F in different protonation states (even simulating one such state of D post-hydrolysis). Our findings show that-as observed experimentally-V72I rigidifies Hsp60 assemblies, significantly impacting internal dynamics. In monomers, V72I introduces a new allosteric route that bypasses the ATP binding site and affects mechanisms driving reactivity. These insights highlight a multiscale complexity of Hsp60 that could inspire the design of experiments to better understand both its WT and V72I variants.

Free energy landscape of the PI3Kα C-terminal activation

The gene PIK3CA, encoding the catalytic subunit p110α of PI3Kα, is the second most frequently mutated gene in cancer, with the highest frequency oncogenic mutants occurring in the C-terminus of the kinase domain. The C-terminus has a dual function in regulating the kinase, playing a putative auto-inhibitory role for kinase activity and being absolutely essential for binding to the cell membrane. However, the molecular mechanisms by which these C-terminal oncogenic mutations cause PI3Kα overactivation remain unclear. To understand how a spectrum of C-terminal mutations of PI3Kα alter kinase activity compared to the WT, we perform unbiased and biased Molecular Dynamics simulations of several C-terminal mutants and report the free energy landscapes for the C-terminal "closed-to-open" transition in the WT, H1047R, G1049R, M1043L and N1068KLKR mutants. Results are consistent with HDX-MS experimental data and provide a molecular explanation why H1047R and G1049R reorient the C-terminus with a different mechanism compared to the WT and M1043L and N1068KLKR mutants. Moreover, we show that in the H1047R mutant, the cavity, where the allosteric ligands STX-478 and RLY-2608 bind, is more accessible contrary to the WT. This study provides insights into the molecular mechanisms underlying activation of oncogenic PI3Kα by C-terminal mutations and represents a valuable resource for continued efforts in the development of mutant selective inhibitors as therapeutics.

Phosphorylation-driven epichaperome assembly is a regulator of cellular adaptability and proliferation

The intricate network of protein-chaperone interactions is crucial for maintaining cellular function. Recent discoveries have unveiled the existence of specialized chaperone assemblies, known as epichaperomes, which serve as scaffolding platforms that orchestrate the reconfiguration of protein-protein interaction networks, thereby enhancing cellular adaptability and proliferation. This study explores the structural and regulatory aspects of epichaperomes, with a particular focus on the role of post-translational modifications (PTMs) in their formation and function. A key finding is the identification of specific PTMs on HSP90, particularly at residues Ser226 and Ser255 within an intrinsically disordered region, as critical determinants of epichaperome assembly. Our data demonstrate that phosphorylation of these serine residues enhances HSP90's interactions with other chaperones and co-chaperones, creating a microenvironment conducive to epichaperome formation. Moreover, we establish a direct link between epichaperome function and cellular physiology, particularly in contexts where robust proliferation and adaptive behavior are essential, such as in cancer and pluripotent stem cell maintenance. These findings not only provide mechanistic insights but also hold promise for the development of novel therapeutic strategies targeting chaperone assemblies in diseases characterized by epichaperome dysregulation, thereby bridging the gap between fundamental research and precision medicine.

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

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.

Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence

Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.

Phosphorylation-Driven Epichaperome Assembly: A Critical Regulator of Cellular Adaptability and Proliferation

The intricate protein-chaperone network is vital for cellular function. Recent discoveries have unveiled the existence of specialized chaperone complexes called epichaperomes, protein assemblies orchestrating the reconfiguration of protein-protein interaction networks, enhancing cellular adaptability and proliferation. This study delves into the structural and regulatory aspects of epichaperomes, with a particular emphasis on the significance of post-translational modifications in shaping their formation and function. A central finding of this investigation is the identification of specific PTMs on HSP90, particularly at residues Ser226 and Ser255 situated within an intrinsically disordered region, as critical determinants in epichaperome assembly. Our data demonstrate that the phosphorylation of these serine residues enhances HSP90's interaction with other chaperones and co-chaperones, creating a microenvironment conducive to epichaperome formation. Furthermore, this study establishes a direct link between epichaperome function and cellular physiology, especially in contexts where robust proliferation and adaptive behavior are essential, such as cancer and stem cell maintenance. These findings not only provide mechanistic insights but also hold promise for the development of novel therapeutic strategies targeting chaperone complexes in diseases characterized by epichaperome dysregulation, bridging the gap between fundamental research and precision medicine.