Dynamic allostery: linkers are not merely flexible

Fold combinations in multi-domain proteins

Domain-domain interactions in multi-domain proteins play an important role in the combined function of individual domains for the overall biological activity of the protein. The functions of the tethered domains are often coupled and hence, limited numbers of domain architectures with defined folds are known in nature. Therefore, it is of interest to document the available fold-fold combinations and their preference in multi-domain proteins. Hence, we analyzed all multi-domain proteins with known structures in the protein databank and observed that only about 860 fold-fold combinations are present among them. Analyses of multi-domain proteins represented in sequence database result in recognition of 29,860 fold-fold combinations and it accounts for only 2.8% of the theoretically possible 1,036,080 (1439C2) fold-fold combinations. The observed preference for fold-fold combinations in multi-domain proteins is interesting in the context of multiple functions through structural adaptation by gene fusion.

Tracing the origin and evolution of pseudokinases across the tree of life

Protein phosphorylation by eukaryotic protein kinases (ePKs) is a fundamental mechanism of cell signaling in all organisms. In model vertebrates, ~10% of ePKs are classified as pseudokinases, which have amino acid changes within the catalytic machinery of the kinase domain that distinguish them from their canonical kinase counterparts. However, pseudokinases still regulate various signaling pathways, usually doing so in the absence of their own catalytic output. To investigate the prevalence, evolutionary relationships, and biological diversity of these pseudoenzymes, we performed a comprehensive analysis of putative pseudokinase sequences in available eukaryotic, bacterial, and archaeal proteomes. We found that pseudokinases are present across all domains of life, and we classified nearly 30,000 eukaryotic, 1500 bacterial, and 20 archaeal pseudokinase sequences into 86 pseudokinase families, including ~30 families that were previously unknown. We uncovered a rich variety of pseudokinases with notable expansions not only in animals but also in plants, fungi, and bacteria, where pseudokinases have previously received cursory attention. These expansions are accompanied by domain shuffling, which suggests roles for pseudokinases in plant innate immunity, plant-fungal interactions, and bacterial signaling. Mechanistically, the ancestral kinase fold has diverged in many distinct ways through the enrichment of unique sequence motifs to generate new families of pseudokinases in which the kinase domain is repurposed for noncanonical nucleotide binding or to stabilize unique, inactive kinase conformations. We further provide a collection of annotated pseudokinase sequences in the Protein Kinase Ontology (ProKinO) as a new mineable resource for the signaling community.

Antigen binding allosterically promotes Fc receptor recognition

A key question in immunology is whether antigen recognition and Fc receptor (FcR) binding are allosterically linked. This question is also relevant for therapeutic antibody design. Antibody Fab and Fc domains are connected by flexible unstructured hinge region. Fc chains have conserved glycosylation sites at Asn297, with each conjugated to a core heptasaccharide and forming biantennary Fc glycan. The glycans modulate the Fc conformations and functions. It is well known that the antibody Fab and Fc domains and glycan affect antibody activity, but whether these elements act independently or synergistically is still uncertain. We simulated four antibody complexes: free antibody, antigen-bound antibody, FcR-bound antibody, and an antigen-antibody-FcR complex. We found that, in the antibody's "T/Y" conformation, the glycans, and the Fc domain all respond to antigen binding, with the antibody population shifting to two dominant clusters, both with the Fc-receptor binding site open. The simulations reveal that the Fc-glycan-receptor complexes also segregate into two conformational clusters, one corresponding to the antigen-free antibody-FcR baseline binding, and the other with an antigen-enhanced antibody-FcR interaction. Our study confirmed allosteric communications in antibody-antigen recognition and following FcR activation. Even though we observed allosteric communications through the IgG domains, the most important mechanism that we observed is the communication via population shift, stimulated by antigen binding and propagating to influence FcR recognition.

The missing links within troponin

The cardiac contraction-relaxation cycle is controlled by a sophisticated set of machinery. Of particular interest, is the revelation that allosteric networks transmit effects of binding at one site to influence troponin complex dynamics and structural-mediated signaling in often distal, functional sites in the myofilament. Our recent observations provide compelling evidence that allostery can explain the function of large-scale macromolecular events. Here we elaborate on our recent findings of interdomain communication within troponin C, using cutting-edge structural biology approaches, and highlight the importance of unveiling the unknown, distant communication networks within this system to obtain more comprehensive knowledge of how allostery impacts cardiac physiology and disease.

Autoinhibition in Ras effectors Raf, PI3Kα, and RASSF5: a comprehensive review underscoring the challenges in pharmacological intervention

Autoinhibition is an effective mechanism that guards proteins against spurious activation. Despite its ubiquity, the distinct organizations of the autoinhibited states and their release mechanisms differ. Signaling is most responsive to the cell environment only if a small shift in the equilibrium is required to switch the system from an inactive (occluded) to an active (exposed) state. Ras signaling follows this paradigm. This underscores the challenge in pharmacological intervention to exploit and enhance autoinhibited states. Here, we review autoinhibition and release mechanisms at the membrane focusing on three representative Ras effectors, Raf protein kinase, PI3Kα lipid kinase, and NORE1A (RASSF5) tumor suppressor, and point to the ramifications to drug discovery. We further touch on Ras upstream and downstream signaling, Ras activation, and the Ras superfamily in this light, altogether providing a broad outlook of the principles and complexities of autoinhibition.

