Protein activity regulation by conformational entropy

Genetically tunable frustration controls allostery in an intrinsically disordered transcription factor

Intrinsically disordered proteins (IDPs) present a functional paradox because they lack stable tertiary structure, but nonetheless play a central role in signaling, utilizing a process known as allostery. Historically, allostery in structured proteins has been interpreted in terms of propagated structural changes that are induced by effector binding. Thus, it is not clear how IDPs, lacking such well-defined structures, can allosterically affect function. Here, we show a mechanism by which an IDP can allosterically control function by simultaneously tuning transcriptional activation and repression, using a novel strategy that relies on the principle of ‘energetic frustration’. We demonstrate that human glucocorticoid receptor tunes this signaling in vivo by producing translational isoforms differing only in the length of the disordered region, which modulates the degree of frustration. We expect this frustration-based model of allostery will prove to be generally important in explaining signaling in other IDPs.

Proteins carry out most of the key tasks inside cells. To perform these roles, proteins must fold up to form complex three-dimensional structures. Researchers used to think that the useful parts of proteins all had set structures. However, we now know that ‘disordered’ proteins with variable structures are common and disordered parts of proteins can have vital roles.

In a process called allosteric regulation, regulator molecules can increase or decrease the activity of a protein by binding to it. This binding was thought to work by changing the structure of the protein, but it was not clear how this works in disordered proteins. To investigate, Li et al. studied a disordered protein called glucocorticoid receptor, and found that disordered regions can have opposing effects on other regions of the protein. This creates a ‘tug-of-war’ that Li et al. term “energetic frustration”, whereby the activity of the protein results from the combination of the opposing interactions.

Further investigation revealed that the glucorticoid receptor produces different versions of itself that have different degrees of energetic frustration, which alters how effectively the proteins perform their tasks. This means that the protein can regulate its own activity even in the absence of binding to regulator molecules.

The concept of energetic frustration could enhance our understanding of the many different proteins that contain disordered regions. Eventually, this knowledge could be used to develop drugs that alter the activity of these proteins and so could form part of treatments for a wide range of conditions including autoimmune diseases (such as rheumatoid arthritis and lupus), cancers, and organ rejection for transplant patients. The results presented by Li et al. suggest where more research is needed to achieve this goal. For example, we need to understand more about the stability of disordered protein regions, and to identify which surfaces of the proteins interact with each other.

Ligand modulation of sidechain dynamics in a wild-type human GPCR

GPCRs regulate all aspects of human physiology, and biophysical studies have deepened our understanding of GPCR conformational regulation by different ligands. Yet there is no experimental evidence for how sidechain dynamics control allosteric transitions between GPCR conformations. To address this deficit, we generated samples of a wild-type GPCR (A_2AR) that are deuterated apart from ¹H/¹³C NMR probes at isoleucine δ1 methyl groups, which facilitated ¹H/¹³C methyl TROSY NMR measurements with opposing ligands. Our data indicate that low [Na⁺] is required to allow large agonist-induced structural changes in A_2AR, and that patterns of sidechain dynamics substantially differ between agonist (NECA) and inverse agonist (ZM241385) bound receptors, with the inverse agonist suppressing fast ps-ns timescale motions at the G protein binding site. Our approach to GPCR NMR creates a framework for exploring how different regions of a receptor respond to different ligands or signaling proteins through modulation of fast ps-ns sidechain dynamics.

A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations

Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.

Atomic view of the energy landscape in the allosteric regulation of Abl kinase

The activity of protein kinases is often regulated in an intramolecular fashion by signaling domains, which feature several phosphorylation or protein-docking sites. How kinases integrate such distinct binding and signaling events to regulate their activities is unclear, especially in quantitative terms. We used NMR spectroscopy to show how structural elements within the Abl regulatory module (RM) synergistically generate a multilayered allosteric mechanism that enables Abl kinase to function as a finely tuned switch. We dissected the structure and energetics of the regulatory mechanism to precisely measure the effects of various activating or inhibiting stimuli on Abl kinase activity. The data provide a mechanistic basis explaining genetic observations and reveal a previously unknown activator region within Abl. Our findings show that drug-resistance mutations in the Abl RM exert their allosteric effect by promoting the activated state of Abl and not by decreasing the drug affinity for the kinase.

Conformational Rigidity and Protein Dynamics at Distinct Timescales Regulate PTP1B Activity and Allostery

