Identification of potential allosteric binding sites in cathepsin K based on intramolecular communication

Network theory methods and molecular dynamics (MD) simulations are accepted tools to study allosteric regulation. Indeed, dynamic networks built upon correlation analysis of MD trajectories provide detailed information about communication paths between distant sites. In this context, we aimed to understand whether the efficiency of intramolecular communication could be used to predict the allosteric potential of a given site. To this end, we performed MD simulations and network theory analyses in cathepsin K (catK), whose allosteric sites are well defined. To obtain a quantitative measure of the efficiency of communication, we designed a new protocol that enables the comparison between properties related to ensembles of communication paths obtained from different sites. Further, we applied our strategy to evaluate the allosteric potential of different catK cavities not yet considered for drug design. Our predictions of the allosteric potential based on intramolecular communication correlate well with previous catK experimental and theoretical data. We also discuss the possibility of applying our approach to other proteins from the same family.     

Nested allosteric interactions in the cytoplasmic chaperonin containing TCP-1

Initial rates of ATP hydrolysis by the chaperonin containing TCP-1 (CCT) from bovine testis were measured as a function of ATP concentration. Two allosteric transitions are observed: one at relatively low concentrations of ATP (<100 microM) and the second at higher concentrations of ATP. The data suggest that CCT has positive intra-ring cooperativity and negative inter-ring cooperativity in ATP hydrolysis, with respect to ATP, as previously observed in the case of GroEL. It is shown that the relatively weak positive intra-ring cooperativity found in the case of CCT may be due to heterogeneity in its subunit composition. Our results suggest that nested allosteric behavior may be common to chaperone double-ring systems.     

Allosteric changes in protein structure computed by a simple mechanical model: hemoglobin T<-->R2 transition

Information on protein dynamics has been usually inferred from spectroscopic studies of parts of the proteins, or indirectly from the comparison of the conformations assumed in the presence of different substrates or ligands. While molecular simulations also provide information on protein dynamics, they usually suffer from incomplete sampling of conformational space, and become prohibitively expensive when exploring the collective dynamics of large macromolecular structures. Here, we explore the dynamics of a well-studied allosteric protein, hemoglobin (Hb), to show that a simple mechanical model based on Gaussian fluctuations of residues can efficiently predict the transition between the tense (T, unliganded) and relaxed (R or R2, O(2) or CO-bound) forms of Hb. The passage from T into R2 is shown to be favored by the global mode of motion, which, in turn is driven by entropic effects. The major difference between the dynamics of the T and R2 forms is the loss of the hinge-bending role of alpha(1)-beta(2) (or alpha(2)-beta(1)) interfacial residues at alpha Phe36-His45 and beta Thr87-Asn102 in the R2 form, which implies a decreased cooperativity in the higher affinity (R2) form of Hb, consistent with many experimental studies. The involvement of the proximal histidine beta His92 in this hinge region suggests that the allosteric propagation of the local structural changes (induced upon O(2) binding) into global ones occur via hinge regions. This is the first demonstration that there is an intrinsic tendency of Hb to undergo T-->R2 transition, induced by purely elastic forces of entropic origin that are uniquely defined for the particular contact topology of the T form.     

Sparse networks of directly coupled,polymorphic, and functional side chains in allosteric proteins

Recent studies have highlighted the role of coupled side‐chain fluctuations alone in the allosteric behavior of proteins. Moreover, examination of X‐ray crystallography data has recently revealed new information about the prevalence of alternate side‐chain conformations (conformational polymorphism), and attempts have been made to uncover the hidden alternate conformations from X‐ray data. Hence, new computational approaches are required that consider the polymorphic nature of the side chains, and incorporate the effects of this phenomenon in the study of information transmission and functional interactions of residues in a molecule. These studies can provide a more accurate understanding of the allosteric behavior. In this article, we first present a novel approach to generate an ensemble of conformations and an efficient computational method to extract direct couplings of side chains in allosteric proteins, and provide sparse network representations of the couplings. We take the side‐chain conformational polymorphism into account, and show that by studying the intrinsic dynamics of an inactive structure, we are able to construct a network of functionally crucial residues. Second, we show that the proposed method is capable of providing a magnified view of the coupled and conformationally polymorphic residues. This model reveals couplings between the alternate conformations of a coupled residue pair. To the best of our knowledge, this is the first computational method for extracting networks of side chains' alternate conformations. Such networks help in providing a detailed image of side‐chain dynamics in functionally important and conformationally polymorphic sites, such as binding and/or allosteric sites. Proteins 2015; 83:497–516. © 2014 Wiley Periodicals, Inc.     

