Modeling signal propagation mechanisms and ligand-based conformational dynamics of the Hsp90 molecular chaperone full-length dimer

Exploring Mutation-Driven Changes in the ATP-ADP Conformational Cycle of Human Hsp70 by All-Atom MD Adaptive Sampling

Hsp70 belongs to a family of molecular chaperones ubiquitous through organisms that assist client protein folding and prevent aggregation. It works through a tightly ATP-regulated allosteric cycle mechanism, which organizes its two NBD and SBD into alternate open and closed arrangements that facilitate loading and unloading of client proteins. The two cytosolic human isoforms Hsc70 and HspA1 are relevant targets for neurodegenerative diseases and cancer. Illuminating the molecular details of Hsp70 functional dynamics is essential to rationalize differences among the well-characterized bacterial homologue DnaK and the less explored human forms and develop subtype- or species-selective allosteric drugs. We present here a molecular dynamics-based analysis of the conformational dynamics of HspA1. By using an "allosterically impaired" mutant for comparison, we can reconstruct the impact of the ADP-ATP swap on interdomain contacts and dynamic coordination in full-length HspA1, supporting previous predictions that were, however, limited to the NBD. We model the initial onset of the conformational cycle by proposing a sequence of structural steps, which reveal the role of a specific human sequence insertion at the linker, and a modulation of the angle formed by the two NBD lobes during the progression of docking. Our findings pinpoint functionally relevant conformations and set the basis for a selective structure-based drug discovery approach targeting allosteric sites in human Hsp70.

Insight into the Nucleotide Based Modulation of the Grp94 Molecular Chaperone Using Multiscale Dynamics

Grp94, an ER-localized molecular chaperone, is required for the folding and activation of many membrane and secretory proteins. Client activation by Grp94 is mediated by nucleotide and conformational changes. In this work, we aim to understand how microscopic changes from nucleotide hydrolysis can potentiate large-scale conformational changes of Grp94. We performed all-atom molecular dynamics simulations on the ATP-hydrolysis competent state of the Grp94 dimer in four different nucleotide bound states. We found that Grp94 was the most rigid when ATP was bound. ATP hydrolysis or nucleotide removal enhanced mobility of the N-terminal domain and ATP lid, resulting in suppression of interdomain communication. In an asymmetric conformation with one hydrolyzed nucleotide, we identified a more compact state, similar to experimental observations. We also identified a potential regulatory role of the flexible linker, as it formed electrostatic interactions with the Grp94 M-domain helix near the region where BiP is known to bind. These studies were complemented with normal-mode analysis of an elastic network model to investigate Grp94's large-scale conformational changes. SPM analysis identified residues that are important in signaling conformational change, many of which have known functional relevance in ATP coordination and catalysis, client binding, and BiP binding. Our findings suggest that ATP hydrolysis in Grp94 alters allosteric wiring and facilitates conformational changes.

Sarcin/Ricin Domain RNA Retains Its Structure Better Than A-RNA in Adaptively Biased Molecular Dynamics Simulations

Less than one in thirty of the RNA sequences transcribed in humans are translated into protein. The noncoding RNA (ncRNA) functions in catalysis, structure, regulation, and more. However, for the most part, these functions are poorly characterized. RNA is modular and described by motifs that include helical A-RNA with canonical Watson-Crick base-pairing as well as structures with only noncanonical base pairs. Understanding the structure and dynamics of motifs will aid in deciphering functions of specific ncRNAs. We present computational studies on a standard sarcin/ricin domain (SRD), citrus bark cracking viroid SRD, as well as A-RNA. We have applied enhanced molecular dynamics techniques that construct an inverse free-energy surface (iFES) determined by collective variables that monitor base-pairing and backbone conformation. Each SRD RNA is flanked on each side by A-RNA, allowing comparison of the behavior of these motifs in the same molecule. The RNA iFESs have single peaks, indicating that the combined motifs should denature as a single cohesive unit, rather than by regional melting. Local root-mean-square deviation (RMSD) analysis and communication propensity (CProp, variance in distances between residue pairs) reveal distinct motif properties. Our analysis indicates that the standard SRD is more stable than the viroid SRD, which is more stable than A-RNA. Base pairs at SRD to A-RNA transitions have limited flexibility. Application of CProp reveals extraordinary stiffness of the SRD, allowing residues on opposite sides of the motif to sense each other's motions.

Mebendazole's Conformational Space and Its Predicted Binding to Human Heat-Shock Protein 90

Recent experimental evidence suggests that mebendazole, a popular antiparasitic drug, binds to heat shock protein 90 (Hsp90) and inhibits acute myeloid leukemia cell growth. In this study we use quantum mechanics (QM), molecular similarity, and molecular dynamics (MD) calculations to predict possible binding poses of mebendazole to the adenosine triphosphate (ATP) binding site of Hsp90. Extensive conformational searches and minimization of the five mebendazole tautomers using the MP2/aug-cc-pVTZ theory level resulted in 152 minima. Mebendazole-Hsp90 complex models were subsequently created using the QM optimized conformations and protein coordinates obtained from experimental crystal structures that were chosen through similarity calculations. Nine different poses were identified from a total of 600 ns of explicit solvent, all-atom MD simulations using two different force fields. All simulations support the hypothesis that mebendazole is able to bind to the ATP binding site of Hsp90.

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.

