The entropic force generated by intrinsically disordered segments tunes protein function

Reinforcement learning to boost molecular docking upon protein conformational ensemble

Intrinsically disordered proteins (IDPs) are widely involved in human diseases and thus are attractive therapeutic targets. In practice, however, it is computationally prohibitive to dock large ligand libraries to thousands and tens of thousands of conformations. Here, we propose a reversible upper confidence bound (UCB) algorithm for the virtual screening of IDPs to address the influence of the conformation ensemble. The docking process is dynamically arranged so that attempts are focused near the boundary to separate top ligands from the bulk accurately. It is demonstrated in the example of transcription factor c-Myc that the average docking number per ligand can be greatly reduced while the performance is merely slightly affected. This study suggests that reinforcement learning is highly efficient in solving the bottleneck of virtual screening due to the conformation ensemble in the rational drug design of IDPs.

Heterogeneous dynamics in partially disordered proteins

Importance of disordered protein regions is increasingly recognized in biology, but their characterization remains challenging due to the lack of suitable experimental and theoretical methods. NMR experiments can detect multiple timescale dynamics and structural details of disordered protein regions, but their detailed interpretation is often difficult. Here we combine protein backbone 15N spin relaxation data with molecular dynamics (MD) simulations to detect not only heterogeneous dynamics of large partially disordered proteins but also their conformational ensembles. We observed that the rotational dynamics of folded regions in partially disordered proteins is dominated by similar rigid body rotation as in globular proteins, thereby being largely independent of flexible disordered linkers. Disordered regions, on the other hand, exhibit complex rotational motions with multiple timescales below ∼30 ns which are difficult to detect from experimental data alone, but can be captured by MD simulations. Combining MD simulations and backbone 15N spin relaxation data, measured applying segmental isotopic labeling with salt-inducible split intein, we resolved the conformational ensemble and dynamics of partially disordered periplasmic domain of TonB protein from Helicobacter pylori containing 250 residues. To demonstrate the universality of our approach, it was applied also to the partially disordered region of chicken Engrailed 2. Our results pave the way in understanding how TonB transfers energy from inner membrane to the outer membrane receptors in Gram-negative bacteria, as well as the function of other proteins with disordered domains.

Intrinsic Disorder in the T Cell Receptor Creates Cooperativity and Controls ZAP70 Binding

Many immunoreceptors have cytoplasmic domains that are intrinsically disordered (i.e., have high configurational entropy), have multiple sites of posttranslational modification (e.g., tyrosine phosphorylation), and participate in nonlinear signaling pathways (e.g., exhibiting switch-like behavior). Several hypotheses to explain the origin of these nonlinearities fall under the broad hypothesis that modification at one site changes the immunoreceptor’s entropy, which in turn changes further modification dynamics. Here, we use coarse-grain simulation to study three scenarios, all related to the chains that constitute the T cell receptor (TCR). We find that first, if phosphorylation induces local changes in the flexibility of the TCR ζ-chain, this naturally leads to rate enhancements and cooperativity. Second, we find that TCR CD3ɛ can provide a switch by modulating its residence in the plasma membrane. By constraining our model to be consistent with the previous observation that both basic residues and phosphorylation control membrane residence, we find that there is only a moderate rate enhancement of 10% between first and subsequent phosphorylation events. Third, we find that volume constraints do not limit the number of ZAP70s that can bind the TCR but that entropic penalties lead to a 200-fold decrease in binding rate by the seventh ZAP70, potentially explaining the observation that each TCR has around six ZAP70 molecules bound after receptor triggering. In all three scenarios, our results demonstrate that phenomena that change an immunoreceptor chain’s entropy (stiffening, confinement to a membrane, and multiple simultaneous binding) can lead to nonlinearities (rate enhancement, switching, and negative cooperativity) in how the receptor participates in signaling. These polymer-entropy-driven nonlinearities may augment the nonlinearities that arise from, e.g., kinetic proofreading and cluster formation. They also suggest different design strategies for engineered receptors, e.g., whether or not to put signaling modules on one chain or multiple clustered chains.

αβDCA method identifies unspecific binding but specific disruption of the group I intron by the StpA chaperone

Chaperone proteins-the most disordered among all protein groups-help RNAs fold into their functional structure by destabilizing misfolded configurations or stabilizing the functional ones. But disentangling the mechanism underlying RNA chaperoning is challenging, mostly because of inherent disorder of the chaperones and the transient nature of their interactions with RNA. In particular, it is unclear how specific the interactions are and what role is played by amino acid charge and polarity patterns. Here, we address these questions in the RNA chaperone StpA. We adapted direct coupling analysis (DCA) into the αβDCA method that can treat in tandem sequences written in two alphabets, nucleotides and amino acids. With αβDCA, we could analyze StpA-RNA interactions and show consistency with a previously proposed two-pronged mechanism: StpA disrupts specific positions in the group I intron while globally and loosely binding to the entire structure. Moreover, the interactions are strongly associated with the charge pattern: Negatively charged regions in the destabilizing StpA amino-terminal affect a few specific positions in the RNA, located in stems and in the pseudoknot. In contrast, positive regions in the carboxy-terminal contain strongly coupled amino acids that promote nonspecific or weakly specific binding to the RNA. The present study opens new avenues to examine the functions of disordered proteins and to design disruptive proteins based on their charge patterns.

