Binding leverage as a molecular basis for allosteric regulation

Predicting allosteric pockets in protein biological assemblages

MOTIVATION

Allostery enables changes to the dynamic behavior of a protein at distant positions induced by binding. Here, we present APOP, a new allosteric pocket prediction method, which perturbs the pockets formed in the structure by stiffening pairwise interactions in the elastic network across the pocket, to emulate ligand binding. Ranking the pockets based on the shifts in the global mode frequencies, as well as their mean local hydrophobicities, leads to high prediction success when tested on a dataset of allosteric proteins, composed of both monomers and multimeric assemblages.

RESULTS

Out of the 104 test cases, APOP predicts known allosteric pockets for 92 within the top 3 rank out of multiple pockets available in the protein. In addition, we demonstrate that APOP can also find new alternative allosteric pockets in proteins. Particularly interesting findings are the discovery of previously overlooked large pockets located in the centers of many protein biological assemblages; binding of ligands at these sites would likely be particularly effective in changing the protein's global dynamics.

AVAILABILITY AND IMPLEMENTATION

APOP is freely available as an open-source code (https://github.com/Ambuj-UF/APOP) and as a web server at https://apop.bb.iastate.edu/.

Learning about allosteric drugs and ways to design them

While the accelerating quest for precision medicine requires new individually targeting and selective drugs, and the ability to work with so-called undruggable targets, the realm of allosteric drugs meeting this need remains largely uncharted. Generalizing the observations on two major drug targets with widely observed inherent allostery, GPCRs and kinases, we describe and discuss basic allosteric modes of action that are universally applicable in all types of structures and functions. Using examples of Class A GPCRs and CMGC protein kinases, we show how Allosteric Signalling and Probing Fingerprints can be used to identify potential allosteric sites and reveal effector-leads that may serve as a starting point for the development of allosteric drugs targeting these regulatory sites. A set of distinct characteristics of allosteric ligands was established, which highlights the versatility of their design and make them advantageous before their orthosteric counterparts in personalized medicine. We argue that rational design of allosteric drugs should begin with the search for latent sites or design of non-natural binding sites followed by fragment-based design of allosteric ligands and by the mutual adjustment of the site-ligand pair in order to achieve required effects. On the basis of the perturbative nature and reversibility of allosteric communication, we propose a generic protocol for computational design of allosteric effectors, enabling also the allosteric tuning of biologics, in obtaining allosteric control over protein functions.

Hidden allosteric sites and De-Novo drug design

INTRODUCTION

Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge.

AREAS COVERED

This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator.

EXPERT OPINION

In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.

Drug discovery from natural products using affinity selection-mass spectrometry

As a starting point for drug discovery, affinity selection-mass spectrometry (AS-MS) is ideal for the discovery of lead compounds from chemically diverse sources such as botanical, fungal and microbial extracts. Based on binding interactions between macromolecular receptors and ligands of low molecular mass, AS-MS enables the rapid isolation of pharmacologically active small molecules from complex mixtures for mass spectrometric characterization and identification. Unlike conventional high-throughput screening, AS-MS requires no radiolabels, no UV or fluorescent chromophores, and is compatible with all classes of receptors, enzymes, incubation buffers, cofactors, and ligands. The most successful types of AS-MS include pulsed ultrafiltration (PUF) AS-MS, size exclusion chromatography (SEC) AS-MS, and magnetic microbead affinity selection screening (MagMASS), which differ in their approaches for separating the ligand-receptor complexes from the non-binding compounds in mixtures. After affinity isolation, the ligand(s) from the mixture are characterized using high resolution UHPLC-MS and tandem mass spectrometry. Based on these elemental composition and structural data, the identities of the lead compounds are determined by searching on-line databases for known natural products and by comparison with standards. The structures of novel natural products are determined using a combination of spectroscopic techniques including two-dimensional NMR and MS.

Mixed-solvent molecular dynamics simulation-based discovery of a putative allosteric site on regulator of G protein signaling 4

Regulator of G protein signaling 4 (RGS4) is an intracellular protein that binds to the G_α subunit ofheterotrimeric G proteins and aids in terminating G protein coupled receptor signaling. RGS4 has been implicated in pain, schizophrenia, and the control of cardiac contractility. Inhibitors of RGS4 have been developed but bind covalently to cysteine residues on the protein. Therefore, we sought to identify alternative druggable sites on RGS4 using mixed-solvent molecular dynamics simulations, which employ low concentrations of organic probes to identify druggable hotspots on the protein. Pseudo-ligands were placed in consensus hotspots, and perturbation with normal mode analysis led to the identification and characterization of a putative allosteric site, which would be invaluable for structure-based drug design of non-covalent, small molecule inhibitors. Future studies on the mechanism of this allostery will aid in the development of novel therapeutics targeting RGS4.

Potential allosteric sites captured in glycolytic enzymes via residue-based network models: Phosphofructokinase, glyceraldehyde-3-phosphate dehydrogenase and pyruvate kinase

Likelihood of new allosteric sites for glycolytic enzymes, phosphofructokinase (PFK), glyceraldehyde-3-phosphate dehydrogenase (GADPH) and pyruvate kinase (PK) was evaluated for bacterial, parasitic and human species. Allosteric effect of a ligand binding at a site was revealed on the basis of low-frequency normal modes via C_α-harmonic residue network model. In bacterial PFK, perturbation of the proposed allosteric site outperformed the known allosteric one, producing a high amount of stabilization or reduced dynamics, on all catalytic regions. Another proposed allosteric spot at the dimer interface in parasitic PFK exhibited major stabilization effect on catalytic regions. In parasitic GADPH, the most desired allosteric response was observed upon perturbation of its tunnel region which incorporated key residues for functional regulation. Proposed allosteric site in bacterial PK produced a satisfactory allosteric response on all catalytic regions, whereas in human and parasitic PKs, a partial inhibition was observed. Residue network model based solely on contact topology identified the 'hub residues' with high betweenness tracing plausible allosteric communication pathways between distant functional sites. For both bacterial PFK and PK, proposed sites accommodated hub residues twice as much as the known allosteric site. Tunnel region in parasitic GADPH with the strongest allosteric effect among species, incorporated the highest number of hub residues. These results clearly suggest a one-to-one correspondence between the degree of allosteric effect and the number of hub residues in that perturbation site, which increases the likelihood of its allosteric nature.

Recent advances in user-friendly computational tools to engineer protein function

Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.

Exploring the Allosteric Territory of Protein Function

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Influence of Disease-Causing Mutations on Protein Structural Networks

The interactions between residues in a protein tertiary structure can be studied effectively using the approach of protein structure network (PSN). A PSN is a node-edge representation of the structure with nodes representing residues and interactions between residues represented by edges. In this study, we have employed weighted PSNs to understand the influence of disease-causing mutations on proteins of known 3D structures. We have used manually curated information on disease mutations from UniProtKB/Swiss-Prot and their corresponding protein structures of wildtype and disease variant from the protein data bank. The PSNs of the wildtype and disease-causing mutant are compared to analyse variation of global and local dissimilarity in the overall network and at specific sites. We study how a mutation at a given site can affect the structural network at a distant site which may be involved in the function of the protein. We have discussed specific examples of the disease cases where the protein structure undergoes limited structural divergence in their backbone but have large dissimilarity in their all atom networks and vice versa, wherein large conformational alterations are observed while retaining overall network. We analyse the effect of variation of network parameters that characterize alteration of function or stability.

Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin

Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.

Allosteric Response of DNA Recognition Helices of Catabolite Activator Protein to cAMP and DNA Binding

The homodimeric catabolite activator protein (CAP) regulates the transcription of several bacterial genes based on the cellular concentration of cyclic adenosine monophosphate (cAMP). The binding of cAMP to CAP triggers allosteric communication between the cAMP binding domains (CBD) and DNA binding domains (DBD) of CAP, which entails repositioning of DNA recognition helices (F-helices) in the DBD to dock favorably to the target DNA. Despite considerable progress, much remains to be understood about the mechanistic details of DNA recognition by CAP and about the map of allosteric pathways involved in CAP-mediated gene transcription. The present study uses molecular dynamics and umbrella sampling simulations to investigate the mechanism of cAMP- and DNA-induced changes in the conformation and energetics of F-helices observed during the allosteric regulation of CAP by cAMP and the subsequent binding to the DNA promoter region. Using novel collective variables, the free energy profiles associated with the orientation and dynamics of F-helices in the unliganded, cAMP-bound, and cAMP-DNA-bound states of CAP are calculated and compared. The binding-induced alterations in the resultant free energy profiles reveal important flexibility constraints imposed on DBD upon cAMP and DNA binding. A comprehensive analysis of residue-wise interaction maps reveals potential allosteric pathways between CBD and DBD that facilitate the allosteric transduction of regulatory signals in CAP. The revelation that the predicted allosteric pathways crisscross the intersubunit interface offers important clues on the microscopic origin of the intersubunit cooperativity and dimer stability of CAP.

AlloSigMA 2: paving the way to designing allosteric effectors and to exploring allosteric effects of mutations

The AlloSigMA 2 server provides an interactive platform for exploring the allosteric signaling caused by ligand binding and/or mutations, for analyzing the allosteric effects of mutations and for detecting potential cancer drivers and pathogenic nsSNPs. It can also be used for searching latent allosteric sites and for computationally designing allosteric effectors for these sites with required agonist/antagonist activity. The server is based on the implementation of the Structure-Based Statistical Mechanical Model of Allostery (SBSMMA), which allows one to evaluate the allosteric free energy as a result of the perturbation at per-residue resolution. The Allosteric Signaling Map (ASM) providing a comprehensive residue-by-residue allosteric control over the protein activity can be obtained for any structure of interest. The Allosteric Probing Map (APM), in turn, allows one to perform the fragment-based-like computational design experiment aimed at finding leads for potential allosteric effectors. The server can be instrumental in elucidating of allosteric mechanisms and actions of allosteric mutations, and in the efforts on design of new elements of allosteric control. The server is freely available at: http://allosigma.bii.a-star.edu.sg.

Disorder driven allosteric control of protein activity

Studies of protein allostery increasingly reveal an involvement of the back and forth order-disorder transitions in this mechanism of protein activity regulation. Here, we investigate the allosteric mechanisms mediated by structural disorder using the structure-based statistical mechanical model of allostery (SBSMMA) that we have previously developed. We show that SBSMMA accounts for the energetics and causality of allosteric communication underlying dimerization of the BirA biotin repressor, activation of the sortase A enzyme, and inhibition of the Rac1 GTPase. Using the SBSMMA, we also show that introducing structural order or disorder in various regions of esterases can originate tunable allosteric modulation of the catalytic triad. On the basis of obtained results, we propose that operating with the order-disorder continuum allows one to establish an allosteric control scale for achieving desired modulation of the protein activity.

•

Back and forth order-disorder transitions can induce allosteric signaling.

•

Allosteric signaling originated by order/disorder follow universal rules.

•

Allosteric control scale facilitates engineering of the protein activity regulation.

Sequence-Based Prediction of Metamorphic Behavior in Proteins

An increasing number of proteins have been demonstrated in recent years to adopt multiple three-dimensional folds with different functions. These metamorphic proteins are characterized by having two or more folds with significant differences in their secondary structure, in which each fold is stabilized by a distinct local environment. So far, ∼90 metamorphic proteins have been identified in the Protein Databank, but we and others hypothesize that a far greater number of metamorphic proteins remain undiscovered. In this work, we introduce a computational model to predict metamorphic behavior in proteins using only knowledge of the sequence. In this model, secondary structure prediction programs are used to calculate diversity indices, which are measures of uncertainty in predicted secondary structure at each position in the sequence; these are then used to assign protein sequences as likely to be metamorphic versus monomorphic (i.e., having just one fold). We constructed a reference data set to train our classification method, which includes a novel compilation of 136 likely monomorphic proteins and a set of 201 metamorphic protein structures taken from the literature. Our model is able to classify proteins as metamorphic versus monomorphic with a Matthews correlation coefficient of ∼0.36 and true positive/true negative rates of ∼65%/80%, suggesting that it is possible to predict metamorphic behavior in proteins using only sequence information.

Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

The zinc transcriptional regulator (ZitR) functions as a metalloregulator that fine tunes transcriptional regulation through zinc-dependent DNA binding. However, the molecular mechanism of zinc-driven allosteric control of the DNA binding to ZitR remains elusive. Here, we performed enhanced sampling accelerated molecular dynamics simulations to figure out the mechanism, revealing the role of protein dynamics in the zinc-induced allosteric control of DNA binding to ZitR. The results suggest that zinc-free ZitR samples distinct conformational states, only a handful of which are compatible with DNA binding. Remarkably, zinc binding reduces the conformational plasticity of the DNA-binding domain of ZitR, promoting the population shift in the ZitR conformational ensemble towards the DNA binding-competent conformation. Further co-binding of DNA to the zinc–ZitR complex stabilizes this competent conformation. These findings suggest that ZitR–DNA interactions are allosterically regulated in a zinc-mediated conformational preselection manner, highlighting the importance of conformational dynamics in the regulation of transcription factor family.

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.

Predicted dynamical couplings of protein residues characterize catalysis, transport and allostery

MOTIVATION

Protein function is intrinsically linked to native dynamics, but the systematic characterization of functionally relevant dynamics remains elusive besides specific examples. Here we exhaustively characterize three types of dynamical couplings between protein residues: co-directionality (moving along collinear directions), coordination (small fluctuations of the interatomic distance) and deformation (the extent by which perturbations applied at one residue modify the local structure of the other one), which we analytically compute through the torsional network model.

RESULTS

We find that ligand binding sites are characterized by large within-site coordination and co-directionality, much larger than expected for generic sets of residues with equivalent sequence distances. In addition, catalytic sites are characterized by high coordination couplings with other residues in the protein, supporting the view that the overall protein structure facilitates the catalytic dynamics. The binding sites of allosteric effectors are characterized by comparably smaller coordination and higher within-site deformation than other ligands, which supports their dynamic nature. Allosteric inhibitors are coupled to the active site more frequently through deformation than through coordination, while the contrary holds for activators. We characterize the dynamical couplings of the sodium-dependent Leucine transporter protein (LeuT). The couplings between and within sites progress consistently along the transport cycle, providing a mechanistic description of the coupling between the uptake and release of ions and substrate, and they highlight qualitative differences between the wild-type and a mutant for which chloride is necessary for transport.

AVAILABILITY AND IMPLEMENTATION

The program tnm is freely available at https://github.com/ugobas/tnm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of Alternative Allosteric Sites in Glycolytic Enzymes for Potential Use as Species-Specific Drug Targets

Three allosteric glycolytic enzymes, phosphofructokinase, glyceraldehyde-3 phosphate dehydrogenase and pyruvate kinase, associated with bacterial, parasitic and human species, were explored to identify potential allosteric sites that would be used as prime targets for species-specific drug design purposes using a newly developed approach which incorporates solvent mapping, elastic network modeling, sequence and structural alignments. The majority of binding sites detected by solvent mapping overlapped with the interface regions connecting the subunits, thus appeared as promising target sites for allosteric regulation. Each binding site was then evaluated by its ability to alter the global dynamics of the receptor defined by the percentage change in the frequencies of the lowest-frequency modes most significantly and as anticipated, the most effective ones were detected in the vicinity of the well-reported catalytic and allosteric sites. Furthermore, some of our proposed regions intersected with experimentally resolved sites which are known to be critical for activity regulation, which further validated our approach. Despite the high degree of structural conservation encountered between bacterial/parasitic and human glycolytic enzymes, the majority of the newly presented allosteric sites exhibited a low degree of sequence conservation which further increased their likelihood to be used as species-specific target regions for drug design studies.

AlloMAPS: allosteric mutation analysis and polymorphism of signaling database

AlloMAPS database provides data on the causality and energetics of allosteric communication obtained with the structure-based statistical mechanical model of allostery (SBSMMA). The database contains data on allosteric signaling in three sets of proteins and protein chains: (i) 46 proteins with comprehensively annotated functional and allosteric sites; (ii) 1908 protein chains from PDBselect set of chains with low (<25%) sequence identity; (iii) 33 proteins with more than 50 known pathological SNPs in each molecule. In addition to energetics of allosteric signaling between known functional and regulatory sites, allosteric modulation caused by the binding to these sites, by SNPs, and by mutations designated by the user can be explored. Allosteric Signaling Maps (ASMs), which are produced via the exhaustive computational scanning for stabilizing and destabilizing mutations and for the modulation range caused by the sequence position are available for each protein/protein chain in the database. We propose to use this database for evaluating the effects of allosteric signaling in the search for latent regulatory sites and in the design of allosteric sites and effectors. The database is freely available at: http://allomaps.bii.a-star.edu.sg.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Leveraging protein dynamics to identify cancer mutational hotspots using 3D structures

Large-scale exome sequencing of tumors has enabled the identification of cancer drivers using recurrence-based approaches. Some of these methods also employ 3D protein structures to identify mutational hotspots in cancer-associated genes. In determining such mutational clusters in structures, existing approaches overlook protein dynamics, despite its essential role in protein function. We present a framework to identify cancer driver genes using a dynamics-based search of mutational hotspot communities. Mutations are mapped to protein structures, which are partitioned into distinct residue communities. These communities are identified in a framework where residue-residue contact edges are weighted by correlated motions (as inferred by dynamics-based models). We then search for signals of positive selection among these residue communities to identify putative driver genes, while applying our method to the TCGA (The Cancer Genome Atlas) PanCancer Atlas missense mutation catalog. Overall, we predict 1 or more mutational hotspots within the resolved structures of proteins encoded by 434 genes. These genes were enriched among biological processes associated with tumor progression. Additionally, a comparison between our approach and existing cancer hotspot detection methods using structural data suggests that including protein dynamics significantly increases the sensitivity of driver detection.

Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks

BACKGROUND

Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks.

METHODS

In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems.

RESULTS

The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators.

CONCLUSIONS

We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases.

GENERAL SIGNIFICANCE

This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.

Toward Comprehensive Allosteric Control over Protein Activity

Universality of allosteric signaling in proteins, molecular machines, and receptors complemented by the great advantages of prospected allosteric drugs in the highly specific, non-competitive, and modulatory nature of their actions calls for deeper theoretical understanding of allosteric communication. We present a computational model that makes it possible to tackle the problem of modulating the energetics of protein allosteric communication. In the context of the energy landscape paradigm, allosteric signaling is always a result of perturbations, such as ligand binding, mutations, and intermolecular interactions. The calculation of local partition functions in the protein harmonic model with perturbations allows us to evaluate the energetics of allosteric communication at the single-residue level. In this framework, Allosteric Signaling Maps are proposed as a tool to exhaustively describe allosteric communication in the protein, to tune already existing signaling, and to design new elements of regulation for taking the protein activity under allosteric control.

Correlation Between Allosteric and Orthosteric Sites

Correlation between an allosteric site and its orthosteric site refers to the phenomenon that perturbations like ligand binding, mutation, or posttranslational modifications at the allosteric site leverage variation in the orthosteric site. Understanding this kind of correlation not only helps to disclose how information is transmitted in allosteric regulation but also provides clues for allosteric drug discovery. This chapter starts with an overview of correlation studies on allosteric and orthosteric sites and then introduces recent progress in evolutionary and simulation-based dynamic studies. Discussions and perspectives on future directions are also given.

Reversing allosteric communication: From detecting allosteric sites to inducing and tuning targeted allosteric response

The omnipresence of allosteric regulation together with the fundamental role of structural dynamics in this phenomenon have initiated a great interest to the detection of regulatory exosites and design of corresponding effectors. However, despite a general consensus on the key role of dynamics most of the earlier efforts on the prediction of allosteric sites are heavily crippled by the static nature of the underlying methods, which are either structure-based approaches seeking for deep surface pockets typical for “traditional” orthosteric drugs or sequence-based techniques exploiting the conservation of protein sequences. Because of the critical role of global protein dynamics in allosteric signaling, we investigate the hypothesis of reversibility in allosteric communication, according to which allosteric sites can be detected via the perturbation of the functional sites. The reversibility is tested here using our structure-based perturbation model of allostery, which allows one to analyze the causality and energetics of allosteric communication. We validate the “reverse perturbation” hypothesis and its predictive power on a set of classical allosteric proteins, then, on the independent extended benchmark set. We also show that, in addition to known allosteric sites, the perturbation of the functional sites unravels rather extended protein regions, which can host latent regulatory exosites. These protein parts that are dynamically coupled with functional sites can also be used for inducing and tuning allosteric communication, and an exhaustive exploration of the per-residue contributions to allosteric effects can eventually lead to the optimal modulation of protein activity. The site-effector interactions necessary for a specific mode and level of allosteric communication can be fine-tuned by adjusting the site’s structure to an available effector molecule and by the design or selection of an appropriate ligand.

Recent advances in the development of allosteric drugs allow one to fully appreciate the sheer power of allosteric effectors in the avoiding toxicity, receptor desensitization and modulatory rather than on/off mode of action, compared to the traditional orthosteric compounds. The detection of allosteric sites is one of the major challenges in the quest for allosteric drugs. This work proposes a “reverse perturbation” approach for identifying allosteric sites as a result of a perturbation applied to the functional ones. We show that according to the traditional Monod-Changeux-Jacob’s definition of allostery, considering non-overlapping regulatory and functional sites is a critical prerequisite for the successful detection of allosteric sites. Using the reverse perturbation method, it is possible to determine wide protein regions with a potential to induce an allosteric response and to adjust its strength. Further studies on inducing and fine-tuning of allosteric signalling seem to be of a great importance for efficient design of non-orthosteric ligands in the development of novel drugs.

Identifying Protein Allosteric Transitions for Drug Discovery with 1D NMR

Allosteric drugs present many advantages over orthosteric drugs and are therefore an attractive approach in drug discovery, despite being highly challenging. First, the binding of ligands in protein allosteric pockets do not ensure an allosteric effect, and second, allosteric ligands can possess diverse modes of pharmacology even within a compound family. Herein we report a new method to: 1) detect allosteric communication between protein binding sites, and 2) compare the effect of allosteric ligands on the allosteric transitions of the protein target. The method, illustrated with glycogen phosphorylase, consists of comparing 1D saturation transfer difference (STD) NMR spectra of a molecular spy (here fragments) in the absence and presence of allosteric ligands. The modification of the STD NMR spectrum of the fragment indicates whether the protein dynamics/conformations have been changed in the presence of the allosteric modulator, thereby highlighting allosteric coupling between the binding pocket of the reference compound (in this case the fragment) and the allosteric pocket.

