Functional Hotspots Research Articles

Identifying Functional Hotspot Residues for Biased Ligand Design in G-Protein-Coupled Receptors.

G-protein-coupled receptors (GPCRs) mediate multiple signaling pathways in the cell, depending on the agonist that activates the receptor and multiple cellular factors. Agonists that show higher potency to specific signaling pathways over others are known as "biased agonists" and have been shown to have better therapeutic index. Although biased agonists are desirable, their design poses several challenges to date. The number of assays to identify biased agonists seems expensive and tedious. Therefore, computational methods that can reliably calculate the possible bias of various ligands ahead of experiments and provide guidance, will be both cost and time effective. In this work, using the mechanism of allosteric communication from the extracellular region to the intracellular transducer protein coupling region in GPCRs, we have developed a computational method to calculate ligand bias ahead of experiments. We have validated the method for several β-arrestin-biased agonists in β2-adrenergic receptor (β2AR), serotonin receptors 5-HT1B and 5-HT2B and for G-protein-biased agonists in the κ-opioid receptor. Using this computational method, we also performed a blind prediction followed by experimental testing and showed that the agonist carmoterol is β-arrestin-biased in β2AR. Additionally, we have identified amino acid residues in the biased agonist binding site in both β2AR and κ-opioid receptors that are involved in potentiating the ligand bias. We call these residues functional hotspots, and they can be used to derive pharmacophores to design biased agonists in GPCRs.

Systematic prioritization of functional hotspot in RIG-1 domains using pattern based conventional molecular dynamic simulation

BackgroundRetinoic acid inducible gene 1 (RIG-1), multi-domain protein has a role-play in detecting viral nucleic acids and stimulates the antiviral response. Dysfunction of this protein due to mutations makes the route vulnerable to viral diseases. AimIdentification of functional hotspots that maintains conformational stability in RIG-1 domains. MethodsIn this study, we employed a systematic in silico strategy on RIG-1 protein to understand the mechanism of structural changes upon mutation. We computationally investigated the protein sequence signature for all the three domains of RIG-1 protein that encloses the mutation within the motif. Further, we carried out a structural comparison between RIG-1 domains with their respective distant orthologs which revealed the minimal number of interactions required to maintain its structural fold. This intra-protein network paved the way to infer hotspot residues crucial for the maintenance of the structural architecture and folding pattern. Key findingsOur analysis revealed about 40 hotspot residues that determine the folding pattern of the RIG-1 domains. Also, conventional molecular dynamic simulation coupled with essential dynamics provides conformational transitions of hot spot residues among native and mutant structures. Structural variations owing to hotspot residues in mutants again confirm the significance of these residues in structural characterization of RIG-1 domains. We believe our results will help the researchers to better comprehend towards regulatory regions and target-binding sites for therapeutic design within the pattern recognition receptor proteins. SignificanceOur protocol employed in this work describes a novel approach in identifying signature residues that would provide structural insights in protein folding.

NHS-Esters As Versatile Reactivity-Based Probes for Mapping Proteome-Wide Ligandable Hotspots.

Most of the proteome is considered undruggable, oftentimes hindering translational efforts for drug discovery. Identifying previously unknown druggable hotspots in proteins would enable strategies for pharmacologically interrogating these sites with small molecules. Activity-based protein profiling (ABPP) has arisen as a powerful chemoproteomic strategy that uses reactivity-based chemical probes to map reactive, functional, and ligandable hotspots in complex proteomes, which has enabled inhibitor discovery against various therapeutic protein targets. Here, we report an alkyne-functionalized N-hydroxysuccinimide-ester (NHS-ester) as a versatile reactivity-based probe for mapping the reactivity of a wide range of nucleophilic ligandable hotspots, including lysines, serines, threonines, and tyrosines, encompassing active sites, allosteric sites, post-translational modification sites, protein interaction sites, and previously uncharacterized potential binding sites. Surprisingly, we also show that fragment-based NHS-ester ligands can be made to confer selectivity for specific lysine hotspots on specific targets including Dpyd, Aldh2, and Gstt1. We thus put forth NHS-esters as promising reactivity-based probes and chemical scaffolds for covalent ligand discovery.

Identification of C3b-Binding Small-Molecule Complement Inhibitors Using Cheminformatics.