The physical origin of rate promoting vibrations in enzymes revealed by structural rigidity

Enzymes are the most efficient catalysts known to date. However, decades of research have failed to fully explain the catalytic power of enzymes, and most of the current attempts to uncloak the details of atomic motions at active sites remain incomplete. Here, a straightforward manner for understanding the interplay between the complex or irregular enzyme topology and dynamical effects at catalytic sites is introduced, by revealing how fast localized vibrations form spontaneously in the stiffest parts of the scaffold. While shedding light on a physical mechanism that allowed the selection of the picosecond (ps) timescale to increase the catalytic proficiency, this approach exposes the functional importance of localized motions as a by-product of the stability-function tradeoff in enzyme evolution. From this framework of analysis—directly accessible from available diffraction data—experimental strategies for engineering the catalytic rate in enzymatic proteins are proposed.

Ensemble-Based Thermodynamics of the Fuzzy Binding between Intrinsically Disordered Proteins and Small-Molecule Ligands

In contrast to the "lock-and-key" model underlying the long-term success of structural biology and rational drug design, intrinsically disordered proteins (IDPs) exist in an ensemble of highly heterogeneous conformations even after binding with small-molecule ligands. It remains controversial how to characterize the thermodynamics of such fuzzy interactions. Here, we derive an ensemble-based thermodynamic framework to analyze the apparent affinity between IDPs and ligands. It is shown that the apparent affinity is related to the interaction free energy between the individual conformation and ligand in a way similar to Jarzynski's equality in nonequilibrium statistics. The oncoprotein c-Myc is adopted as an example to demonstrate the related properties, for example, the distribution of conformation-ligand interaction free energy, the entropic contribution from the ensemble, the conformation shift under ligand binding, and how to control the error under a limited number of sampled conformations.

Splice-Junction-Based Mapping of Alternative Isoforms in the Human Proteome

The protein-level translational status and function of many alternative splicing events remain poorly understood. We use an RNA sequencing (RNA-seq)-guided proteomics method to identify protein alternative splicing isoforms in the human proteome by constructing tissue-specific protein databases that prioritize transcript splice junction pairs with high translational potential. Using the custom databases to reanalyze ~80 million mass spectra in public proteomics datasets, we identify more than 1,500 noncanonical protein isoforms across 12 human tissues, including ~400 sequences undocumented on TrEMBL and RefSeq databases. We apply the method to original quantitative mass spectrometry experiments and observe widespread isoform regulation during human induced pluripotent stem cell cardiomyocyte differentiation. On a proteome scale, alternative isoform regions overlap frequently with disordered sequences and post-translational modification sites, suggesting that alternative splicing may regulate protein function through modulating intrinsically disordered regions. The described approach may help elucidate functional consequences of alternative splicing and expand the scope of proteomics investigations in various systems.

The translation and function of many alternative splicing events await confirmation at the protein level. Lau et al. use an integrated proteotranscriptomics approach to identify non-canonical and undocumented isoforms from 12 organs in the human proteome. Alternative isoforms interfere with functional sequence features and are differentially regulated during iPSC cardiomyocyte differentiation.

Disorder in a two-domain neuronal Ca2+-binding protein regulates domain stability and dynamics using ligand mimicry

Understanding the interplay between sequence, structure and function of proteins has been complicated in recent years by the discovery of intrinsically disordered proteins (IDPs), which perform biological functions in the absence of a well-defined three-dimensional fold. Disordered protein sequences account for roughly 30% of the human proteome and in many proteins, disordered and ordered domains coexist. However, few studies have assessed how either feature affects the properties of the other. In this study, we examine the role of a disordered tail in the overall properties of the two-domain, calcium-sensing protein neuronal calcium sensor 1 (NCS-1). We show that loss of just six of the 190 residues at the flexible C-terminus is sufficient to severely affect stability, dynamics, and folding behavior of both ordered domains. We identify specific hydrophobic contacts mediated by the disordered tail that may be responsible for stabilizing the distal N-terminal domain. Moreover, sequence analyses indicate the presence of an LSL-motif in the tail that acts as a mimic of native ligands critical to the observed order-disorder communication. Removing the disordered tail leads to a shorter life-time of the ligand-bound complex likely originating from the observed destabilization. This close relationship between order and disorder may have important implications for how investigations into mixed systems are designed and opens up a novel avenue of drug targeting exploiting this type of behavior.