Induced fit versus conformational selection: From rate constants to fluxes… and back to rate constants

Induced fit- (IF) and conformational selection (CS) binding mechanisms have long been regarded to be mutually exclusive. Yet, they are now increasingly considered to produce the final ligand-target complex alongside within a thermodynamic cycle. This viewpoint benefited from the introduction of binding fluxes as a tool for analyzing the overall behavior of such cycle. This study aims to provide more vivid and applicable insights into this emerging field. In this respect, combining differential equation- based simulations and hitherto little explored alternative modes of calculation provide concordant information about the intricate workings of such cycle. In line with previous reports, we observe that the relative contribution of IF increases with the ligand concentration at equilibrium. Yet the baseline contribution may vary from one case to another and simulations as well as calculations show that this parameter is essentially regulated by the dissociation rate of both pathways. Closer attention should be paid to how the contributions of IF and CS compare at physiologically relevant drug/ligand concentrations. To this end, a simple equation discloses how changing a limited set of "microscopic" rate constants can extend the concentration range at which CS contributes most effectively. Finally, it could also be beneficial to extend the utilization of flux- based approaches to more physiologically relevant time scales and alternative binding models.

Targeting the mitochondrial chaperone TRAP1: strategies and therapeutic perspectives

TRAP1, the mitochondrial isoform of heat shock protein (Hsp)90 chaperones, is a key regulator of metabolism and organelle homeostasis in diverse pathological states. While selective TRAP1 targeting is an attractive goal, classical active-site-directed strategies have proved difficult, due to high active site conservation among Hsp90 paralogs. Here, we discuss advances in developing TRAP1-directed strategies, from lead modification with mitochondria delivery groups to the computational discovery of allosteric sites and ligands. Specifically, we address the unique opportunities that targeting TRAP1 opens up in tackling fundamental questions on its biology and in unveiling new therapeutic approaches. Finally, we show how crucial to this endeavor is our ability to predict the activities of TRAP1-selective allosteric ligands and to optimize target engagement to avoid side effects.

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.

Discovery of a covalent inhibitor of heat shock protein 90 with antitumor activity that blocks the co-chaperone binding via C-terminal modification

Heat shock protein (Hsp90), a critical molecular chaperone that regulates the maturation of a large number of oncogenic client proteins, plays an essential role in the growth of neoplastic cells. Herein, DDO-6600 is identified to covalent modification of Cys598 on Hsp90 from in silico study and is verified by a series of biological assays. We demonstrated that DDO-6600 covalently bound to Cys598 on the Hsp90 C terminus and exhibited antiproliferative activities against multiple tumor cells without inhibiting ATPase activity. Further studies showed that DDO-6600 disrupted the interaction between Hsp90 and Cdc37, which induced the degradation of kinase client proteins in multiple tumor cell lines, promoted apoptosis, and inhibited cell motility. Our findings offer mechanic insights into the covalent modification of Hsp90 and provide an alternative strategy for the development of Hsp90 covalent regulators or chemical probes to explore the therapeutical potential of Hsp90.

Structural Communication between the E. coli Chaperones DnaK and Hsp90

The 70 kDa and 90 kDa heat shock proteins Hsp70 and Hsp90 are two abundant and highly conserved ATP-dependent molecular chaperones that participate in the maintenance of cellular homeostasis. In Escherichia coli, Hsp90 (Hsp90Ec) and Hsp70 (DnaK) directly interact and collaborate in protein remodeling. Previous work has produced a model of the direct interaction of both chaperones. The locations of the residues involved have been confirmed and the model has been validated. In this study, we investigate the allosteric communication between Hsp90Ec and DnaK and how the chaperones couple their conformational cycles. Using elastic network models (ENM), normal mode analysis (NMA), and a structural perturbation method (SPM) of asymmetric and symmetric DnaK-Hsp90Ec, we extract biologically relevant vibrations and identify residues involved in allosteric signaling. When one DnaK is bound, the dominant normal modes favor biological motions that orient a substrate protein bound to DnaK within the substrate/client binding site of Hsp90Ec and release the substrate from the DnaK substrate binding domain. The presence of one DnaK molecule stabilizes the entire Hsp90Ec protomer to which it is bound. Conversely, the symmetric model of DnaK binding results in steric clashes of DnaK molecules and suggests that the Hsp90Ec and DnaK chaperone cycles operate independently. Together, this data supports an asymmetric binding of DnaK to Hsp90Ec.

Machine Learning of Allosteric Effects: The Analysis of Ligand-Induced Dynamics to Predict Functional Effects in TRAP1