The complement system is an elegantly regulated biochemical cascade formed by the collective molecular recognition properties and proteolytic activities of more than two dozen membrane-bound or serum proteins. Complement plays diverse roles in human physiology, such as acting as a sentry against invading microorganisms, priming of the adaptive immune response, and removal of immune complexes. However, dysregulation of complement can serve as a trigger for a wide range of human diseases, which include autoimmune, inflammatory, and degenerative conditions. Despite several potential advantages of modulating complement with small-molecule inhibitors, small-molecule drugs are highly underrepresented in the current complement-directed therapeutics pipeline. In this study, we have employed a cheminformatics drug discovery approach based on the extensive structural and functional knowledge available for the central proteolytic fragment of the cascade, C3b. Using parallel in silico screening methodologies, we identified 45 small molecules that putatively bind C3b near ligand-guided functional hot spots. Surface plasmon resonance experiments resulted in the validation of seven dose-dependent C3b-binding compounds. Competition-based biochemical assays demonstrated the ability of several C3b-binding compounds to interfere with binding of the original C3b ligand that guided their discovery. In vitro assays of complement function identified a single complement inhibitory compound, termed cmp-5, and mechanistic studies of the cmp-5 inhibitory mode revealed it acts at the level of C5 activation. This study has led to the identification of a promising new class of C3b-binding small-molecule complement inhibitors and, to our knowledge, provides the first demonstration of cheminformatics-based, complement-directed drug discovery.

Oxy-haemoglobin protein engineering: An automated design for hotspots stability, site-specific mutations and smart libraries by using HotSpot Wizard 2.0 software

Oxy-haemoglobin is an important metallo-protein, which helps in oxygen binding and transporting to the tissues. The objective of the present study was to detect hot spots and design of smart libraries for engineering protein stability, substrate specificity, tunnels and cavities as well as suitable mutability position of oxy-haemoglobin protein by using a software, HotSpot Wizard, version 2.0., is a free online software. The prediction results were obtained in output interface for functional hot spots, stability hot spots (structural flexibility), correlated hot spots and stability hot spots (sequence consensus). In conclusion, pocket identification and mutability prediction of oxy-haemoglobin can lead to detect structural alternation mainly in disease diagnosis and space for ligand binding pocket in new drug development for disease therapy. This computational prediction is suggesting to compare with experimental hotspots for oxy-haemoglobin in relation to therapeutic efficacies and druggability assessment . Keywords: oxyhaemoglobin; metallo-protein; protein engineering; HotSpot Wizard; computational biology; druggability

Computational prediction and analysis of deleterious cancer associated missense mutations in DYNC1H1

Dynein, cytoplasmic 1, heavy chain 1 (DYNC1H1) gene encodes a subunit of human cytoplasmic dynein complex, which has several crucial functions in the cell, such as intracellular transport of DNA damage proteins and mitotic spindle positioning. Recent studies reported the altered expression of DYNC1H1 in different cancers and DYNC1H1 was suggested to be potential biomarker in colorectal cancers. Previously, DYNC1H1 mutations have been associated with neurodegenerative diseases, however mutations of DYNC1H1 have not been fully investigated in cancers except for different types of pancreatic cancers. In this study we aimed to identify the cancer related mutations in DYNC1H1, which are deleterious for the DYNC1H1 structure and/or function. We investigated 523 cancer related missense mutations in DYNC1H1, which were collected from COSMIC database, to predict the effect of mutations on DYNC1H1 function. Of the 523 mutations, we identified 28 amino acid substitutions, which were predicted to be deleterious by PredictSNP1.0. When we searched for the effect of 28 deleterious mutations on protein stability by MUpro and I-Mutant2.0, we observed that most of the mutations decrease the protein stability. We analyzed the localization of deleterious mutations on primary protein structure and identified that predicted deleterious mutations were mainly located in the motor domain, which is crucial for the DYNC1H1 function. In addition, we detected close positioning of mutated residues in AAA + regions on 3D structure by STRUM and UCSF Chimera. When we searched the mutations in COSMIC database, we observed the occurrence of the mutations in different cancers, which might show the importance of these regions in corresponding cancers. Therefore, our findings provide potential structural and functional mutations and hotspots for DYNC1H1.

Changes in the location of biodiversity-ecosystem function hot spots across the seafloor landscape with increasing sediment nutrient loading.