Local and global anatomy of antibody-protein antigen recognition

Deciphering antibody-protein antigen recognition is of fundamental and practical significance. We constructed an antibody structural dataset, partitioned it into human and murine subgroups, and compared it with nonantibody protein-protein complexes. We investigated the physicochemical properties of regions on and away from the antibody-antigen interfaces, including net charge, overall antibody charge distributions, and their potential role in antigen interaction. We observed that amino acid preference in antibody-protein antigen recognition is entropy driven, with residues having low side-chain entropy appearing to compensate for the high backbone entropy in interaction with protein antigens. Antibodies prefer charged and polar antigen residues and bridging water molecules. They also prefer positive net charge, presumably to promote interaction with negatively charged protein antigens, which are common in proteomes. Antibody-antigen interfaces have large percentages of Tyr, Ser, and Asp, but little Lys. Electrostatic and hydrophobic interactions in the Ag binding sites might be coupled with Fab domains through organized charge and residue distributions away from the binding interfaces. Here we describe some features of antibody-antigen interfaces and of Fab domains as compared with nonantibody protein-protein interactions. The distributions of interface residues in human and murine antibodies do not differ significantly. Overall, our results provide not only a local but also a global anatomy of antibody structures.

Intramolecular Fuzzy Interactions Involving Intrinsically Disordered Domains

Structural disorder is an essential ingredient for function in many proteins and protein complexes. Fuzzy complexes describe the many instances where disorder is maintained as a critical element of protein interactions. In this minireview we discuss how intramolecular fuzzy interactions function in signaling complexes. Focussing on the Src family of kinases, we argue that the intrinsically disordered domains that are unique for each of the family members and display a clear fingerprint of long range interactions in Src, might have critical roles as functional sensor or effectors and mediate allosteric communication via fuzzy interactions.

Exploring modular allostery via interchangeable regulatory domains

Most proteins comprise two or more domains from a limited suite of protein families. These domains are often rearranged in various combinations through gene fusion events to evolve new protein functions, including the acquisition of protein allostery through the incorporation of regulatory domains. The enzyme 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAH7PS) is the first enzyme of aromatic amino acid biosynthesis and displays a diverse range of allosteric mechanisms. DAH7PSs adopt a common architecture with a shared (β/α)₈ catalytic domain which can be attached to an ACT-like or a chorismate mutase regulatory domain that operates via distinct mechanisms. These respective domains confer allosteric regulation by controlling DAH7PS function in response to ligand Tyr or prephenate. Starting with contemporary DAH7PS proteins, two protein chimeras were created, with interchanged regulatory domains. Both engineered proteins were catalytically active and delivered new functional allostery with switched ligand specificity and allosteric mechanisms delivered by their nonhomologous regulatory domains. This interchangeability of protein domains represents an efficient method not only to engineer allostery in multidomain proteins but to create a new bifunctional enzyme.

Coupled regulation by the juxtamembrane and sterile α motif (SAM) linker is a hallmark of ephrin tyrosine kinase evolution

Ephrin (Eph) receptor tyrosine kinases have evolutionarily diverged from other tyrosine kinases to respond to specific activation and regulatory signals that require close coupling of kinase catalytic and regulatory functions. However, the evolutionary basis for such functional coupling is not fully understood. We employed an evolutionary systems approach involving statistical mining of large sequence and structural data sets to define the hallmarks of Eph kinase evolution and functional specialization. We found that some of the most distinguishing Eph-specific residues structurally tether the flanking juxtamembrane and sterile α motif (SAM) linker regions to the kinase domain, and substitutions of these residues in EphA3 resulted in faster kinase activation. We report for the first time that the SAM domain linker is functionally coupled to the juxtamembrane through co-conserved residues in the kinase domain and that together these residues provide a structural framework for coupling catalytic and regulatory functions. The unique organization of Eph-specific tethering networks and the identification of other Eph-specific sequence features of unknown functions provide new hypotheses for future functional studies and new clues to disease mutations altering Eph kinase-specific functions.

Same but not alike: Structure, flexibility and energetics of domains in multi-domain proteins are influenced by the presence of other domains