Predicting Protein Dynamics and Allostery Using Multi-Protein Atomic Distance Constraints

The related concepts of protein dynamics, conformational ensembles and allostery are often difficult to study with molecular dynamics (MD) due to the timescales involved. We present ExProSE (Exploration of Protein Structural Ensembles), a distance geometry-based method that generates an ensemble of protein structures from two input structures. ExProSE provides a unified framework for the exploration of protein structure and dynamics in a fast and accessible way. Using a dataset of apo/holo pairs it is shown that existing coarse-grained methods often cannot span large conformational changes. For T4-lysozyme, ExProSE is able to generate ensembles that are more native-like than tCONCOORD and NMSim, and comparable with targeted MD. By adding additional constraints representing potential modulators, ExProSE can predict allosteric sites. ExProSE ranks an allosteric pocket first or second for 27 out of 58 allosteric proteins, which is similar and complementary to existing methods. The ExProSE source code is freely available.

Proteome-Scale Investigation of Protein Allosteric Regulation Perturbed by Somatic Mutations in 7,000 Cancer Genomes

The allosteric regulation triggering the protein's functional activity via conformational changes is an intrinsic function of protein under many physiological and pathological conditions, including cancer. Identification of the biological effects of specific somatic variants on allosteric proteins and the phenotypes that they alter during tumor initiation and progression is a central challenge for cancer genomes in the post-genomic era. Here, we mapped more than 47,000 somatic missense mutations observed in approximately 7,000 tumor-normal matched samples across 33 cancer types into protein allosteric sites to prioritize the mutated allosteric proteins and we tested our prediction in cancer cell lines. We found that the deleterious mutations identified in cancer genomes were more significantly enriched at protein allosteric sites than tolerated mutations, suggesting a critical role for protein allosteric variants in cancer. Next, we developed a statistical approach, namely AlloDriver, and further identified 15 potential mutated allosteric proteins during pan-cancer and individual cancer-type analyses. More importantly, we experimentally confirmed that p.Pro360Ala on PDE10A played a potential oncogenic role in mediating tumorigenesis in non-small cell lung cancer (NSCLC). In summary, these findings shed light on the role of allosteric regulation during tumorigenesis and provide a useful tool for the timely development of targeted cancer therapies.