Declining biodiversity and loss of ecosystem function threatens the ability of habitats to contribute ecosystem services. However, the form of the relationship between biodiversity and ecosystem function (BEF) and how relationships change with environmental change is poorly understood. This limits our ability to predict the consequences of biodiversity loss on ecosystem function, particularly in real-world marine ecosystems that are species rich, and where multiple ecosystem functions are represented by multiple indicators. We investigated spatial variation in BEF relationships across a 300 000 m2 intertidal sandflat by nesting experimental manipulations of sediment pore water nitrogen concentration into sites with contrasting macrobenthic community composition. Our results highlight the significance of many different elements of biodiversity associated with environmental characteristics, community structure, functional diversity, ecological traits or particular species (ecosystem engineers) to important functions of coastal marine sediments (benthic oxygen consumption, ammonium pore water concentrations and flux across the sediment-water interface). Using the BEF relationships developed from our experiment, we demonstrate patchiness across a landscape in functional performance and the potential for changes in the location of functional hot and cold spots with increasing nutrient loading that have important implications for mapping and predicating change in functionality and the concomitant delivery of ecosystem services.

Identification and analysis of mutational hotspots in oncogenes and tumour suppressors.

BackgroundThe key to interpreting the contribution of a disease-associated mutation in the development and progression of cancer is an understanding of the consequences of that mutation both on the function of the affected protein and on the pathways in which that protein is involved. Protein domains encapsulate function and position-specific domain based analysis of mutations have been shown to help elucidate their phenotypes.ResultsIn this paper we examine the domain biases in oncogenes and tumour suppressors, and find that their domain compositions substantially differ. Using data from over 30 different cancers from whole-exome sequencing cancer genomic projects we mapped over one million mutations to their respective Pfam domains to identify which domains are enriched in any of three different classes of mutation; missense, indels or truncations. Next, we identified the mutational hotspots within domain families by mapping small mutations to equivalent positions in multiple sequence alignments of protein domains We find that gain of function mutations from oncogenes and loss of function mutations from tumour suppressors are normally found in different domain families and when observed in the same domain families, hotspot mutations are located at different positions within the multiple sequence alignment of the domain.ConclusionsBy considering hotspots in tumour suppressors and oncogenes independently, we find that there are different specific positions within domain families that are particularly suited to accommodate either a loss or a gain of function mutation. The position is also dependent on the class of mutation. We find rare mutations co-located with well-known functional mutation hotspots, in members of homologous domain superfamilies, and we detect novel mutation hotspots in domain families previously unconnected with cancer. The results of this analysis can be accessed through the MOKCa database (http://strubiol.icr.ac.uk/extra/MOKCa).

Unifying view of mechanical and functional hotspots across class A GPCRs.

G protein-coupled receptors (GPCRs) are the largest superfamily of signaling proteins. Their activation process is accompanied by conformational changes that have not yet been fully uncovered. Here, we carry out a novel comparative analysis of internal structural fluctuations across a variety of receptors from class A GPCRs, which currently has the richest structural coverage. We infer the local mechanical couplings underpinning the receptors’ functional dynamics and finally identify those amino acids whose virtual deletion causes a significant softening of the mechanical network. The relevance of these amino acids is demonstrated by their overlap with those known to be crucial for GPCR function, based on static structural criteria. The differences with the latter set allow us to identify those sites whose functional role is more clearly detected by considering dynamical and mechanical properties. Of these sites with a genuine mechanical/dynamical character, the top ranking is amino acid 7x52, a previously unexplored, and experimentally verifiable key site for GPCR conformational response to ligand binding.

Spatial congruence between taxonomic, phylogenetic and functional hotspots: true pattern or methodological artefact?

AbstractAimTo assess the spatial congruence between hotspots based on taxonomic, phylogenetic and functional diversity, after accounting for the correlation between diversity metrics, and the spatial scale and sampling completeness of data.LocationThe Ordesa and Monte Perdido National Park (Central Pyrenees, Spain), a species‐rich area subjected to intensive botanical sampling.MethodsWe selected hotspots using different diversity metrics and two different data sources (~49,000 occurrence records of 1379 vascular plants in 1 × 1 km grid cells and 1218 inventories of plant communities containing a total of 859 taxa) and compared their spatial congruence. The effect of sampling completeness of data was explicitly assessed. Phylogenetic diversity and functional diversity (measured with richness‐dependent and richness‐independent metrics) were based on a molecular phylogeny, and a functional dendrogram, respectively. The effectiveness of different types of hotspots in representing other diversity components was tested with permutation tests.ResultsWe found that spurious correlations between diversity metrics explained the congruence between taxonomic, phylogenetic and functional hotspots. When richness‐independent metrics were used, diversity hotspots were no longer congruent regardless of the source of data. Hotspots were biased towards intensively sampled grid cells, and the amount of diversity they captured was exaggerated due to the coarse spatial scale of species‐occurrence data. The efficiency of hotspots in terms of integrating different diversity components was lower at community scale and not significantly higher than expected at random, regardless of the sampling completeness.Main conclusionsOur results stress that the arbitrary use of diversity metrics and the scale of analyses along with the sampling bias in data can distort the true location of hotspots, and exaggerate their spatial congruence. After accounting for such methodological issues, we found a clear mismatch between diversity components that questions the utility of hotspots as a conservation tool of multiple diversity components.

Microsecond molecular dynamics simulations provide insight into the ATP-competitive inhibitor-induced allosteric protection of Akt kinase phosphorylation.

Akt is a serine/threonine protein kinase, a critical mediator of growth factor-induced survival in key cellular pathways. Allosteric signaling between protein intramolecular domains requires long-range communication mediated by hotspot residues, often triggered by ligand binding. Here, based on extensive 3μs explicit solvent molecular dynamics (MD) simulations of Akt1 kinase domain in the unbound (apo) and ATP-competitive inhibitor, GDC-0068-bound states, we propose a molecular mechanism for allosteric regulation of Akt1 kinase phosphorylation by GDC-0068 binding to the ATP-binding site. MD simulations revealed that the apo Akt1 is flexible with two disengaged N- and C-lobes, equilibrated between the open and closed conformations. GDC-0068 occupancy of the ATP-binding site shifts the conformational equilibrium of Akt1 from the open conformation toward the closed conformation and stabilizes the closed state. This effect enables allosteric signal propagation from the GDC-0068 to the phosphorylated T308 (pT308) in the activation loop and restrains phosphatase access to pT308, thereby protecting the pT308 in the GDC-0068-bound Akt1. Importantly, functional hotspots involved in the allosteric communication from the GDC-0068 to the pT308 are identified. Our analysis of GDC-0068-induced allosteric protection of Akt kinase phosphorylation yields important new insights into the molecular mechanism of allosteric regulation of Akt kinase activity.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus).

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake’s body length and nestle within the rib cavity. Despite the snake’s notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