Protein function originates from a cooperation of structural rigidity, dynamics at different timescales, and allostery. However, how these three pillars of protein function are integrated is still only poorly understood. Here we show how these pillars are connected in Protein Tyrosine Phosphatase 1B (PTP1B), a drug target for diabetes and cancer that catalyzes the dephosphorylation of numerous substrates in essential signaling pathways. By combining new experimental and computational data on WT-PTP1B and ≥10 PTP1B variants in multiple states, we discovered a fundamental and evolutionarily conserved CH/π switch that is critical for positioning the catalytically important WPD loop. Furthermore, our data show that PTP1B uses conformational and dynamic allostery to regulate its activity. This shows that both conformational rigidity and dynamics are essential for controlling protein activity. This connection between rigidity and dynamics at different timescales is likely a hallmark of all enzyme function.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

New insight on the S100A1–STIP1 complex highlights the important relationship between allostery and entropy in protein function

Calcium signaling serves as a nexus of many vital cellular processes. Of particular importance is the role the calcium signaling plays in the prevention of protein misfolding, and the S100 family of calcium-binding proteins is a key player in this pathway. While the S100 proteins carry out a range of roles, the interaction of S100A1 and the stress-inducible phosphoprotein 1 (STIP1) has been shown to be particularly important. A recent study by Maciejewski et al. in Biochemical Journal (Biochemical Journal (2017) 474, 1853–1866) revealed new insights into the nature of the S100A1–STIP1 interaction. Not only did the present paper indicate the stoichiometry of binding for this interaction (three S100A1 dimers : one STIP1), it also demonstrated that the binding interaction is highly co-operative and that each S100A1–STIP1-binding interaction is entropically driven. The findings presented raise important new questions regarding the relationship between entropy and allostery in protein function. Recently, the dynamical underpinnings of allostery in protein function have become a topic of increased interest. A broad range of investigations have demonstrated that allostery can be mediated by entropic processes such as changes in the flexibility of the protein backbone and in the range of motions explored by side chains. The S100A1–STIP1 complex as described by Maciejewski et al. suggests a new system in which an allosteric-binding interaction driven by entropic processes may be systematically dissected in the future.

Thermodynamic and Structural Adaptation Differences between the Mesophilic and Psychrophilic Lactate Dehydrogenases

The thermodynamics of substrate binding and enzymatic activity of a glycolytic enzyme, lactate dehydrogenase (LDH), from both porcine heart, phLDH (Sus scrofa; a mesophile), and mackerel icefish, cgLDH (Chamapsocephalus gunnari; a psychrophile), were investigated. Using a novel and quite sensitive fluorescence assay that can distinguish protein conformational changes close to and distal from the substrate binding pocket, a reversible global protein structural transition preceding the high-temperature transition (denaturation) was surprisingly found to coincide with a marked change in enzymatic activity for both LDHs. A similar reversible structural transition of the active site structure was observed for phLDH but not for cgLDH. An observed lower substrate binding affinity for cgLDH compared to that for phLDH was accompanied by a larger contribution of entropy to ΔG, which reflects a higher functional plasticity of the psychrophilic cgLDH compared to that of the mesophilic phLDH. The natural osmolyte, trimethylamine N-oxide (TMAO), increases stability and shifts all structural transitions to higher temperatures for both orthologs while simultaneously reducing catalytic activity. The presence of TMAO causes cgLDH to adopt catalytic parameters like those of phLDH in the absence of the osmolyte. Our results are most naturally understood within a model of enzyme dynamics whereby different conformations of the enzyme that have varied catalytic parameters (i.e., binding and catalytic proclivity) and whose population profiles are temperature-dependent and influenced by osmolytes interconvert among themselves. Our results also show that adaptation can be achieved by means other than gene mutations and complements the synchronic evolution of the cellular milieu.

Biased Allostery

G-protein-coupled receptors (GPCRs) constitute a large group of integral membrane proteins that transduce extracellular signals from a wide range of agonists into targeted intracellular responses. Although the responses can vary depending on the category of G-proteins activated by a particular receptor, responses were also found to be triggered by interactions of the receptor with β-arrestins. It was subsequently discovered that for the same receptor molecule (e.g., the β-adrenergic receptor), some agonists have a propensity to specifically favor responses by G-proteins, others by β-arrestins, as has now been extensively studied. This feature of the GPCR system is known as biased agonism and is subject to various interpretations, including agonist-induced conformational change versus selective stabilization of preexisting active conformations. Here, we explore a complete allosteric framework for biased agonism based on alternative preexisting conformations that bind more strongly, but nonexclusively, either G-proteins or β-arrestins. The framework incorporates reciprocal effects among all interacting molecules. As a result, G-proteins and β-arrestins are in steric competition for binding to the cytoplasmic surface of either the G-protein-favoring or β-arrestin-favoring GPCR conformation. Moreover, through linkage relations, the strength of the interactions of G-proteins or β-arrestins with the corresponding active conformation potentiates the apparent affinity for the agonist, effectively equating these two proteins to allosteric modulators. The balance between response alternatives can also be influenced by the physiological concentrations of either G-proteins or β-arrestins, as well as by phosphorylation or interactions with positive or negative allosteric modulators. The nature of the interactions in the simulations presented suggests novel experimental tests to distinguish more fully among alternative mechanisms.