Structural analyses of 2015-updated drug-resistant mutations in HIV-1 protease: an implication of protease inhibitor cross-resistance

Background

Strategies to control HIV for improving the quality of patient lives have been aided by the Highly Active Anti-Retroviral Therapy (HAART), which consists of a cocktail of inhibitors targeting key viral enzymes. Numerous new drugs have been developed over the past few decades but viral resistances to these drugs in the targeted viral enzymes are increasingly reported. Nonetheless the acquired mutations often reduce viral fitness and infectivity. Viral compensatory secondary-line mutations mitigate this loss of fitness, equipping the virus with a broad spectrum of resistance against these drugs. While structural understanding of the viral protease and its drug resistance mutations have been well established, the interconnectivity and development of structural cross-resistance remain unclear. This paper reports the structural analyses of recent clinical mutations on the drug cross-resistance effects from various protease and protease inhibitors (PIs) complexes.

Methods

Using the 2015 updated clinical HIV protease mutations, we constructed a structure-based correlation network and a minimum-spanning tree (MST) based on the following features: (i) topology of the PI-binding pocket, (ii) allosteric effects of the mutations, and (iii) protease structural stability.

Results and conclusion

Analyis of the network and the MST of dominant mutations conferring resistance to the seven PIs (Atazanavir-ATV, Darunavir-DRV, Indinavir-IDV, Lopinavir-LPV, Nelfinavir-NFV, Saquinavir-SQV, and Tipranavir-TPV) showed that cross-resistance can develop easily across NFV, SQV, LPV, IDV, and DRV, but not for ATV or TPV. Through estimation of the changes in vibrational entropies caused by each reported mutation, some secondary mutations were found to destabilize protease structure. Our findings provide an insight into the mechanism of PI cross-resistance and may also be useful in guiding the selection of PI in clinical treatment to delay the onset of cross drug resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1372-3) contains supplementary material, which is available to authorized users.