Allosteric molecules provide a powerful means to modulate protein function. However, the effect of such ligands on distal orthosteric sites cannot be easily described by classical docking methods. Here, we applied machine learning (ML) approaches to expose the links between local dynamic patterns and different degrees of allosteric inhibition of the ATPase function in the molecular chaperone TRAP1. We focused on 11 novel allosteric modulators with similar affinities to the target but with inhibitory efficacy between the 26.3 and 76%. Using a set of experimentally related local descriptors, ML enabled us to connect the molecular dynamics (MD) accessible to ligand-bound (perturbed) and unbound (unperturbed) systems to the degree of ATPase allosteric inhibition. The ML analysis of the comparative perturbed ensembles revealed a redistribution of dynamic states in the inhibitor-bound versus inhibitor-free systems following allosteric binding. Linear regression models were built to quantify the percentage of experimental variance explained by the predicted inhibitor-bound TRAP1 states. Our strategy provides a comparative MD–ML framework to infer allosteric ligand functionality. Alleviating the time scale issues which prevent the routine use of MD, a combination of MD and ML represents a promising strategy to support in silico mechanistic studies and drug design.

Dissecting Molecular Principles of the Hsp90 Chaperone Regulation by Allosteric Modulators Using a Hierarchical Simulation Approach and Network Modeling of Allosteric Interactions: Conformational Selection Dictates the Diversity of Protein Responses and Ligand-Specific Functional Mechanisms

Conformational plasticity of the Hsp90 molecular chaperones underlies the diversity of functional mechanisms that these versatile molecular machines employ to coordinate their vast protein clientele in the cellular environment. Despite a steady progress in studies of the Hsp90 machinery, a great deal remains unknown about molecular principles and ligand-specific functional mechanisms of the Hsp90 regulation by allosteric modulators that attracted significant attention because of their therapeutic potential. Due to structural complexity and dynamic nature of the Hsp90 responses to allosteric modulators, the atomistic details about the mode of action of these small molecules continue to be fairly scarce and controversial. In this work, we employ an integrative strategy that encompassed atomistic simulations of the Hsp90 proteins and hierarchical modeling of Hsp90-ligand binding with network analysis to explore functional mechanisms of the Hsp90 regulation by a panel of allosteric modulators (novobiocin, KU-135, KU-174, and KU-32) with different models of action. The results show that functional mechanisms of allosteric modulation in the Hsp90 proteins may be driven by conformational selection principles in which ligands elicit pre-existing states of the unbound chaperone to drive ligand-specific protein responses and distinct scenarios of Hsp90 regulation. We found that novobiocin can selectively sequester an ensemble of open chaperone conformations and inhibit the progression of the functional cycle through a cascade of cumulative dynamic changes. In contrast, KU-32 displayed unique preferences toward partially closed dynamic states, inducing robust allosteric signaling and stimulation of the ATPase cycle. The proposed model of the Hsp90 regulation by allosteric modulators reconciled diverse experimental data and showed that allosteric modulators may operate via targeted exploitation of dynamic landscapes eliciting vastly different protein responses and diverse mechanisms of action.

Dynamically Shaping Chaperones. Allosteric Modulators of HSP90 Family as Regulatory Tools of Cell Metabolism in Neoplastic Progression

Molecular chaperones have recently emerged as fundamental regulators of salient biological routines, including metabolic adaptations to environmental changes. Yet, many of the molecular mechanisms at the basis of their functions are still unknown or at least uncertain. This is in part due to the lack of chemical tools that can interact with the chaperones to induce measurable functional perturbations. In this context, the use of small molecules as modulators of protein functions has proven relevant for the investigation of a number of biomolecular systems. Herein, we focus on the functions, interactions and signaling pathways of the HSP90 family of molecular chaperones as possible targets for the discovery of new molecular entities aimed at tuning their activity and interactions. HSP90 and its mitochondrial paralog, TRAP1, regulate the activity of crucial metabolic circuitries, making cells capable of efficiently using available energy sources, with relevant implications both in healthy conditions and in a variety of disease states and especially cancer. The design of small-molecules targeting the chaperone cycle of HSP90 and able to inhibit or stimulate the activity of the protein can provide opportunities to finely dissect their biochemical activities and to obtain lead compounds to develop novel, mechanism-based drugs.

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.

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

Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.

Ligand Binding, Unbinding, and Allosteric Effects: Deciphering Small-Molecule Modulation of HSP90

The molecular chaperone HSP90 oversees the functional activation of a large number of client proteins. Because of its role in multiple pathways linked to cancer and neurodegeneration, drug discovery targeting HSP90 has been actively pursued. Yet, a number of inhibitors failed to meet expectations due to induced toxicity problems. In this context, allosteric perturbation has emerged as an alternative strategy for the pharmacological modulation of HSP90 functions. Specifically, novel allosteric stimulators showed the interesting capability of accelerating HSP90 closure dynamics and ATPase activities while inducing tumor cell death. Here, we gain atomistic insight into the mechanisms of allosteric ligand recognition and their consequences on the functional dynamics of HSP90, starting from the fully unbound state. We integrate advanced computational sampling methods based on FunnelMetadynamics, with the analysis of internal dynamics of the structural ensembles visited during the simulations. We observe several binding/unbinding events, and from these, we derive an accurate estimation of the absolute binding free energy. Importantly, we show that different binding poses induce different dynamics states. Our work for the first time explicitly correlates HSP90 responses to binding/unbinding of an allosteric ligand to the modulation of functionally oriented protein motions.

