Coupling of global and local vibrational modes in dynamic allostery of proteins

Interaction of Clostridium perfringens Epsilon Toxin with the Plasma Membrane: The Role of Amino Acids Y42, Y43 and H162

Clostridium perfringens epsilon toxin (Etx) is a pore forming toxin that causes enterotoxaemia in ruminants and may be a cause of multiple sclerosis in humans. To date, most in vitro studies of Etx have used the Madin-Darby canine kidney (MDCK) cell line. However, studies using Chinese hamster ovary (CHO) cells engineered to express the putative Etx receptor, myelin and lymphocyte protein (MAL), suggest that amino acids important for Etx activity differ between species. In this study, we investigated the role of amino acids Y42, Y43 and H162, previously identified as important in Etx activity towards MDCK cells, in Etx activity towards CHO-human MAL (CHO-hMAL) cells, human red blood cells (hRBCs) and synthetic bilayers using site-directed mutants of Etx. We show that in CHO-hMAL cells Y42 is critical for Etx binding and not Y43 as in MDCK cells, indicating that surface exposed tyrosine residues in the receptor binding domain of Etx impact efficiency of cell binding to MAL-expressing cells in a species-specific manner. We also show that Etx mutant H162A was unable to lyse CHO-hMAL cells, lysed hRBCs, whilst it was able to form pores in synthetic bilayers, providing evidence of the complexity of Etx pore formation in different lipid environments.

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combining steady state and temperature jump IR spectroscopy to investigate the allosteric effects of ligand binding to dsDNA

Changes in the structural dynamics of double stranded (ds)DNA upon ligand binding have been linked to the mechanism of allostery without conformational change, but direct experimental evidence remains elusive. To address this, a combination of steady state infrared (IR) absorption spectroscopy and ultrafast temperature jump IR absorption measurements has been used to quantify the extent of fast (∼100 ns) fluctuations in (ds)DNA·Hoechst 33258 complexes at a range of temperatures. Exploiting the direct link between vibrational band intensities and base stacking shows that the absolute magnitude of the change in absorbance caused by fast structural fluctuations following the temperature jump is only weakly dependent on the starting temperature of the sample. The observed fast dynamics are some two orders of magnitude faster than strand separation and associated with all points along the 10-base pair duplex d(GCATATATCC). Binding the Hoechst 33258 ligand causes a small but consistent reduction in the extent of these fast fluctuations of base pairs located outside of the ligand binding region. These observations point to a ligand-induced reduction in the flexibility of the dsDNA near the binding site, consistent with an estimated allosteric propagation length of 15 Å, about 5 base pairs, which agrees well with both molecular simulation and coarse-grained statistical mechanics models of allostery leading to cooperative ligand binding.

Computational analysis of dynamic allostery and control in the SARS-CoV-2 main protease

The COVID-19 pandemic caused by the novel coronavirus SARS-CoV-2 has no publicly available vaccine or antiviral drugs at the time of writing. An attractive coronavirus drug target is the main protease (Mpro, also known as 3CLpro) because of its vital role in the viral cycle. A significant body of work has been focused on finding inhibitors which bind and block the active site of the main protease, but little has been done to address potential non-competitive inhibition, targeting regions other than the active site, partly because the fundamental biophysics of such allosteric control is still poorly understood. In this work, we construct an elastic network model (ENM) of the SARS-CoV-2 Mpro homodimer protein and analyse its dynamics and thermodynamics. We found a rich and heterogeneous dynamical structure, including allosterically correlated motions between the homodimeric protease's active sites. Exhaustive 1-point and 2-point mutation scans of the ENM and their effect on fluctuation free energies confirm previously experimentally identified bioactive residues, but also suggest several new candidate regions that are distant from the active site, yet control the protease function. Our results suggest new dynamically driven control regions as possible candidates for non-competitive inhibiting binding sites in the protease, which may assist the development of current fragment-based binding screens. The results also provide new insights into the active biophysical research field of protein fluctuation allostery and its underpinning dynamical structure.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

Transient association between proteins elicits alteration of dynamics at sites far away from interfaces

Proteins are known to undergo structural changes upon binding to partner proteins. However, the prevalence, extent, location, and function of change in protein dynamics due to transient protein-protein interactions is not well documented. Here, we have analyzed a dataset of 58 protein-protein complexes of known three-dimensional structure and structures of their corresponding unbound forms to evaluate dynamics changes induced by binding. Fifty-five percent of cases showed significant dynamics change away from the interfaces. This change is not always accompanied by an observed structural change. Binding of protein partner is found to alter inter-residue communication within the tertiary structure in about 90% of cases. Also, residue motions accessible to proteins in unbound form were not always maintained in the bound form. Further analyses revealed functional roles for the distant site where dynamics change was observed. Overall, the results presented here strongly suggest that alteration of protein dynamics due to binding of a partner protein commonly occurs.

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.

How proteins' negative cooperativity emerges from entropic optimisation of versatile collective fluctuations

The fact that allostery, a nonlocal signaling between distant binding sites, can arise mainly from the entropy balance of collective thermal modes, without conformational changes, is by now well known. However, the propensity to generate negative cooperativity is still unclear. Starting from an elastic-network picture of small protein complexes, in which effector binding is modeled by locally altering interaction strengths in lieu of adding a node-spring pair, we elucidate mechanisms particularly for such negative cooperativity. The approach via a few coupled harmonic oscillators with internal elastic strengths allows us to trace individual eigenmodes, their frequencies, and their statistical weights through successive bindings. We find that the alteration of the oscillators' couplings is paramount to covering both signs of allostery. Binding-modified couplings create a rich set of eigenmodes individually for each binding state, modes inaccessible to an ensemble of noninteracting units. The associated shifts of collective-mode frequencies, nonuniform with respect to modes and binding states, result in an enhanced optimizability, reflected by a subtle phase map of allosteric behaviors.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