Toward Allosterically Increased Catalytic Activity of Insulin-Degrading Enzyme against Amyloid Peptides

The physiological role of insulin-degrading enzyme (IDE) in the intracytosolic clearance of amyloid β (Aβ) and other amyloid-like peptides supports a hypothesis that human IDE hyperactivation could be therapeutically beneficial for the treatment of late-onset Alzheimer's disease (AD). The major challenge standing in the way of this goal is increasing the specific catalytic activity of IDE against the Aβ substrate. There were previous indications that the allosteric mode of IDE activity regulation could potentially provide a highly specific path toward degradation of amyloid-like peptides, while not dramatically affecting activity against other substrates. Recently developed theoretical concepts are used here to explore potential allosteric modulation of the IDE activity as a result of single-residue mutations. Five candidates are selected for experimental follow-up and allosteric free energy calculations: Ser137Ala, Lys396Ala, Asp426Ala, Phe807Ala, and Lys898Ala. Our experiments show that three mutations (Ser137Ala, Phe807Ala, and Lys898Ala) decrease the K_m of the Aβ substrate. Mutation Lys898Ala results in increased catalytic activity of IDE; on the other hand, Lys364Ala does not change the activity and Asp426Ala diminishes it. Quantifying effects of mutations in terms of allosteric free energy, we show that favorable mutations lead to stabilization of the catalytic sites and other function-relevant distal sites as well as increased dynamics of the IDE-N and IDE-C halves that allow efficient substrate entrance and cleavage. A possibility for intramolecular upregulation of IDE activity against amyloid peptides via allosteric mutations calls for further investigations in this direction. Ultimately, we are hopeful it will lead to the development of IDE-based drugs for the treatment of the late-onset form of AD characterized by an overall impairment of Aβ clearance.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Protein function machinery: from basic structural units to modulation of activity