Structure, Function, and Regulation of the Hsp90 Machinery

Heat shock protein 90 (Hsp90) is a molecular chaperone involved in the maturation of a plethora of substrates ("clients"), including protein kinases, transcription factors, and E3 ubiquitin ligases, positioning Hsp90 as a central regulator of cellular proteostasis. Hsp90 undergoes large conformational changes during its ATPase cycle. The processing of clients by cytosolic Hsp90 is assisted by a cohort of cochaperones that affect client recruitment, Hsp90 ATPase function or conformational rearrangements in Hsp90. Because of the importance of Hsp90 in regulating central cellular pathways, strategies for the pharmacological inhibition of the Hsp90 machinery in diseases such as cancer and neurodegeneration are being developed. In this review, we summarize recent structural and mechanistic progress in defining the function of organelle-specific and cytosolic Hsp90, including the impact of individual cochaperones on the maturation of specific clients and complexes with clients as well as ways of exploiting Hsp90 as a drug target.

Autoencoder-Based Detection of Dynamic Allostery Triggered by Ligand Binding Based on Molecular Dynamics

Dynamic allostery on proteins, triggered by regulator binding or chemical modifications, transmits information from the binding site to distant regions, dramatically altering protein function. It is accompanied by subtle changes in side-chain conformations of the protein, indicating that the changes in dynamics, and not rigid or large conformational changes, are essential to understand regulation of protein function. Although a lot of experimental and theoretical studies have been dedicated to investigate this issue, the regulation mechanism of protein function is still being debated. Here, we propose an autoencoder-based method that can detect dynamic allostery. The method is based on the comparison of time fluctuations of protein structures, in the form of distance matrices, obtained from molecular dynamics simulations in ligand-bound and -unbound forms. Our method detected that the changes in dynamics by ligand binding in the PDZ2 domain led to the reorganization of correlative fluctuation motions among residue pairs, which revealed a different view of the correlated motions from the PCA and DCCM. In addition, other correlative motions were also found as a result of the dynamic perturbation from the ligand binding, which may lead to dynamic allostery. This autoencoder-based method would be usefully applied to the signal transduction and mutagenesis systems involved in protein functions and severe diseases.

In silico survey of the central conserved regions in viroids of the Pospiviroidae family for conserved asymmetric loop structures

Viroids are the smallest replicative pathogens, consisting of RNA circles (∼300 nucleotides) that require host machinery to replicate. Structural RNA elements recruit these host factors. Currently, many of these structural elements and the nature of their interactions are unknown. All Pospiviroidae have homology in the central conserved region (CCR). The CCR of potato spindle tuber viroid (PSTVd) contains a sarcin/ricin domain (SRD), the only viroid structural element with an unequivocal replication role. We assumed that every member of this family uses this region to recruit host factors, and that each CCR has an SRD-like asymmetric loop within it. Potential SRD or SRD-like motifs were sought in the CCR of each Pospiviroidae member as follows. Motif location in each CCR was predicted with MUSCLE alignment and Vienna RNAfold. Viroid-specific models of SRD-like motifs were built by superimposing noncanonical base pairs and nucleotides on a model of an SRD. The RNA geometry search engine FR3D was then used to find nucleotide groups close to the geometry suggested by this superimposition. Atomic resolution structures were assembled using the molecular visualization program Chimera, and the stability of each motif was assessed with molecular dynamics (MD). Some models required a protonated cytosine. To be stable within a cell, the pK_a of that cytosine must be shifted up. Constant pH-replica exchange MD analysis showed such a shift in the proposed structures. These data show that every Pospiviroidae member could form a motif that resembles an SRD in its CCR, and imply there could be undiscovered mimics of other RNA domains.

Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example

Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Modulation of Human Hsp90α Conformational Dynamics by Allosteric Ligand Interaction at the C-Terminal Domain

Recent years have seen heat shock protein 90 kDa (Hsp90) attract significant interest as a viable drug target, particularly for cancer. To date, designed inhibitors that target the ATPase domain demonstrate potent anti-proliferative effects, but have failed clinical trials due to high levels of associated toxicity. To circumvent this, the focus has shifted away from the ATPase domain. One option involves modulation of the protein through allosteric activation/inhibition. Here, we propose a novel approach: we use previously obtained information via residue perturbation scanning coupled with dynamic residue network analysis to identify allosteric drug targeting sites for inhibitor docking. We probe the open conformation of human Hsp90α for druggable sites that overlap with these allosteric control elements, and identify three putative natural compound allosteric modulators: Cephalostatin 17, 20(29)-Lupene-3β-isoferulate and 3'-Bromorubrolide F. We assess the allosteric potential of these ligands by examining their effect on the conformational dynamics of the protein. We find evidence for the selective allosteric activation and inhibition of Hsp90's conformational transition toward the closed state in response to ligand binding and shed valuable insight to further the understanding of allosteric drug design and Hsp90's complex allosteric mechanism of action.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.

Photoprogramming Allostery in Human Serum Albumin