A Biophysical Perspective on Enzyme Catalysis

Even after a century of investigation, our understanding of how enzymes work remains far from complete. In particular, several factors that enable enzymes to achieve high catalytic efficiencies remain only poorly understood. A number of theories have been developed, which propose or reaffirm that enzymes work as structural scaffolds, serving to bring together and properly orient the participants so that the reaction can proceed; therefore, leading to enzymes being viewed as only passive participants in the catalyzed reaction. A growing body of evidence shows that enzymes are not rigid structures but are constantly undergoing a wide range of internal motions and conformational fluctuations. In this Perspective, on the basis of studies from our group, we discuss the emerging biophysical model of enzyme catalysis that provides a detailed understanding of the interconnection among internal protein motions, conformational substates, enzyme mechanisms, and the catalytic efficiency of enzymes. For a number of enzymes, networks of conserved residues that extend from the surface of the enzyme all the way to the active site have been discovered. These networks are hypothesized to serve as pathways of energy transfer that enables thermodynamical coupling of the surrounding solvent with enzyme catalysis and play a role in promoting enzyme function. Additionally, the role of enzyme structure and electrostatic effects has been well acknowledged for quite some time. Collectively, the recent knowledge gained about enzyme mechanisms suggests that the conventional paradigm of enzyme structure encoding function is incomplete and needs to be extended to structure encodes dynamics, and together these enzyme features encode function including catalytic rate acceleration.

Ligand-regulated oligomerisation of allosterically interacting proteins

The binding of ligands to distinct sites at proteins or at protein clusters is often cooperative or anti-cooperative due to allosteric signalling between those sites. The allostery is usually attributed to a configurational change of the proteins from a relaxed to a configurationally different tense state. Alternatively, as originally proposed by Cooper and Dryden, a tense state may be achieved by merely restricting the thermal vibrations of the protein around its mean configuration. In this work, we provide theoretical tools to investigate fluctuation allostery using cooling and titration experiments in which ligands regulate dimerisation, or ring or chain formation. We discuss in detail how ligands may regulate the supramolecular (co)polymerisation of liganded and unliganded proteins.

Site-specific ion occupation in the selectivity filter causes voltage-dependent gating in a viral K+ channel

Many potassium channels show voltage-dependent gating without a dedicated voltage sensor domain. This is not fully understood yet, but often explained by voltage-induced changes of ion occupation in the five distinct K+ binding sites in the selectivity filter. To better understand this mechanism of filter gating we measured the single-channel current and the rate constant of sub-millisecond channel closure of the viral K+ channel KcvNTS for a wide range of voltages and symmetric and asymmetric K+ concentrations in planar lipid membranes. A model-based analysis employed a global fit of all experimental data, i.e., using a common set of parameters for current and channel closure under all conditions. Three different established models of ion permeation and various relationships between ion occupation and gating were tested. Only one of the models described the data adequately. It revealed that the most extracellular binding site (S0) in the selectivity filter functions as the voltage sensor for the rate constant of channel closure. The ion occupation outside of S0 modulates its dependence on K+ concentration. The analysis uncovers an important role of changes in protein flexibility in mediating the effect from the sensor to the gate.

The 'allosteron' model for entropic allostery of self-assembly

Using the simple 'allosteron' model, we show that it is possible, in principle, to elicit pathways by which fluctuation allostery affects self-assembly of protein complexes. We treat the cases of (i) protein fibrils and nucleation, (ii) n-mer protein complexes, and (iii) weakly attractive allosteric interactions in protein-like soft nanoscale objects that can be tuned to define exclusive self-associating families.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Evolutionary conservation of the allosteric activation of factor VIIa by tissue factor in lamprey

Essentials Tissue factor (TF) enhances factor VIIa (FVIIa) activity through structural and dynamic changes. We analyzed conservation of TF-activated FVIIa allosteric networks in extant vertebrate lamprey. Lamprey Tf/FVIIa molecular dynamics show conserved Tf-induced structural/dynamic FVIIa changes. Lamprey Tf activation of FVIIa allosteric networks follows molecular pathways similar to human.

SUMMARY

Background Previous studies have provided insight into the molecular basis of human tissue factor (TF) activation of activated factor VII (FVIIa). TF-induced allosteric networks of FVIIa activation have been rationalized through analysis of the dynamic changes and residue connectivities in the human soluble TF (sTF)/FVIIa complex structure during molecular dynamics (MD) simulation. Evolutionary conservation of the molecular mechanisms for TF-induced allosteric FVIIa activation between humans and extant vertebrate jawless fish (lampreys), where blood coagulation emerged more than 500 million years ago, is unknown and of considerable interest. Objective To model the sTf/FVIIa complex from cloned Petromyzon marinus lamprey sequences, and with comparisons to human sTF/FVlla investigate conservation of allosteric mechanisms of FVIIa activity enhancement by soluble TF using MD simulations. Methods Full-length cDNAs of lamprey tf and f7 were cloned and characterized. Comparative models of lamprey sTf/FVIIa complex and free FVIIa were determined based on constructed human sTF/FVIIa complex and free FVIIa models, used in full-atomic MD simulations, and characterized using dynamic network analysis approaches. Results Allosteric paths of correlated motion from Tf contact points in lamprey sTf/FVIIa to the FVIIa active site were determined and quantified, and were found to encompass residue-residue interactions along significantly similar paths compared with human. Conclusions Despite low conservation of residues between lamprey and human proteins, 30% TF and 39% FVII, the structural and protein dynamic effects of TF activation of FVIIa appear conserved and, moreover, present in extant vertebrate proteins from 500 million years ago when TF/FVIIa-initiated extrinsic pathway blood coagulation emerged.