Mutations in Antibody Fragments Modulate Allosteric Response Via Hydrogen-Bond Network Fluctuations

A mechanical perturbation method that locally restricts conformational entropy along the protein backbone is used to identify putative allosteric sites in a series of antibody fragments. The method is based on a distance constraint model that integrates mechanical and thermodynamic viewpoints of protein structure wherein mechanical clamps that mimic substrate or cosolute binding are introduced. Across a set of six single chain-Fv fragments of the anti-lymphotoxin-β receptor antibody, statistically significant responses are obtained by averaging over 10 representative structures sampled from a molecular dynamics simulation. As expected, the introduced clamps locally rigidify the protein, but long-ranged increases in both rigidity and flexibility are also frequently observed. Expanding our analysis to every molecular dynamics frame demonstrates that the allosteric responses are modulated by fluctuations within the hydrogen-bond network where the native ensemble is comprised of conformations that both are, and are not, affected by the perturbation in question. Population shifts induced by the mutations alter the allosteric response by adjusting which hydrogen-bond networks are the most probable. These effects are compared using response maps that track changes across each single chain-Fv fragment, thus providing valuable insight into how sensitive allosteric mechanisms are to mutations.

Entropy in molecular recognition by proteins

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

Elucidating the energetic contributions to the binding free energy

New exact equations are derived for the terms contributing to the binding free energy (ΔG⁰) of a ligand-receptor pair using our recently introduced formalism which we here call perturbation-divergence formalism (PDF). Specifically, ΔG⁰ equals the sum of the average of the perturbation (pertaining to new interactions) and additional dissipative terms. The average of the perturbation includes the sum of the average receptor-ligand interactions and the average of the change of solvation energies upon association. The Kullback-Leibler (KL) divergence quantifies the energetically dissipative terms, which are due to the configurational changes and, using the chain rule of KL divergence, can be decomposed into (i) dissipation due to limiting the external liberation (translation and rotation) of the ligand relative to the receptor and (ii) dissipation due to conformational (internal) changes inside the receptor and the ligand. We also identify all exactly canceling energetic terms which do not contribute to ΔG⁰. Furthermore, the PDF provides a new approach towards dimensionality reduction in the representation of the association process and towards relating the dynamic (high dimensional) with the thermodynamic (one-dimensional) changes.

The Dynameomics Entropy Dictionary: A Large-Scale Assessment of Conformational Entropy across Protein Fold Space

Molecular dynamics (MD) simulations contain considerable information with regard to the motions and fluctuations of a protein, the magnitude of which can be used to estimate conformational entropy. Here we survey conformational entropy across protein fold space using the Dynameomics database, which represents the largest existing data set of protein MD simulations for representatives of essentially all known protein folds. We provide an overview of MD-derived entropies accounting for all possible degrees of dihedral freedom on an unprecedented scale. Although different side chains might be expected to impose varying restrictions on the conformational space that the backbone can sample, we found that the backbone entropy and side chain size are not strictly coupled. An outcome of these analyses is the Dynameomics Entropy Dictionary, the contents of which have been compared with entropies derived by other theoretical approaches and experiment. As might be expected, the conformational entropies scale linearly with the number of residues, demonstrating that conformational entropy is an extensive property of proteins. The calculated conformational entropies of folding agree well with previous estimates. Detailed analysis of specific cases identifies deviations in conformational entropy from the average values that highlight how conformational entropy varies with sequence, secondary structure, and tertiary fold. Notably, α-helices have lower entropy on average than do β-sheets, and both are lower than coil regions.

Quantifying Allosteric Communication via Both Concerted Structural Changes and Conformational Disorder with CARDS

Allosteric (i.e., long-range) communication within proteins is crucial for many biological processes, such as the activation of signaling cascades in response to specific stimuli. However, the physical basis for this communication remains unclear. Existing computational methods for identifying allostery focus on the role of concerted structural changes, but recent experimental work demonstrates that disorder is also an important factor. Here, we introduce the Correlation of All Rotameric and Dynamical States (CARDS) framework for quantifying correlations between both the structure and disorder of different regions of a protein. To quantify disorder, we draw inspiration from methods for quantifying "dynamic heterogeneity" from chemical physics to classify segments of a dihedral's time evolution as being in either ordered or disordered regimes. To demonstrate the utility of this approach, we apply CARDS to the Catabolite Activator Protein (CAP), a transcriptional activator that is regulated by Cyclic Adenosine MonoPhosphate (cAMP) binding. We find that CARDS captures allosteric communication between the two cAMP-Binding Domains (CBDs). Importantly, CARDS reveals that this coupling is dominated by disorder-mediated correlations, consistent with NMR experiments that establish allosteric coupling between the CBDs occurs without a concerted structural change. CARDS also recapitulates an enhanced role for disorder in the communication between the DNA-Binding Domains (DBDs) and CBDs in the S62F variant of CAP. Finally, we demonstrate that using CARDS to find communication hotspots identifies regions of CAP that are in allosteric communication without foreknowledge of their identities. Therefore, we expect CARDS to be of great utility for both understanding and predicting allostery.