Identification of core allosteric sites through temperature- and nucleus-invariant chemical shift covariance

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or in the design of allosteric drugs. Here, we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. T- and CLASS-CHESCAs, as well as complete-linkage CHESCA, were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, T-, CLASS-, and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.     

Perturbing dimer interactions and allosteric communication modulates the immunosuppressive activity of human galectin-7

The design of allosteric modulators to control protein function is a key objective in drug discovery programs. Altering functionally essential allosteric residue networks provides unique protein family subtype specificity, minimizes unwanted off-target effects, and helps avert resistance acquisition typically plaguing drugs that target orthosteric sites. In this work, we used protein engineering and dimer interface mutations to positively and negatively modulate the immunosuppressive activity of the proapoptotic human galectin-7 (GAL-7). Using the PoPMuSiC and BeAtMuSiC algorithms, mutational sites and residue identity were computationally probed and predicted to either alter or stabilize the GAL-7 dimer interface. By designing a covalent disulfide bridge between protomers to control homodimer strength and stability, we demonstrate the importance of dimer interface perturbations on the allosteric network bridging the two opposite glycan-binding sites on GAL-7, resulting in control of induced apoptosis in Jurkat T cells. Molecular investigation of G16X GAL-7 variants using X-ray crystallography, biophysical, and computational characterization illuminates residues involved in dimer stability and allosteric communication, along with discrete long-range dynamic behaviors involving loops 1, 3, and 5. We show that perturbing the protein–protein interface between GAL-7 protomers can modulate its biological function, even when the overall structure and ligand-binding affinity remains unaltered. This study highlights new avenues for the design of galectin-specific modulators influencing both glycan-dependent and glycan-independent interactions.     

Detection of multiple network-based allosteric interactions between peptides and arrays of DNA binding sites

Quantitative DNase I footprinting shows that three designed peptides containing N-methylpyrrole (Py) moieties display different types of network-based allosteric communication in binding to DNA: circuit type, incomplete-circuit type, and non-circuit type characterized by interstrand bidentate interactions. Positive cooperative binding of all three peptides to individual DNA binding sites is commonly observed. CD spectral characterization of the interaction between peptides and model undecanucleotide duplexes is consistent with the footprinting results and supports the allosteric model. This study provides insights relating to the interaction network nature of allostery in complex DNA–small molecule interactions.     

Human melanocortin 1 receptor-mediated ubiquitination of nonvisual arrestins. Role of Mahogunin Ring Finger 1 E3 ligase

Signaling from the melanocortin 1 receptor (MC1R), a Gs protein-coupled receptor (GPCR) crucial for melanocyte proliferation and differentiation, is regulated by cytosolic β-arrestins (ARRBs). MC1R signaling is also negatively modulated by the E3-ubiquitin ligase Mahogunin Ring Finger-1 (MGRN1), whose mutation causes hyperpigmentation, congenital heart defects and neurodegeneration in mice. We showed previously that although MC1R interacts stably with human ARRB1 or ARRB2, only ARRB2 mediates receptor desensitization and internalization. We analyzed MC1R-dependent ARRB ubiquitination, and the possible role of MGRN1. ARRB1 expressed in heterologous cells or human melanoma cells migrated in SDS-PAGE as a 55 kDa protein whereas ARRB2 migrated as two major bands of apparent molecular weight near 45 and 55 kDa, with an intermediate mobility band occasionally detected. These forms were related by post-translational modification rather than by proteolysis. Presence of MC1R favored expression of the 45 kDa protein, the form that interacted preferentially with MC1R. MC1R also mediated poly- or multimonoubiquitination of ARRB2. Ubiquitination was agonist-independent, but required a native MC1R conformation and/or normal receptor trafficking to the plasma membrane, as it was not observed for loss-of-function MC1R variants. In a heterologous expression system, MC1R-dependent ARRB ubiquitination was enhanced by overexpression of MGRN1 and was impaired by siRNA-mediated MGRN1 knockdown thus pointing to MGRN1 as the responsible E3-ligase. Co-immunoprecipitation experiments demonstrated interaction of MGRN1 and ARRBs in the presence of MC1R, suggesting a scaffolding role for the GPCR that may determine the selectivity of E3-ubiquitin ligase engagement and the functional outcome of ARRB ubiquitination.     