Chlorophyll-Derivative Modulation of Rhodopsin Signaling Properties through Evolutionarily Conserved Interaction Pathways

Retinal is the light-absorbing chromophore that is responsible for the activation of visual pigments and light-driven ion pumps. Evolutionary changes in the intermolecular interactions of the retinal with specific amino acids allow for adaptation of the spectral characteristics, referred to as spectral tuning. However, it has been proposed that a specific species of dragon fish has bypassed the adaptive evolutionary process of spectral tuning and replaced it with a single evolutionary event: photosensitization of rhodopsin by chlorophyll derivatives. Here, by using a combination of experimental measurements and computational modeling to probe retinal-receptor interactions in rhodopsin, we show how the binding of the chlorophyll derivative, chlorin-e6 (Ce6) in the intracellular domain (ICD) of the receptor allosterically excites G-protein coupled receptor class A (GPCR-A) conserved long-range correlated fluctuations that connect distant parts of the receptor. These long-range correlated motions are associated with regulating the dynamics and intermolecular interactions of specific amino acids in the retinal ligand-binding pocket that have been associated with shifts in the absorbance peak maximum (λ_max) and hence, spectral sensitivity of the visual system. Moreover, the binding of Ce6 affects the overall global properties of the receptor. Specifically, we find that Ce6-induced dynamics alter the thermal stability of rhodopsin by adjusting hydrogen-bonding interactions near the receptor active-site that consequently also influences the intrinsic conformational equilibrium of the receptor. Due to the conservation of the ICD residues amongst different receptors in this class and the fact that all GPCR-A receptors share a common mechanism of activation, it is possible that the allosteric associations excited in rhodopsin with Ce6 binding are a common feature in all class A GPCRs.

Vibrational resonance, allostery, and activation in rhodopsin-like G protein-coupled receptors

G protein-coupled receptors are a large family of membrane proteins activated by a variety of structurally diverse ligands making them highly adaptable signaling molecules. Despite recent advances in the structural biology of this protein family, the mechanism by which ligands induce allosteric changes in protein structure and dynamics for its signaling function remains a mystery. Here, we propose the use of terahertz spectroscopy combined with molecular dynamics simulation and protein evolutionary network modeling to address the mechanism of activation by directly probing the concerted fluctuations of retinal ligand and transmembrane helices in rhodopsin. This approach allows us to examine the role of conformational heterogeneity in the selection and stabilization of specific signaling pathways in the photo-activation of the receptor. We demonstrate that ligand-induced shifts in the conformational equilibrium prompt vibrational resonances in the protein structure that link the dynamics of conserved interactions with fluctuations of the active-state ligand. The connection of vibrational modes creates an allosteric association of coupled fluctuations that forms a coherent signaling pathway from the receptor ligand-binding pocket to the G-protein activation region. Our evolutionary analysis of rhodopsin-like GPCRs suggest that specific allosteric sites play a pivotal role in activating structural fluctuations that allosterically modulate functional signals.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

Function changing mutations in glucocorticoid receptor evolution correlate with their relevance to mode coupling

Nonlinear effects in protein dynamics are expected to play role in function, particularly of allosteric nature, by facilitating energy transfer between vibrational modes. A recently proposed method focusing on the non-Gaussian shape of the configurational population near equilibrium projects this information onto real space in order to identify the aminoacids relevant to function. We here apply this method to three ancestral proteins in glucocorticoid receptor (GR) family and show that the mutations that restrict functional activity during GR evolution correlate significantly with locations that are highlighted by the nonlinear contribution to the near-native configurational distribution. Our findings demonstrate that the analysis of nonlinear effects in protein dynamics can be harnessed into a predictive tool for functional site determination.

Dynamic Transmission of Protein Allostery without Structural Change: Spatial Pathways or Global Modes?

We examine the contrast between mechanisms for allosteric signaling that involve structural change, and those that do not, from the perspective of allosteric pathways. In particular we treat in detail the case of fluctuation-allostery by which amplitude modulation of the thermal fluctuations of the elastic normal modes conveys the allosteric signal, and address the question of what an allosteric pathway means in this case. We find that a perturbation theory of thermal elastic solids and nonperturbative approach (by super-coarse-graining elasticity into internal bending modes) have opposite signatures in their structure of correlated pathways. We illustrate the effect from analysis of previous results from GlxR of Corynebacterium glutamicum, an example of the CRP/FNR transcription family of allosteric homodimers. We find that the visibility of both correlated pathways and disconnected sites of correlated motion in this protein suggests that mechanisms of local elastic stretch and bend are recruited for the purpose of creating and controlling allosteric cooperativity.

On the possible role of protein vibrations in information processing in the brain: three Russian dolls

