The origin of allosteric functional modulation: multiple pre-existing pathways

Allostery in disease and in drug discovery

Allostery is largely associated with conformational and functional transitions in individual proteins. This concept can be extended to consider the impact of conformational perturbations on cellular function and disease states. Here, we clarify the concept of allostery and how it controls physiological activities. We focus on the challenging questions of how allostery can both cause disease and contribute to development of new therapeutics. We aim to increase the awareness of the linkage between disease symptoms on the cellular level and specific aberrant allosteric actions on the molecular level and to emphasize the potential of allosteric drugs in innovative therapies.

Multiple pathways promote dynamical coupling between catalytic domains in Escherichia coli prolyl-tRNA synthetase

Aminoacyl-tRNA synthetases are multidomain enzymes that catalyze covalent attachment of amino acids to their cognate tRNA. Cross-talk between functional domains is a prerequisite for this process. In this study, we investigate the molecular mechanism of site-to-site communication in Escherichia coli prolyl-tRNA synthetase (Ec ProRS). Earlier studies have demonstrated that evolutionarily conserved and/or co-evolved residues that are engaged in correlated motion are critical for the propagation of functional conformational changes from one site to another in modular proteins. Here, molecular simulation and bioinformatics-based analysis were performed to identify dynamically coupled and evolutionarily constrained residues that form contiguous pathways of residue-residue interactions between the aminoacylation and editing domains of Ec ProRS. The results of this study suggest that multiple pathways exist between these two domains to maintain the dynamic coupling essential for enzyme function. Moreover, residues in these interaction networks are generally highly conserved. Site-directed changes of on-pathway residues have a significant impact on enzyme function and dynamics, suggesting that any perturbation along these pathways disrupts the native residue-residue interactions that are required for effective communication between the two functional domains. Free energy analysis revealed that communication between residues within a pathway and cross-talk between pathways are important for coordinating functions of different domains of Ec ProRS for efficient catalysis.

MCPath: Monte Carlo path generation approach to predict likely allosteric pathways and functional residues

Allosteric mechanism of proteins is essential in biomolecular signaling. An important aspect underlying this mechanism is the communication pathways connecting functional residues. Here, a Monte Carlo (MC) path generation approach is proposed and implemented to define likely allosteric pathways through generating an ensemble of maximum probability paths. The protein structure is considered as a network of amino acid residues, and inter-residue interactions are described by an atomistic potential function. PDZ domain structures are presented as case studies. The analysis for bovine rhodopsin and three myosin structures are also provided as supplementary case studies. The suggested pathways and the residues constituting the pathways are maximally probable and mostly agree with the previous studies. Overall, it is demonstrated that the communication pathways could be multiple and intrinsically disposed, and the MC path generation approach provides an effective tool for the prediction of key residues that mediate the allosteric communication in an ensemble of pathways and functionally plausible residues. The MCPath server is available at http://safir.prc.boun.edu.tr/clbet_server.

A mechanistic understanding of allosteric immune escape pathways in the HIV-1 envelope glycoprotein

The HIV-1 envelope (Env) spike, which consists of a compact, heterodimeric trimer of the glycoproteins gp120 and gp41, is the target of neutralizing antibodies. However, the high mutation rate of HIV-1 and plasticity of Env facilitates viral evasion from neutralizing antibodies through various mechanisms. Mutations that are distant from the antibody binding site can lead to escape, probably by changing the conformation or dynamics of Env; however, these changes are difficult to identify and define mechanistically. Here we describe a network analysis-based approach to identify potential allosteric immune evasion mechanisms using three known HIV-1 Env gp120 protein structures from two different clades, B and C. First, correlation and principal component analyses of molecular dynamics (MD) simulations identified a high degree of long-distance coupled motions that exist between functionally distant regions within the intrinsic dynamics of the gp120 core, supporting the presence of long-distance communication in the protein. Then, by integrating MD simulations with network theory, we identified the optimal and suboptimal communication pathways and modules within the gp120 core. The results unveil both strain-dependent and -independent characteristics of the communication pathways in gp120. We show that within the context of three structurally homologous gp120 cores, the optimal pathway for communication is sequence sensitive, i.e. a suboptimal pathway in one strain becomes the optimal pathway in another strain. Yet the identification of conserved elements within these communication pathways, termed inter-modular hotspots, could present a new opportunity for immunogen design, as this could be an additional mechanism that HIV-1 uses to shield vulnerable antibody targets in Env that induce neutralizing antibody breadth.

Phosphorylation of the retinoic acid receptor alpha induces a mechanical allosteric regulation and changes in internal dynamics

Nuclear receptor proteins constitute a superfamily of proteins that function as ligand dependent transcription factors. They are implicated in the transcriptional cascades underlying many physiological phenomena, such as embryogenesis, cell growth and differentiation, and apoptosis, making them one of the major signal transduction paradigms in metazoans. Regulation of these receptors occurs through the binding of hormones, and in the case of the retinoic acid receptor (RAR), through the binding of retinoic acid (RA). In addition to this canonical scenario of RAR activity, recent discoveries have shown that RAR regulation also occurs as a result of phosphorylation. In fact, RA induces non-genomic effects, such as the activation of kinase signaling pathways, resulting in the phosphorylation of several targets including RARs themselves. In the case of RARα, phosphorylation of Ser369 located in loop L9–10 of the ligand-binding domain leads to an increase in the affinity for the protein cyclin H, which is part of the Cdk-activating kinase complex of the general transcription factor TFIIH. The cyclin H binding site in RARα is situated more than 40 Å from the phosphorylated serine. Using molecular dynamics simulations of the unphosphorylated and phosphorylated forms of the receptor RARα, we analyzed the structural implications of receptor phosphorylation, which led to the identification of a structural mechanism for the allosteric coupling between the two remote sites of interest. The results show that phosphorylation leads to a reorganization of a local salt bridge network, which induces changes in helix extension and orientation that affects the cyclin H binding site. This results in changes in conformation and flexibility of the latter. The high conservation of the residues implicated in this signal transduction suggests a mechanism that could be applied to other nuclear receptor proteins.