Discovery of 3α,7α,11β-Trihydroxy-6α-ethyl-5β-cholan-24-oic Acid (TC-100), a Novel Bile Acid as Potent and Highly Selective FXR Agonist for Enterohepatic Disorders.

As a continuation of previous efforts in mapping functional hot spots on the bile acid scaffold, we here demonstrate that the introduction of a hydroxy group at the C11β position affords high selectivity for FXR. In particular, the synthesis and FXR/TGR5 activity of novel bile acids bearing different hydroxylation patterns at the C ring are reported and discussed from a structure-activity standpoint. The results obtained led us to discover the first bile acid derivative endowed with high potency and selectivity at the FXR receptor, 3α,7α,11β-trihydroxy-6α-ethyl-5β-cholan-24-oic acid (TC-100, 7) which also shows a remarkable physicochemical and pharmacological profile. Compound 7 combines the excellent physicochemical properties of hydrophilic bile acids such as ursodeoxycholic acid, with the distinct ability to specifically bind and regulate FXR activity in vivo, thus providing a bona fide novel therapeutic agent to treat enterohepatic disorders such as cholestasis, NASH, and inflammatory bowel disease.

The Generation and Exploitation of Protein Mutability Landscapes for Enzyme Engineering.

The increasing number of enzyme applications in chemical synthesis calls for new engineering methods to develop the biocatalysts of the future. An interesting concept in enzyme engineering is the generation of large‐scale mutational data in order to chart protein mutability landscapes. These landscapes allow the important discrimination between beneficial mutations and those that are neutral or detrimental, thus providing detailed insight into sequence–function relationships. As such, mutability landscapes are a powerful tool with which to identify functional hotspots at any place in the amino acid sequence of an enzyme. These hotspots can be used as targets for combinatorial mutagenesis to yield superior enzymes with improved catalytic properties, stability, or even new enzymatic activities. The generation of mutability landscapes for multiple properties of one enzyme provides the exciting opportunity to select mutations that are beneficial either for one or for several of these properties. This review presents an overview of the recent advances in the construction of mutability landscapes and discusses their importance for enzyme engineering.