Until recently it was held that the neurocomputations conducted by the brain involved only whole neurons as the operating units. This may however represent only a part of the mechanism. This theoretical and academic position article reviews the considerable evidence that allosteric interactions between proteins (as extensively described by Fuxe et al., 2014), and in particular protein vibrations in neurons, form small scale codes that are involved as parts of the complex information processing systems of the brain. The argument is then developed to suggest that the protein allosteric and vibration codes (that operate at the molecular level) are nested within a medium scale coding system whose computational units are organelles (such as microtubules). This medium scale code is nested in turn inside a large scale coding system, whose computational units are individual neurons. The hypothesis suggests that these three levels interact vertically in both directions thus materially increasing the computational capacity of the brain. The whole hierarchy is thus similar to three nested Russian dolls. This theoretical development may be of use in the design of experiments to test it.

The Role of Protein-Ligand Contacts in Allosteric Regulation of the Escherichia coli Catabolite Activator Protein

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. Both experimental and theoretical evidence demonstrate that allostery can be communicated through altered slow relaxation protein dynamics without conformational change. The catabolite activator protein (CAP) of Escherichia coli is an exemplar for the analysis of such entropically driven allostery. Negative allostery in CAP occurs between identical cAMP binding sites. Changes to the cAMP-binding pocket can therefore impact the allosteric properties of CAP. Here we demonstrate, through a combination of coarse-grained modeling, isothermal calorimetry, and structural analysis, that decreasing the affinity of CAP for cAMP enhances negative cooperativity through an entropic penalty for ligand binding. The use of variant cAMP ligands indicates the data are not explained by structural heterogeneity between protein mutants. We observe computationally that altered interaction strength between CAP and cAMP variously modifies the change in allosteric cooperativity due to second site CAP mutations. As the degree of correlated motion between the cAMP-contacting site and a second site on CAP increases, there is a tendency for computed double mutations at these sites to drive CAP toward noncooperativity. Naturally occurring pairs of covarying residues in CAP do not display this tendency, suggesting a selection pressure to fine tune allostery on changes to the CAP ligand-binding pocket without a drive to a noncooperative state. In general, we hypothesize an evolutionary selection pressure to retain slow relaxation dynamics-induced allostery in proteins in which evolution of the ligand-binding site is occurring.

Global low-frequency motions in protein allostery: CAP as a model system

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. There is considerable evidence that allosteric cooperativity can be communicated by the modulation of protein dynamics without conformational change. The Catabolite Activator Protein (CAP) of Escherichia coli is an important experimental exemplar for entropically driven allostery. Here we discuss recent experimentally supported theoretical analysis that highlights the role of global low-frequency dynamics in allostery in CAP and identify how allostery arises as a natural consequence of changes in global low-frequency protein fluctuations on ligand binding.

Nonlinear excitations match correlated motions unveiled by NMR in proteins: a new perspective on allosteric cross-talk

In this paper we propose a novel theoretical framework for interpreting long-range dynamical correlations unveiled in proteins through NMR measurements. The theoretical rationale relies on the hypothesis that correlated motions in proteins may be reconstructed as large-scale, collective modes sustained by long-lived nonlinear vibrations known as discrete breathers (DB) localized at key, hot-spot sites. DBs are spatially localized modes, whose nonlinear nature hinders resonant coupling with the normal modes, thus conferring them long lifetimes as compared to normal modes. DBs have been predicted to exist in proteins, localized at few hot-spot residues typically within the stiffest portions of the structure. We compute DB modes analytically in the framework of the nonlinear network model, showing that the displacement patterns of many DBs localized at key sites match to a remarkable extent the experimentally uncovered correlation blueprint. The computed dispersion relations prove that it is physically possible for some of these DBs to be excited out of thermal fluctuations at room temperature. Based on our calculations, we speculate that transient energy redistribution among the vibrational modes in a protein might favor the emergence of DB-like bursts of long-lived energy at hot-spot sites with lifetimes in the ns range, able to sustain critical, function-encoding correlated motions. More generally, our calculations provide a novel quantitative tool to predict fold-spanning dynamical pathways of correlated residues that may be central to allosteric cross-talk in proteins.

Epitope fluctuations in the human papillomavirus are under dynamic allosteric control: a computational evaluation of a new vaccine design strategy

The dynamic properties of the capsid of the human papillomavirus (HPV) type 16 were examined using classical molecular dynamics simulations. By systematically comparing the structural fluctuations of the capsid protein, a strong dynamic allosteric connection between the epitope containing loops and the h4 helix located more than 50 Å away is identified, which was not recognized thus far. Computer simulations show that restricting the structural fluctuations of the h4 helix is key to rigidifying the epitopes, which is thought to be required for eliciting a proper immune response. The allostery identified in the components of the HPV is nonclassical because the mean structure of the epitope carrying loops remains unchanged, but as a result of allosteric effect the structural fluctuations are altered significantly, which in turn changes the biochemical reactivity profile of the epitopes. Exploiting this novel insight, a new vaccine design strategy is proposed wherein a relatively small virus capsid fragment is deposited on a silica nanoparticle in such a way that the fluctuations of the h4 helix are suppressed. The structural and dynamic properties of the epitope carrying loops on this hybrid nanoparticle match the characteristics of epitopes found on the full virus-like particle precisely, suggesting that these nanoparticles may serve as potent, cost-effective, and safe alternatives to traditionally developed vaccines. The structural and dynamic properties of the hybrid nanoparticle are examined in detail to establish the general concepts of the proposed new design.

Allostery without conformation change: modelling protein dynamics at multiple scales