The majority of the proteins encoded in the genomes of eukaryotes contain more than one domain. Reasons for high prevalence of multi-domain proteins in various organisms have been attributed to higher stability and functional and folding advantages over single-domain proteins. Despite these advantages, many proteins are composed of only one domain while their homologous domains are part of multi-domain proteins. In the study presented here, differences in the properties of protein domains in single-domain and multi-domain systems and their influence on functions are discussed. We studied 20 pairs of identical protein domains, which were crystallized in two forms (a) tethered to other proteins domains and (b) tethered to fewer protein domains than (a) or not tethered to any protein domain. Results suggest that tethering of domains in multi-domain proteins influences the structural, dynamic and energetic properties of the constituent protein domains. 50% of the protein domain pairs show significant structural deviations while 90% of the protein domain pairs show differences in dynamics and 12% of the residues show differences in the energetics. To gain further insights on the influence of tethering on the function of the domains, 4 pairs of homologous protein domains, where one of them is a full-length single-domain protein and the other protein domain is a part of a multi-domain protein, were studied. Analyses showed that identical and structurally equivalent functional residues show differential dynamics in homologous protein domains; though comparable dynamics between in-silico generated chimera protein and multi-domain proteins were observed. From these observations, the differences observed in the functions of homologous proteins could be attributed to the presence of tethered domain. Overall, we conclude that tethered domains in multi-domain proteins not only provide stability or folding advantages but also influence pathways resulting in differences in function or regulatory properties.

High prevalence of multi-domain proteins in proteomes has been attributed to higher stability and functional and folding advantages of the multi-domain proteins. Influence of tethering of domains on the overall properties of proteins has been well studied but its influence on the properties of the constituent domains is largely unaddressed. Here, we investigate the influence of tethering of domains in multi-domain proteins on the structural, dynamics and energetics properties of the constituent domains and its implications on the functions of proteins. To this end, comparative analyses were carried out for identical protein domains crystallized in tethered and untethered forms. Also, comparative analyses of single-domain proteins and their homologous multi-domain proteins were performed. The analyses suggest that tethering influences the structural, dynamic and energetic properties of constituent protein domains. Our observations hint at regulation of protein domains by tethered domains in multi-domain systems, which may manifest at the differential function observed between single-domain and homologous multi-domain proteins.

Combinatorial Assembly of Lumitoxins

Ion channels are among the most important proteins in neuroscience and serve as drug targets for many brain disorders. During development, learning, disease progression, and other processes, the activity levels of specific ion channels are tuned in a cell-type specific manner. However, it is difficult to assess how cell-specific changes in ion channel activity alter emergent brain functions. We have developed a protein architecture for fully genetically encoded light-activated modulation of endogenous ion channel activity. Fusing a genetically encoded photoswitch and an ion channel-modulating peptide toxin in a computationally designed fashion, this reagent, which we call Lumitoxins, can mediate light-modulation of specific endogenous ion channel activities in targeted cells. The modular lumitoxin architecture may be useful in a diversity of neuroscience tools. Here, we delineate how to construct lumitoxin genes from synthesized components, and provide a general outline for how to test their function in mammalian cell culture.

Bioinformatic Analysis Reveals Conservation of Intrinsic Disorder in the Linker Sequences of Prokaryotic Dual-family Immunophilin Chaperones

The two classical immunophilin families, found essentially in all living cells, are: cyclophilin (CYN) and FK506-binding protein (FKBP). We previously reported a novel class of immunophilins that are natural chimera of these two, which we named dual-family immunophilin (DFI). The DFIs were found in either of two conformations: CYN-linker-FKBP (CFBP) or FKBP-3TPR-CYN (FCBP). While the 3TPR domain can serve as a flexible linker between the FKBP and CYN modules in the FCBP-type DFI, the linker sequences in the CFBP-type DFIs are relatively short, diverse in sequence, and contain no discernible motif or signature. Here, I present several lines of computational evidence that, regardless of their primary structure, these CFBP linkers are intrinsically disordered. This report provides the first molecular foundation for the model that the CFBP linker acts as an unstructured, flexible loop, allowing the two flanking chaperone modules function independently while linked in cis, likely to assist in the folding of multisubunit client complexes.

•

Dual-family immunophilins (DFIs) are recently discovered chimeric chaperones.

•

Bacterial DFIs possess the structure CYN-linker-FKBP, abbreviated as CFBP.

•

The linker sequences in CFBP are diverse but possess intrinsic disorder (ID).

•

The large CFBP linker with ID is a novel discovery in prokaryotes.

•

ID befits the linker's role as a flexible connector between two chaperone modules.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Intrinsic protein disorder in oncogenic KRAS signaling

How Ras, and in particular its most abundant oncogenic isoform K-Ras4B, is activated and signals in proliferating cells, poses some of the most challenging questions in cancer cell biology. In this paper, we ask how intrinsically disordered regions in K-Ras4B and its effectors help promote proliferative signaling. Conformational disorder allows spanning long distances, supports hinge motions, promotes anchoring in membranes, permits segments to fulfil multiple roles, and broadly is crucial for activation mechanisms and intensified oncogenic signaling. Here, we provide an overview illustrating some of the key mechanisms through which conformational disorder can promote oncogenesis, with K-Ras4B signaling serving as an example. We discuss (1) GTP-bound KRas4B activation through membrane attachment; (2) how farnesylation and palmitoylation can promote isoform functional specificity; (3) calmodulin binding and PI3K activation; (4) how Ras activates its RASSF5 cofactor, thereby stimulating signaling of the Hippo pathway and repressing proliferation; and (5) how intrinsically disordered segments in Raf help its attachment to the membrane and activation. Collectively, we provide the first inclusive review of the roles of intrinsic protein disorder in oncogenic Ras-driven signaling. We believe that a broad picture helps to grasp and formulate key mechanisms in Ras cancer biology and assists in therapeutic intervention.