Discovering hotspots in functional genomic data superposed on 3D chromatin configuration reconstructions.

The spatial organization of the genome influences cellular function, notably gene regulation. Recent studies have assessed the three-dimensional (3D) co-localization of functional annotations (e.g. centromeres, long terminal repeats) using 3D genome reconstructions from Hi-C (genome-wide chromosome conformation capture) data; however, corresponding assessments for continuous functional genomic data (e.g. chromatin immunoprecipitation-sequencing (ChIP-seq) peak height) are lacking. Here, we demonstrate that applying bump hunting via the patient rule induction method (PRIM) to ChIP-seq data superposed on a Saccharomyces cerevisiae 3D genome reconstruction can discover ‘functional 3D hotspots’, regions in 3-space for which the mean ChIP-seq peak height is significantly elevated. For the transcription factor Swi6, the top hotspot by P-value contains MSB2 and ERG11 – known Swi6 target genes on different chromosomes. We verify this finding in a number of ways. First, this top hotspot is relatively stable under PRIM across parameter settings. Second, this hotspot is among the top hotspots by mean outcome identified by an alternative algorithm, k-Nearest Neighbor (k-NN) regression. Third, the distance between MSB2 and ERG11 is smaller than expected (by resampling) in two other 3D reconstructions generated via different normalization and reconstruction algorithms. This analytic approach can discover functional 3D hotspots and potentially reveal novel regulatory interactions.

Interactions of noradrenaline and cortisol and the induction of indelible memories.

The glutamate amplifies noradrenergic effects (GANE) model emphasizes the role of focal glutamate-noradrenaline interactions in creating functional hotspots for prioritized processing of salient stimuli. Here, we briefly outline current evidence that synergistic action of noradrenaline and cortisol enables emotional stimuli to gain privileged access to amygdala-hippocampus circuits, eventually resulting in the formation of indelible memories and posttraumatic stress disorder (PTSD).

"Phylogenetic and evolutionary analysis of functional divergence among Gamma glutamyl transpeptidase (GGT) subfamilies".

Backgroundγ-glutamyltranspeptidase (GGT) is a bi-substrate enzyme conserved in all three domains of life. It catalyzes the cleavage and transfer of γ-glutamyl moiety of glutathione to either water (hydrolysis) or substrates like peptides (transpeptidation). GGTs exhibit great variability in their enzyme kinetics although the mechanism of catalysis is conserved. Recently, GGT has been shown to be a virulence factor in microbes like Helicobacter pylori and Bacillus anthracis. In mammalian cells also, GGT inhibition prior to chemotherapy has been shown to sensitize tumors to the therapy. Therefore, lately both bacterial and eukaryotic GGTs have emerged as potential drug targets, but the efforts directed towards finding suitable inhibitors have not yielded any significant results yet. We propose that delineating the residues responsible for the functional diversity associated with these proteins could help in design of species/clade specific inhibitors.ResultsIn the present study, we have carried out phylogenetic analysis on a set of 47 GGT-like proteins to address the functional diversity. These proteins segregate into various subfamilies, forming separate clades on the tree. Sequence conservation and motif prediction studies show that even though most of the highly conserved residues have been characterized biochemically in previous studies, a significant number of novel putative sites and motifs are discovered that vary in a clade specific manner. Many of the putative sites predicted during the functional divergence type I and type II analysis, lie close to the known catalytic residues and line the walls of the substrate binding cavity, reinforcing their role in modulating the substrate specificity, catalytic rates and stability of this protein.ConclusionThe study offers interesting insights into the evolution of GGT-like proteins in pathogenic vs. non-pathogenic bacteria, archaea and eukaryotes. Our analysis delineates residues that are highly specific to each GGT subfamily. We propose that these sites not only explain the differences in stability and catalytic variability of various GGTs but can also aid in design of specific inhibitors against particular GGTs. Thus, apart from the commonly used in-silico inhibitor screening approaches, evolutionary analysis identifying the functional divergence hotspots in GGT proteins could augment the structure based drug design approaches.ReviewersThis article was reviewed by Andrei Osterman, Christine Orengo, and Srikrishna Subramanian. For complete reports, see the Reviewers’ reports sectionElectronic supplementary materialThe online version of this article (doi:10.1186/s13062-015-0080-7) contains supplementary material, which is available to authorized users.