Re-engineering Electrochemical Aptamer-Based Biosensors to Tune Their Useful Dynamic Range via Distal-Site Mutation and Allosteric Inhibition

Electrochemical aptamer-based (E-AB) sensors, exploiting binding-induced changes in biomolecular conformation, are rapid, specific, and selective and perform well even in a complex matrix, such as directly in whole blood and even in vivo. However, like all sensors employing biomolecular recognitions, E-AB sensors suffer from an inherent limitation of single-site binding, i.e., its fixed dose-response curve. To circumvent this, we employ here distal-site mutation and allosteric inhibition to rationally tune the dynamic range of E-AB sensors, achieving sets of sensors with a significantly varied target affinity (∼3 orders of magnitude). Using their combination, we recreate several approaches to narrow (down to 5-fold) or extend (up to 2000-fold) the dynamic range of biological receptors. The thermodynamic consequences of aptamer-surface interactions are estimated via the free-energy difference in solution-phase and surface-bound biosensors employing the same aptamer as a recognition element, revealing that an allostery strategy provides a more predictable and efficient means to finely control the target affinity and dynamic range. Such an ability to rationally modulate the affinity of biomolecule receptors would open the door to applications including cancer therapy, bioelectronics, and many other fields employing biomolecule recognition.

Integration of Sugar Metabolism and Proteoglycan Synthesis by UDP-glucose Dehydrogenase

Regulation of proteoglycan and glycosaminoglycan synthesis is critical throughout development, and to maintain normal adult functions in wound healing and the immune system, among others. It has become increasingly clear that these processes are also under tight metabolic control and that availability of carbohydrate and amino acid metabolite precursors has a role in the control of proteoglycan and glycosaminoglycan turnover. The enzyme uridine diphosphate (UDP)-glucose dehydrogenase (UGDH) produces UDP-glucuronate, an essential precursor for new glycosaminoglycan synthesis that is tightly controlled at multiple levels. Here, we review the cellular mechanisms that regulate UGDH expression, discuss the structural features of the enzyme, and use the structures to provide a context for recent studies that link post-translational modifications and allosteric modulators of UGDH to its function in downstream pathways.

Dehydration entropy drives liquid-liquid phase separation by molecular crowding

Complex coacervation driven liquid-liquid phase separation (LLPS) of biopolymers has been attracting attention as a novel phase in living cells. Studies of LLPS in this context are typically of proteins harboring chemical and structural complexity, leaving unclear which properties are fundamental to complex coacervation versus protein-specific. This study focuses on the role of polyethylene glycol (PEG)-a widely used molecular crowder-in LLPS. Significantly, entropy-driven LLPS is recapitulated with charged polymers lacking hydrophobicity and sequence complexity, and its propensity dramatically enhanced by PEG. Experimental and field-theoretic simulation results are consistent with PEG driving LLPS by dehydration of polymers, and show that PEG exerts its effect without partitioning into the dense coacervate phase. It is then up to biology to impose additional variations of functional significance to the LLPS of biological systems.

Physical Chemistry of the Protein Backbone: Enabling the Mechanisms of Intrinsic Protein Disorder

Over the last two decades it has become clear that well-defined structure is not a requisite for proteins to properly function. Rather, spectra of functionally competent, structurally disordered states have been uncovered requiring canonical paradigms in molecular biology to be revisited or reimagined. It is enticing and oftentimes practical to divide the proteome into structured and unstructured, or disordered, proteins. While function, composition, and structural properties largely differ, these two classes of protein are built upon the same scaffold, namely, the protein backbone. The versatile physicochemical properties of the protein backbone must accommodate structural disorder, order, and transitions between these states. In this review, we survey these properties through the conceptual lenses of solubility and conformational populations and in the context of protein-disorder mediated phenomena (e.g., phase separation, order-disorder transitions, allostery). Particular attention is paid to the results of computational studies, which, through thermodynamic decomposition and dissection of molecular interactions, can provide valuable mechanistic insight and testable hypotheses to guide further solution experiments. Lastly, we discuss changes in the dynamics of side chains and order-disorder transitions of the protein backbone as two modes or realizations of "entropic reservoirs" capable of tuning coupled thermodynamic processes.

DisProt: intrinsic protein disorder annotation in 2020

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the ‘dark’ proteome.

Flexibility of intrinsically disordered degrons in AUX/IAA proteins reinforces auxin co-receptor assemblies