Developing strategies to interfere with allosteric interactions in proteins not only promises to deepen our understanding of vital cellular processes but also allows their regulation using external triggers. Light is particularly attractive as a trigger being spatiotemporally selective and compatible with the physiological environment. Here, we engineered a hybrid protein in which irradiation with light opens a new allosteric communication route that is not inherent to the natural system. We select human serum albumin, a promiscuous protein responsible for transporting a variety of ligands in plasma, and show that by covalently incorporating a synthetic photoswitch to subdomain IA we achieve optical control of the ligand binding in subdomain IB. Molecular dynamics simulations confirm the allosteric nature of the interactions between IA and IB in the engineered protein. Specifically, upon illumination, photoconversion of the switch is found to correlate with a less-coordinated motion of the two subdomains and an increased flexibility of the binding pocket in subdomain IB, whose fluctuations are cooperatively enhanced by the presence of ligands, ultimately facilitating their release. Our combined experimental and computational work demonstrates how harnessing artificial molecular switches enables photoprogramming the allosteric regulation of binding activities in such a prominent protein.

Dual Roles for Yeast Sti1/Hop in Regulating the Hsp90 Chaperone Cycle

The Hsp90 chaperone is regulated by many cochaperones that tune its activities, but how they act to coordinate various steps in the reaction cycle is unclear. The primary role of Saccharomyces cerevisiae Hsp70/Hsp90 cochaperone Sti1 (Hop in mammals) is to bridge Hsp70 and Hsp90 to facilitate client transfer. Sti1 is not essential, so Hsp90 can interact with Hsp70 in vivo without Sti1. Nevertheless, many Hsp90 mutations make Sti1 necessary. We noted that Sti1-dependent mutations cluster in regions proximal to N-terminal domains (SdN) or C-terminal domains (SdC), which are known to be important for interaction with Hsp70 or clients, respectively. To uncover mechanistic details of Sti1-Hsp90 cooperation, we identified intramolecular suppressors of the Hsp90 mutants and assessed their physical, functional, and genetic interactions with Hsp70, Sti1, and other cochaperones. Our findings suggest Hsp90 SdN and SdC mutants depend on the same interaction with Sti1, but for different reasons. Sti1 promoted an essential Hsp70 interaction in the SdN region and supported SdC-region function by establishing an Hsp90 conformation crucial for capturing clients and progressing through the reaction cycle. We find the Hsp70 interaction and relationship with Sti1/Hop is conserved in the human Hsp90 system. Our work consolidates and clarifies much structural, biochemical, and computational data to define in vivo roles of Sti1/Hop in coordinating Hsp70 binding and client transfer with progression of the Hsp90 reaction cycle.

Dynamics-based community analysis and perturbation response scanning of allosteric interaction networks in the TRAP1 chaperone structures dissect molecular linkage between conformational asymmetry and sequential ATP hydrolysis

Allosteric interactions of the Hsp90 chaperones with cochaperones and diverse protein clients can often exhibit distinct asymmetric features that determine regulatory mechanisms and cellular functions in many signaling networks. The recent crystal structures of the mitochondrial Hsp90 isoform TRAP1 in complexes with ATP analogs have provided first evidence of significant asymmetry in the closed dimerized state that triggers independent activity of the chaperone protomers, whereby preferential hydrolysis of the buckled protomer is followed by conformational flipping between protomers and hydrolysis of the second protomer. Despite significant insights in structural characterizations of the TRAP1 chaperone, the atomistic details and mechanics of allosteric interactions that couple sequential ATP hydrolysis with asymmetric conformational switching in the TRAP1 protomers remain largely unknown. In this work, we explored atomistic and coarse-grained simulations of the TRAP1 dimer structures in combination with the ensemble-based network modeling and perturbation response scanning of residue interaction networks to probe salient features underlying allosteric signaling mechanism. This study has revealed that key effector sites that orchestrate allosteric interactions occupy the ATP binding region and N-terminal interface of the buckled protomer, whereas the main sensors of allosteric signals that drive functional conformational changes during ATPase cycle are consolidated near the client binding region of the straight protomer, channeling the energy of ATP hydrolysis for client remodeling. The community decomposition analysis of the interaction networks and reconstruction of allosteric communication pathways in the TRAP1 structures have quantified mechanism of allosteric regulation, revealing control points and interactions that coordinate asymmetric switching during ATP hydrolysis.

A switch point in the molecular chaperone Hsp90 responding to client interaction

Heat shock protein 90 (Hsp90) is a dimeric molecular chaperone that undergoes large conformational changes during its functional cycle. It has been established that conformational switch points exist in the N-terminal (Hsp90-N) and C-terminal (Hsp90-C) domains of Hsp90, however information for switch points in the large middle-domain (Hsp90-M) is scarce. Here we report on a tryptophan residue in Hsp90-M as a new type of switch point. Our study shows that this conserved tryptophan senses the interaction of Hsp90 with a stringent client protein and transfers this information via a cation–π interaction with a neighboring lysine. Mutations at this position hamper the communication between domains and the ability of a client protein to affect the Hsp90 cycle. The residue thus allows Hsp90 to transmit information on the binding of a client from Hsp90-M to Hsp90-N which is important for progression of the conformational cycle and the efficient processing of client proteins.