Mapping the intramolecular signal transduction of G‐protein coupled receptors

G‐protein coupled receptors (GPCRs), a major gatekeeper of extracellular signals on plasma membrane, are unarguably one of the most important therapeutic targets. Given the recent discoveries of allosteric modulations, an allosteric wiring diagram of intramolecular signal transductions would be of great use to glean the mechanism of receptor regulation. Here, by evaluating betweenness centrality (C_B) of each residue, we calculate maps of information flow in GPCRs and identify key residues for signal transductions and their pathways. Compared with preexisting approaches, the allosteric hotspots that our C_B‐based analysis detects for A_2A adenosine receptor (A_2AAR) and bovine rhodopsin are better correlated with biochemical data. In particular, our analysis outperforms other methods in locating the rotameric microswitches, which are generally deemed critical for mediating orthosteric signaling in class A GPCRs. For A_2AAR, the inter‐residue cross‐correlation map, calculated using equilibrium structural ensemble from molecular dynamics simulations, reveals that strong signals of long‐range transmembrane communications exist only in the agonist‐bound state. A seemingly subtle variation in structure, found in different GPCR subtypes or imparted by agonist bindings or a point mutation at an allosteric site, can lead to a drastic difference in the map of signaling pathways and protein activity. The signaling map of GPCRs provides valuable insights into allosteric modulations as well as reliable identifications of orthosteric signaling pathways. Proteins 2014; 82:727–743. © 2013 Wiley Periodicals, Inc.     

Coupling dynamics and evolutionary information with structure to identify protein regulatory and functional binding sites

Binding sites in proteins can be either specifically functional binding sites (active sites) that bind specific substrates with high affinity or regulatory binding sites (allosteric sites), that modulate the activity of functional binding sites through effector molecules. Owing to their significance in determining protein function, the identification of protein functional and regulatory binding sites is widely acknowledged as an important biological problem. In this work, we present a novel binding site prediction method, Active and Regulatory site Prediction (AR-Pred), which supplements protein geometry, evolutionary, and physicochemical features with information about protein dynamics to predict putative active and allosteric site residues. As the intrinsic dynamics of globular proteins plays an essential role in controlling binding events, we find it to be an important feature for the identification of protein binding sites. We train and validate our predictive models on multiple balanced training and validation sets with random forest machine learning and obtain an ensemble of discrete models for each prediction type. Our models for active site prediction yield a median area under the curve (AUC) of 91% and Matthews correlation coefficient (MCC) of 0.68, whereas the less well-defined allosteric sites are predicted at a lower level with a median AUC of 80% and MCC of 0.48. When tested on an independent set of proteins, our models for active site prediction show comparable performance to two existing methods and gains compared to two others, while the allosteric site models show gains when tested against three existing prediction methods. AR-Pred is available as a free downloadable package at https://github.com/sambitmishra0628/AR-PRED_source .     

A novel homology model of TRPC3 reveals allosteric coupling between gate and selectivity filter

Utilizing a novel molecular model of TRPC3, based on the voltage-gated sodium channel from Arcobacter butzleri (Na_VAB) as template, we performed structure-guided mutagenesis experiments to identify amino acid residues involved in divalent permeation and gating. Substituted cysteine accessibility screening within the predicted selectivity filter uncovered amino acids 629–631 as the narrowest part of the permeation pathway with an estimated pore diameter of <5.8 Å. E630 was found to govern not only divalent permeability but also sensitivity of the channel to block by ruthenium red. Mutations in a hydrophobic cluster at the cytosolic termini of transmembrane segment 6, corresponding to the S6 bundle crossing structure in Na_VAB, distorted channel gating. Removal of a large hydrophobic residue (I667A or I667E) generated channels with approximately 60% constitutive activity, suggesting I667 as part of the dynamic structure occluding the permeation path. Destabilization of the gate was associated with reduced Ca²⁺ permeability, altered cysteine cross-linking in the selectivity filter and promoted channel block by ruthenium red. Collectively, we present a structural model of the TRPC3 permeation pathway and localize the channel's selectivity filter and the occluding gate. Moreover, we provide evidence for allosteric coupling between the gate and the selectivity filter in TRPC3.     