Cullin RING-type E3 ubiquitin ligases SCF^TIR1/AFB1-5 and their AUX/IAA targets perceive the phytohormone auxin. The F-box protein TIR1 binds a surface-exposed degron in AUX/IAAs promoting their ubiquitylation and rapid auxin-regulated proteasomal degradation. Here, by adopting biochemical, structural proteomics and in vivo approaches we unveil how flexibility in AUX/IAAs and regions in TIR1 affect their conformational ensemble allowing surface accessibility of degrons. We resolve TIR1·auxin·IAA7 and TIR1·auxin·IAA12 complex topology, and show that flexible intrinsically disordered regions (IDRs) in the degron’s vicinity, cooperatively position AUX/IAAs on TIR1. We identify essential residues at the TIR1 N- and C-termini, which provide non-native interaction interfaces with IDRs and the folded PB1 domain of AUX/IAAs. We thereby establish a role for IDRs in modulating auxin receptor assemblies. By securing AUX/IAAs on two opposite surfaces of TIR1, IDR diversity supports locally tailored positioning for targeted ubiquitylation, and might provide conformational flexibility for a multiplicity of functional states.

Auxin-mediated recruitment of AUX/IAAs by the F-box protein TIR1 prompts rapid AUX/IAA ubiquitylation and degradation. By resolving auxin receptor topology, the authors show that intrinsically disordered regions near the degrons of two Aux/IAA proteins reinforce complex assembly and position Aux/IAAs for ubiquitylation.

NMR Structure and Dynamics of TonB Investigated by Scar-Less Segmental Isotopic Labeling Using a Salt-Inducible Split Intein

The growing understanding of partially unfolded proteins increasingly points to their biological relevance in allosteric regulation, complex formation, and protein design. However, the structural characterization of disordered proteins remains challenging. NMR methods can access both the dynamics and structures of such proteins, yet suffering from a high degeneracy of NMR signals. Here, we overcame this bottleneck utilizing a salt-inducible split intein to produce segmentally isotope-labeled samples with the native sequence, including the ligation junction. With this technique, we investigated the NMR structure and conformational dynamics of TonB from Helicobacter pylori in the presence of a proline-rich low complexity region. Spin relaxation experiments suggest that the several nano-second time scale dynamics of the C-terminal domain (CTD) is almost independent of the faster pico-to-nanosecond dynamics of the low complexity central region (LCCR). Our results demonstrate the utility of segmental isotopic labeling for proteins with heterogenous dynamics such as TonB and could advance NMR studies of other partially unfolded proteins.

Allostery of multidomain proteins with disordered linkers

Intrinsically disordered regions are often involved in allosteric regulation of multidomain proteins. They can act as disordered linkers to connect and interact with domains, resulting in rather complex allosteric mechanism and novel protein behavior. Therefore, it is necessary to analyze the diverse functions of disordered linkers in order to better understand allostery and relevant regulation process. Here we summarize recent advances in understanding the function of linkers and the advantages of adopting mutlidomain architecture with disorder linkers. It was shown that linkers between domains enhance the local domain concentration and make the allosteric regulation of weakly interacting partners possible, while linkers with only one tethered end cause an entropy effect to reduce binding affinity and prevent aggregation.

The mechanism of Hsp90-induced oligomerizaton of Tau

Aggregation of the microtubule-associated protein Tau is a hallmark of Alzheimer's disease with Tau oligomers suspected as the most toxic agent. Tau is a client of the molecular chaperone Hsp90, although it is unclear whether and how the chaperone massages the structure of intrinsically disordered Tau. Using electron paramagnetic resonance, we extract structural information from the very broad conformational ensemble of Tau: Tau in solution is highly dynamic and polymorphic, although "paper clip"-shaped by long-range contacts. Interaction with Hsp90 promotes an open Tau conformation, which we identify as the molecular basis for the formation of small Tau oligomers by exposure of the aggregation-prone repeat domain to other Tau molecules. At the same time, formation of Tau fibrils is inhibited. We therefore provide the nanometer-scale zoom into chaperoning an amyloid client, highlighting formation of oligomers as the consequence of this biologically relevant interaction.

Flexibility and intrinsic disorder are conserved features of hepatitis C virus E2 glycoprotein

The glycoproteins of hepatitis C virus, E1E2, are unlike any other viral fusion machinery yet described, and are the current focus of immunogen design in HCV vaccine development; thus, making E1E2 both scientifically and medically important. We used pre-existing, but fragmentary, structures to model a complete ectodomain of the major glycoprotein E2 from three strains of HCV. We then performed molecular dynamic simulations to explore the conformational landscape of E2, revealing a number of important features. Despite high sequence divergence, and subtle differences in the models, E2 from different strains behave similarly, possessing a stable core flanked by highly flexible regions, some of which perform essential functions such as receptor binding. Comparison with sequence data suggest that this consistent behaviour is conferred by a network of conserved residues that act as hinge and anchor points throughout E2. The variable regions (HVR-1, HVR-2 and VR-3) exhibit particularly high flexibility, and bioinformatic analysis suggests that HVR-1 is a putative intrinsically disordered protein region. Dynamic cross-correlation analyses demonstrate intramolecular communication and suggest that specific regions, such as HVR-1, can exert influence throughout E2. To support our computational approach we performed small-angle X-ray scattering with purified E2 ectodomain; this data was consistent with our MD experiments, suggesting a compact globular core with peripheral flexible regions. This work captures the dynamic behaviour of E2 and has direct relevance to the interaction of HCV with cell-surface receptors and neutralising antibodies.