The heat shock protein 90 (Hsp90) chaperone undergoes large conformational changes during its functional cycle. Here the authors combine in vivo, biochemical, biophysical and computational approaches and provide insights into the allosteric regulation of Hsp90 by identifying and characterizing a switch point in the Hsp90 middle domain.

Regulation of the Hsp90 system

Hsp90 is a highly conserved and abundant chaperone. It participates in essential cellular activities by supporting the maturation process of its client proteins, many of which are protein kinases and steroid receptors. Client processing is achieved via extensive conformational changes within the dimeric chaperone. This requires an ATP hydrolysis activity that is controlled by auto-inhibitory mechanisms and several structurally diverse cofactors. Especially the client-specificity of Hsp90 depends on client-specific cofactors, which can adapt Hsp90's activities to the client requirements at different conditions and in different cell types. Additionally, post-translational modifications can influence almost every aspect of Hsp90's interactions and activities. In this review, we present these regulatory principles, discuss the factors that have an impact on Hsp90's function and elaborate the mechanisms that are responsible for regulating the Hsp90 machinery.

Allosteric Modulation of Human Hsp90α Conformational Dynamics

Central to Hsp90's biological function is its ability to interconvert between various conformational states. Drug targeting of Hsp90's regulatory mechanisms, including its modulation by cochaperone association, presents as an attractive therapeutic strategy for Hsp90 associated pathologies. In this study, we utilized homology modeling techniques to calculate full-length structures of human Hsp90α in closed and partially open conformations and used these structures as a basis for several molecular dynamics based analyses aimed at elucidating allosteric mechanisms and modulation sites in human Hsp90α. Atomistic simulations demonstrated that bound adenosine triphosphate (ATP) stabilizes the dimer by "tensing" each protomer, while adenosine diphosphate (ADP) and apo configurations "relax" the complex by increasing global flexibility, the former case resulting in a fully open "v-like" conformation. Dynamic residue network analysis revealed regions of the protein involved in intraprotein communication and identified several key communication hubs that correlate with known functional sites. Pairwise comparison of betweenness centrality, shortest path, and residue fluctuations revealed that a proportional relationship exists between the latter two measurables and an inverse relationship between these two and betweenness centrality. This analysis showed how protein flexibility, degree of compactness, and the distance cutoff used for network construction influence the correlations between these metrics. These findings are novel and suggest shortest path and betweenness centrality to be more relevant quantities to follow for detecting functional residues in proteins compared to residue fluctuations. Perturbation response scanning analysis identified several potential residue sites capable of modulating conformational change in favor of interstate conversion. For the ATP-bound open conformation, these sites were found to overlap with known Aha1 and client binding sites, demonstrating how naturally occurring forces associated with cofactor binding could allosterically modulate conformational dynamics.

Phosphorylation induced cochaperone unfolding promotes kinase recruitment and client class-specific Hsp90 phosphorylation

During the Hsp90-mediated chaperoning of protein kinases, the core components of the machinery, Hsp90 and the cochaperone Cdc37, recycle between different phosphorylation states that regulate progression of the chaperone cycle. We show that Cdc37 phosphorylation at Y298 results in partial unfolding of the C-terminal domain and the population of folding intermediates. Unfolding facilitates Hsp90 phosphorylation at Y197 by unmasking a phosphopeptide sequence, which serves as a docking site to recruit non-receptor tyrosine kinases to the chaperone complex via their SH2 domains. In turn, Hsp90 phosphorylation at Y197 specifically regulates its interaction with Cdc37 and thus affects the chaperoning of only protein kinase clients. In summary, we find that by providing client class specificity, Hsp90 cochaperones such as Cdc37 do not merely assist in client recruitment but also shape the post-translational modification landscape of Hsp90 in a client class-specific manner.

The Hsp90 chaperone cycle is influenced by multiple phosphorylation events but their regulatory functions are poorly understood. Here, the authors show that phosphorylation and unfolding of cochaperone Cdc37 tailors the Hsp90 chaperone cycle by recruiting kinases that promote distinct phosphorylation patterns.

Dynamics revelation of conformational changes and binding modes of heat shock protein 90 induced by inhibitor associations

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment. The current work integrates molecular dynamics (MD) simulations, calculations of binding free energy, and principal component (PC) analysis with scanning of inhibitor–residue interaction to probe the binding modes of inhibitors YK9, YKJ and YKI to Hsp90 and identify the hot spot of the inhibitor–Hsp90 binding. The results suggest that the introductions of two groups G1 and G2 into YKJ and YKI strengthen the binding ability of YKJ and YKI to Hsp90 compared to YK9. PC analysis based MD trajectories prove that inhibitor bindings exert significant effects on the conformational changes, internal dynamics and motion modes of Hsp90, especially for the helix α2 and the loops L1 and L2. The calculations of residue-based free energy decomposition and scanning of the inhibitor–Hsp90 interaction suggest that six residues L107, G108, F138, Y139, W162 and F170 construct the common hot spot of the inhibitor–residue interactions. Moreover the substitutions of the groups G1 and G2 in YKJ and YKI lead to two additional hydrogen bonding interactions and multiple hydrophobic interactions for bindings of YKJ and YKI to Hsp90. This work is also expected to contribute theoretical hints for the design of potent inhibitors toward Hsp90.

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment.