Contemporary protein structure is a result of the trade off between the laws of physics and the evolutionary selection. The polymer nature of proteins played a decisive role in establishing the basic structural and functional units of soluble proteins. We discuss how these elementary building blocks work in the hierarchy of protein domain structure, co-translational folding, as well as in enzymatic activity and molecular interactions. Next, we consider modulators of the protein function, such as intermolecular interactions, disorder-to-order transitions, and allosteric signaling, acting via interference with the protein's structural dynamics. We also discuss the post-translational modifications, which is a complementary intricate mechanism evolved for regulation of protein functions and interactions. In conclusion, we assess an anticipated contribution of discussed topics to the future advancements in the field.

Nonlinear backbone torsional pair correlations in proteins

Protein allostery requires dynamical structural correlations. Physical origin of which, however, remain elusive despite intensive studies during last two and half decades. Based on analysis of molecular dynamics (MD) simulation trajectories for ten proteins with different sizes and folds, we found that nonlinear backbone torsional pair (BTP) correlations, which are mainly spatially long-ranged and are dominantly executed by loop residues, exist extensively in most analyzed proteins. Examination of torsional motion for correlated BTPs suggested that such nonlinear correlations are mainly associated aharmonic torsional state transitions and in some cases strongly anisotropic local torsional motion of participating torsions, and occur on widely different and relatively longer time scales. In contrast, correlations between backbone torsions in stable α helices and β strands are mainly linear and spatially short-ranged, and are more likely to associate with harmonic local torsional motion. Further analysis revealed that the direct cause of nonlinear contributions are heterogeneous linear correlations. These findings implicate a general search strategy for novel allosteric modulation sites of protein activities.