Allosteric regulation of proteins is critically important in many biological processes. Here we focused on the allosteric pathway of communication within a ligand-regulated transcription factor, the Retinoic Acid Receptor (RAR). Recent experimental studies performed with the RARα subtype have shown that phosphorylation of a residue located at one extremity of an α-helix in RAR, leads to a changes in binding affinity at the other extremity of the same helix for cyclin H, a binding partner that is necessary for gene transcription activation. The purpose of our study was to understand the conformational changes occurring within the receptor upon phosphorylation. Molecular dynamics simulations are well suited for this sort of study. Through this approach, we were able to show that although the overall structure of the phosphorylated RAR shows no distinct difference from the unphosphorylated form, evidence is provided for an allosteric regulation pathway that implicates more subtle changes, such as changes in side chain orientations, which affect the internal dynamics of the receptor.

Tertiary and quaternary effects in the allosteric regulation of animal hemoglobins

In the last decade, protein allostery has experienced a major resurgence, boosted by the extension of the concept to systems of increasing complexity and by its exploitation for the development of drugs. Expansion of the field into new directions has not diminished the key role of hemoglobin as a test molecule for theory and experimental validation of allosteric models. Indeed, the diffusion of hemoglobins in all kingdoms of life and the variety of functions and of quaternary assemblies based on a common tertiary fold indicate that this superfamily of proteins is ideally suited for investigating the physical and molecular basis of allostery and firmly maintains its role as a main player in the field. This review is an attempt to briefly recollect common and different strategies adopted by metazoan hemoglobins, from monomeric molecules to giant complexes, exploiting homotropic and heterotropic allostery to increase their functional dynamic range. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

The underappreciated role of allostery in the cellular network

Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.

Characterizing conformation changes in proteins through the torsional elastic response

The relationship between functional conformation changes and thermal dynamics of proteins is investigated with the help of the torsional network model (TNM), an elastic network model in torsion angle space that we recently introduced. We propose and test a null-model of "random" conformation changes that assumes that the contributions of normal modes to conformation changes are proportional to their contributions to thermal fluctuations. Deviations from this null model are generally small. When they are large and significant, they consist in conformation changes that are represented by very few low frequency normal modes and overcome small energy barriers. We interpret these features as the result of natural selection favoring the intrinsic protein dynamics consistent with functional conformation changes. These "selected" conformation changes are more frequently associated to ligand binding, and in particular phosphorylation, than to pairs of conformations with the same ligands. This deep relationship between the thermal dynamics of a protein, represented by its normal modes, and its functional dynamics can reconcile in a unique framework the two models of conformation changes, conformational selection and induced fit. The program TNM that computes torsional normal modes and analyzes conformation changes is available upon request. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

An interdomain energetic tug-of-war creates the allosterically active state in Hsp70 molecular chaperones

The allosteric mechanism of Hsp70 molecular chaperones enables ATP binding to the N-terminal nucleotide-binding domain (NBD) to alter substrate affinity to the C-terminal substrate-binding domain (SBD) and substrate binding to enhance ATP hydrolysis. Cycling between ATP-bound and ADP/substrate-bound states requires Hsp70s to visit a state with high ATPase activity and fast on/off kinetics of substrate binding. We have trapped this "allosterically active" state for the E. coli Hsp70, DnaK, and identified how interactions among the NBD, the β subdomain of the SBD, the SBD α-helical lid, and the conserved hydrophobic interdomain linker enable allosteric signal transmission between ligand-binding sites. Allostery in Hsp70s results from an energetic tug-of-war between domain conformations and formation of two orthogonal interfaces: between the NBD and SBD, and between the helical lid and the β subdomain of the SBD. The resulting energetic tension underlies Hsp70 functional properties and enables them to be modulated by ligands and cochaperones and "tuned" through evolution.

Substrate and inhibitor-induced dimerization and cooperativity in caspase-1 but not caspase-3

Caspases are intracellular cysteine-class proteases with aspartate specificity that is critical for driving processes as diverse as the innate immune response and apoptosis, exemplified by caspase-1 and caspase-3, respectively. Interestingly, caspase-1 cleaves far fewer cellular substrates than caspase-3 and also shows strong positive cooperativity between the two active sites of the homodimer, unlike caspase-3. Biophysical and kinetic studies here present a molecular basis for this difference. Analytical ultracentrifugation experiments show that mature caspase-1 exists predominantly as a monomer under physiological concentrations that undergoes dimerization in the presence of substrate; specifically, substrate binding shifts the KD for dimerization by 20-fold. We have created a hemi-active site-labeled dimer of caspase-1, where one site is blocked with the covalent active site inhibitor, benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone. This hemi-labeled enzyme is about 9-fold more active than the apo-dimer of caspase-1. These studies suggest that substrate not only drives dimerization but also, once bound to one site in the dimer, promotes an active conformation in the other monomer. Steady-state kinetic analysis and modeling independently support this model, where binding of one substrate molecule not only increases substrate binding in preformed dimers but also drives the formation of heterodimers. Thus, the cooperativity in caspase-1 is driven both by substrate-induced dimerization as well as substrate-induced activation. Substrate-induced dimerization and activation seen in caspase-1 and not in caspase-3 may reflect their biological roles. Whereas caspase-1 cleaves a dramatically smaller number of cellular substrates that need to be concentrated near inflammasomes, caspase-3 is a constitutively active dimer that cleaves many more substrates located diffusely throughout the cell.