Computational Modeling of the Hsp90 Interactions with Cochaperones and Small-Molecule Inhibitors

Allosteric interactions of the molecular chaperone Hsp90 with a diverse array of cochaperones and client proteins, such as protein kinases and transcription factors, allow for efficient molecular communication in signal transduction networks. Deregulation of pathways involving these proteins is commonly associated with cancer pathologies and allosteric inhibition of oncogenic clients by targeting Hsp90 provides a powerful therapeutic strategy in cancer research. We review several validated computational approaches and tools used in the studies of the Hsp90 interactions with proteins and small molecules. These methods include experimentally guided docking to predict Hs90-protein interactions, molecular and binding free energy simulations to analyze Hsp90 binding with small molecules, and structure-based network modeling to evaluate allosteric interactions and communications in the Hsp90 regulatory complexes. Through the lens of allosteric-centric view on Hsp90 function and regulation, we discuss newly emerging computational tools that link protein structure modeling with biophysical simulations and network-based systems biology approaches.

Protein conformational flexibility modulates kinetics and thermodynamics of drug binding

Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.

Atomistic simulations and network-based modeling of the Hsp90-Cdc37 chaperone binding with Cdk4 client protein: A mechanism of chaperoning kinase clients by exploiting weak spots of intrinsically dynamic kinase domains

A fundamental role of the Hsp90 and Cdc37 chaperones in mediating conformational development and activation of diverse protein kinase clients is essential in signal transduction. There has been increasing evidence that the Hsp90-Cdc37 system executes its chaperoning duties by recognizing conformational instability of kinase clients and modulating their folding landscapes. The recent cryo-electron microscopy structure of the Hsp90-Cdc37-Cdk4 kinase complex has provided a framework for dissecting regulatory principles underlying differentiation and recruitment of protein kinase clients to the chaperone machinery. In this work, we have combined atomistic simulations with protein stability and network-based rigidity decomposition analyses to characterize dynamic factors underlying allosteric mechanism of the chaperone-kinase cycle and identify regulatory hotspots that control client recognition. Through comprehensive characterization of conformational dynamics and systematic identification of stabilization centers in the unbound and client- bound Hsp90 forms, we have simulated key stages of the allosteric mechanism, in which Hsp90 binding can induce instability and partial unfolding of Cdk4 client. Conformational landscapes of the Hsp90 and Cdk4 structures suggested that client binding can trigger coordinated dynamic changes and induce global rigidification of the Hsp90 inter-domain regions that is coupled with a concomitant increase in conformational flexibility of the kinase client. This process is allosteric in nature and can involve reciprocal dynamic exchanges that exert global effect on stability of the Hsp90 dimer, while promoting client instability. The network-based rigidity analysis and emulation of thermal unfolding of the Cdk4-cyclin D complex and Hsp90-Cdc37-Cdk4 complex revealed weak spots of kinase instability that are present in the native Cdk4 structure and are targeted by the chaperone during client recruitment. Our findings suggested that this mechanism may be exploited by the Hsp90-Cdc37 chaperone to recruit and protect intrinsically dynamic kinase clients from degradation. The results of this investigation are discussed and interpreted in the context of diverse experimental data, offering new insights into mechanisms of chaperone regulation and binding.

Cooperative Nucleotide Binding in Hsp90 and Its Regulation by Aha1

The function of the molecular chaperone Hsp90 depends on large conformational changes, the rearrangement of local motifs, and the binding and hydrolysis of ATP. The size and complexity of the Hsp90 system impedes the detailed investigation of their interplay using standard methods. To overcome this limitation, we developed a three-color single-molecule FRET assay to study the interaction of Hsp90 with a fluorescently labeled reporter nucleotide in detail. It allows us to directly observe the cooperativity between the two nucleotide binding pockets in the protein dimer. Furthermore, our approach disentangles the protein conformation and the nucleotide binding state of Hsp90 and extracts the kinetics of the state transitions. Thereby, we can identify the kinetic causes mediating the cooperativity. We find that the presence of the first nucleotide prolongs the binding of the second nucleotide to Hsp90. In addition, we observe changes in the kinetics for both the open and the closed conformation of Hsp90 in dependence on the number of occupied nucleotide binding sites. Our analysis also reveals how the co-chaperone Aha1, known to accelerate Hsp90's ATPase activity, affects those transitions in a nucleotide-dependent and independent manner, thereby adding another layer of regulation to Hsp90.

Patterns of coevolving amino acids unveil structural and dynamical domains