The involvement of cannabinoids and mTOR in the reconsolidation of an emotional memory in the hippocampal-amygdala-insular circuit

Memory reconsolidation is the process in which reactivated long-term memory becomes transiently sensitive to amnesic agents. We evaluated the ability of post reactivation administration of the mTOR inhibitor rapamycin, separately and in combination with the cannabinoid CB1/2 receptor agonist WIN55,212-2 (WIN), given systemically or specifically into the hippocampal CA1 area, basolateral amygdala (BLA) or insular cortex (IC), to reduce inhibitory avoidance fear in rats. Systemic administration of rapamycin after reactivation of fear memory impaired reconsolidation and facilitated extinction. A combined treatment with WIN and rapamycin resulted in similar effects. WIN injected systemically facilitated extinction, with no effect on reconsolidation. WIN alone and with rapamycin also decreased anxiety-like behavior. Further, when spontaneous recovery was tested, the WIN+rapamycin group did not demonstrate recovery of fear which can occur spontaneously after the passage of time. Rapamycin and WIN had differential effects on reconsolidation and extinction when microinjected into the CA1, BLA and IC. Furthermore, exposure to shock increased p70s6K activation in the BLA, indicating activation by mTOR. Treatment with rapamycin, WIN or WIN+rapamycin decreased activation and there was a strong positive correlation between fear retrieval and p70s6K activation in the BLA, suggesting that enhanced fear retrieval is associated with enhanced p70s6K activation. Taken together, the results suggest that rapamycin or a combined treatment that involves blocking mTOR and activating cannabinoids may be a promising pharmacological approach for the attenuation of reactivated emotional memories, and thus, it could represent a potential treatment strategy for disorders associated with traumatic memories.

Entropy redistribution controls allostery in a metalloregulatory protein

Allosteric communication between two ligand-binding sites in a protein is a central aspect of biological regulation that remains mechanistically unclear. Here we show that perturbations in equilibrium picosecond-nanosecond motions impact zinc (Zn)-induced allosteric inhibition of DNA binding by the Zn efflux repressor CzrA (chromosomal zinc-regulated repressor). DNA binding leads to an unanticipated increase in methyl side-chain flexibility and thus stabilizes the complex entropically; Zn binding redistributes these motions, inhibiting formation of the DNA complex by restricting coupled fast motions and concerted slower motions. Allosterically impaired CzrA mutants are characterized by distinct nonnative fast internal dynamics "fingerprints" upon Zn binding, and DNA binding is weakly regulated. We demonstrate the predictive power of the wild-type dynamics fingerprint to identify key residues in dynamics-driven allostery. We propose that driving forces arising from dynamics can be harnessed by nature to evolve new allosteric ligand specificities in a compact molecular scaffold.

Site-to-site interdomain communication may mediate different loss-of-function mechanisms in a cancer-associated NQO1 polymorphism

Disease associated genetic variations often cause intracellular enzyme inactivation, dysregulation and instability. However, allosteric communication of mutational effects to distant functional sites leading to loss-of-function remains poorly understood. We characterize here interdomain site-to-site communication by which a common cancer-associated single nucleotide polymorphism (c.C609T/p.P187S) reduces the activity and stability in vivo of NAD(P)H:quinone oxidoreductase 1 (NQO1). NQO1 is a FAD-dependent, two-domain multifunctional stress protein acting as a Phase II enzyme, activating cancer pro-drugs and stabilizing p53 and p73α oncosuppressors. We show that p.P187S causes structural and dynamic changes communicated to functional sites far from the mutated site, affecting the FAD binding site located at the N-terminal domain (NTD) and accelerating proteasomal degradation through dynamic effects on the C-terminal domain (CTD). Structural protein:protein interaction studies reveal that the cancer-associated polymorphism does not abolish the interaction with p73α, indicating that oncosuppressor destabilization largely mirrors the low intracellular stability of p.P187S. In conclusion, we show how a single disease associated amino acid change may allosterically perturb several functional sites in an oligomeric and multidomain protein. These results have important implications for the understanding of loss-of-function genetic diseases and the identification of novel structural hot spots as targets for pharmacological intervention.