Thermodynamics of allostery paves a way to allosteric drugs

We overview here our recent work on the thermodynamic view of allosteric regulation and communication. Starting from the geometry-based prediction of regulatory binding sites in a static structure, we move on to exploring a connection between ligand binding and the intrinsic dynamics of the protein molecule. We describe here two recently introduced measures, binding leverage and leverage coupling, which allow one to analyze the molecular basis of allosteric regulation. We discuss the advantages of these measures and show that they work universally in proteins of different sizes, oligomeric states, and functions. We also point the problems that have to be solved before completing an atomic level description of allostery, and briefly discuss ideas for computational design of allosteric drugs. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Allosteric inhibition of a zinc-sensing transcriptional repressor: insights into the arsenic repressor (ArsR) family

The molecular basis of allosteric regulation remains a subject of intense interest. Staphylococcus aureus CzrA is a member of the ubiquitous arsenic repressor (ArsR) family of bacterial homodimeric metal-sensing proteins and has emerged as a model system for understanding allosteric regulation of operator DNA binding by transition metal ions. Using unnatural amino acid substitution and a standard linkage analysis, we show that a His97' NH(ε2)...O=C His67 quaternary structural hydrogen bond is an energetically significant contributor to the magnitude of the allosteric coupling free energy, ∆Gc. A "cavity" introduced just beneath this hydrogen bond in V66A/L68V CzrA results in a significant reduction in regulation by Zn(II) despite adopting a wild-type global structure and Zn(II) binding and DNA binding affinities only minimally affected from wild type. The energetics of Zn(II) binding and heterotropic coupling free energies (∆Hc, -T∆Sc) of the double mutant are also radically altered and suggest that increased internal dynamics leads to poorer allosteric negative regulation in V66A/L68V CzrA. A statistical coupling analysis of 3000 ArsR proteins reveals a sector that links the DNA-binding determinants and the α5 Zn(II)-sensing sites through V66/L68 in CzrA. We propose that distinct regulatory sites uniquely characteristic of individual ArsR proteins result from evolution of distinct connectivities to this sector, each capable of driving the same biological outcome, transcriptional derepression.

Diversity of allosteric regulation in proteases

Allostery is a fundamental regulatory mechanism that is based on a functional modulation of a site by a distant site. Allosteric regulation can be triggered by binding of diverse allosteric effectors, ranging from small molecules to macromolecules, and is therefore offering promising opportunities for functional modulation in a wide range of applications including the development of chemical probes or drug discovery. Here, we provide an overview of key classes of allosteric protease effectors, their corresponding molecular mechanisms, and their practical implications.

Connecting protein conformational dynamics with catalytic function as illustrated in dihydrofolate reductase

Combined quantum mechanics/molecular mechanics molecular dynamics simulations reveal that the M20 loop conformational dynamics of dihydrofolate reductase (DHFR) is severely restricted at the transition state of the hydride transfer as a result of the M42W/G121V double mutation. Consequently, the double-mutant enzyme has a reduced entropy of activation, i.e., increased entropic barrier, and altered temperature dependence of kinetic isotope effects in comparison with those of wild-type DHFR. Interestingly, in both wild-type DHFR and the double mutant, the average donor-acceptor distances are essentially the same in the Michaelis complex state (~3.5 Å) and the transition state (2.7 Å). It was found that an additional hydrogen bond is formed to stabilize the M20 loop in the closed conformation in the M42W/G121V double mutant. The computational results reflect a similar aim designed to knock out precisely the dynamic flexibility of the M20 loop in a different double mutant, N23PP/S148A.

Analysis of 15N-1H NMR relaxation in proteins by a combined experimental and molecular dynamics simulation approach: picosecond-nanosecond dynamics of the Rho GTPase binding domain of plexin-B1 in the dimeric state indicates allosteric pathways

We investigate picosecond–nanosecond dynamics of the Rho-GTPase Binding Domain (RBD) of plexin-B1, which plays a key role in plexin-mediated cell signaling. Backbone 15N relaxation data of the dimeric RBD are analyzed with the model-free (MF) method, and with the slowly relaxing local structure/molecular dynamics (SRLS-MD) approach. Independent analysis of the MD trajectories, based on the MF paradigm, is also carried out. MF is a widely popular and simple method, SRLS is a general approach, and SRLS-MD is an integrated approach we developed recently. Corresponding parameters from the RBD dimer, a previously studied RBD monomer mutant, and the previously studied complex of the latter with the GTPase Rac1, are compared. The L2, L3, and L4 loops of the plexin-B1 RBD are involved in interactions with other plexin domains, GTPase binding, and RBD dimerization, respectively. Peptide groups in the loops of both the monomeric and dimeric RBD are found to experience weak and moderately asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes, and nanosecond backbone motion. Peptide groups in the α-helices and the β-strands of the dimer (the β-strands of the monomer) experience strong and highly asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes (N–H bonds). N–H fluctuations occur on the picosecond time scale. An allosteric pathway for GTPase binding, providing new insights into plexin function, is delineated.

A broad view of scaffolding suggests that scaffolding proteins can actively control regulation and signaling of multienzyme complexes through allostery