Pan-Cancer Analysis of Mutation Hotspots in Protein Domains

In cancer genomics, recurrence of mutations in independent tumor samples is a strong indicator of functional impact. However, rare functional mutations can escape detection by recurrence analysis owing to lack of statistical power. We enhance statistical power by extending the notion of recurrence of mutations from single genes to gene families that share homologous protein domains. Domain mutation analysis also sharpens the functional interpretation of the impact of mutations, as domains more succinctly embody function than entire genes. By mapping mutations in 22 different tumor types to equivalent positions in multiple sequence alignments of domains, we confirm well-known functional mutation hotspots, identify uncharacterized rare variants in one gene that are equivalent to well-characterized mutations in another gene, detect previously unknown mutation hotspots, and provide hypotheses about molecular mechanisms and downstream effects of domain mutations. With the rapid expansion of cancer genomics projects, protein domain hotspot analysis will likely provide many more leads linking mutations in proteins to the cancer phenotype.

Subgroup-specific intrinsic disorder profiles of arabidopsis NAC transcription factors: Identification of functional hotspots

Protein intrinsic disorder (ID), referring to the lack of a fixed tertiary structure, is significant in signaling and transcription. We recently characterized ID in 6 phylogenetically representative Arabidopsis thaliana NAC transcription factors. Their transcription regulatory domains are mostly disordered but contain short, functionally important regions with structure propensities known as molecular recognition features. Here, we analyze for NAC subgroup-specific ID patterns. Some subgroups, such as the VND subgroup implicated in secondary cell wall biosynthesis, and the NAP/SHYG subgroup have highly conserved ID profiles. For the stress-associated ATAF1 subgroup and the CUC/ORE1 subgroup involved in development, only sub clades have similar ID patterns. For similar ID profiles, conserved molecular recognition features and sequence motifs represent likely functional determinants of e.g. transcriptional activation and interactions. Based on our analysis, we suggest that ID profiling of regulatory proteins in general can be used to guide identification of interaction partners of network proteins.

Allosteric Regulation of the Hsp90 Dynamics and Stability by Client Recruiter Cochaperones: Protein Structure Network Modeling

The fundamental role of the Hsp90 chaperone in supporting functional activity of diverse protein clients is anchored by specific cochaperones. A family of immune sensing client proteins is delivered to the Hsp90 system with the aid of cochaperones Sgt1 and Rar1 that act cooperatively with Hsp90 to form allosterically regulated dynamic complexes. In this work, functional dynamics and protein structure network modeling are combined to dissect molecular mechanisms of Hsp90 regulation by the client recruiter cochaperones. Dynamic signatures of the Hsp90-cochaperone complexes are manifested in differential modulation of the conformational mobility in the Hsp90 lid motif. Consistent with the experiments, we have determined that targeted reorganization of the lid dynamics is a unifying characteristic of the client recruiter cochaperones. Protein network analysis of the essential conformational space of the Hsp90-cochaperone motions has identified structurally stable interaction communities, interfacial hubs and key mediating residues of allosteric communication pathways that act concertedly with the shifts in conformational equilibrium. The results have shown that client recruiter cochaperones can orchestrate global changes in the dynamics and stability of the interaction networks that could enhance the ATPase activity and assist in the client recruitment. The network analysis has recapitulated a broad range of structural and mutagenesis experiments, particularly clarifying the elusive role of Rar1 as a regulator of the Hsp90 interactions and a stability enhancer of the Hsp90-cochaperone complexes. Small-world organization of the interaction networks in the Hsp90 regulatory complexes gives rise to a strong correspondence between highly connected local interfacial hubs, global mediator residues of allosteric interactions and key functional hot spots of the Hsp90 activity. We have found that cochaperone-induced conformational changes in Hsp90 may be determined by specific interaction networks that can inhibit or promote progression of the ATPase cycle and thus control the recruitment of client proteins.