Fast kinetic studies on the allosteric interactions between acetylcholine receptor and local anesthetic binding sites.

Preincubation of receptor-rich membrane fragments from Torpedo marmorata with tertiary amine local anesthetics and several toxins such as histrionicotoxin, crotoxin and cerulotoxin, modifies the amplitude and time course of the relaxation processes monitored upon rapid mixing of the membrane fragments with the fluorescent agonist, Dns-C6-Cho. In particular, the amplitude of the rapid relaxation process, which is proportional to the fraction of acetylcholine receptor sites in a high-affinity state, increases; accordingly, the rate constant of the 'slow' and 'intermediate' relaxation processes also increases up to ten times (except with histrionicotoxin) whereas in a higher range of local anesthetic concentrations the rate constant of the 'rapid' relaxation process decreases. The data are accounted for by a two-state model of the acetylcholine regulator, assuming distinct binding sites for cholinergic agonists and local anesthetics and allosteric interactions between these two classes of sites; local anesthetics stabilize the regulator in a high-affinity state for agonists even in the absence of agonist, and modify the rate constants for th interconversions between the low-affinity and high-affinity states. The model accounts for the 'slow' fluorescence increase monitored upon addition of local anesthetics to a suspension of receptor-rich membranes supplemented with trace amounts of Dns-C6-Cho. The effect of local anesthetics on the apparent rate constant of the 'rapid' relaxation process can be accounted for on the basis of an additional low-affinity binding of local anesthetics to the acetylcholine receptor site. Finally the increase of the apparent rate constant of the 'intermediate' relaxation process can be simply accounted for by assuming the existence of a third state, corresponding to the 'active' state, to which local anesthetics bind and block ionic transport.     

A sphingosine 1-phosphate receptor 2 selective allosteric agonist

Molecular probe tool compounds for the Sphingosine 1-phosphate receptor 2 (S1PR2) are important for investigating the multiple biological processes in which the S1PR2 receptor has been implicated. Amongst these are NF-κB-mediated tumor cell survival and fibroblast chemotaxis to fibronectin. Here we report our efforts to identify selective chemical probes for S1PR2 and their characterization. We employed high throughput screening to identify two compounds which activate the S1PR2 receptor. SAR optimization led to compounds with high nanomolar potency. These compounds, XAX-162 and CYM-5520, are highly selective and do not activate other S1P receptors. Binding of CYM-5520 is not competitive with the antagonist JTE-013. Mutation of receptor residues responsible for binding to the zwitterionic headgroup of sphingosine 1-phosphate (S1P) abolishes S1P activation of the receptor, but not activation by CYM-5520. Competitive binding experiments with radiolabeled S1P demonstrate that CYM-5520 is an allosteric agonist and does not displace the native ligand. Computational modeling suggests that CYM-5520 binds lower in the orthosteric binding pocket, and that co-binding with S1P is energetically well tolerated. In summary, we have identified an allosteric S1PR2 selective agonist compound.     

Highly multiplexed bioactivity screening reveals human and microbiota metabolome-GPCRome interactions

1. Download : Download high-res image (221KB)
2. Download : Download full-size image

     

V1b vasopressin receptor trafficking and signaling: Role of arrestins,G proteins and Src kinase