Enzymes often work sequentially in pathways; and consecutive reaction steps are typically carried out by molecules associated in the same multienzyme complex. Localization confines the enzymes; anchors them; increases the effective concentration of substrates and products; and shortens pathway timescales; however, it does not explain enzyme coordination or pathway branching. Here, we distinguish between metabolic and signaling multienzyme complexes. We argue for a central role of scaffolding proteins in regulating multienzyme complexes signaling and suggest that metabolic multienzyme complexes are less dependent on scaffolding because they undergo conformational control through direct subunit-subunit contacts. In particular, we propose that scaffolding proteins have an essential function in controlling branching in signaling pathways. This new broadened definition of scaffolding proteins goes beyond cases such as the classic yeast mitogen-activated protein kinase Ste5 and encompasses proteins such as E3 ligases which lack active sites and work via allostery. With this definition, we classify the mechanisms of multienzyme complexes based on whether the substrates are transferred through the involvement of scaffolding proteins, and outline the functional merits to metabolic or signaling pathways. Overall, while co-localization topography helps multistep pathways non-specifically, allosteric regulation requires precise multienzyme organization and interactions and works via population shift, either through direct enzyme subunit-subunit interactions or through active involvement of scaffolding proteins. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Structured crowding and its effects on enzyme catalysis

Macromolecular crowding decreases the diffusion rate, shifts the equilibrium of protein-protein and protein-substrate interactions, and changes protein conformational dynamics. Collectively, these effects contribute to enzyme catalysis. Here we describe how crowding may bias the conformational change and dynamics of enzyme populations and in this way affect catalysis. Crowding effects have been studied using artificial crowding agents and in vivo-like environments. These studies revealed a correlation between protein dynamics and function in the crowded environment. We suggest that crowded environments be classified into uniform crowding and structured crowding. Uniform crowding represents random crowding conditions created by synthetic particles with a narrow size distribution. Structured crowding refers to the highly coordinated cellular environment, where proteins and other macromolecules are clustered and organized. In structured crowded environments the perturbation of protein thermal stability may be lower; however, it may still be able to modulate functions effectively and dynamically. Dynamic, allosteric enzymes could be more sensitive to cellular perturbations if their free energy landscape is flatter around the native state; on the other hand, if their free energy landscape is rougher, with high kinetic barriers separating deep minima, they could be more robust. Above all, cells are structured; and this holds both for the cytosol and for the membrane environment. The crowded environment is organized, which limits the search, and the crowders are not necessarily inert. More likely, they too transmit allosteric effects, and as such play important functional roles. Overall, structured cellular crowding may lead to higher enzyme efficiency and specificity.

Using distal-site mutations and allosteric inhibition to tune, extend, and narrow the useful dynamic range of aptamer-based sensors

Here we demonstrate multiple, complementary approaches by which to tune, extend, or narrow the dynamic range of aptamer-based sensors. Specifically, we employ both distal-site mutations and allosteric control to tune the affinity and dynamic range of a fluorescent aptamer beacon. We show that allosteric control, achieved by using a set of easily designed oligonucleotide inhibitors that competes against the folding of the aptamer, allows rational fine-tuning of the affinity of our model aptamer across 3 orders of magnitude of target concentration with greater precision than that achieved using mutational approaches. Using these methods, we generate sets of aptamers varying significantly in target affinity and then combine them to recreate several of the mechanisms employed by nature to narrow or broaden the dynamic range of biological receptors. Such ability to finely control the affinity and dynamic range of aptamers may find many applications in synthetic biology, drug delivery, and targeted therapies, fields in which aptamers are of rapidly growing importance.

Impact of mutations on the allosteric conformational equilibrium