Knowledge-based entropies improve the identification of native protein structures

Evaluating protein structures requires reliable free energies with good estimates of both potential energies and entropies. Although there are many demonstrated successes from using knowledge-based potential energies, computing entropies of proteins has lagged far behind. Here we take an entirely different approach and evaluate knowledge-based conformational entropies of proteins based on the observed frequencies of contact changes between amino acids in a set of 167 diverse proteins, each of which has two alternative structures. The results show that charged and polar interactions break more often than hydrophobic pairs. This pattern correlates strongly with the average solvent exposure of amino acids in globular proteins, as well as with polarity indices and the sizes of the amino acids. Knowledge-based entropies are derived by using the inverse Boltzmann relationship, in a manner analogous to the way that knowledge-based potentials have been extracted. Including these new knowledge-based entropies almost doubles the performance of knowledge-based potentials in selecting the native protein structures from decoy sets. Beyond the overall energy-entropy compensation, a similar compensation is seen for individual pairs of interacting amino acids. The entropies in this report have immediate applications for 3D structure prediction, protein model assessment, and protein engineering and design.

Water-Restructuring Mutations Can Reverse the Thermodynamic Signature of Ligand Binding to Human Carbonic Anhydrase

This study uses mutants of human carbonic anhydrase (HCAII) to examine how changes in the organization of water within a binding pocket can alter the thermodynamics of protein-ligand association. Results from calorimetric, crystallographic, and theoretical analyses suggest that most mutations strengthen networks of water-mediated hydrogen bonds and reduce binding affinity by increasing the enthalpic cost and, to a lesser extent, the entropic benefit of rearranging those networks during binding. The organization of water within a binding pocket can thus determine whether the hydrophobic interactions in which it engages are enthalpy-driven or entropy-driven. Our findings highlight a possible asymmetry in protein-ligand association by suggesting that, within the confines of the binding pocket of HCAII, binding events associated with enthalpically favorable rearrangements of water are stronger than those associated with entropically favorable ones.

Recent advances in Major Histocompatibility Complex (MHC) class I antigen presentation: Plastic MHC molecules and TAPBPR-mediated quality control

We have known since the late 1980s that the function of classical major histocompatibility complex (MHC) class I molecules is to bind peptides and display them at the cell surface to cytotoxic T cells. Recognition by these sentinels of the immune system can lead to the destruction of the presenting cell, thus protecting the host from pathogens and cancer. Classical MHC class I molecules (MHC I hereafter) are co-dominantly expressed, polygenic, and exceptionally polymorphic and have significant sequence diversity. Thus, in most species, there are many different MHC I allotypes expressed, each with different peptide-binding specificity, which can have a dramatic effect on disease outcome. Although MHC allotypes vary in their primary sequence, they share common tertiary and quaternary structures. Here, we review the evidence that, despite this commonality, polymorphic amino acid differences between allotypes alter the ability of MHC I molecules to change shape (that is, their conformational plasticity). We discuss how the peptide loading co-factor tapasin might modify this plasticity to augment peptide loading. Lastly, we consider recent findings concerning the functions of the non-classical MHC I molecule HLA-E as well as the tapasin-related protein TAPBPR (transporter associated with antigen presentation binding protein-related), which has been shown to act as a second quality-control stage in MHC I antigen presentation.

Asymmetric configurations in a reengineered homodimer reveal multiple subunit communication pathways in protein allostery

Many allosteric proteins form homo-oligomeric complexes to regulate a biological function. In homo-oligomers, subunits establish communication pathways that are modulated by external stimuli like ligand binding. A challenge for dissecting the communication mechanisms in homo-oligomers is identifying intermediate liganded states, which are typically transiently populated. However, their identities provide the most mechanistic information on how ligand-induced signals propagate from bound to empty subunits. Here, we dissected the directionality and magnitude of subunit communication in a reengineered single-chain version of the homodimeric transcription factor cAMP receptor protein. By combining wild-type and mutant subunits in various asymmetric configurations, we revealed a linear relationship between the magnitude of cooperative effects and the number of mutant subunits. We found that a single mutation is sufficient to change the global allosteric behavior of the dimer even when one subunit was wild type. Dimers harboring two mutations with opposite cooperative effects had different allosteric properties depending on the arrangement of the mutations. When the two mutations were placed in the same subunit, the resulting cooperativity was neutral. In contrast, when placed in different subunits, the observed cooperativity was dominated by the mutation with strongest effects over cAMP affinity relative to wild type. These results highlight the distinct roles of intrasubunit interactions and intersubunit communication in allostery. Finally, dimers bound to either one or two cAMP molecules had similar DNA affinities, indicating that both asymmetric and symmetric liganded states activate DNA interactions. These studies have revealed the multiple communication pathways that homo-oligomers employ to transduce signals.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Mutation of a kinase allosteric node uncouples dynamics linked to phosphotransfer