The signaling pathway of G protein‐coupled receptors is strongly linked to their trafficking profile. Little is known about the molecular mechanisms involved in the vasopressin receptor V_1b subtype (V_1bR) trafficking and its impact on receptor signaling and regulation. For this purpose, we investigated the role of β‐arrestins in receptor desensitization, internalization and recycling and attempted to dissect the V_1bR‐mediated MAP kinase pathway. Using MEF cells Knocked‐out for β‐arrestins 1 and 2, we demonstrated that both β‐arrestins 1 and 2 play a fundamental role in internalization and recycling of V_1bR with a rapid and transient V_1bR‐β‐arrestin interaction in contrast to a slow and long‐lasting β‐arrestin recruitment of the V₂ vasopressin receptor subtype (V₂R). Using V_1bR‐V₂R chimeras and V_1bR C‐terminus truncations, we demonstrated the critical role of the V_1bR C‐terminus in its interaction with β‐arrestins thereby regulating the receptor internalization and recycling kinetics in a phosphorylation‐independent manner. In parallel, V_1bR MAP kinase activation was dependent on arrestins and Src‐kinase but independent on G proteins. Interestingly, Src interacted with hV_1bR at basal state and dissociated when receptor internalization occurred. Altogether, our data describe for the first time the trafficking profile and MAP kinase pathway of V_1bR involving both arrestins and Src kinase family.      

Structural analysis of protein tyrosine phosphatase 1B reveals potentially druggable allosteric binding sites

Catalytic proteins such as human protein tyrosine phosphatase 1B (PTP1B), with conserved and highly polar active sites, warrant the discovery of druggable nonactive sites, such as allosteric sites, and potentially, therapeutic small molecules that can bind to these sites. Catalyzing the dephosphorylation of numerous substrates, PTP1B is physiologically important in intracellular signal transduction pathways in diverse cell types and tissues. Aberrant PTP1B is associated with obesity, diabetes, cancers, and neurodegenerative disorders. Utilizing clustering methods (based on root mean square deviation, principal component analysis, nonnegative matrix factorization, and independent component analysis), we have examined multiple PTP1B structures. Using the resulting representative structures in different conformational states, we determined consensus clustroids and used them to identify both known and novel binding sites, some of which are potentially allosteric. We report several lead compounds that could potentially bind to the novel PTP1B binding sites and can be further optimized. Considering the possibility for drug repurposing, we discovered homologous binding sites in other proteins, with ligands that could potentially bind to the novel PTP1B binding sites.     

In silico mutational studies of Hsp70 disclose sites with distinct functional attributes

The Mutation‐Minimization Method (MuMi) to study the local response of proteins to point mutations has been introduced here. The heat shock protein Hsp70 as the test system since it displays features that have been studied in great detail has been used here. It has many conserved residues, serves several different functions on each of its domains, and displays interdomain allostery. For the analysis of spatial arrangement of residues within the protein, the network properties of the wild type (WT) protein as well as its all single alanine residue mutants using MuMi has been investigated. The measures to express the amount of change from the WT structure upon mutation and compare these deviations to find potential critical sites have been proposed. The functional significance of the potential sites to the parameter that uncovers them has been mapped. It was found that sites directly involved in binding were sensitive to mutations and were characterized by large displacements. On the other hand, sites that steer large conformational changes typically had increased reachability upon alanine mutations occurring elsewhere in the protein. Finally, residues that control communication within and between domains reside on the largest number of paths connecting pairs of residues in the protein. Proteins 2015; 83:2077–2090. © 2015 Wiley Periodicals, Inc.     

Conformational changes and allosteric communications in human serum albumin due to ligand binding

It is well recognized that knowledge of structure alone is not sufficient to understand the fundamental mechanism of biomolecular recognition. Information of dynamics is necessary to describe motions involving relevant conformational states of functional importance. We carried out principal component analysis (PCA) of structural ensemble, derived from 84 crystal structures of human serum albumin (HSA) with different ligands and/or different conditions, to identify the functionally important collective motions, and compared with the motions along the low-frequency modes obtained from normal mode analysis of the elastic network model (ENM) of unliganded HSA. Significant overlap is observed in the collective motions derived from PCA and ENM. PCA and ENM analysis revealed that ligand selects the most favored conformation from accessible equilibrium structures of unliganded HSA. Further, we analyzed dynamic network obtained from molecular dynamics simulations of unliganded HSA and fatty acids- bound HSA. Our results show that fatty acids-bound HSA has more robust community network with several routes to communicate among different parts of the protein. Critical nodes (residues) identified from dynamic network analysis are in good agreement with allosteric residues obtained from sequence-based statistical coupling analysis method. This work underscores the importance of intrinsic structural dynamics of proteins in ligand recognition and can be utilized for the development of novel drugs with optimum activity.     

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.