Computational approaches to investigating allostery

Allosteric regulation plays a key role in many biological processes, such as signal transduction, transcriptional regulation, and many more. It is rooted in fundamental thermodynamic and dynamic properties of macromolecular systems that are still poorly understood and are moreover modulated by the cellular context. Here we review the computational approaches used in the investigation of allosteric processes in protein systems. We outline how the models of allostery have evolved from their initial formulation in the sixties to the current views, which more fully account for the roles of the thermodynamic and dynamic properties of the system. We then describe the major classes of computational approaches employed to elucidate the mechanisms of allostery, the insights they have provided, as well as their limitations. We complement this analysis by highlighting the role of computational approaches in promising practical applications, such as the engineering of regulatory modules and identifying allosteric binding sites.

Regulatory network of the allosteric ATP inhibition of E. coli phosphofructokinase-2 studied by hybrid dimers

We have proposed an allosteric ATP inhibition mechanism of Pfk-2 determining the structure of different forms of the enzyme together with a kinetic enzyme analysis. Here we complement the mechanism by using hybrid oligomers of the homodimeric enzyme to get insights about the allosteric communication pathways between the same sites or different ones located in different subunits. Kinetic analysis of the hybrid enzymes indicate that homotropic interactions between allosteric sites for ATP or between substrate sites for fructose-6-P have a minor effect on the enzymatic inhibition induced by ATP. In fact, the sigmoid response for fructose-6-P observed at elevated ATP concentrations can be eliminated even though the enzymatic inhibition is still operative. Nevertheless, leverage coupling analysis supports heterotropic interactions between the allosteric ATP and fructose-6-P binding occurring between and within each subunit.