The expertise of protein kinases lies in their dynamic structure, wherein they are able to modulate cellular signaling by their phosphotransferase activity. Only a few hundreds of protein kinases regulate key processes in human cells, and protein kinases play a pivotal role in health and disease. The present study dwells on understanding the working of the protein kinase-molecular switch as an allosteric network of "communities" composed of congruently dynamic residues that make up the protein kinase core. Girvan-Newman algorithm-based community maps of the kinase domain of cAMP-dependent protein kinase A allow for a molecular explanation for the role of protein conformational entropy in its catalytic cycle. The community map of a mutant, Y204A, is analyzed vis-à-vis the wild-type protein to study the perturbations in its dynamic profile such that it interferes with transfer of the γ-phosphate to a protein substrate. Conventional biochemical measurements are used to ascertain the effect of these dynamic perturbations on the kinetic profiles of both proteins. These studies pave the way for understanding how mutations far from the kinase active site can alter its dynamic properties and catalytic function even when major structural perturbations are not obvious from static crystal structures.

Solvent effects on ligand binding to a serine protease

ITC experiments and MD simulations reveal the mechanism behind enthalpy/entropy compensation upon trypsin-benzamidine binding at different solvation conditions.

Structural biology of supramolecular assemblies by magic-angle spinning NMR spectroscopy

In recent years, exciting developments in instrument technology and experimental methodology have advanced the field of magic-angle spinning (MAS) nuclear magnetic resonance (NMR) to new heights. Contemporary MAS NMR yields atomic-level insights into structure and dynamics of an astounding range of biological systems, many of which cannot be studied by other methods. With the advent of fast MAS, proton detection, and novel pulse sequences, large supramolecular assemblies, such as cytoskeletal proteins and intact viruses, are now accessible for detailed analysis. In this review, we will discuss the current MAS NMR methodologies that enable characterization of complex biomolecular systems and will present examples of applications to several classes of assemblies comprising bacterial and mammalian cytoskeleton as well as human immunodeficiency virus 1 and bacteriophage viruses. The body of work reviewed herein is representative of the recent advancements in the field, with respect to the complexity of the systems studied, the quality of the data, and the significance to the biology.

Photoactivation Reduces Side-Chain Dynamics of a LOV Photoreceptor

We used neutron-scattering experiments to probe the conformational dynamics of the light, oxygen, voltage (LOV) photoreceptor PpSB1-LOV from Pseudomonas putida in both the dark and light states. Global protein diffusion and internal macromolecular dynamics were measured using incoherent neutron time-of-flight and backscattering spectroscopy on the picosecond to nanosecond timescales. Global protein diffusion of PpSB1-LOV is not influenced by photoactivation. Observation-time-dependent global diffusion coefficients were found, which converge on the nanosecond timescale toward diffusion coefficients determined by dynamic light scattering. Mean-square displacements of localized internal motions and effective force constants, <k'>, describing the resilience of the proteins were determined on the respective timescales. Photoactivation significantly modifies the flexibility and the resilience of PpSB1-LOV. On the fast, picosecond timescale, small changes in the mean-square displacement and <k'> are observed, which are enhanced on the slower, nanosecond timescale. Photoactivation results in a slightly larger resilience of the photoreceptor on the fast, picosecond timescale, whereas in the nanosecond range, a significantly less resilient structure of the light-state protein is observed. For a residue-resolved interpretation of the experimental neutron-scattering data, we analyzed molecular dynamics simulations of the PpSB1-LOV X-ray structure. Based on these data, it is tempting to speculate that light-induced changes in the protein result in altered side-chain mobility mostly for residues on the protruding Jα helix and on the LOV-LOV dimer interface. Our results provide strong experimental evidence that side-chain dynamics play a crucial role in photoactivation and signaling of PpSB1-LOV via modulation of conformational entropy.

Exploring Molecular Mechanisms of Paradoxical Activation in the BRAF Kinase Dimers: Atomistic Simulations of Conformational Dynamics and Modeling of Allosteric Communication Networks and Signaling Pathways