Stiffness of the C-terminal disordered linker affects the geometry of the active site in endoglucanase Cel8A

Cellulosomes are complex multi-enzyme machineries which efficiently degrade plant cell-wall polysaccharides. The multiple domains of the cellulosome proteins are often tethered together by intrinsically disordered regions. The properties and functions of these disordered linkers are not well understood. In this work, we study endoglucanase Cel8A, which is a relevant enzymatic component of the cellulosomes of Clostridium thermocellum. We use both all-atom and coarse-grained simulations to investigate how the conformations of the catalytic domain of Cel8A are affected by the disordered linker at its C terminus. We find that when the endoglucanase is bound to its substrate, the effective stiffness of the linker can influence the distances between groups of amino-acid residues throughout the entire enzymatic domain. In particular, variations in the linker stiffness can lead to small changes in the geometry of the active-site cleft. We suggest that such geometrical changes may have an effect on the catalytic activity of the enzyme.

Influence of Linker Flexibility on the Binding Affinity of Bidentate Binders

The design of responsive nanosensors typically relies on the availability of probes capable of capturing their target with high affinity and specificity. This can be achieved by coupling two or more binding units through a linker. In this work, we study the dependence on the binder architecture of the binding affinity between a target molecule and a semirigid bidentate binder. Using two different binder architectures, central-rigid and extreme-rigid, and modifying the length and the flexibility degree of the linker we generated 153 different architectures. We computed their dissociation free energies by means of Monte Carlo simulations and thermodynamic integration. We found that central-rigid bidentate binders are a poor choice, as they dissociate more easily than analogous fully flexible bidentate binders. On the other hand, molecular architectures presenting extreme-rigid units were shown effective for a wide range of set-ups.

Cullin neddylation may allosterically tune polyubiquitin chain length and topology

Conjugation of Nedd8 (neddylation) to Cullins (Cul) in Cul-RING E3 ligases (CRLs) stimulates ubiquitination and polyubiquitination of protein substrates. CRL is made up of two Cul-flanked arms: one consists of the substrate-binding and adaptor proteins and the other consists of E2 and Ring-box protein (Rbx). Polyubiquitin chain length and topology determine the substrate fate. Here, we ask how polyubiquitin chains are accommodated in the limited space available between the two arms and what determines the polyubiquitin linkage topology. We focus on Cul5 and Rbx1 in three states: before Cul5 neddylation (closed state), after neddylation (open state), and after deneddylation, exploiting molecular dynamics simulations and the Gaussian Network Model. We observe that regulation of substrate ubiquitination and polyubiquitination takes place through Rbx1 rotations, which are controlled by Nedd8-Rbx1 allosteric communication. Allosteric propagation proceeds from Nedd8 via Cul5 dynamic hinges and hydrogen bonds between the C-terminal domain of Cul5 (Cul5CTD) and Rbx1 (Cul5CTD residues R538/R569 and Rbx1 residue E67, or Cul5CTD E474/E478/N491 and Rbx1 K105). Importantly, at each ubiquitination step (homogeneous or heterogeneous, linear or branched), the polyubiquitin linkages fit into the distances between the two arms, and these match the inherent CRL conformational tendencies. Hinge sites may constitute drug targets.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Allosteric therapies for lung cancer

Allostery is a regulation at a distance by conveying information from one site to another and an intrinsic property of dynamic proteins. Allostery plays an essential role in receptor trafficking, signal transmission, controlled catalysis, gene turn on/off, or cell apoptosis. Allosteric mutations are considered as one of causes responsible for cancer development, leading to "allosteric diseases" by stabilizing an active or inactive conformation or changing the dynamic distribution of preexisting propagation pathways. The present article mainly focuses on the potential of allosteric therapies for lung cancer. Allosteric drugs may have several advantages over traditional drugs. The epidermal growth factor receptor mutations and signaling pathways downstream (such as PI3K/AKT/mTOR and RAS/RAF/MEK/ERK pathways) were suggested to play a key role in lung cancer and considered as targets of allosteric therapy. Some allosteric inhibitors for lung cancer-specific targets and a series of preclinical trials of allosteric inhibitors for lung cancer have been developed and reported. We expect that allosteric therapies will gain more attentions to develop combinatorial strategies for lung cancer and metastasis.

Nonlinear backbone torsional pair correlations in proteins

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

Oncogenic Mutations Differentially Affect Bax Monomer, Dimer, and Oligomeric Pore Formation in the Membrane