The original ideas of Cooper and Dryden, that allosteric signalling can be induced between distant binding sites on proteins without any change in mean structural conformation, has proved to be a remarkably prescient insight into the rich structure of protein dynamics. It represents an alternative to the celebrated Monod-Wyman-Changeux mechanism and proposes that modulation of the amplitude of thermal fluctuations around a mean structure, rather than shifts in the structure itself, give rise to allostery in ligand binding. In a complementary approach to experiments on real proteins, here we take a theoretical route to identify the necessary structural components of this mechanism. By reviewing and extending an approach that moves from very coarse-grained to more detailed models, we show that, a fundamental requirement for a body supporting fluctuation-induced allostery is a strongly inhomogeneous elastic modulus. This requirement is reflected in many real proteins, where a good approximation of the elastic structure maps strongly coherent domains onto rigid blocks connected by more flexible interface regions.

Modulation of global low-frequency motions underlies allosteric regulation: demonstration in CRP/FNR family transcription factors

Allostery in bacterial transcription factors arises from changes in global low-frequency protein dynamics. Amino acids that regulate low-frequency dynamics are identified and seen to be evolutionarily conserved.

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distinct site. There is growing evidence that allosteric cooperativity can be communicated by modulation of protein dynamics without conformational change. The mechanisms, however, for communicating dynamic fluctuations between sites are debated. We provide a foundational theory for how allostery can occur as a function of low-frequency dynamics without a change in structure. We have generated coarse-grained models that describe the protein backbone motions of the CRP/FNR family transcription factors, CAP of Escherichia coli and GlxR of Corynebacterium glutamicum. The latter we demonstrate as a new exemplar for allostery without conformation change. We observe that binding the first molecule of cAMP ligand is correlated with modulation of the global normal modes and negative cooperativity for binding the second cAMP ligand without a change in mean structure. The theory makes key experimental predictions that are tested through an analysis of variant proteins by structural biology and isothermal calorimetry. Quantifying allostery as a free energy landscape revealed a protein “design space” that identified the inter- and intramolecular regulatory parameters that frame CRP/FNR family allostery. Furthermore, through analyzing CAP variants from diverse species, we demonstrate an evolutionary selection pressure to conserve residues crucial for allosteric control. This finding provides a link between the position of CRP/FNR transcription factors within the allosteric free energy landscapes and evolutionary selection pressures. Our study therefore reveals significant features of the mechanistic basis for allostery. Changes in low-frequency dynamics correlate with allosteric effects on ligand binding without the requirement for a defined spatial pathway. In addition to evolving suitable three-dimensional structures, CRP/FNR family transcription factors have been selected to occupy a dynamic space that fine-tunes biological activity and thus establishes the means to engineer allosteric mechanisms driven by low-frequency dynamics.

Allostery is a process by which a molecule binding to one site of a protein alters the activity of the protein at another site. Allostery is typically thought to occur through a change in protein structure, but there is now clear evidence that the dynamic properties of a protein can also regulate allostery without a change in overall conformation. Here we examine two members of a large family of bacterial transcription factors and provide a mechanism to describe the allosteric binding of their activating ligands. We demonstrate, in these systems, that allostery arises as a natural consequence of changes in global low-frequency protein fluctuations on ligand binding. We further demonstrate that the higher dimensional parameter space that describes all potential variant transcription factors can be reduced to a two-dimensional free energy landscape that determines the key molecular parameters that predominantly regulate allostery. We additionally show that the amino acids we determine as contributing sensitively to allosteric control tend to be conserved in diverse bacteria; thus we identify a link between residues that contribute to low-frequency fluctuations and evolutionary selection pressures.

How does hemoglobin generate such diverse functionality of physiological relevance?

The absolute values of the O2-affinities (P50, Klow, and Khigh) of hemoglobin (Hb) are regulated neither by changes in the static T-/R-quaternary and associated tertiary structures nor the ligation states. They are pre-determined and regulated by the extrinsic environmental factors such as pH, buffers, and heterotropic effectors. The effect and role of O2 on Hb are reversibly to drive the structural allosteric equilibrium between the T(deoxy)- and R(oxy)-Hb toward R(oxy)-Hb (the structural allostery). R(oxy)-Hb has a higher O2-affinity (Khigh) relative to that (Klow) of the T(deoxy)-Hb (Khigh>Klow) under any fixed environmental conditions. The apparent O2-affinity of Hb is high, as the globin matrix interferes with the dissociation process of O2, forcing the dissociated O2 geminately to re-bind to the heme Fe. This artificially increases [oxy-Hb] and concomitantly decreases [deoxy-Hb], leading to the apparent increases of the O2-affinity of Hb. The effector-linked high-frequency thermal fluctuations of the globin matrix act as a gating mechanism to modulate such physical, energetic, and kinetic barriers to enhance the dissociation process of O2, resulted in increases in [deoxy-Hb] and concomitant decrease in [oxy-Hb], leading to apparent reductions of the O2-affinity of Hb (the entropic allostery). The heme in Hb is simply a low-affinity O2-trap, the coordination structure of which is not altered by static T-/R-quaternary and associated tertiary structural changes of Hb. Thus, heterotrophic effectors are the signal molecule, which acts as a functional link between these two allosteries and generates the diverse functionality of Hb of physiological relevance. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Local motifs in proteins combine to generate global functional moves