The recent studies have revealed that most BRAF inhibitors can paradoxically induce kinase activation by promoting dimerization and enzyme transactivation. Despite rapidly growing number of structural and functional studies about the BRAF dimer complexes, the molecular basis of paradoxical activation phenomenon is poorly understood and remains largely hypothetical. In this work, we have explored the relationships between inhibitor binding, protein dynamics and allosteric signaling in the BRAF dimers using a network-centric approach. Using this theoretical framework, we have combined molecular dynamics simulations with coevolutionary analysis and modeling of the residue interaction networks to determine molecular determinants of paradoxical activation. We have investigated functional effects produced by paradox inducer inhibitors PLX4720, Dabrafenib, Vemurafenib and a paradox breaker inhibitor PLX7904. Functional dynamics and binding free energy analyses of the BRAF dimer complexes have suggested that negative cooperativity effect and dimer-promoting potential of the inhibitors could be important drivers of paradoxical activation. We have introduced a protein structure network model in which coevolutionary residue dependencies and dynamic maps of residue correlations are integrated in the construction and analysis of the residue interaction networks. The results have shown that coevolutionary residues in the BRAF structures could assemble into independent structural modules and form a global interaction network that may promote dimerization. We have also found that BRAF inhibitors could modulate centrality and communication propensities of global mediating centers in the residue interaction networks. By simulating allosteric communication pathways in the BRAF structures, we have determined that paradox inducer and breaker inhibitors may activate specific signaling routes that correlate with the extent of paradoxical activation. While paradox inducer inhibitors may facilitate a rapid and efficient communication via an optimal single pathway, the paradox breaker may induce a broader ensemble of suboptimal and less efficient communication routes. The central finding of our study is that paradox breaker PLX7904 could mimic structural, dynamic and network features of the inactive BRAF-WT monomer that may be required for evading paradoxical activation. The results of this study rationalize the existing structure-functional experiments by offering a network-centric rationale of the paradoxical activation phenomenon. We argue that BRAF inhibitors that amplify dynamic features of the inactive BRAF-WT monomer and intervene with the allosteric interaction networks may serve as effective paradox breakers in cellular environment.

Enthalpy-entropy compensation: the role of solvation

Structural modifications to interacting systems frequently lead to changes in both the enthalpy (heat) and entropy of the process that compensate each other, so that the Gibbs free energy is little changed: a major barrier to the development of lead compounds in drug discovery. The conventional explanation for such enthalpy–entropy compensation (EEC) is that tighter contacts lead to a more negative enthalpy but increased molecular constraints, i.e., a compensating conformational entropy reduction. Changes in solvation can also contribute to EEC but this contribution is infrequently discussed. We review long-established and recent cases of EEC and conclude that the large fluctuations in enthalpy and entropy observed are too great to be a result of only conformational changes and must result, to a considerable degree, from variations in the amounts of water immobilized or released on forming complexes. Two systems exhibiting EEC show a correlation between calorimetric entropies and local mobilities, interpreted to mean conformational control of the binding entropy/free energy. However, a substantial contribution from solvation gives the same effect, as a consequence of a structural link between the amount of bound water and the protein flexibility. Only by assuming substantial changes in solvation—an intrinsically compensatory process—can a more complete understanding of EEC be obtained. Faced with such large, and compensating, changes in the enthalpies and entropies of binding, the best approach to engineering elevated affinities must be through the addition of ionic links, as they generate increased entropy without affecting the enthalpy.

Metal Ion Binding at the Catalytic Site Induces Widely Distributed Changes in a Sequence Specific Protein-DNA Complex

Metal ion cofactors can alter the energetics and specificity of sequence specific protein–DNA interactions, but it is unknown if the underlying effects on structure and dynamics are local or dispersed throughout the protein–DNA complex. This work uses EcoRV endonuclease as a model, and catalytically inactive lanthanide ions, which replace the Mg²⁺ cofactor. Nuclear magnetic resonance (NMR) titrations indicate that four Lu³⁺ or two La³⁺ cations bind, and two new crystal structures confirm that Lu³⁺ binding is confined to the active sites. NMR spectra show that the metal-free EcoRV complex with cognate (GATATC) DNA is structurally distinct from the nonspecific complex, and that metal ion binding sites are not assembled in the nonspecific complex. NMR chemical shift perturbations were determined for ¹H–¹⁵N amide resonances, for ¹H–¹³C Ile-δ-CH₃ resonances, and for stereospecifically assigned Leu-δ-CH₃ and Val-γ-CH₃ resonances. Many chemical shifts throughout the cognate complex are unperturbed, so metal binding does not induce major conformational changes. However, some large perturbations of amide and side chain methyl resonances occur as far as 34 Å from the metal ions. Concerted changes in specific residues imply that local effects of metal binding are propagated via a β-sheet and an α-helix. Both amide and methyl resonance perturbations indicate changes in the interface between subunits of the EcoRV homodimer. Bound metal ions also affect amide hydrogen exchange rates for distant residues, including a distant subdomain that contacts DNA phosphates and promotes DNA bending, showing that metal ions in the active sites, which relieve electrostatic repulsion between protein and DNA, cause changes in slow dynamics throughout the complex.