Dysfunction of Bax, a pro-apoptotic regulator of cellular metabolism is implicated in neurodegenerative diseases and cancer. We have constructed the first atomistic models of the Bax oligomeric pore consisting with experimental residue-residue distances. The models are stable, capturing well double electron-electron resonance (DEER) spectroscopy measurements and provide structural details in line with the DEER data. Comparison with the latest experimental results revealed that our models agree well with both Bax and Bak pores, pointed to a converged structural arrangement for Bax and Bak pore formation. Using multi-scale molecular dynamics simulations, we probed mutational effects on Bax transformation from monomer → dimer → membrane pore formation at atomic resolution. We observe that two cancer-related mutations, G40E and S118I, allosterically destabilize the monomer and stabilize an off-pathway swapped dimer, preventing productive pore formation. This observation suggests a mechanism whereby the mutations may work mainly by over-stabilizing the monomer → dimer transformation toward an unproductive off-pathway swapped-dimer state. Our observations point to misfolded Bax states, shedding light on the molecular mechanism of Bax mutation-elicited cancer. Most importantly, the structure of the Bax pore facilitates future study of releases cytochrome C in atomic detail.

Investigating the Role of Large-Scale Domain Dynamics in Protein-Protein Interactions

Intrinsically disordered linkers provide multi-domain proteins with degrees of conformational freedom that are often essential for function. These highly dynamic assemblies represent a significant fraction of all proteomes, and deciphering the physical basis of their interactions represents a considerable challenge. Here we describe the difficulties associated with mapping the large-scale domain dynamics and describe two recent examples where solution state methods, in particular NMR spectroscopy, are used to investigate conformational exchange on very different timescales.

Computational approaches to investigating allostery

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

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

A review of multi-domain and flexible molecular chaperones studies by small-angle X-ray scattering

Intrinsic flexibility is closely related to protein function, and a plethora of important regulatory proteins have been found to be flexible, multi-domain or even intrinsically disordered. On the one hand, understanding such systems depends on how these proteins behave in solution. On the other, small-angle X-ray scattering (SAXS) is a technique that fulfills the requirements to study protein structure and dynamics relatively quickly with few experimental limitations. Molecular chaperones from Hsp70 and Hsp90 families are multi-domain proteins containing flexible and/or disordered regions that play central roles in cellular proteostasis. Here, we review the structure and function of these proteins by SAXS. Our general approach includes the use of SAXS data to determine size and shape parameters, as well as protein shape reconstruction and their validation by using accessory biophysical tools. Some remarkable examples are presented that exemplify the potential of the SAXS technique. Protein structure can be determined in solution even at limiting protein concentrations (for example, human mortalin, a mitochondrial Hsp70 chaperone). The protein organization, flexibility and function (for example, the J-protein co-chaperones), oligomeric status, domain organization, and flexibility (for the Hsp90 chaperone and the Hip and Hep1 co-chaperones) may also be determined. Lastly, the shape, structural conservation, and protein dynamics (for the Hsp90 chaperone and both p23 and Aha1 co-chaperones) may be studied by SAXS. We believe this review will enhance the application of the SAXS technique to the study of the molecular chaperones.

Evolution of the concept of conformational dynamics of enzyme functions over half of a century: A personal view

To most physicists, it was always evident that conformational fluctuation is an inherent property of all molecules. Its existence in proteins was mentioned first by Linderström-Lang and Schellman in 1959 based on their hydrogen-deuterium exchange experiments. The "induced fit" mechanism to explain ligand-induced conformational changes was suggested by Koshland in 1958. Straub combined these two concepts in his "fluctuation fit" theory in 1964. The era of protein X-ray crystallography imposed a static view of protein structures. With proteins becoming accessible to NMR analysis, conformational dynamics could be mapped, and a new wave of dynamic interpretation of enzymatic catalysis and molecular recognition appeared. Energy landscapes, energy funnels, conformational selection, conformational distribution shifts are now frequent terms in interpreting biomolecular recognition and enzymatic catalysis. All these interpretations are based on the concept that evolution uses the conformational fluctuations of enzymes to develop efficient and dynamic catalytic machines. In a resurrection of the original "fluctuation fit" concept, it is generally recognized now that spatial and temporal events of catalysis are equally important to describe its mechanism. This special issue, dedicated to the memory of Henryk Eisenberg, prompted us to look back at the last 50 years of development of a concept that-like other important concepts-appeared, evolved and became accepted during the period covered by the scientific lifespan of Henryk.

Structure and Function of Fusicoccadiene Synthase, a Hexameric Bifunctional Diterpene Synthase