Literature on the topological properties of folded proteins that has emerged as a field in its own right in the past decade is reviewed. Physics-based construction of coarse-grained models of proteins from knowledge of all-atom coordinates of the average structure is discussed. Once network is thus obtained with the node and link information, local motifs provide plethora of information on protein function. The hierarchical structure of the proteins manifested in the interrelations of local motifs is emphasized. Motifs are also related to modularity of the structure, and they quantify shifts in the landscapes upon conformational changes induced by, e.g. ligand binding. Redundancy emerges as a balance between local and global network descriptors and is related to the collectivity of the protein motions. Introducing weight on links followed by sequential removal of least cohesive contacts allows interactions in proteins to be represented as the superposition of essential and redundant sets. Lack of the former makes the network non-functional, while the latter ensures robust functioning under a wide range of perturbation scenarios.

Crystal structure of unliganded TRAP: implications for dynamic allostery

Allostery is vital to the function of many proteins. In some cases, rather than a direct steric effect, mutual modulation of ligand binding at spatially separated sites may be achieved through a change in protein dynamics. Thus changes in vibrational modes of the protein, rather than conformational changes, allow different ligand sites to communicate. Evidence for such an effect has been found in TRAP (trp RNA-binding attenuation protein), a regulatory protein found in species of Bacillus. TRAP is part of a feedback system to modulate expression of the trp operon, which carries genes involved in tryptophan synthesis. Negative feedback is thought to depend on binding of tryptophan-bound, but not unbound, TRAP to a specific mRNA leader sequence. We find that, contrary to expectations, at low temperatures TRAP is able to bind RNA in the absence of tryptophan, and that this effect is particularly strong in the case of Bacillus stearothermophilus TRAP. We have solved the crystal structure of this protein with no tryptophan bound, and find that much of the structure shows little deviation from the tryptophan-bound form. These data support the idea that tryptophan may exert its effect on RNA binding by TRAP through dynamic and not structural changes, and that tryptophan binding may be mimicked by low temperature.

Regulation of catch binding by allosteric transitions

An allosteric model is used to describe changes in lifetimes of biological receptor-ligand bonds subjected to an external force. Force-induced transitions between the two states of the allosteric site lead to changes in the receptor conformation. The ligand bound to the receptor fluctuates between two different potentials formed by the two receptor conformations. The effect of the force on the receptor-ligand interaction potential is described by the Bell mechanism. The probability of detecting the ligand in the bound state is found to depend on the relaxation times of both ligand and allosteric sites. An analytic expression for the bond lifetime is derived as a function of force. The formal theoretical results are used to explain the anomalous force and time dependences of the integrin-fibronectin bond lifetimes measured by atomic force microscopy (Kong, F.; et al J. Cell Biol. 2009, 185, 1275-1284). The analytic expression and model parameters describe very well all anomalous dependences identified in the experiments.

Substrate-modulated thermal fluctuations affect long-range allosteric signaling in protein homodimers: exemplified in CAP

The role of conformational dynamics in allosteric signaling of proteins is increasingly recognized as an important and subtle aspect of this ubiquitous phenomenon. Cooperative binding is commonly observed in proteins with twofold symmetry that bind two identical ligands. We construct a coarse-grained model of an allosteric coupled dimer and show how the signal can be propagated between the distant binding sites via change in slow global vibrational modes alone. We demonstrate that modulation on substrate binding of as few as 5-10 slow modes can give rise to cooperativity observed in biological systems and that the type of cooperativity is given by change of interaction between the two monomers upon ligand binding. To illustrate the application of the model, we apply it to a challenging test case: the catabolite activator protein (CAP). CAP displays negative cooperativity upon association with two identical ligands. The conformation of CAP is not affected by the binding, but its vibrational spectrum undergoes a strong modification. Intriguingly, the first binding enhances thermal fluctuations, yet the second quenches them. We show that this counterintuitive behavior is, in fact, necessary for an optimal anticooperative system, and captured within a well-defined region of the model's parameter space. From analyzing the experimental results, we conclude that fast local modes take an active part in the allostery of CAP, coupled to the more-global slow modes. By including them into the model, we elucidate the role of the modes on different timescales. We conclude that such dynamic control of allostery in homodimers may be a general phenomenon and that our model framework can be used for extended interpretation of thermodynamic parameters in other systems.

All-atom normal-mode analysis reveals an RNA-induced allostery in a bacteriophage coat protein

Assembly of the T=3 bacteriophage MS2 is initiated by the binding of a 19 nucleotide RNA stem loop from within the phage genome to a symmetric coat protein dimer. This binding event effects a folding of the FG loop in one of the protein subunits of the dimer and results in the formation of an asymmetric dimer. Since both the symmetric and asymmetric forms of the dimer are needed for the assembly of the protein container, this allosteric switch plays an important role in the life cycle of the phage. We provide here details of an all-atom normal-mode analysis of this allosteric effect. The results suggest that asymmetric contacts between the A -duplex RNA phosphodiester backbone of the stem loop with the EF loop in one coat protein subunit results in an increased dynamic behavior of its FG loop. The four lowest-frequency modes, which encompass motions predominantly on the FG loops, account for over 90% of the increased dynamic behavior due to a localization of the vibrational pattern on a single FG loop. Finally, we show that an analysis of the allosteric effect using an elastic network model fails to predict this localization effect, highlighting the importance of using an all-atom full force field method for this problem.

Circadian KaiC phosphorylation: a multi-layer network