Distinct Roles for Conformational Dynamics in Protein-Ligand Interactions

Conformational dynamics has an established role in enzyme catalysis, but its contribution to ligand binding and specificity is largely unexplored. Here we used the Tiam1 PDZ domain and an engineered variant (QM PDZ) with broadened specificity to investigate the role of structure and conformational dynamics in molecular recognition. Crystal structures of the QM PDZ domain both free and bound to ligands showed structural features central to binding (enthalpy), while nuclear-magnetic-resonance-based methyl relaxation experiments and isothermal titration calorimetry revealed that conformational entropy contributes to affinity. In addition to motions relevant to thermodynamics, slower microsecond to millisecond switching was prevalent in the QM PDZ ligand-binding site consistent with a role in ligand specificity. Our data indicate that conformational dynamics plays distinct and fundamental roles in tuning the affinity (conformational entropy) and specificity (excited-state conformations) of molecular interactions. More broadly, our results have important implications for the evolution, regulation, and design of protein-ligand interactions.

The flexibility and dynamics of protein disulfide isomerase

We have studied the mobility of the multidomain folding catalyst, protein disulfide isomerase (PDI), by a coarse‐graining approach based on flexibility. We analyze our simulations of yeast PDI (yPDI) using measures of backbone movement, relative positions and orientations of domains, and distances between functional sites. We find that there is interdomain flexibility at every interdomain junction but these show very different characteristics. The extent of interdomain flexibility is such that yPDI's two active sites can approach much more closely than is found in crystal structures—and indeed hinge motion to bring these sites into proximity is the lowest energy normal mode of motion of the protein. The flexibility predicted for yPDI (based on one structure) includes the other known conformation of yPDI and is consistent with (i) the mobility observed experimentally for mammalian PDI and (ii) molecular dynamics. We also observe intradomain flexibility and clear differences between the domains in their propensity for internal motion. Our results suggest that PDI flexibility enables it to interact with many different partner molecules of widely different sizes and shapes, and highlights considerable similarities of yPDI and mammalian PDI. Proteins 2016; 84:1776–1785. © 2016 Wiley Periodicals, Inc.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°_bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S²_axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°_conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S²_axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°_conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

Evolution of a Protein Interaction Domain Family by Tuning Conformational Flexibility

Conformational flexibility allows proteins to adopt multiple functionally important conformations but can also lead to nonfunctional structures. We analyzed the dynamic behavior of the enzyme guanylate kinase as it evolved into the GK protein interaction domain (GK_PID) to investigate the role of flexibility in the evolution of new protein functions. We found that the ancestral enzyme is very flexible, allowing it to adopt open conformations that can bind nucleotide and closed ones that enable catalysis of phosphotransfer from ATP to GMP. Historical mutations that converted the GK from an enzyme to a protein interaction domain dramatically reduce flexibility, predominantly by inhibiting rotations of the protein backbone that are coupled to the global closing motion. Removing flexibility prevents adoption of conformations that cannot fit the protein partner in the binding site. Our results highlight the importance of mutations that optimize protein conformational flexibility with function during evolution.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Accurate determination of rates from non-uniformly sampled relaxation data

The application of non-uniform sampling (NUS) to relaxation experiments traditionally used to characterize the fast internal motion of proteins is quantitatively examined. Experimentally acquired Poisson-gap sampled data reconstructed with iterative soft thresholding are compared to regular sequentially sampled (RSS) data. Using ubiquitin as a model system, it is shown that 25 % sampling is sufficient for the determination of quantitatively accurate relaxation rates. When the sampling density is fixed at 25 %, the accuracy of rates is shown to increase sharply with the total number of sampled points until eventually converging near the inherent reproducibility of the experiment. Perhaps contrary to some expectations, it is found that accurate peak height reconstruction is not required for the determination of accurate rates. Instead, inaccuracies in rates arise from inconsistencies in reconstruction across the relaxation series that primarily manifest as a non-linearity in the recovered peak height. This indicates that the performance of an NUS relaxation experiment cannot be predicted from comparison of peak heights using a single RSS reference spectrum. The generality of these findings was assessed using three alternative reconstruction algorithms, eight different relaxation measurements, and three additional proteins that exhibit varying degrees of spectral complexity. From these data, it is revealed that non-linearity in peak height reconstruction across the relaxation series is strongly correlated with errors in NUS-derived relaxation rates. Importantly, it is shown that this correlation can be exploited to reliably predict the performance of an NUS-relaxation experiment by using three or more RSS reference planes from the relaxation series. The RSS reference time points can also serve to provide estimates of the uncertainty of the sampled intensity, which for a typical relaxation times series incurs no penalty in total acquisition time.