Allostery in a protein involves effector binding at an allosteric site that changes the structure and/or dynamics at a distant, functional site. In addition to the chemical equilibrium of ligand binding, allostery involves a conformational equilibrium between one protein substate that binds the effector and a second substate that less strongly binds the effector. We run molecular dynamics simulations using simple, smooth energy landscapes to sample specific ligand-induced conformational transitions, as defined by the effector-bound and effector-unbound protein structures. These simulations can be performed using our web server (http://salilab.org/allosmod/). We then develop a set of features to analyze the simulations and capture the relevant thermodynamic properties of the allosteric conformational equilibrium. These features are based on molecular mechanics energy functions, stereochemical effects, and structural/dynamic coupling between sites. Using a machine-learning algorithm on a data set of 10 proteins and 179 mutations, we predict both the magnitude and the sign of the allosteric conformational equilibrium shift by the mutation; the impact of a large identifiable fraction of the mutations can be predicted with an average unsigned error of 1k(B)T. With similar accuracy, we predict the mutation effects for an 11th protein that was omitted from the initial training and testing of the machine-learning algorithm. We also assess which calculated thermodynamic properties contribute most to the accuracy of the prediction.

Allosteric regulation of DegS protease subunits through a shared energy landscape

The PDZ domains of the trimeric DegS protease bind unassembled outer-membrane proteins (OMPs) that accumulate in the E. coli periplasm. This cooperative binding reaction triggers a proteolytic cascade that activates a transcriptional stress response. To dissect the mechanism of allosteric activation, we generated hybrid DegS trimers with different numbers of PDZ domains and/or protease-domain mutations. By studying the chemical reactivity and enzymatic properties of these hybrids, we show that all subunits experience a strongly coupled energetic landscape. For example, OMP-peptide binding to a single PDZ domain stimulates active-site chemical modification and proteolytic cleavage in the attached and neighboring protease domains. OMP-peptide binding relieves inhibitory PDZ interactions, whereas the interfaces between protease domains in the trimeric DegS core mediate positively cooperative activation driven both by substrate binding and inhibition relief.

Comparing proteins by their internal dynamics: exploring structure-function relationships beyond static structural alignments

The growing interest for comparing protein internal dynamics owes much to the realisation that protein function can be accompanied or assisted by structural fluctuations and conformational changes. Analogously to the case of functional structural elements, those aspects of protein flexibility and dynamics that are functionally oriented should be subject to evolutionary conservation. Accordingly, dynamics-based protein comparisons or alignments could be used to detect protein relationships that are more elusive to sequence and structural alignments. Here we provide an account of the progress that has been made in recent years towards developing and applying general methods for comparing proteins in terms of their internal dynamics and advance the understanding of the structure-function relationship.

Exploiting protein flexibility to predict the location of allosteric sites

BACKGROUND

Allostery is one of the most powerful and common ways of regulation of protein activity. However, for most allosteric proteins identified to date the mechanistic details of allosteric modulation are not yet well understood. Uncovering common mechanistic patterns underlying allostery would allow not only a better academic understanding of the phenomena, but it would also streamline the design of novel therapeutic solutions. This relatively unexplored therapeutic potential and the putative advantages of allosteric drugs over classical active-site inhibitors fuel the attention allosteric-drug research is receiving at present. A first step to harness the regulatory potential and versatility of allosteric sites, in the context of drug-discovery and design, would be to detect or predict their presence and location. In this article, we describe a simple computational approach, based on the effect allosteric ligands exert on protein flexibility upon binding, to predict the existence and position of allosteric sites on a given protein structure.

RESULTS

By querying the literature and a recently available database of allosteric sites, we gathered 213 allosteric proteins with structural information that we further filtered into a non-redundant set of 91 proteins. We performed normal-mode analysis and observed significant changes in protein flexibility upon allosteric-ligand binding in 70% of the cases. These results agree with the current view that allosteric mechanisms are in many cases governed by changes in protein dynamics caused by ligand binding. Furthermore, we implemented an approach that achieves 65% positive predictive value in identifying allosteric sites within the set of predicted cavities of a protein (stricter parameters set, 0.22 sensitivity), by combining the current analysis on dynamics with previous results on structural conservation of allosteric sites. We also analyzed four biological examples in detail, revealing that this simple coarse-grained methodology is able to capture the effects triggered by allosteric ligands already described in the literature.

CONCLUSIONS

We introduce a simple computational approach to predict the presence and position of allosteric sites in a protein based on the analysis of changes in protein normal modes upon the binding of a coarse-grained ligand at predicted cavities. Its performance has been demonstrated using a newly curated non-redundant set of 91 proteins with reported allosteric properties. The software developed in this work is available upon request from the authors.

Allosteric communication across the native and mutated KIT receptor tyrosine kinase

A fundamental goal in cellular signaling is to understand allosteric communication, the process by which signals originated at one site in a protein propagate dependably to affect remote functional sites. Here, we describe the allosteric regulation of the receptor tyrosine kinase KIT. Our analysis evidenced that communication routes established between the activation loop (A-loop) and the distant juxtamembrane region (JMR) in the native protein were disrupted by the oncogenic mutation D816V positioned in the A-loop. In silico mutagenesis provided a plausible way of restoring the protein communication detected in the native KIT by introducing a counter-balancing second mutation D792E. The communication patterns observed in the native and mutated KIT correlate perfectly with the structural and dynamical features of these proteins. Particularly, a long-distance effect of the D816V mutation manifested as an important structural re-organization of the JMR in the oncogenic mutant was completely vanished in the double mutant D816V/D792E. This detailed characterization of the allosteric communication in the different forms of KIT, native and mutants, was performed by using a modular network representation composed of communication pathways and independent dynamic segments. Such representation permits to enrich a purely mechanistic interaction-based model of protein communication by the introduction of concerted local atomic fluctuations. This method, validated on KIT receptor, may guide a rational modulation of the physiopathological activities of other receptor tyrosine kinases.

The majority of functionally important biological processes are regulated by allosteric communication within individual proteins and across protein complexes. Receptor tyrosine kinases (RTKs) control signal transduction pathways and consequently represent a typical paradigm. The mutation-induced deregulation of RTK activity impairs crucial cellular physiological functions and causes serious human diseases. The present study focuses on the allosteric communication across the three-dimensional structure of the RTK KIT cytoplasmic region. Combining a mechanistic model of information transmission with the analysis of concerted local atomic fluctuations we examined and compared the communication profiles in the native and D816V-mutated proteins. This approach permitted to localize and visualize communication routes in the native KIT and revealed that these routes were disrupted in the mutant D816V. We proposed in silico mutagenesis as a mean to restore the communication detected in the native KIT. Our work sheds light on the allosteric communication in RTKs, a phenomenon playing an essential role in signaling pathways albeit experiments do not provide the atomic details of the path followed in going from one structural element to the other. A rational understanding of the molecular determinants underlying the effects of disease-related kinase mutations may contribute to the improvement of targeted therapies.

Function of hyperekplexia-causing α1R271Q/L glycine receptors is restored by shifting the affected residue out of the allosteric signalling pathway

BACKGROUND AND PURPOSE Glycine receptor α1 subunit R271Q and R271L (α1R271Q/L) mutations cause the neuromotor disorder, hereditary hyperekplexia. Studies suggest that the 271 residue is located within the allosteric signalling pathway linking the agonist binding site to the channel gate. The present study aimed to investigate a possible mechanism for restoring the function of the α1R271Q/L glycine receptor. EXPERIMENTAL APPROACH A 12-amino-acid segment incorporating the 271 residue on the glycine receptor α1271Q/L subunit was replaced by the homologous segment from the glycine receptor β subunit (α1(Ch) 271Q/L). The function of the α1(Ch) 271Q/L glycine receptor was examined by whole-cell patch-clamp recording and voltage-clamp fluorometry techniques. KEY RESULTS The function of the α1(Ch) 271Q/L glycine receptor was restored to the level of the wild-type (WT) α1 glycine receptor. Moreover, in the α1(Ch) glycine receptor, in contrast to the α1 glycine receptor, the channel function was not sensitive to various substitutions of the 271 residue, and the conformational change in the vicinity of the 271 residue was uncoupled from the channel gating. CONCLUSIONS AND IMPLICATIONS The 271 residue is shifted out of the allosteric signalling pathway in the α1(Ch) glycine receptor. We propose that this mechanism provides a novel drug design strategy not only for glycine receptor α1R271Q/L-caused hereditary hyperekplexia, but also for any pathological condition that is caused by missense mutation- or covalent modification-induced disorders involving residues in allosteric signalling pathways. Such a strategy makes it possible to design an ideal drug, which only corrects the function of the mutant or modified protein without affecting the WT or naive protein.

A dynamic model of long-range conformational adaptations triggered by nucleotide binding in GroEL-GroES

The molecular chaperone, GroEL, essential for correct protein folding in E. coli, is composed of 14 identical subunits organized in two interacting rings, each providing a folding chamber for non-native substrate proteins. The oligomeric assembly shows positive cooperativity within each ring and negative cooperativity between the rings. Although it is well known that ATP and long-range allosteric interactions drive the functional cycle of GroEL, an atomic resolution view of how ligand binding modulates conformational adaptations over long distances remains a major challenge. Moreover, little is known on the relation between equilibrium dynamics at physiological temperatures and the allosteric transitions in GroEL. Here we present multiple all-atom molecular dynamics simulations of the GroEL-GroES assemblies at different stages of the functional cycle. Combined with an extensive analysis of the complete set of experimentally available structures, principal component analysis and conformer plots, we provide an explicit evaluation of the accessible conformational space of unliganded GroEL. Our results suggest the presence of pre-existing conformers at the equatorial domain level, and a shift of the conformational ensemble upon ATP-binding. At the inter-ring interface the simulations capture a remarkable offset motion of helix D triggered by ATP-binding to the folding active ring. The reorientation of helix D, previously only observed upon GroES association, correlates with a change of the internal dynamics in the opposite ring. This work contributes to the understanding of the molecular mechanisms in GroEL and highlights the ability of all-atom MD simulations to model long-range structural changes and allosteric events in large systems.

The contribution of RCK domains to human BK channel allosteric activation

Large conductance voltage- and Ca(2+)-activated K(+) (BK) channels are potent regulators of cellular processes including neuronal firing, synaptic transmission, cochlear hair cell tuning, insulin release, and smooth muscle tone. Their unique activation pathway relies on structurally distinct regulatory domains including one transmembrane voltage-sensing domain (VSD) and two intracellular high affinity Ca(2+)-sensing sites per subunit (located in the RCK1 and RCK2 domains). Four pairs of RCK1 and RCK2 domains form a Ca(2+)-sensing apparatus known as the "gating ring." The allosteric interplay between voltage- and Ca(2+)-sensing apparati is a fundamental mechanism of BK channel function. Using voltage-clamp fluorometry and UV photolysis of intracellular caged Ca(2+), we optically resolved VSD activation prompted by Ca(2+) binding to the gating ring. The sudden increase of intracellular Ca(2+) concentration ([Ca(2+)](i)) induced a hyperpolarizing shift in the voltage dependence of both channel opening and VSD activation, reported by a fluorophore labeling position 202, located in the upper side of the S4 transmembrane segment. The neutralization of the Ca(2+) sensor located in the RCK2 domain abolished the effect of [Ca(2+)](i) increase on the VSD rearrangements. On the other hand, the mutation of RCK1 residues involved in Ca(2+) sensing did not prevent the effect of Ca(2+) release on the VSD, revealing a functionally distinct interaction between RCK1 and RCK2 and the VSD. A statistical-mechanical model quantifies the complex thermodynamics interplay between Ca(2+) association in two distinct sites, voltage sensor activation, and BK channel opening.

Structural and Entropic Allosteric Signal Transduction Strength via Correlated Motions

Allosteric signal transduction in biomacromolecules can play an essential role in their function. Internal motional correlations in proteins provide a possible communication mechanism, but the quantitative relationship between statistical correlations and allostery is unknown. Quantitative relationships between internal motional correlations and the efficiency of propagation of allosteric structural and entropic effects are introduced and validated against conformational ensembles obtained from molecular dynamics simulations. This framework can explain a range of phenomena, such as the occurrence of an allosteric entropy change in the absence of any noticeable structural change.

On the effect of protein conformation diversity in discriminating among neutral and disease related single amino acid substitutions

BACKGROUND

Non-synonymous coding SNPs (nsSNPs) that are associated to disease can also be related with alterations in protein stability. Computational methods are available to predict the effect of single amino acid substitutions (SASs) on protein stability based on a single folded structure. However, the native state of a protein is not unique and it is better represented by the ensemble of its conformers in dynamic equilibrium. The maintenance of the ensemble is essential for protein function. In this work we investigated how protein conformational diversity can affect the discrimination of neutral and disease related SASs based on protein stability estimations. For this purpose, we used 119 proteins with 803 associated SASs, 60% of which are disease related. Each protein was associated with its corresponding set of available conformers as found in the Protein Conformational Database (PCDB). Our dataset contains proteins with different extensions of conformational diversity summing up a total number of 1023 conformers.

RESULTS

The existence of different conformers for a given protein introduces great variability in the estimation of the protein stability (ΔΔG) after a single amino acid substitution (SAS) as computed with FoldX. Indeed, in 35% of our protein set at least one SAS can be described as stabilizing, destabilizing or neutral when a cutoff value of ±2 kcal/mol is adopted for discriminating neutral from perturbing SASs. However, when the ΔΔG variability among conformers is taken into account, the correlation among the perturbation of protein stability and the corresponding disease or neutral phenotype increases as compared with the same analysis on single protein structures. At the conformer level, we also found that the different conformers correlate in a different way to the corresponding phenotype.

CONCLUSIONS

Our results suggest that the consideration of conformational diversity can improve the discrimination of neutral and disease related protein SASs based on the evaluation of the corresponding Gibbs free energy change.

Allosterically gated enzyme dynamics in the cysteine synthase complex regulate cysteine biosynthesis in Arabidopsis thaliana

Plants and bacteria assimilate sulfur into cysteine. Cysteine biosynthesis involves a bienzyme complex, the cysteine synthase complex (CSC), which consists of serine-acetyl-transferase (SAT) and O-acetyl-serine-(thiol)-lyase (OAS-TL) enzymes. The activity of OAS-TL is reduced by formation of the CSC. Although this reduction is an inherent part of the self-regulation cycle of cysteine biosynthesis, there has until now been no explanation as to how OAS-TL loses activity in plants. Complexation of SAT and OAS-TL involves binding of the C-terminal tail of SAT in one of the active sites of the homodimeric OAS-TL. We here explore the flexibility of the unoccupied active site in Arabidopsis thaliana cytosolic and mitochondrial OAS-TLs. Our results reveal two gates in the OAS-TL active site that define its accessibility. The observed dynamics of the gates show allosteric closure of the unoccupied active site of OAS-TL in the CSC, which can hinder substrate binding, abolishing its turnover to cysteine.

Identification of the bile salt binding site on IpaD from Shigella flexneri and the influence of ligand binding on IpaD structure

Type III secretion (TTS) is an essential virulence factor for Shigella flexneri, the causative agent of shigellosis. The Shigella TTS apparatus (TTSA) is an elegant nanomachine that is composed of a basal body, an external needle to deliver effectors into human cells, and a needle tip complex that controls secretion activation. IpaD is at the tip of the nascent TTSA needle where it controls the first step of TTS activation. The bile salt deoxycholate (DOC) binds to IpaD to induce recruitment of the translocator protein IpaB into the maturing tip complex. We recently used spectroscopic analyses to show that IpaD undergoes a structural rearrangement that accompanies binding to DOC. Here, we report a crystal structure of IpaD with DOC bound and test the importance of the residues that make up the DOC binding pocket on IpaD function. IpaD binds DOC at the interface between helices α3 and α7, with concomitant movement in the orientation of helix α7 relative to its position in unbound IpaD. When the IpaD residues involved in DOC binding are mutated, some are found to lead to altered invasion and secretion phenotypes. These findings suggest that adoption of a DOC bound structural state for IpaD primes the Shigella TTSA for contact with host cells. The data presented here and in the studies leading up to this work provide the foundation for developing a model of the first step in Shigella TTS activation.

Probing molecular mechanisms of the Hsp90 chaperone: biophysical modeling identifies key regulators of functional dynamics

Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based "conformational selection" of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients.

Insights into the fold organization of TIM barrel from interaction energy based structure networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or “sequence conservation” as the basis for their understanding. Recently “interaction energy” based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the “interaction conservation” viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Proteins are polymers of amino-acids that fold into unique three-dimensional structures to perform cellular functions. This structure formation has been shown to depend on the amino-acid sequences. But examples of proteins with diverse sequences retaining a similar structural fold are quite substantial that we can no longer consider such phenomenon as exceptions. Therefore, this non-canonical relationship has been studied extensively mostly by studying the remote sequence similarities between proteins. Here we have attempted to address the above-mentioned problem by analyzing the similarities in the spatial interactions among amino-acids. Since the protein structure is a resultant of different interactions, we have considered the proteins as networks of interacting amino-acids to derive the common interactions within a popular structural fold called the TIM barrel fold. We were able to find common interactions among different families of the TIM fold and generalize the patterns of interactions by which the fold is being maintained despite sequence diversity. The results substantiate our hypothesis that interaction conservation might by a driving factor in fold formation and this new outlook can be used extensively in engineering proteins with better biophysical characteristics.

Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins

BACKGROUND

Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question "What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?"

RESULTS

A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex.Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors.

CONCLUSIONS

The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.

Allosteric modulators can restore function in an amino acid neurotransmitter receptor by slightly altering intra-molecular communication pathways

Mutations, even if not directly in the ligand binding sites of proteins, can lead to disease. In cell surface receptors, this can happen if they uncouple conformational changes that take place upon agonist (or antagonist) binding to the extracellular domain and the intracellular response. Uncoupling can take place by disrupting a major allosteric propagation pathway between the extra- and intracellular domains. Here I provide a mechanistic explanation: I first describe how propagation takes place; second, what can happen in the presence of a disease-related mutation which is distant from the binding site; and finally, how drugs may overcome this disruption and rescue function. The mutations in the glycine receptor α1 subunit (α1R271Q/L) which cause the neuromotor disorder hyperekplexia are on example of such allosteric mutations. In this issue of the BJP, Shan et al. show that normal function was restored to these mutant receptors by substitution of the segment which contained the mutated position, by a homologous one. An allosteric drug could mimic the effects of such substitution. Within this framework, I highlight the advantages of allosteric drugs and the challenges in their design.

Conformational control of the binding of the transactivation domain of the MLL protein and c-Myb to the KIX domain of CREB

The KIX domain of CBP is a transcriptional coactivator. Concomitant binding to the activation domain of proto-oncogene protein c-Myb and the transactivation domain of the trithorax group protein mixed lineage leukemia (MLL) transcription factor lead to the biologically active ternary MLL∶KIX∶c-Myb complex which plays a role in Pol II-mediated transcription. The binding of the activation domain of MLL to KIX enhances c-Myb binding. Here we carried out molecular dynamics (MD) simulations for the MLL∶KIX∶c-Myb ternary complex, its binary components and KIX with the goal of providing a mechanistic explanation for the experimental observations. The dynamic behavior revealed that the MLL binding site is allosterically coupled to the c-Myb binding site. MLL binding redistributes the conformational ensemble of KIX, leading to higher populations of states which favor c-Myb binding. The key element in the allosteric communication pathways is the KIX loop, which acts as a control mechanism to enhance subsequent binding events. We tested this conclusion by in silico mutations of loop residues in the KIX∶MLL complex and by comparing wild type and mutant dynamics through MD simulations. The loop assumed MLL binding conformation similar to that observed in the KIX∶c-Myb state which disfavors the allosteric network. The coupling with c-Myb binding site faded, abolishing the positive cooperativity observed in the presence of MLL. Our major conclusion is that by eliciting a loop-mediated allosteric switch between the different states following the binding events, transcriptional activation can be regulated. The KIX system presents an example how nature makes use of conformational control in higher level regulation of transcriptional activity and thus cellular events.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.

The role of conserved waters in conformational transitions of Q61H K-ras

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 µs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the C_α atoms of G60 at switch 2 and G10 at the P-loop and the N-C_α-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.

K-ras belongs to the Ras family of G-proteins that regulate cell proliferation and development. To execute its function, K-ras adopts different conformational states when it is active and inactive. In addition to these two states, it samples many transient intermediate conformations as it makes the transition from one state to the other. Mutations that affect the population of these states can cause cancer or developmental disorder. Using simulation approaches, here we show that a number of water molecules buried within the structure of an oncogenic K-ras protein modulate the distribution of its conformational states. Moreover, a detailed analysis based on two novel structural parameters revealed the existence of long-range water-mediated interactions that facilitate a dynamic coupling between the two lobes of the protein. These findings pave the way for a dynamics-guided strategy to inhibit abnormal Ras signaling.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.

Protein dynamics and conformational selection in bidirectional signal transduction

Protein conformational dynamics simultaneously allow promiscuity and specificity in binding. The multiple conformations of the free EphA4 ligand-binding domain observed in two new EphA4 crystal structures provide a unique insight into the conformational dynamics of EphA4 and its signaling pathways. The heterogeneous ensemble and loop dynamics explain how the EphA4 receptor is able to bind multiple A- and B-ephrin ligands and small molecules via conformational selection, which helps to fine-tune cellular signal response in both receptor and ligand cells.

See research article http://www.biomedcentral.com/2046-1682/5/2

A critical evaluation of correlated mutation algorithms and coevolution within allosteric mechanisms

The notion of using the evolutionary history encoded within multiple sequence alignments to predict allosteric mechanisms is appealing. In this approach, correlated mutations are expected to reflect coordinated changes that maintain intramolecular coupling between residue pairs. Despite much early fanfare, the general suitability of correlated mutations to predict allosteric couplings has not yet been established. Lack of progress along these lines has been hindered by several algorithmic limitations including phylogenetic artifacts within alignments masking true covariance and the computational intractability of consideration of more than two correlated residues at a time. Recent progress in algorithm development, however, has been substantial with a new generation of correlated mutation algorithms that have made fundamental progress toward solving these difficult problems. Despite these encouraging results, there remains little evidence to suggest that the evolutionary constraints acting on allosteric couplings are sufficient to be recovered from multiple sequence alignments. In this review, we argue that due to the exquisite sensitivity of protein dynamics, and hence that of allosteric mechanisms, the latter vary widely within protein families. If it turns out to be generally true that even very similar homologs display a wide divergence of allosteric mechanisms, then even a perfect correlated mutation algorithm could not be reliably used as a general mechanism for discovery of allosteric pathways.

How do dynamic cellular signals travel long distances?

Communication is essential. It is vital between cells in multi-cellular organisms, and within cells. A signaling molecule binds to a receptor protein, and initiates a cascade of dynamic events. Signaling is a multistep pathway, which allows signal amplification: if some of the molecules in a pathway transmit the signal to multiple molecules, the result can be a large number of activated molecules across the cell and multiple reactions. That is how a small number of extracellular signaling molecules can produce a major cellular response. The pathway can relay signals from the extracellular space to the nucleus. How do signals travel efficiently over long-distances across the cell? Here we argue that evolution has utilized three properties: a modular functional organization of the cellular network; sequences in some key regions of proteins, such as linkers or loops, which were pre-encoded by evolution to facilitate signaling among domains; and compact interactions between proteins which is achieved via conformational disorder.