Circadian KaiC phosphorylation in cyanobacteria reconstituted in vitro recently initiates a series of studies experimentally and theoretically to explore its mechanism. In this paper, we report a dynamic diversity in hexameric KaiC phosphoforms using a multi-layer reaction network based on the nonequivalence of the dual phosphorylation sites (S431 and T432) in each KaiC subunit. These diverse oscillatory profiles can generate a kaleidoscopic phase modulation pattern probably responsible for the genome-wide transcription rhythms directly and/or indirectly in cyanobacteria. Particularly, our model reveals that a single KaiC hexamer is an energy-based, phosphorylation-dependent and self-regulated circadian oscillator modulated by KaiA and KaiB. We suggest that T432 is the main regulator for the oscillation amplitude, while S431 is the major phase regulator. S431 and T432 coordinately control the phosphorylation period. Robustness of the Kai network was examined by mixing samples in different phases, and varying protein concentrations and temperature. Similar results were obtained regardless of the deterministic or stochastic method employed. Therefore, the dynamic diversities and robustness of Kai oscillator make it a qualified core pacemaker that controls the cellular processes in cyanobacteria pervasively and accurately.

Circadian clocks are endogenous timing mechanisms that allow living organisms to coordinate their activities with daily environmental fluctuations. In cyanobacteria, almost all the genes are rhythmically expressed with the same ∼24 h period yet exhibit a variety of phase relationships and waveforms. Remarkably, the core pacemaker ticks robustly via simple biochemical reactions carried out by three Kai proteins: KaiC undergoes circadian phosphorylation in the presence of KaiA, KaiB and ATP. In this work, we propose a reaction network modeling the Kai oscillator based on the differentiation of dual phosphorylation sites. We found a dynamic diversity in KaiC phosphorylation which may serve as a potential regulatory mechanism related to the diverse-phased genome-wide expressions in cyanobacteria. In addition, we deduce that each KaiC hexamer is a single oscillator in regulating its own phosphorylation and interactions with KaiA or/and KaiB. In complex organisms, a number of key clock components possess similar activities (e.g., phosphorylation) with multiple nonequivalent active sites, and they may also show some unusual dynamic features that are embedded in the proteins' own reaction networks. We hope our work could be helpful to study the correlations between gene expressions and circadian rhythm in prokaryotic cells, even in eukaryotic cells.

Dynamic allostery in the methionine repressor revealed by force distribution analysis

Many fundamental cellular processes such as gene expression are tightly regulated by protein allostery. Allosteric signal propagation from the regulatory to the active site requires long-range communication, the molecular mechanism of which remains a matter of debate. A classical example for long-range allostery is the activation of the methionine repressor MetJ, a transcription factor. Binding of its co-repressor SAM increases its affinity for DNA several-fold, but has no visible conformational effect on its DNA binding interface. Our molecular dynamics simulations indicate correlated domain motions within MetJ, and quenching of these dynamics upon SAM binding entropically favors DNA binding. From monitoring conformational fluctuations alone, it is not obvious how the presence of SAM is communicated through the largely rigid core of MetJ and how SAM thereby is able to regulate MetJ dynamics. We here directly monitored the propagation of internal forces through the MetJ structure, instead of relying on conformational changes as conventionally done. Our force distribution analysis successfully revealed the molecular network for strain propagation, which connects collective domain motions through the protein core. Parts of the network are directly affected by SAM binding, giving rise to the observed quenching of fluctuations. Our results are in good agreement with experimental data. The force distribution analysis suggests itself as a valuable tool to gain insight into the molecular function of a whole class of allosteric proteins.

Evolutionary basis for the coupled-domain motions in Thermus thermophilus leucyl-tRNA synthetase

Aminoacyl-tRNA synthetases are multidomain proteins that catalyze the covalent attachment of amino acids to their cognate transfer RNA. Various domains of an aminoacyl-tRNA synthetase perform their specific functions in a highly coordinated manner to maintain high accuracy in protein synthesis in cells. The coordination of their function, therefore, requires communication between domains. In this study we explored the relevance of enzyme motion in domain-domain communications. Specifically, we attempted to probe whether the communication between distantly located domains of a multidomain protein is accomplished through a coordinated movement of structural elements. We investigated the collective motion in Thermus thermophilus leucyl-tRNA synthetase by studying the low frequency normal modes. We identified the mode that best described the experimentally observed conformational changes of T. thermophilus leucyl-tRNA synthetase upon substrate binding and analyzed the correlated and anticorrelated motions between different domains. Furthermore, we used statistical coupling analysis to explore if the amino acid pairs and/or clusters whose motions are thermally coupled have also coevolved. Our study demonstrates that a small number of residues belong to the category whose coupled thermal motions correspond to evolutionary coupling as well. These residue clusters constitute a distinguished set of interacting networks that are sparsely distributed in the domain interface. Residues of these networking clusters are within van der Waals contact, and we suggest that they are critical in the propagation of long range mechanochemical motions in T. thermophilus leucyl-tRNA synthetase.

Allostery: absence of a change in shape does not imply that allostery is not at play

Allostery is essential for controlled catalysis, signal transmission, receptor trafficking, turning genes on and off, and apoptosis. It governs the organism's response to environmental and metabolic cues, dictating transient partner interactions in the cellular network. Textbooks taught us that allostery is a change of shape at one site on the protein surface brought about by ligand binding to another. For several years, it has been broadly accepted that the change of shape is not induced; rather, it is observed simply because a larger protein population presents it. Current data indicate that while side chains can reorient and rewire, allostery may not even involve a change of (backbone) shape. Assuming that the enthalpy change does not reverse the free-energy change due to the change in entropy, entropy is mainly responsible for binding.

Contributions to catalysis and potential interactions of the three catalytic domains in a contiguous trimeric creatine kinase