Identifying Allosteric Hotspots with Dynamics: Application to Inter- and Intra-species Conservation

The rapidly growing volume of data being produced by next-generation sequencing initiatives is enabling more in-depth analyses of conservation than previously possible. Deep sequencing is uncovering disease loci and regions under selective constraint, despite the fact that intuitive biophysical reasons for such constraint are sometimes absent. Allostery may often provide the missing explanatory link. We use models of protein conformational change to identify allosteric residues by finding essential surface pockets and information-flow bottlenecks, and we develop a software tool that enables users to perform this analysis on their own proteins of interest. Though fundamentally 3D-structural in nature, our analysis is computationally fast, thereby allowing us to run it across the PDB and to evaluate general properties of predicted allosteric residues. We find that these tend to be conserved over diverse evolutionary time scales. Finally, we highlight examples of allosteric residues that help explain poorly understood disease-associated variants.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

Allosteric drug development holds promise for delivering medicines that are more selective and less toxic than those that target orthosteric sites. To date, the discovery of allosteric binding sites and lead compounds has been mostly serendipitous, achieved through high-throughput screening. Over the past decade, structural data has become more readily available for larger protein systems and more membrane protein classes (e.g., GPCRs and ion channels), which are common allosteric drug targets. In parallel, improved simulation methods now provide better atomistic understanding of the protein dynamics and cooperative motions that are critical to allosteric mechanisms. As a result of these advances, the field of predictive allosteric drug development is now on the cusp of a new era of rational structure-based computational methods. Here, we review algorithms that predict allosteric sites based on sequence data and molecular dynamics simulations, describe tools that assess the druggability of these pockets, and discuss how Markov state models and topology analyses provide insight into the relationship between protein dynamics and allosteric drug binding. In each section, we first provide an overview of the various method classes before describing relevant algorithms and software packages.