Patterns of interacting amino acids are so preserved within protein families that the sole analysis of evolutionary comutations can identify pairs of contacting residues. It is also known that evolution conserves functional dynamics, i.e., the concerted motion or displacement of large protein regions or domains. Is it, therefore, possible to use a pure sequence-based analysis to identify these dynamical domains? To address this question, we introduce here a general coevolutionary coupling analysis strategy and apply it to a curated sequence database of hundreds of protein families. For most families, the sequence-based method partitions amino acids into a few clusters. When viewed in the context of the native structure, these clusters have the signature characteristics of viable protein domains: They are spatially separated but individually compact. They have a direct functional bearing too, as shown for various reference cases. We conclude that even large-scale structural and functionally related properties can be recovered from inference methods applied to evolutionary-related sequences. The method introduced here is available as a software package and web server (spectrus.sissa.it/spectrus-evo_webserver).

Role of Terahertz (THz) Fluctuations in the Allosteric Properties of the PDZ Domains

With the aim of investigating the relationship between the fast fluctuations of proteins and their allosteric behavior, we perform molecular dynamics simulations of two model PDZ domains with differential allosteric responses. We focus on protein dynamics in the THz regime (0.1–3 THz) as opposed to lower frequencies. By characterizing the dynamic modulation of the protein backbone induced by ligand binding in terms of single residue and pairwise distance fluctuations, we identify a response nucleus modulated by the ligand that is visible only at THz frequencies. The residues of this nucleus undergo a significant stiffening and an increase in mutual coordination upon binding. Additionally, we find that the dynamic modulation is significantly more intense for the side chains, where it is also redistributed to distal regions not immediately in contact with the ligand allowing us to better define the response nucleus at THz frequencies. The overlap between the known allosterically responding residues of the investigated PDZ domains and the modulated region highlighted here suggests that fast THz dynamics could play a role in allosteric mechanisms.

A protease-resistant Escherichia coli asparaginase with outstanding stability and enhanced anti-leukaemic activity in vitro

L-Asparaginases (ASNases) have been used as first line drugs for paediatric Acute Lymphoblastic Leukaemia (ALL) treatment for more than 40 years. Both the Escherichia coli (EcAII) and Erwinia chrysanthemi (ErAII) type II ASNases currently used in the clinics are characterized by high in vivo instability, short half-life and the requirement of several administrations to obtain a pharmacologically active concentration. Moreover, they are sensitive to proteases (cathepsin B and asparagine endopeptidase) that are over-expressed by resistant leukaemia lymphoblasts, thereby impairing drug activity and pharmacokinetics. Herein, we present the biochemical, structural and in vitro antiproliferative characterization of a new EcAII variant, N24S. The mutant shows completely preserved asparaginase and glutaminase activities, long-term storage stability, improved thermal parameters, and outstanding resistance to proteases derived from leukaemia cells. Structural analysis demonstrates a modification in the hydrogen bond network related to residue 24, while Normal Mode-based geometric Simulation and Molecular Dynamics predict a general rigidification of the monomer as compared to wild-type. These improved features render N24S a potential alternative treatment to reduce the number of drug administrations in vivo and to successfully address one of the major current challenges of ALL treatment: spontaneous, protease-dependent and immunological inactivation of ASNase.

Allosteric Communication Networks in Proteins Revealed through Pocket Crosstalk Analysis

The detection and characterization of binding pockets and allosteric communication in proteins is crucial for studying biological regulation and performing drug design. Nowadays, ever-longer molecular dynamics (MD) simulations are routinely used to investigate the spatiotemporal evolution of proteins. Yet, there is no computational tool that can automatically detect all the pockets and potential allosteric communication networks along these extended MD simulations. Here, we use a novel and fully automated algorithm that examines pocket formation, dynamics, and allosteric communication embedded in microsecond-long MD simulations of three pharmaceutically relevant proteins, namely, PNP, A2A, and Abl kinase. This dynamic analysis uses pocket crosstalk, defined as the temporal exchange of atoms between adjacent pockets, along the MD trajectories as a fingerprint of hidden allosteric communication networks. Importantly, this study indicates that dynamic pocket crosstalk analysis provides new mechanistic understandings on allosteric communication networks, enriching the available experimental data. Thus, our results suggest the prospective use of this unprecedented dynamic analysis to characterize transient binding pockets for structure-based drug design.

Design of Elastic Networks with Evolutionary Optimized Long-Range Communication as Mechanical Models of Allosteric Proteins

Allosteric effects often underlie the activity of proteins, and elucidating generic design aspects and functional principles unique to allosteric phenomena represent a major challenge. Here an approach consisting of the in silico design of synthetic structures, which, as the principal element of allostery, encode dynamical long-range coupling among two sites, is presented. The structures are represented by elastic networks, similar to coarse-grained models of real proteins. A strategy of evolutionary optimization was implemented to iteratively improve allosteric coupling. In the designed structures, allosteric interactions were analyzed in terms of strain propagation, and simple pathways that emerged during evolution were identified as signatures through which long-range communication was established. Moreover, robustness of allosteric performance with respect to mutations was demonstrated. As it turned out, the designed prototype structures reveal dynamical properties resembling those found in real allosteric proteins. Hence, they may serve as toy models of complex allosteric systems, such as proteins. Application of the developed modeling scheme to the allosteric transition in the myosin V molecular motor was also demonstrated.