Fusicoccin A is a diterpene glucoside phytotoxin generated by the fungal pathogen Phomopsis amygdali that causes the plant disease constriction canker, first discovered in New Jersey peach orchards in the 1930s. Fusicoccin A is also an emerging new lead in cancer chemotherapy. The hydrocarbon precursor of fusicoccin A is the tricyclic diterpene fusicoccadiene, which is generated by a bifunctional terpenoid synthase. Here, we report X-ray crystal structures of the individual catalytic domains of fusicoccadiene synthase: the C-terminal domain is a chain elongation enzyme that generates geranylgeranyl diphosphate, and the N-terminal domain catalyzes the cyclization of geranylgeranyl diphosphate to form fusicoccadiene. Crystal structures of each domain complexed with bisphosphonate substrate analogues suggest that three metal ions and three positively charged amino acid side chains trigger substrate ionization in each active site. While in vitro incubations reveal that the cyclase domain can utilize farnesyl diphosphate and geranyl diphosphate as surrogate substrates, these shorter isoprenoid diphosphates are mainly converted into acyclic alcohol or hydrocarbon products. Gel filtration chromatography and analytical ultracentrifugation experiments indicate that full-length fusicoccadiene synthase adopts hexameric quaternary structure, and small-angle X-ray scattering data yield a well-defined molecular envelope illustrating a plausible model for hexamer assembly.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

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

Allosteric Partial Inhibition of Monomeric Proteases. Sulfated Coumarins Induce Regulation, not just Inhibition, of Thrombin

Allosteric partial inhibition of soluble, monomeric proteases can offer major regulatory advantages, but remains a concept on paper to date; although it has been routinely documented for receptors and oligomeric proteins. Thrombin, a key protease of the coagulation cascade, displays significant conformational plasticity, which presents an attractive opportunity to discover small molecule probes that induce sub-maximal allosteric inhibition. We synthesized a focused library of some 36 sulfated coumarins to discover two agents that display sub-maximal efficacy (~50%), high potency (<500 nM) and high selectivity for thrombin (>150-fold). Michaelis-Menten, competitive inhibition, and site-directed mutagenesis studies identified exosite 2 as the site of binding for the most potent sulfated coumarin. Stern-Volmer quenching of active site-labeled fluorophore suggested that the allosteric regulators induce intermediate structural changes in the active site as compared to those that display ~80–100% efficacy. Antithrombin inactivation of thrombin was impaired in the presence of the sulfated coumarins suggesting that allosteric partial inhibition arises from catalytic dysfunction of the active site. Overall, sulfated coumarins represent first-in-class, sub-maximal inhibitors of thrombin. The probes establish the concept of allosteric partial inhibition of soluble, monomeric proteins. This concept may lead to a new class of anticoagulants that are completely devoid of bleeding.

An atomistic view of Hsp70 allosteric crosstalk: from the nucleotide to the substrate binding domain and back

The Hsp70 is an allosterically regulated family of molecular chaperones. They consist of two structural domains, NBD and SBD, connected by a flexible linker. ATP hydrolysis at the NBD modulates substrate recognition at the SBD, while peptide binding at the SBD enhances ATP hydrolysis. In this study we apply Molecular Dynamics (MD) to elucidate the molecular determinants underlying the allosteric communication from the NBD to the SBD and back. We observe that local structural and dynamical modulation can be coupled to large-scale rearrangements, and that different combinations of ligands at NBD and SBD differently affect the SBD domain mobility. Substituting ADP with ATP in the NBD induces specific structural changes involving the linker and the two NBD lobes. Also, a SBD-bound peptide drives the linker docking by increasing the local dynamical coordination of its C-terminal end: a partially docked DnaK structure is achieved by combining ATP in the NBD and peptide in the SBD. We propose that the MD-based analysis of the inter domain dynamics and structure modulation could be used as a tool to computationally predict the allosteric behaviour and functional response of Hsp70 upon introducing mutations or binding small molecules, with potential applications for drug discovery.

Functional Diversification after Gene Duplication: Paralog Specific Regions of Structural Disorder and Phosphorylation in p53, p63, and p73

Conformational and functional flexibility promote protein evolvability. High evolvability allows related proteins to functionally diverge and perhaps to neostructuralize. p53 is a multifunctional protein frequently referred to as the Guardian of the Genome–a hub for e.g. incoming and outgoing signals in apoptosis and DNA repair. p53 has been found to be structurally disordered, an extreme form of conformational flexibility. Here, p53, and its paralogs p63 and p73, were studied for further insights into the evolutionary dynamics of structural disorder, secondary structure, and phosphorylation. This study is focused on the post gene duplication phase for the p53 family in vertebrates, but also visits the origin of the protein family and the early domain loss and gain events. Functional divergence, measured by rapid evolutionary dynamics of protein domains, structural properties, and phosphorylation propensity, is inferred across vertebrate p53 proteins, in p63 and p73 from fish, and between the three paralogs. In particular, structurally disordered regions are redistributed among paralogs, but within clades redistribution of structural disorder also appears to be an ongoing process. Despite its deemed importance as the Guardian of the Genome, p53 is indeed a protein with high evolvability as seen not only in rearranged structural disorder, but also in fluctuating domain sequence signatures among lineages.

Discovery of Novel Allosteric Eg5 Inhibitors Through Structure-Based Virtual Screening