Discovery of a small-molecule binder of the oncoprotein gankyrin that modulates gankyrin activity in the cell

Gankyrin is an ankyrin-repeat oncoprotein whose overexpression has been implicated in the development of many cancer types. Elevated gankyrin levels are linked to aberrant cellular events including enhanced degradation of tumour suppressor protein p53, and inhibition of gankyrin activity has therefore been identified as an attractive anticancer strategy. Gankyrin interacts with several partner proteins, and a number of these protein-protein interactions (PPIs) are of relevance to cancer. Thus, molecules that bind the PPI interface of gankyrin and interrupt these interactions are of considerable interest. Herein, we report the discovery of a small molecule termed cjoc42 that is capable of binding to gankyrin. Cell-based experiments demonstrate that cjoc42 can inhibit gankyrin activity in a dose-dependent manner: cjoc42 prevents the decrease in p53 protein levels normally associated with high amounts of gankyrin, and it restores p53-dependent transcription and sensitivity to DNA damage. The results represent the first evidence that gankyrin is a "druggable" target with small molecules.

Structure-Based Statistical Mechanical Model Accounts for the Causality and Energetics of Allosteric Communication

Allostery is one of the pervasive mechanisms through which proteins in living systems carry out enzymatic activity, cell signaling, and metabolism control. Effective modeling of the protein function regulation requires a synthesis of the thermodynamic and structural views of allostery. We present here a structure-based statistical mechanical model of allostery, allowing one to observe causality of communication between regulatory and functional sites, and to estimate per residue free energy changes. Based on the consideration of ligand free and ligand bound systems in the context of a harmonic model, corresponding sets of characteristic normal modes are obtained and used as inputs for an allosteric potential. This potential quantifies the mean work exerted on a residue due to the local motion of its neighbors. Subsequently, in a statistical mechanical framework the entropic contribution to allosteric free energy of a residue is directly calculated from the comparison of conformational ensembles in the ligand free and ligand bound systems. As a result, this method provides a systematic approach for analyzing the energetics of allosteric communication based on a single structure. The feasibility of the approach was tested on a variety of allosteric proteins, heterogeneous in terms of size, topology and degree of oligomerization. The allosteric free energy calculations show the diversity of ways and complexity of scenarios existing in the phenomenology of allosteric causality and communication. The presented model is a step forward in developing the computational techniques aimed at detecting allosteric sites and obtaining the discriminative power between agonistic and antagonistic effectors, which are among the major goals in allosteric drug design.

The 50th anniversary of Monod-Changeux-Jacob seminal paper “Allosteric proteins and cellular control systems” became the hallmark of a new wave in the allostery studies and the turning point in our vision of allostery and its implications in protein engineering and drug design. Recent experimental and theoretical works clearly show relevance of allosteric phenomenon to drug design, unraveling advantages of allosteric drugs in comparison to traditional orthosteric compounds. Remarkable simplicity of allosteric effectors and, at the same time, their potentially high specificity is one of the most important traits. The non conserved nature of allosteric ligands is a basis for avoiding drug resistance, and existence of latent regulatory sites make them attractive drug targets. The model presented in this work provides a theoretical framework for the quantification of the causality and energetics of allosteric regulation, which is a prerequisite for design of effector molecules with required characteristics. The synthesis between the thermodynamics of allostery and the intrinsic atomic nature of proteins and their interactions with the allosteric effectors accomplished in this work is a small initial step in the long endeavor towards future allosteric drugs.