Three separate creatine kinase (CK) isoform families exist in animals. Two of these (cytoplasmic and mitochondrial) are obligate oligomers. A third, flagellar, is monomeric but contains the residues for three complete CK domains. It is not known whether the active sites in each of the contiguous flagellar domains are catalytically competent, and, if so, whether they are capable of acting independently. Here we have utilized site-directed mutagenesis to selectively disable individual active sites and all possible combinations thereof. Kinetic studies showed that these mutations had minimal impact on substrate binding and synergism. Interestingly, the active sites were not catalytically equivalent, and were in fact interdependent, a phenomenon that has previously been reported only in the oligomeric CK isoforms.

Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation

In recent years, there has been a surge in the number of studies exploring the relationship between proteins' equilibrium dynamics and structural changes involved in function. An emerging concept, supported by both theory and experiments, is that under native state conditions proteins have an intrinsic ability to sample conformations that meet functional requirements. A typical example is the ability of enzymes to sample open and closed forms, irrespective of substrate, succeeded by the stabilization of one form (usually closed) upon substrate binding. This ability is structure-encoded, and plays a key role in facilitating allosteric regulation, which suggests complementing the sequence-encodes-structure paradigm of protein science by structure-encodes-dynamics-encodes-function. The emerging connection implies an evolutionary role in selecting/conserving structures based on their ability to achieve functional dynamics, and in turn, selecting sequences that fold into such 'apt' structures.

Nonequilibrium dynamic mechanism for allosteric effect

Allosteric regulation is often viewed as thermodynamic in nature. However, protein internal motions during an enzymatic reaction cycle can slow the hoping processes over numerous potential barriers. We propose that regulating molecules may function by modifying the nonequilibrium protein dynamics. The theory predicts that an enzyme under the new mechanism has a different temperature dependence, waiting time distribution of the turnover cycle, and dynamic fluctuation patterns with and without an effector. Experimental tests of the theory are proposed.

A simple physical description of DNA dynamics: quasi-harmonic analysis as a route to the configurational entropy

It has become increasingly apparent that the dynamic as well as the structural properties of biological macromolecules are important to their function. However, information concerning molecular flexibility can be difficult to obtain experimentally at the atomic level. Computer modelling techniques such as molecular dynamics (MD) have therefore proved invaluable in advancing our understanding of biomolecular flexibility. This paper describes how a combination of atomistic MD simulations and quasi-harmonic analysis can be used to describe the dynamics of duplex DNA, with a particular emphasis on methods for calculating differences in configurational entropies. We demonstrate that DNA possesses remarkably simple mechanical properties relative to globular proteins, making it an ideal system for exploring biomolecular flexibility in general. Our results also highlight the importance of solvent viscosity in determining the dynamic behaviour of DNA in aqueous solution.

Modular architecture of protein structures and allosteric communications: potential implications for signaling proteins and regulatory linkages

BACKGROUND

Allosteric communications are vital for cellular signaling. Here we explore a relationship between protein architectural organization and shortcuts in signaling pathways.

RESULTS

We show that protein domains consist of modules interconnected by residues that mediate signaling through the shortest pathways. These mediating residues tend to be located at the inter-modular boundaries, which are more rigid and display a larger number of long-range interactions than intra-modular regions. The inter-modular boundaries contain most of the residues centrally conserved in the protein fold, which may be crucial for information transfer between amino acids. Our approach to modular decomposition relies on a representation of protein structures as residue-interacting networks, and removal of the most central residue contacts, which are assumed to be crucial for allosteric communications. The modular decomposition of 100 multi-domain protein structures indicates that modules constitute the building blocks of domains. The analysis of 13 allosteric proteins revealed that modules characterize experimentally identified functional regions. Based on the study of an additional functionally annotated dataset of 115 proteins, we propose that high-modularity modules include functional sites and are the basic functional units. We provide examples (the Galphas subunit and P450 cytochromes) to illustrate that the modular architecture of active sites is linked to their functional specialization.

CONCLUSION

Our method decomposes protein structures into modules, allowing the study of signal transmission between functional sites. A modular configuration might be advantageous: it allows signaling proteins to expand their regulatory linkages and may elicit a broader range of control mechanisms either via modular combinations or through modulation of inter-modular linkages.

A mutation in the S-switch region of the Runt domain alters the dynamics of an allosteric network responsible for CBFbeta regulation

The Runt domain is the DNA binding domain of the core binding factor (CBF) Runx subunits. The CBFs are transcription factors that play critical roles in hematopoiesis, bone, and neuron development in mammals. A common non-DNA binding CBFbeta subunit heterodimerizes with the Runt domain of the Runx proteins and allosterically regulates its affinity for DNA. Previous NMR dynamics studies suggested a model whereby CBFbeta allosterically regulates DNA binding by quenching conformational exchange in the Runt domain, particularly in the S-switch region and the betaE'-F loop. We sought to test this model, and to this end introduced all possible single amino acid substitutions into the S-switch region and the betaE'-F loop, and screened for mutations that enhanced DNA-binding. We demonstrate that one Runt domain mutant, R164N, binds both DNA and CBFbeta with higher affinity, but it is less sensitive to allosteric regulation by CBFbeta. Analysis of NMR relaxation data shows that the chemical exchange exhibited by the wild-type Runt domain is largely quenched by the R164N substitution. These data support a model in which the dynamic behavior of a network of residues connecting the CBFbeta and DNA binding sites on the Runt domain plays a critical role in the mechanism of allosteric regulation. This study provides an important functional link between dynamic behavior and protein allosteric function, consistent with results on other allosterically regulated proteins.