Mitotic kinesin Eg5 is an attractive anticancer drug target. Discovery of Eg5 inhibitors has been focused on targeting the 'monastrol-binding site'. However, acquired drug resistance has been reported for such inhibitors. Therefore, identifying new Eg5 inhibitors which function through a different mechanism(s) could complement current drug candidates and improve drug efficacy. In this study, we explored a novel allosteric site of Eg5 and identified new Eg5 inhibitors through structure-based virtual screening. Experiments with the saturation-transfer difference NMR demonstrated that the identified Eg5 inhibitor SRI35566 binds directly to Eg5 without involving microtubules. Moreover, SRI35566 and its two analogs significantly induced monopolar spindle formation in colorectal cancer HCT116 cells and suppressed cancer cell viability and colony formation. Together, our findings reveal a new allosteric regulation mechanism of Eg5 and a novel drug targeting site for cancer therapy.

Conformational Dynamics and Allostery in Pyruvate Kinase

Pyruvate kinase catalyzes the final step in glycolysis and is allosterically regulated to control flux through the pathway. Two models are proposed to explain how Escherichia coli pyruvate kinase type 1 is allosterically regulated: the "domain rotation model" suggests that both the domains within the monomer and the monomers within the tetramer reorient with respect to one another; the "rigid body reorientation model" proposes only a reorientation of the monomers within the tetramer causing rigidification of the active site. To test these hypotheses and elucidate the conformational and dynamic changes that drive allostery, we performed time-resolved electrospray ionization mass spectrometry coupled to hydrogen-deuterium exchange studies followed by mutagenic analysis to test the activation mechanism. Global exchange experiments, supported by thermostability studies, demonstrate that fructose 1,6-bisphosphate binding to the allosteric domain causes a shift toward a globally more dynamic ensemble of conformations. Mapping deuterium exchange to peptides within the enzyme highlight site-specific regions with altered conformational dynamics, many of which increase in conformational flexibility. Based upon these and mutagenic studies, we propose an allosteric mechanism whereby the binding of fructose 1,6-bisphosphate destabilizes an α-helix that bridges the allosteric and active site domains within the monomeric unit. This destabilizes the β-strands within the (β/α)8-barrel domain and the linked active site loops that are responsible for substrate binding. Our data are consistent with the domain rotation model but inconsistent with the rigid body reorientation model given the increased flexibility at the interdomain interface, and we can for the first time explain how fructose 1,6-bisphosphate affects the active site.

Conserved linker regions and their regulation determine multiple chromatin-binding modes of UHRF1

Ubiquitin-like with PHD and RING Finger domains 1 (UHRF1) is an important nuclear protein that is mutated and aberrantly expressed in many tumors. The protein integrates different chromatin modifications and is essential for their maintenance throughout the cell cycle. Separate chromatin-binding modules of UHRF1 have been studied on a functional and structural level. The unmodified N-terminus of histone H3 is recognized by a PHD domain, while a TTD domain specifically interacts with histone H3 Lysine 9 trimethylation. A SRA region binds hemimethylatd DNA. Emerging evidence indicates that the modules of UHRF1 do not act independently of each other but establish complex modes of interaction with patterns of chromatin modifications. This multivalent readout is regulated by allosteric binding of phosphatidylinositol 5-phosphate to a region outside the PHD, TTD and SRA domains as well as by phosphorylation of one of the linker regions connecting these modules. Here, we summarize the current knowledge on UHRF1 chromatin interaction and introduce a novel model of conformational transitions of the protein that are directed by the flexible and highly charged linker regions. We propose that these are essential in setting up defined structural states of the protein where different domains or combinations thereof are available for binding chromatin modifications or are prevented from doing so. Lastly, we suggest that controlled tuning of intramolecular linker interactions by ligands and posttranslational modifications establishes a rational framework for comprehending UHRF1 regulation and putatively the working mode of other chromatin factors in different physiological contexts.

Tandem Domains with Tuned Interactions Are a Powerful Biological Design Principle

Allosteric effects of mutations, ligand binding, or post-translational modifications on protein function occur through changes to the protein's shape, or conformation. In a cell, there are many copies of the same protein, all experiencing these perturbations in a dynamic fashion and fluctuating through different conformations and activity states. According to the "conformational selection and population shift" theory, ligand binding selects a particular conformation. This perturbs the ensemble and induces a population shift. In a new PLOS Biology paper, Melacini and colleagues describe a novel model of protein regulation, the "Double-Conformational Selection Model", which demonstrates how two tandem ligand-binding domains interact to regulate protein function. Here we explain how tandem domains with tuned interactions-but not single domains-can provide a blueprint for sensitive activation sensors within a narrow window of ligand concentration, thereby promoting signaling control.

Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins

Hsp70 and Hsp110 chaperones play an important role in regulating cellular processes that involve protein folding and stabilization, which are essential for the integrity of signaling networks. Although many aspects of allosteric regulatory mechanisms in Hsp70 and Hsp110 chaperones have been extensively studied and significantly advanced in recent experimental studies, the atomistic picture of signal propagation and energetics of dynamics-based communication still remain unresolved. In this work, we have combined molecular dynamics simulations and protein stability analysis of the chaperone structures with the network modeling of residue interaction networks to characterize molecular determinants of allosteric mechanisms. We have shown that allosteric mechanisms of Hsp70 and Hsp110 chaperones may be primarily determined by nucleotide-induced redistribution of local conformational ensembles in the inter-domain regions and the substrate binding domain. Conformational dynamics and energetics of the peptide substrate binding with the Hsp70 structures has been analyzed using free energy calculations, revealing allosteric hotspots that control negative cooperativity between regulatory sites. The results have indicated that cooperative interactions may promote a population-shift mechanism in Hsp70, in which functional residues are organized in a broad and robust allosteric network that can link the nucleotide-binding site and the substrate-binding regions. A smaller allosteric network in Hsp110 structures may elicit an entropy-driven allostery that occurs in the absence of global structural changes. We have found that global mediating residues with high network centrality may be organized in stable local communities that are indispensable for structural stability and efficient allosteric communications. The network-centric analysis of allosteric interactions has also established that centrality of functional residues could correlate with their sensitivity to mutations across diverse chaperone functions. This study reconciles a wide spectrum of structural and functional experiments by demonstrating how integration of molecular simulations and network-centric modeling may explain thermodynamic and mechanistic aspects of allosteric regulation in chaperones.

Dynamic allostery governs cyclophilin A-HIV capsid interplay

Host factor protein Cyclophilin A (CypA) regulates HIV-1 viral infectivity through direct interactions with the viral capsid, by an unknown mechanism. CypA can either promote or inhibit viral infection, depending on host cell type and HIV-1 capsid (CA) protein sequence. We have examined the role of conformational dynamics on the nanosecond to millisecond timescale in HIV-1 CA assemblies in the escape from CypA dependence, by magic-angle spinning (MAS) NMR and molecular dynamics (MD). Through the analysis of backbone (1)H-(15)N and (1)H-(13)C dipolar tensors and peak intensities from 3D MAS NMR spectra of wild-type and the A92E and G94D CypA escape mutants, we demonstrate that assembled CA is dynamic, particularly in loop regions. The CypA loop in assembled wild-type CA from two strains exhibits unprecedented mobility on the nanosecond to microsecond timescales, and the experimental NMR dipolar order parameters are in quantitative agreement with those calculated from MD trajectories. Remarkably, the CypA loop dynamics of wild-type CA HXB2 assembly is significantly attenuated upon CypA binding, and the dynamics profiles of the A92E and G94D CypA escape mutants closely resemble that of wild-type CA assembly in complex with CypA. These results suggest that CypA loop dynamics is a determining factor in HIV-1's escape from CypA dependence.

Reassembly and co-crystallization of a family 9 processive endoglucanase from its component parts: structural and functional significance of the intermodular linker

Non-cellulosomal processive endoglucanase 9I (Cel9I) from Clostridium thermocellum is a modular protein, consisting of a family-9 glycoside hydrolase (GH9) catalytic module and two family-3 carbohydrate-binding modules (CBM3c and CBM3b), separated by linker regions. GH9 does not show cellulase activity when expressed without CBM3c and CBM3b and the presence of the CBM3c was previously shown to be essential for endoglucanase activity. Physical reassociation of independently expressed GH9 and CBM3c modules (containing linker sequences) restored 60-70% of the intact Cel9I endocellulase activity. However, the mechanism responsible for recovery of activity remained unclear. In this work we independently expressed recombinant GH9 and CBM3c with and without their interconnecting linker in Escherichia coli. We crystallized and determined the molecular structure of the GH9/linker-CBM3c heterodimer at a resolution of 1.68 Å to understand the functional and structural importance of the mutual spatial orientation of the modules and the role of the interconnecting linker during their re-association. Enzyme activity assays and isothermal titration calorimetry were performed to study and compare the effect of the linker on the re-association. The results indicated that reassembly of the modules could also occur without the linker, albeit with only very low recovery of endoglucanase activity. We propose that the linker regions in the GH9/CBM3c endoglucanases are important for spatial organization and fixation of the modules into functional enzymes.

Calcium-controlled conformational choreography in the N-terminal half of adseverin

Adseverin is a member of the calcium-regulated gelsolin superfamily of actin-binding proteins. Here we report the crystal structure of the calcium-free N-terminal half of adseverin (iA1-A3) and the Ca(2+)-bound structure of A3, which reveal structural similarities and differences with gelsolin. Solution small-angle X-ray scattering combined with ensemble optimization revealed a dynamic Ca(2+)-dependent equilibrium between inactive, intermediate and active conformations. Increasing calcium concentrations progressively shift this equilibrium from a main population of inactive conformation to the active form. Molecular dynamics simulations of iA1-A3 provided insights into Ca(2+)-induced destabilization, implicating a critical role for the A2 type II calcium-binding site and the A2A3 linker in the activation process. Finally, mutations that disrupt the A1/A3 interface increase Ca(2+)-independent F-actin severing by A1-A3, albeit at a lower efficiency than observed for gelsolin domains G1-G3. Together, these data address the calcium dependency of A1-A3 activity in relation to the calcium-independent activity of G1-G3.