Analyzing the Sequences of Intrinsically Disordered Regions with CIDER and localCIDER

Intrinsically disordered proteins and protein regions are ubiquitous across eukaryotic proteomes where they play a range of functional roles. Unlike folded proteins, IDRs lack a well-defined native state but exist in heterogeneous ensembles of conformations. In the absence of a defined native state, structure-guided mutations to test specific mechanistic hypotheses are generally not possible. Despite this, the use of mutations to alter sequence properties has become a relatively common approach for teasing out the relationship between sequence, ensemble, and function. A key step in designing informative mutants is the ability to identify specific sequence features that may reveal an interpretable response if perturbed. Here, we provide guidance on using the CIDER and localCIDER tools for amino acid sequence analysis, with a focus on building intuition with respect to the most commonly described features.

Linker Domains: Why ABC Transporters 'Live in Fragments no Longer'

ATP-binding cassette (ABC) transporters are membrane proteins present in all kingdoms of life. We have considered the disordered region that connects the N- and C-terminal halves in many eukaryotic ABC transporters, allowing all four consensus functional domains to be linked. The recent availability of structures of ABC transporters containing linker regions has allowed us to identify the start and end points of the connectors as well as hinting at their localisation. We address questions such as: Where did the linker regions come from? Why do some ABC transporters have connectors and others not? What are the rules and roles of the linker regions? What are the consequences of mutations in these connector regions for disease in humans?

Ribosome Evolution and Structural Capacitance

In addition to the canonical loss-of-function mutations, mutations in proteins may additionally result in gain-of-function through the binary activation of cryptic "structural capacitance elements." Our previous bioinformatic analysis allowed us to propose a new mechanism of protein evolution - structural capacitance - that arises via the generation of new elements of microstructure upon mutations that cause a disorder-to-order (D→O) transition in previously disordered regions of proteins. Here we propose that the D→O transition is a necessary follow-on from expected early codon-anticodon and tRNA acceptor stem-amino acid usage, via the accumulation of structural capacitance elements - reservoirs of disorder in proteins. We develop this argument further to posit that structural capacitance is an inherent consequence of the evolution of the genetic code.

The Plasmodium falciparum Hsp70-x chaperone assists the heat stress response of the malaria parasite

Plasmodium falciparum is the most lethal of human-infective malaria parasites. A hallmark of P. falciparum malaria is extensive remodeling of host erythrocytes by the parasite, which facilitates the development of virulence properties such as host cell adhesion to the endothelial lining of the microvasculature. Host remodeling is mediated by a large complement of parasite proteins exported to the erythrocyte; among them is a single heat shock protein (Hsp)70-class protein chaperone, P. falciparum Hsp70-x (PfHsp70-x). PfHsp70-x was previously shown to assist the development of virulent cytoadherence characteristics. Here, we show that PfHsp70-x also supports parasite growth under elevated temperature conditions that simulate febrile episodes, especially at the beginning of the parasite life cycle when most of host cell remodeling takes place. Biochemical and biophysical analyses of PfHsp70-x, including crystallographic structures of its catalytic domain and the J-domain of its stimulatory Hsp40 cochaperone, suggest that PfHsp70-x is highly similar to human Hsp70 chaperones endogenous to the erythrocyte. Nevertheless, our results indicate that selective inhibition of PfHsp70-x function using small molecules may be possible and highlight specific sites of its catalytic domain as potentially of high interest. We discuss the likely roles of PfHsp70-x and human chaperones in P. falciparum biology and how specific inhibitors may assist us in disentangling their relative contributions.-Day, J., Passecker, A., Beck, H.-P., Vakonakis, I. The Plasmodium falciparum Hsp70-x chaperone assists the heat stress response of the malaria parasite.

Mechanochemistry in Translation

As the influence of translation rates on protein folding and function has come to light, the mechanisms by which translation speed is modulated have become an important issue. One mechanism entails the generation of force by the nascent protein. Cotranslational processes, such as nascent protein folding, the emergence of unfolded nascent chain segments from the ribosome's exit tunnel, and insertion of the nascent chain into or translocation of the nascent chain through membranes, can generate forces that are transmitted back to the peptidyl transferase center and affect translation rates. In this Perspective, we examine the processes that generate these forces, the mechanisms of transmission along the ribosomal exit tunnel to the peptidyl transferase center, and the effects of force on the ribosome's catalytic cycle. We also discuss the physical models that have been developed to predict and explain force generation for individual processes and speculate about other processes that may generate forces that have yet to be tested.

MglA functions as a three-state GTPase to control movement reversals of Myxococcus xanthus

In Myxococcus xanthus, directed movement is controlled by pole-to-pole oscillations of the small GTPase MglA and its GAP MglB. Direction reversals require that MglA is inactivated by MglB, yet paradoxically MglA and MglB are located at opposite poles at reversal initiation. Here we report the complete MglA/MglB structural cycle combined to GAP kinetics and in vivo motility assays, which uncovers that MglA is a three-state GTPase and suggests a molecular mechanism for concerted MglA/MglB relocalizations. We show that MglA has an atypical GTP-bound state (MglA-GTP*) that is refractory to MglB and is re-sensitized by a feedback mechanism operated by MglA-GDP. By identifying and mutating the pole-binding region of MglB, we then provide evidence that the MglA-GTP* state exists in vivo. These data support a model in which MglA-GDP acts as a soluble messenger to convert polar MglA-GTP* into a diffusible MglA-GTP species that re-localizes to the opposite pole during reversals.

In Myxococcus xanthus, directed movement is controlled by pole-to-pole oscillations of the small GTPase MglA and its GAP MglB. Here authors report the complete MglA/MglB structural cycle and uncover that MglA is a three-state GTPase that adopts an atypical GTP-bound state that is refractory to inactivation by MglB.

The Hypervariable Region of K-Ras4B Governs Molecular Recognition and Function

The flexible C-terminal hypervariable region distinguishes K-Ras4B, an important proto-oncogenic GTPase, from other Ras GTPases. This unique lysine-rich portion of the protein harbors sites for post-translational modification, including cysteine prenylation, carboxymethylation, phosphorylation, and likely many others. The functions of the hypervariable region are diverse, ranging from anchoring K-Ras4B at the plasma membrane to sampling potentially auto-inhibitory binding sites in its GTPase domain and participating in isoform-specific protein-protein interactions and signaling. Despite much research, there are still many questions about the hypervariable region of K-Ras4B. For example, mechanistic details of its interaction with plasma membrane lipids and with the GTPase domain require further clarification. The roles of the hypervariable region in K-Ras4B-specific protein-protein interactions and signaling are incompletely defined. It is also unclear why post-translational modifications frequently found in protein polylysine domains, such as acetylation, glycation, and carbamoylation, have not been observed in K-Ras4B. Expanding knowledge of the hypervariable region will likely drive the development of novel highly-efficient and selective inhibitors of K-Ras4B that are urgently needed by cancer patients.

Conservation of Atypical Allostery in C. elegans UDP-Glucose Dehydrogenase

Human UDP-glucose dehydrogenase (hUGDH) oxidizes uridine diphosphate (UDP)-glucose to UDP-glucuronic acid, an essential substrate in the phase II metabolism of drugs. The activity of hUGDH is controlled by an atypical allosteric mechanism in which the feedback inhibitor UDP-xylose competes with the substrate for the active site and triggers a buried allosteric switch to produce an inactive complex (E^Ω). Previous comparisons with a nonallosteric UGDH identified six large-to-small substitutions that produce packing defects in the protein core and provide the conformational flexibility necessary for the allosteric transition. Here, we test the hypothesis that these large-to-small substitutions form a motif that can be used to identify allosteric UGDHs. Caenorhabditis elegans UGDH (cUGDH) conserves this motif with the exception of an Ala-to-Pro substitution in position 109. The crystal structures of unliganded and UDP-xylose bound cUGDH show that the A109P substitution is accommodated by an Asn-to-Ser substitution at position 290. Steady-state analysis and sedimentation velocity studies show that the allosteric transition is conserved in cUGDH. The enzyme also exhibits hysteresis in progress curves and negative cooperativity with respect to NAD⁺ binding. Both of these phenomena are conserved in the human enzyme, which is strong evidence that these represent fundamental features of atypical allostery in UGDH. A phylogenetic analysis of UGDH shows that the atypical allostery motif is ancient and identifies a potential transition point in the evolution of the UGDH family.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

Spontaneous Rupture and Entanglement of Human Neuronal Tau Protein Induced by Piconewton Compressive Force

Mechanical force vector fluctuations in living cells can have a significant impact on protein behavior and functions. Here we report that a human tau protein tertiary structure can abruptly and spontaneously rupture, like a balloon, under biologically available piconewton compressive force, using a home-modified atomic force microscopy single-molecule manipulation. The rupture behavior is dependent on the physiological level of presence of ions, such as K⁺ and Mg²⁺. We observed rupture events in the presence of K⁺ but not in the presence of Mg²⁺ ions. We have also explored the entangled protein state formed following the events of the multiple and simultaneous protein ruptures under crowding. Crowded proteins simultaneously rupture and then spontaneously refold to an entangled folding state, different from either folded and unfolded states of the tau protein, which can be a plausible pathway for the tau protein aggregation that is related to a number of neurodegenerative diseases.

Allosteric Regulation of HIV-1 Capsid Structure for Gag Assembly, Virion Production, and Viral Infectivity by a Disordered Interdomain Linker

The retroviral Gag capsid (Gag-CA) interdomain linker is an unstructured peptide segment connecting structured N-terminal and C-terminal domains. Although the region is reported to play roles in virion morphogenesis and infectivity, underlying molecular mechanisms remain unexplored. To address this issue, we determined biological and molecular phenotypes of HIV-1 CA linker mutants by experimental and in silico approaches. Among the nine linker mutants tested, eight exhibited attenuation of viral particle production to various extents mostly in parallel with a reduction in viral infectivity. Sucrose density gradient, confocal microscopy, and live-cell protein interaction analyses indicated that the defect is accompanied by attenuation of Gag-Gag interactions following Gag plasma membrane targeting in the cells. In silico analyses revealed distinct distributions of interaction-prone hydrophobic patches between immature and mature CA proteins. Molecular dynamics simulations predicted that the linker mutations can allosterically alter structural fluctuations, including the interaction surfaces apart from the mutation sites in both the immature and mature CA proteins. These results suggest that the HIV-1 CA interdomain linker is a cis-modulator of the CA interaction surfaces to optimize efficiency of Gag assembly, virion production, and viral infectivity.IMPORTANCE HIV-1 particle production and infection are highly ordered processes. Viral Gag proteins play a central role in the assembly and disassembly of viral molecules. Of these, capsid protein (CA) is a major contributor to the Gag-Gag interactions. CA consists of two structured domains, i.e., N-terminal (NTD) and C-terminal (CTD) domains, connected by an unstructured domain named the interdomain linker. While multiple regions in the NTD and CTD are reported to play roles in virion morphogenesis and infectivity, the roles of the linker region in Gag assembly and virus particle formation remain elusive. In this study, we showed by biological and molecular analyses that the linker region functions as an intramolecular modulator to tune Gag assembly, virion production, and viral infectivity. Our study thus illustrates a hitherto-unrecognized mechanism, an allosteric regulation of CA structure by the disordered protein element, for HIV-1 replication.

A Note on the Entropy Force in Kinetic Theory and Black Holes

The entropy force is the collective effect of inhomogeneity in disorder in a statistical many particle system. We demonstrate its presumable effect on one particular astrophysical object, the black hole. We then derive the kinetic equations of a large system of particles including the entropy force. It adds a collective therefore integral term to the Klimontovich equation for the evolution of the one-particle distribution function. Its integral character transforms the basic one particle kinetic equation into an integro-differential equation already on the elementary level, showing that not only the microscopic forces but the hole system reacts to its evolution of its probability distribution in a holistic way. It also causes a collisionless dissipative term which however is small in the inverse particle number and thus negligible. However it contributes an entropic collisional dissipation term. The latter is defined via the particle correlations but lacks any singularities and thus is large scale. It allows also for the derivation of a kinetic equation for the entropy density in phase space. This turns out to be of same structure as the equation for the phase space density. The entropy density determines itself holistically via the integral entropy force thus providing a self-controlled evolution of entropy in phase space.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein's structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15-150), small (ca. 4 Hz) 3 JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy ( D∥/ D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Entropy-Driven Mechanisms between Disulfide-Bond Formation Protein A (DsbA) and B (DsbB) in Escherichia coli

A disulfide-bond formation system for nascent proteins in the Escherichia coli periplasm contains efficient electron transfer systems for the catalysis of oxidation. This electrochemical system has interesting implications in vivo. Disulfide bonds are formed by disulfide-bond formation protein A (DsbA), which contains two reactive cysteines. DsbA is reoxidized by a membrane protein, disulfide-bond formation protein B (DsbB), which has four catalytic cysteines. The oxidation of DsbA by DsbB seems energetically unfavorable on the basis of the redox potential. The oxidizing power of ubiquinone (UQ), which endogenously binds with DsbB, is believed to promote this reaction. However, using UQ-deficient DsbB, it was found that the oxidation of DsbA by DsbB proceeds independently of UQ. Thus, the reaction mechanism of DsbA oxidation by DsbB is under debate. In this study, we used the quartz crystal microbalance technique, which detects the intermediate complex between DsbA and DsbB during DsbA oxidation as a change in mass, to obtain kinetic parameters of DsbA oxidation under both the oxidized and reduced states of UQ at acidic and basic pH. In addition, we utilized sodium dodecyl sulfate polyacrylamide gel electrophoresis mobility shift assay technique to determine the pK _a of the cysteine thiol groups in DsbA and DsbB. We found that DsbA oxidation proceeded independently of UQ and was greatly affected in kinetics by the shuffling of electrons among the four cysteine residues in DsbB, regardless of pH. These results suggest that DsbA oxidation is driven in an entropy-dependent manner, in which the electron-delocalized intermediate complex is stabilized by preventing a reverse reaction. These findings could contribute to the design of bio-inspired electrochemical systems for industrial applications.

Regulation of measles virus gene expression by P protein coiled-coil properties

The polymerase of negative-stranded RNA viruses consists of the large protein (L) and the phosphoprotein (P), the latter serving both as a chaperon and a cofactor for L. We mapped within measles virus (MeV) P the regions responsible for binding and stabilizing L and showed that the coiled-coil multimerization domain (MD) of P is required for gene expression. MeV MD is kinked as a result of the presence of a stammer. Both restoration of the heptad regularity and displacement of the stammer strongly decrease or abrogate activity in a minigenome assay. By contrast, P activity is rather tolerant of substitutions within the stammer. Single substitutions at the "a" or "d" hydrophobic anchor positions with residues of variable hydrophobicity revealed that P functionality requires a narrow range of cohesiveness of its MD. Results collectively indicate that, beyond merely ensuring P oligomerization, the MD finely tunes viral gene expression through its cohesiveness.

Distinct USP25 and USP28 Oligomerization States Regulate Deubiquitinating Activity

The evolutionarily related deubiquitinating enzymes (DUBs) USP25 and USP28 comprise an identical overall domain architecture but are functionally non-redundant: USP28 stabilizes c-MYC and other nuclear proteins, and USP25 regulates inflammatory TRAF signaling. We here compare molecular features of USP25 and USP28. Active enzymes form distinctively shaped dimers, with a dimerizing insertion spatially separating independently active catalytic domains. In USP25, but not USP28, two dimers can form an autoinhibited tetramer, where a USP25-specific, conserved insertion sequence blocks ubiquitin binding. In full-length enzymes, a C-terminal domain with a previously unknown fold has no impact on oligomerization, but N-terminal regions affect the dimer-tetramer equilibrium in vitro. We confirm oligomeric states of USP25 and USP28 in cells and show that modulating oligomerization affects substrate stabilization in accordance with in vitro activity data. Our work highlights how regions outside of the catalytic domain enable a conceptually intriguing interplay of DUB oligomerization and activity.

Understanding Intramolecular Crosstalk in an Intrinsically Disordered Protein

The interaction between N_TAIL and XD from the measles virus represents a paradigmatic example of molecular recognition between an intrinsically disordered protein and a folded partner. By binding to XD, a small portion of N_TAIL (classically denoted as MoRE) undergoes a disorder-to-order transition, populating an α-helical structure, while the reminder of the protein remains disordered. Here, we demonstrate an unexpected crosstalk between such a disordered region and the adjacent molecular recognition element (MoRE). This result was obtained by producing a series of truncation and site-directed variants of N_TAIL while measuring the effects on the kinetics of folding and binding. We show that the disordered region of N_TAIL exerts its inhibitory role by slowing the folding step of the MoRE, thereby tuning the affinity of the interaction.

Omnipresent Maxwell's demons orchestrate information management in living cells

The development of synthetic biology calls for accurate understanding of the critical functions that allow construction and operation of a living cell. Besides coding for ubiquitous structures, minimal genomes encode a wealth of functions that dissipate energy in an unanticipated way. Analysis of these functions shows that they are meant to manage information under conditions when discrimination of substrates in a noisy background is preferred over a simple recognition process. We show here that many of these functions, including transporters and the ribosome construction machinery, behave as would behave a material implementation of the information-managing agent theorized by Maxwell almost 150 years ago and commonly known as Maxwell's demon (MxD). A core gene set encoding these functions belongs to the minimal genome required to allow the construction of an autonomous cell. These MxDs allow the cell to perform computations in an energy-efficient way that is vastly better than our contemporary computers.

Correlation Between Allosteric and Orthosteric Sites

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

Modulation of Measles Virus NTAIL Interactions through Fuzziness and Sequence Features of Disordered Binding Sites

In this paper we review our recent findings on the different interaction mechanisms of the C-terminal domain of the nucleoprotein (N) of measles virus (MeV) N_TAIL, a model viral intrinsically disordered protein (IDP), with two of its known binding partners, i.e., the C-terminal X domain of the phosphoprotein of MeV XD (a globular viral protein) and the heat-shock protein 70 hsp70 (a globular cellular protein). The N_TAIL binds both XD and hsp70 via a molecular recognition element (MoRE) that is flanked by two fuzzy regions. The long (85 residues) N-terminal fuzzy region is a natural dampener of the interaction with both XD and hsp70. In the case of binding to XD, the N-terminal fuzzy appendage of N_TAIL reduces the rate of α-helical folding of the MoRE. The dampening effect of the fuzzy appendage on XD and hsp70 binding depends on the length and fuzziness of the N-terminal region. Despite this similarity, N_TAIL binding to XD and hsp70 appears to rely on completely different requirements. Almost any mutation within the MoRE decreases XD binding, whereas many of them increase the binding to hsp70. In addition, XD binding is very sensitive to the α-helical state of the MoRE, whereas hsp70 is not. Thus, contrary to hsp70, XD binding appears to be strictly dependent on the wild-type primary and secondary structure of the MoRE.