Protein dynamics and conformational disorder in molecular recognition

Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

Pathogen elicitor peptide (pep), systemin, and their receptors in tomato: sequence analysis sheds light on standing disagreements about biotic stress signaling components

Peps are endogenous damage-associated polypeptides that evoke defense responses in plants. Like other damage-associated molecular patterns, Pep signals are transduced by receptors. PEPRs are the receptors that transduce Pep danger signals. This paper identifies new putative Peps in the Solanaceae (including Solanum spp., Nicotiana spp., and Petunia spp.) and Coffea and explores their properties. Using these newly identified Peps we derive sequence logos that present a refinement of the current understanding of the importance of specific residues in the Pep signaling molecules in Solanaceae, including several arginines, prolines that restrict peptide's conformations, and C-terminal asparagine. We examine the degree of disorder in Pep, which is likely important to the mechanism of Pep perception. This work also calls into question some of the evolutionary relationships between Peps in Solanaceae and specific Arabidopsis Peps published in previous literature, culminating in a conclusion that SlPep should not be named SlPep6 due to the lack of conservation of protein sequences in AtPROPEP6 and SlPROPEP, and that SlPep probably does not have two receptors in tomato, based on phylogenetic analysis. Our analyses advance understanding of the Pep signaling system in Solanaceae.

Regulation of the DLC3 tumor suppressor by a novel phosphoswitch

Deleted in liver cancer 3 (DLC3) is a Rho GTPase-activating protein (RhoGAP) that plays a crucial role in maintaining adherens junction integrity and coordinating polarized vesicle transport by modulating Rho activity at the plasma membrane and endomembranes. By employing bioinformatical sequence analysis, in vitro experiments, and in cellulo assays we here identified a polybasic region (PBR) in DLC3 that facilitates the association of the protein with cellular membranes. Within the PBR, we mapped two serines whose phosphorylation can alter the electrostatic character of the region. Consequently, phosphomimetic mutations of these sites impaired the membrane association of DLC3. Furthermore, we found a new PBR-dependent localization of DLC3 at the midbody region, where the protein locally controlled Rho activity. Here, the phosphorylation-dependent regulation of DLC3 appeared to be required for proper cytokinesis. Our work thus provides a novel mechanism for spatiotemporal termination of Rho signaling by the RhoGAP protein DLC3.

Membrane mimetic-dependence of GPCR energy landscapes

We leveraged variable-temperature 19F-NMR spectroscopy to compare the conformational equilibria of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), across a range of temperatures ranging from lower temperatures typically employed in 19F-NMR experiments to physiological temperature. A2AAR complexes with partial agonists and full agonists showed large increases in the population of a fully active conformation with increasing temperature. NMR data measured at physiological temperature were more in line with functional data. This was pronounced for complexes with partial agonists, where the population of active A2AAR was nearly undetectable at lower temperature but became evident at physiological temperature. Temperature-dependent behavior of complexes with either full or partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. Cellular signaling experiments correlated with the temperature-dependent conformational equilibria of A2AAR in lipid nanodiscs but not in some detergents, underscoring the importance of the membrane environment in studies of GPCR function.

Exploring NAD+ metabolism and NNAT: Insights from structure, function, and computational modeling

Nicotinamide Adenine Dinucleotide (NAD+), a coenzyme, is ubiquitously distributed and serves crucial functions in diverse biological processes, encompassing redox reactions, energy metabolism, and cellular signalling. This review article explores the intricate realm of NAD + metabolism, with a particular emphasis on the complex relationship between its structure, function, and the pivotal enzyme, Nicotinate Nucleotide Adenylyltransferase (NNAT), also known as nicotinate mononucleotide adenylyltransferase (NaMNAT), in the process of its biosynthesis. Our findings indicate that NAD + biosynthesis in humans and bacteria occurs via the same de novo synthesis route and the pyridine ring salvage pathway. Maintaining NAD homeostasis in bacteria is imperative, as most bacterial species cannot get NAD+ from their surroundings. However, due to lower sequence identity and structurally distant relationship of bacteria, including E. faecium and K. pneumonia, to its human counterpart, inhibiting NNAT, an indispensable enzyme implicated in NAD + biosynthesis, is a viable alternative in curtailing infections orchestrated by E. faecium and K. pneumonia. By merging empirical and computational discoveries and connecting the intricate NAD + metabolism network with NNAT's crucial role, it becomes clear that the synergistic effect of these insights may lead to a more profound understanding of the coenzyme's function and its potential applications in the fields of therapeutics and biotechnology.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

The Functionality of Membrane-Inserting Proteins and Peptides: Curvature Sensing, Generation, and Pore Formation

Proteins and peptides with hydrophobic and amphiphilic segments are responsible for many biological functions. The sensing and generation of membrane curvature are the functions of several protein domains or motifs. While some specific membrane proteins play an essential role in controlling the curvature of distinct intracellular membranes, others participate in various cellular processes such as clathrin-mediated endocytosis, where several proteins sort themselves at the neck of the membrane bud. A few membrane-inserting proteins form nanopores that permeate selective ions and water to cross the membrane. In addition, many natural and synthetic small peptides and protein toxins disrupt the membrane by inducing nonspecific pores in the membrane. The pore formation causes cell death through the uncontrolled exchange between interior and exterior cellular contents. In this article, we discuss the insertion depth and orientation of protein/peptide helices, and their role as a sensor and inducer of membrane curvature as well as a pore former in the membrane. We anticipate that this extensive review will assist biophysicists to gain insight into curvature sensing, generation, and pore formation by membrane insertion.

Functional unfoldomics: Roles of intrinsic disorder in protein (multi)functionality

Intrinsically disordered proteins (IDPs), which are functional proteins without stable tertiary structure, and hybrid proteins containing ordered domains and intrinsically disordered regions (IDRs) constitute prominent parts of all proteomes collectively known as unfoldomes. IDPs/IDRs exist as highly dynamic structural ensembles of rapidly interconverting conformations and are characterized by the exceptional structural heterogeneity, where their different parts are (dis)ordered to different degree, and their overall structure represents a complex mosaic of foldons, inducible foldons, inducible morphing foldons, non-foldons, semifoldons, and even unfoldons. Despite their lack of unique 3D structures, IDPs/IDRs play crucial roles in the control of various biological processes and the regulation of different cellular pathways and are commonly involved in recognition and signaling, indicating that the disorder-based functional repertoire is complementary to the functions of ordered proteins. Furthermore, IDPs/IDRs are frequently multifunctional, and this multifunctionality is defined by their structural flexibility and heterogeneity. Intrinsic disorder phenomenon is at the roots of the structure-function continuum model, where the structure continuum is defined by the presence of differently (dis)ordered regions, and the function continuum arises from the ability of all these differently (dis)ordered parts to have different functions. In their everyday life, IDPs/IDRs utilize a broad spectrum of interaction mechanisms thereby acting as interaction specialists. They are crucial for the biogenesis of numerous proteinaceous membrane-less organelles driven by the liquid-liquid phase separation. This review introduces functional unfoldomics by representing some aspects of the intrinsic disorder-based functionality.

Recognition and coacervation of G-quadruplexes by a multifunctional disordered region in RECQ4 helicase

Biomolecular polyelectrolyte complexes can be formed between oppositely charged intrinsically disordered regions (IDRs) of proteins or between IDRs and nucleic acids. Highly charged IDRs are abundant in the nucleus, yet few have been functionally characterized. Here, we show that a positively charged IDR within the human ATP-dependent DNA helicase Q4 (RECQ4) forms coacervates with G-quadruplexes (G4s). We describe a three-step model of charge-driven coacervation by integrating equilibrium and kinetic binding data in a global numerical model. The oppositely charged IDR and G4 molecules form a complex in the solution that follows a rapid nucleation-growth mechanism leading to a dynamic equilibrium between dilute and condensed phases. We also discover a physical interaction with Replication Protein A (RPA) and demonstrate that the IDR can switch between the two extremes of the structural continuum of complexes. The structural, kinetic, and thermodynamic profile of its interactions revealed a dynamic disordered complex with nucleic acids and a static ordered complex with RPA protein. The two mutually exclusive binding modes suggest a regulatory role for the IDR in RECQ4 function by enabling molecular handoffs. Our study extends the functional repertoire of IDRs and demonstrates a role of polyelectrolyte complexes involved in G4 binding.

Membrane Mimetic-Dependence of GPCR Energy Landscapes

Protein function strongly depends on temperature, which is related to temperature-dependent changes in the equilibria of protein conformational states. We leveraged variable-temperature 19F-NMR spectroscopy to interrogate the temperature dependence of the conformational landscape of the human A2A adenosine receptor (A2AAR), a class A GPCR. Temperature-induced changes in the conformational equilibria of A2AAR in lipid nanodiscs were markedly dependent on the efficacy of bound drugs. While antagonist complexes displayed only modest changes as the temperature rose, both full and partial agonist complexes exhibited substantial increases in the active state population. Importantly, the temperature-dependent response of complexes with both full and partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. In striking contrast to observations within lipid nanodiscs, in detergent micelles the active state population exhibited different behavior for A2AAR complexes with both full and partial agonists. This underscores the importance of the protein environment in understanding the thermodynamics of GPCR activation.

Driving forces of the complex formation between highly charged disordered proteins

Highly disordered complexes between oppositely charged intrinsically disordered proteins present a new paradigm of biomolecular interactions. Here, we investigate the driving forces of such interactions for the example of the highly positively charged linker histone H1 and its highly negatively charged chaperone, prothymosin α (ProTα). Temperature-dependent single-molecule Förster resonance energy transfer (FRET) experiments and isothermal titration calorimetry reveal ProTα-H1 binding to be enthalpically unfavorable, and salt-dependent affinity measurements suggest counterion release entropy to be an important thermodynamic driving force. Using single-molecule FRET, we also identify ternary complexes between ProTα and H1 in addition to the heterodimer at equilibrium and show how they contribute to the thermodynamics observed in ensemble experiments. Finally, we explain the observed thermodynamics quantitatively with a mean-field polyelectrolyte theory that treats counterion release explicitly. ProTα-H1 complex formation resembles the interactions between synthetic polyelectrolytes, and the underlying principles are likely to be of broad relevance for interactions between charged biomolecules in general.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

Mutational scan inferred binding energetics and structure in intrinsically disordered protein CcdA

Unlike globular proteins, mutational effects on the function of Intrinsically Disordered Proteins (IDPs) are not well-studied. Deep Mutational Scanning of a yeast surface displayed mutant library yields insights into sequence-function relationships in the CcdA IDP. The approach enables facile prediction of interface residues and local structural signatures of the bound conformation. In contrast to previous titration-based approaches which use a number of ligand concentrations, we show that use of a single rationally chosen ligand concentration can provide quantitative estimates of relative binding constants for large numbers of protein variants. This is because the extended interface of IDP ensures that energetic effects of point mutations are spread over a much smaller range than for globular proteins. Our data also provides insights into the much-debated role of helicity and disorder in partner binding of IDPs. Based on this exhaustive mutational sensitivity dataset, a rudimentary model was developed in an attempt to predict mutational effects on binding affinity of IDPs that form alpha-helical structures upon binding.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Enrichment patterns of intrinsic disorder in proteins

Intrinsically disordered regions in proteins have been shown to be important in protein function. However, not all proteins contain the same amount of intrinsic disorder. The variation in the levels of intrinsic disorder in different types of proteins has been extensively studied over the last two decades. It is now known that the levels of intrinsic disorder vary in proteins across organisms, functions, diseases, and cellular locations. This review consolidates the known trends in the abundance of intrinsic disorder identified in groups of proteins across varying conditions and functions. It also presents new data towards the understanding of intrinsic disorder in cell type-specific proteins.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01016-7.

A Mechanistic Model for Cell Cycle Control in Which CDKs Act as Switches of Disordered Protein Phase Separation

Cyclin-dependent kinases (CDKs) are presumed to control the cell cycle by phosphorylating a large number of proteins involved in S-phase and mitosis, two mechanistically disparate biological processes. While the traditional qualitative model of CDK-mediated cell cycle control relies on differences in inherent substrate specificity between distinct CDK-cyclin complexes, they are largely dispensable according to the opposing quantitative model, which states that changes in the overall CDK activity level promote orderly progression through S-phase and mitosis. However, a mechanistic explanation for how such an activity can simultaneously regulate many distinct proteins is lacking. New evidence suggests that the CDK-dependent phosphorylation of ostensibly very diverse proteins might be achieved due to underlying similarity of phosphorylation sites and of the biochemical effects of their phosphorylation: they are preferentially located within intrinsically disordered regions of proteins that are components of membraneless organelles, and they regulate phase separation. Here, we review this evidence and suggest a mechanism for how a single enzyme’s activity can generate the dynamics required to remodel the cell at mitosis.

NMR Reveals Specific Tracts within the Intrinsically Disordered Regions of the SARS-CoV-2 Nucleocapsid Protein Involved in RNA Encountering

The SARS-CoV-2 nucleocapsid (N) protein is crucial for the highly organized packaging and transcription of the genomic RNA. Studying atomic details of the role of its intrinsically disordered regions (IDRs) in RNA recognition is challenging due to the absence of structure and to the repetitive nature of their primary sequence. IDRs are known to act in concert with the folded domains of N and here we use NMR spectroscopy to identify the priming events of N interacting with a regulatory SARS-CoV-2 RNA element. ¹³C-detected NMR experiments, acquired simultaneously to ¹H detected ones, provide information on the two IDRs flanking the N-terminal RNA binding domain (NTD) within the N-terminal region of the protein (NTR, 1–248). We identify specific tracts of the IDRs that most rapidly sense and engage with RNA, and thus provide an atom-resolved picture of the interplay between the folded and disordered regions of N during RNA interaction.

Multivalent Interaction of Beta-Catenin With its Intrinsically Disordered Binding Partner Adenomatous Polyposis Coli

The Wnt signalling pathway plays key roles in cell proliferation, differentiation and fate decisions in embryonic development and maintenance of adult tissues, and the twelve Armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here we investigate the interaction between β-catenin's ARM repeat domain and the intrinsically disordered protein adenomatous polyposis coli (APC). APC is a giant multivalent scaffold that brings together the different components of the so-called "β-catenin destruction complex", which drives β-catenin degradation via the ubiquitin-proteasome pathway. Mutations and truncations in APC, resulting in loss of APC function and hence elevated β-catenin levels and upregulation of Wnt signalling, are associated with numerous cancers including colorectal carcinomas. APC has a long intrinsically disordered region (IDR) that contains a series of 15-residue and 20-residue binding regions for β-catenin. Here we explore the multivalent nature of the interaction of β-catenin with the highest affinity APC repeat, both at equilibrium and under kinetic conditions. We use a combination of single-site substitutions, deletions and insertions to dissect the mechanism of molecular recognition and the roles of the three β-catenin-binding subdomains of APC.

Physiological functions and roles in cancer of the proliferation marker Ki-67

What do we know about Ki-67, apart from its usefulness as a cell proliferation biomarker in histopathology? Discovered in 1983, the protein and its regulation of expression and localisation throughout the cell cycle have been well characterised. However, its function and molecular mechanisms have received little attention and few answers. Although Ki-67 has long been thought to be required for cell proliferation, recent genetic studies have conclusively demonstrated that this is not the case, as loss of Ki-67 has little or no impact on cell proliferation. In contrast, Ki-67 is important for localising nucleolar material to the mitotic chromosome periphery and for structuring perinucleolar heterochromatin, and emerging data indicate that it also has critical roles in cancer development. However, its mechanisms of action have not yet been fully identified. Here, we review recent findings and propose the hypothesis that Ki-67 is involved in structuring cellular sub-compartments that assemble by liquid-liquid phase separation. At the heterochromatin boundary, this may control access of chromatin regulators, with knock-on effects on gene expression programmes. These changes allow adaptation of the cell to its environment, which, for cancer cells, is a hostile one. We discuss unresolved questions and possible avenues for future exploration.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Intrinsically Disordered Proteins: Critical Components of the Wetware

Despite the wealth of knowledge gained about intrinsically disordered proteins (IDPs) since their discovery, there are several aspects that remain unexplored and, hence, poorly understood. A living cell is a complex adaptive system that can be described as a wetware─a metaphor used to describe the cell as a computer comprising both hardware and software and attuned to logic gates─capable of "making" decisions. In this focused Review, we discuss how IDPs, as critical components of the wetware, influence cell-fate decisions by wiring protein interaction networks to keep them minimally frustrated. Because IDPs lie between order and chaos, we explore the possibility that they can be modeled as attractors. Further, we discuss how the conformational dynamics of IDPs manifests itself as conformational noise, which can potentially amplify transcriptional noise to stochastically switch cellular phenotypes. Finally, we explore the potential role of IDPs in prebiotic evolution, in forming proteinaceous membrane-less organelles, in the origin of multicellularity, and in protein conformation-based transgenerational inheritance of acquired characteristics. Together, these ideas provide a new conceptual framework to discern how IDPs may perform critical biological functions despite their lack of structure.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Probing allosteric interactions in homo-oligomeric molecular machines using solution NMR spectroscopy

Developments in solution NMR spectroscopy have significantly impacted the biological questions that can now be addressed by this methodology. By means of illustration, we present here a perspective focusing on studies of a number of molecular machines that are critical for cellular homeostasis. The role of NMR in elucidating the structural dynamics of these important molecules is emphasized, focusing specifically on intersubunit allosteric communication in homo-oligomers. In many biophysical studies of oligomers, allostery is inferred by showing that models specifically including intersubunit communication best fit the data of interest. Ideally, however, experimental studies focusing on one subunit of a multisubunit system would be performed as an important complement to the more traditional bulk measurements in which signals from all components are measured simultaneously. Using an approach whereby asymmetric molecules are prepared in concert with NMR experiments focusing on the structural dynamics of individual protomers, we present examples of how intersubunit allostery can be directly observed in high-molecular-weight protein systems. These examples highlight some of the unique roles of solution NMR spectroscopy in studies of complex biomolecules and emphasize the important synergy between NMR and other atomic resolution biophysical methods.

Folded domain charge properties influence the conformational behavior of disordered tails

Intrinsically disordered proteins and protein regions (IDRs) make up around 30% of the human proteome where they play essential roles in dictating and regulating many core biological processes. While IDRs are often studied as isolated domains, in naturally occurring proteins most IDRs are found adjacent to folded domains, where they exist as either N- or C-terminal tails or as linkers connecting two folded domains. Prior work has shown that charge properties of IDRs can influence their conformational behavior, both in isolation and in the context of folded domains. In contrast, the converse scenario is less well-explored: how do the charge properties of folded domains influence IDR conformational behavior? To answer this question, we combined a large-scale structural bioinformatics analysis with all-atom implicit solvent simulations of both rationally designed and naturally occurring proteins. Our results reveal three key takeaways. Firstly, the relative position and accessibility of charged residues across the surface of a folded domain can dictate IDR conformational behavior, overriding expectations based on net surface charge properties. Secondly, naturally occurring proteins possess multiple charge patches that are physically accessible to local IDRs. Finally, even modest changes in the local electrostatic environment of a folded domain can substantially modulate IDR-folded domain interactions. Taken together, our results suggest that folded domain surfaces can act as local determinants of IDR conformational behavior.

•

Intrinsically disordered regions (IDRs) are mostly found adjacent to folded domains.

•

Here we propose that the folded domain surface properties influence IDR behavior.

•

We combine all-atom simulations and sequence design of IDRs and folded domains.

•

IDR conformational behavior is determined by a complex combination of factors.

•

Folded domains can substantially alter IDR conformational biases.

Subtractive Proteomics and Immuno-informatics Approaches for Multi-peptide Vaccine Prediction Against Klebsiella oxytoca and Validation Through In Silico Expression

Klebsiella oxytoca is a gram-negative bacterium. It is opportunistic in nature and causes hospital acquired infections. Subtractive proteomics and reverse vaccinology approaches were employed to screen out the best proteins for vaccine designing. Whole proteome of K. oxytoca strain ATCC 8724, consisting of 5483 proteins, was used for designing the vaccine. Total 1670 cytotoxic T lymphocyte (CTL) epitope were predicted through NetCTL while 1270 helper T lymphocyte (HTL) epitopes were predicted through IEDB server. The epitopes were screened for non-toxicity, allergenicity, antigenicity and water solubility. After epitope screening 300 CTL and 250 HTL epitopes were submitted to IFN-γ epitope server to predict their Interferon-γ induction response. The selected IFN-γ positive epitopes were tested for their binding affinity with MHCI-DRB1 by MHCPred. The 15 CTL and 13 HTL epitopes were joined by linkers AAY and GPGPG respectively in vaccine construct. Chain C of Pam3CSK4 (PDB ID; 2Z7X) was linked to the vaccine construct as an adjuvant. A 450aa long vaccine construct was submitted to I-TASSER server for 3D structure prediction. Thirteen Linear B cells were predicted by ABCPred server and 10 sets of discontinues epitopes for 3D vaccine structure were predicted by DiscoTope server. The modeled 3D vaccine construct was docked with human Toll-like receptor 2 (PDB ID: 6NIG) by PatchDock. The docked complexes were refined by FireDock. The selected docked complex showed five hydrogen bonds and one salt bridge. The vaccine sequence was reverse transcribed to get nucleotide sequence for In silico cloning. The reverse transcribed sequence strand was cloned in pET28a(+) expression vector. A clone containing 6586 bp was constructed including the 450 bp of query gene sequence.

Intrinsically disordered protein regions and phase separation: sequence determinants of assembly or lack thereof

Intrinsically disordered protein regions (IDRs) - regions that do not fold into a fixed three-dimensional structure but instead exist in a heterogeneous ensemble of conformations - have recently entered mainstream cell biology in the context of liquid-liquid phase separation (LLPS). IDRs are frequently found to be enriched in phase-separated compartments. Due to this observation, the presence of an IDR in a protein is frequently assumed to be diagnostic of its ability to phase separate. In this review, we clarify the role of IDRs in biological assembly and explore the physical principles through which amino acids can confer the attractive molecular interactions that underlie phase separation. While some disordered regions will robustly drive phase separation, many others will not. We emphasize that rather than 'disorder' driving phase separation, multivalency drives phase separation. As such, whether or not a disordered region is capable of driving phase separation will depend on the physical chemistry encoded within its amino acid sequence. Consequently, an in-depth understanding of that physical chemistry is a prerequisite to make informed inferences on how and why an IDR may be involved in phase separation or, more generally, in protein-mediated intermolecular interactions.

Parallel and Sequential Pathways of Molecular Recognition of a Tandem-Repeat Protein and Its Intrinsically Disordered Binding Partner

The Wnt signalling pathway plays an important role in cell proliferation, differentiation, and fate decisions in embryonic development and the maintenance of adult tissues. The twelve armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here, we investigated the interaction between β-catenin and the intrinsically disordered transcription factor TCF7L2, comprising a very long nanomolar-affinity interface of approximately 4800 Å2 that spans ten of the twelve ARM repeats of β-catenin. First, a fluorescence reporter system for the interaction was engineered and used to determine the kinetic rate constants for the association and dissociation. The association kinetics of TCF7L2 and β-catenin were monophasic and rapid (7.3 ± 0.1 × 107 M-1·s-1), whereas dissociation was biphasic and slow (5.7 ± 0.4 × 10-4 s-1, 15.2 ± 2.8 × 10-4 s-1). This reporter system was then combined with site-directed mutagenesis to investigate the striking variability in the conformation adopted by TCF7L2 in the three different crystal structures of the TCF7L2-β-catenin complex. We found that the mutation had very little effect on the association kinetics, indicating that most interactions form after the rate-limiting barrier for association. Mutations of the N- and C-terminal subdomains of TCF7L2 that adopt relatively fixed conformations in the crystal structures had large effects on the dissociation kinetics, whereas the mutation of the labile sub-domain connecting them had negligible effect. These results point to a two-site avidity mechanism of binding with the linker region forming a "fuzzy" complex involving transient contacts that are not site-specific. Strikingly, the two mutations in the N-terminal subdomain that had the largest effects on the dissociation kinetics showed two additional phases, indicating partial flux through an alternative dissociation pathway that is inaccessible to the wild type. The results presented here provide insights into the kinetics of the molecular recognition of a long intrinsically disordered region with an elongated repeat-protein surface, a process found to involve parallel routes with sequential steps in each.

PhosIDP: a web tool to visualize the location of phosphorylation sites in disordered regions

Charge is a key determinant of intrinsically disordered protein (IDP) and intrinsically disordered region (IDR) properties. IDPs and IDRs are enriched in sites of phosphorylation, which alters charge. Visualizing the degree to which phosphorylation modulates the charge profile of a sequence would assist in the functional interpretation of IDPs and IDRs. PhosIDP is a web tool that shows variation of charge and fold propensity upon phosphorylation. In combination with the displayed location of protein domains, the information provided by the web tool can lead to functional inferences for the consequences of phosphorylation. IDRs are components of many proteins that form biological condensates. It is shown that IDR charge, and its modulation by phosphorylation, is more tightly controlled for proteins that are essential for condensate formation than for those present in condensates but inessential.

Novel RyR2 Mutation (G3118R) Is Associated With Autosomal Recessive Ventricular Fibrillation and Sudden Death: Clinical, Functional, and Computational Analysis

Background

The cardiac ryanodine receptor type 2 (RyR2) is a large homotetramer, located in the sarcoplasmic reticulum (SR), which releases Ca²⁺ from the SR during systole. The molecular mechanism underlying Ca²⁺ sensing and gating of the RyR2 channel in health and disease is only partially elucidated. Catecholaminergic Polymorphic Ventricular Tachycardia (CPVT1) is the most prevalent syndrome caused by RyR2 mutations.

Methods and Results

This study involves investigation of a family with 4 cases of ventricular fibrillation and sudden death and physiological tests in HEK 293 cells and normal mode analysis (NMA) computation. We found 4 clinically affected members who were homozygous for a novel RyR2 mutation, G3118R, whereas their heterozygous relatives are asymptomatic. G3118R is located in the periphery of the protein, far from the mutation hotspot regions. HEK293 cells harboring G3118R mutation inhibited Ca²⁺ release in response to increasing doses of caffeine, but decreased the termination threshold for store‐overload‐induced Ca²⁺ release, thus increasing the fractional Ca²⁺ release in response to increasing extracellular Ca²⁺. NMA showed that G3118 affects RyR2 tetramer in a dose‐dependent manner, whereas in the model of homozygous mutant RyR2, the highest entropic values are assigned to the pore and the central regions of the protein.

Conclusions

RyR2 G3118R is related to ventricular fibrillation and sudden death in recessive mode of inheritance and has an effect of gain of function on the protein. Despite a peripheral location, it has an allosteric effect on the stability of central and pore regions in a dose‐effect manner.

The complexity of protein interactions unravelled from structural disorder

The importance of unstructured biology has quickly grown during the last decades accompanying the explosion of the number of experimentally resolved protein structures. The idea that structural disorder might be a novel mechanism of protein interaction is widespread in the literature, although the number of statistically significant structural studies supporting this idea is surprisingly low. At variance with previous works, our conclusions rely exclusively on a large-scale analysis of all the 134337 X-ray crystallographic structures of the Protein Data Bank averaged over clusters of almost identical protein sequences. In this work, we explore the complexity of the organisation of all the interaction interfaces observed when a protein lies in alternative complexes, showing that interfaces progressively add up in a hierarchical way, which is reflected in a logarithmic law for the size of the union of the interface regions on the number of distinct interfaces. We further investigate the connection of this complexity with different measures of structural disorder: the standard missing residues and a new definition, called "soft disorder", that covers all the flexible and structurally amorphous residues of a protein. We show evidences that both the interaction interfaces and the soft disordered regions tend to involve roughly the same amino-acids of the protein, and preliminary results suggesting that soft disorder spots those surface regions where new interfaces are progressively accommodated by complex formation. In fact, our results suggest that structurally disordered regions not only carry crucial information about the location of alternative interfaces within complexes, but also about the order of the assembly. We verify these hypotheses in several examples, such as the DNA binding domains of P53 and P73, the C3 exoenzyme, and two known biological orders of assembly. We finally compare our measures of structural disorder with several disorder bioinformatics predictors, showing that these latter are optimised to predict the residues that are missing in all the alternative structures of a protein and they are not able to catch the progressive evolution of the disordered regions upon complex formation. Yet, the predicted residues, when not missing, tend to be characterised as soft disordered regions.

An analytical theory to describe sequence-specific inter-residue distance profiles for polyampholytes and intrinsically disordered proteins

Intrinsically Disordered Proteins (IDPs), unlike folded proteins, lack a unique folded structure and rapidly interconvert among ensembles of disordered states. However, they have specific conformational properties when averaged over their ensembles of disordered states. It is critical to develop a theoretical formalism to predict these ensemble average conformational properties that are encoded in the IDP sequence (the specific order in which amino acids/residues are linked). We present a general heteropolymer theory that analytically computes the ensemble average distance profiles (⟨R_ij ²⟩) between any two (i, j) monomers (amino acids for IDPs) as a function of the sequence. Information rich distance profiles provide a detailed description of the IDP in contrast to typical metrics such as scaling exponents, radius of gyration, or end-to-end distance. This generalized formalism supersedes homopolymer-like models or models that are built only on the composition of amino acids but ignore sequence details. The prediction of these distance profiles for highly charged polyampholytes and naturally occurring IDPs unmasks salient features that are hidden in the sequence. Moreover, the model reveals strategies to modulate the entire distance map to achieve local or global swelling/compaction by subtle changes/modifications-such as phosphorylation, a biologically relevant process-in specific hotspots in the sequence. Sequence-specific distance profiles and their modulation have been benchmarked against all-atom simulations. Our new formalism also predicts residue-pair specific coil-globule transitions. The analytical nature of the theory will facilitate design of new sequences to achieve specific target distance profiles with broad applications in synthetic biology and polymer science.

Polyelectrolyte interactions enable rapid association and dissociation in high-affinity disordered protein complexes

Highly charged intrinsically disordered proteins can form complexes with very high affinity in which both binding partners fully retain their disorder and dynamics, exemplified by the positively charged linker histone H1.0 and its chaperone, the negatively charged prothymosin α. Their interaction exhibits another surprising feature: The association/dissociation kinetics switch from slow two-state-like exchange at low protein concentrations to fast exchange at higher, physiologically relevant concentrations. Here we show that this change in mechanism can be explained by the formation of transient ternary complexes favored at high protein concentrations that accelerate the exchange between bound and unbound populations by orders of magnitude. Molecular simulations show how the extreme disorder in such polyelectrolyte complexes facilitates (i) diffusion-limited binding, (ii) transient ternary complex formation, and (iii) fast exchange of monomers by competitive substitution, which together enable rapid kinetics. Biological polyelectrolytes thus have the potential to keep regulatory networks highly responsive even for interactions with extremely high affinities.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Selective Modulation of Dynamic Protein Complexes

Dynamic proteins perform critical roles in cellular machines, including those that control proteostasis, transcription, translation, and signaling. Thus, dynamic proteins are prime candidates for chemical probe and drug discovery but difficult targets because they do not conform to classical rules of design and screening. Selectivity is pivotal for candidate probe molecules due to the extensive interaction network of these dynamic hubs. Recognition that the traditional rules of probe discovery are not necessarily applicable to dynamic proteins and their complexes, as well as technological advances in screening, have produced remarkable results in the last 2-4 years. Particularly notable are the improvements in target selectivity for small-molecule modulators of dynamic proteins, especially with techniques that increase the discovery likelihood of allosteric regulatory mechanisms. We focus on approaches to small-molecule screening that appear to be more suitable for highly dynamic targets and have the potential to streamline identification of selective modulators.

RORγ Structural Plasticity and Druggability

Retinoic acid receptor-related orphan receptor γ (RORγ) is a transcription factor regulating the expression of the pro-inflammatory cytokine IL-17 in human T helper 17 (Th17) cells. Activating RORγ can induce multiple IL-17-mediated autoimmune diseases but may also be useful for anticancer therapy. Its deep immunological functions make RORɣ an attractive drug target. Over 100 crystal structures have been published describing atomic interactions between RORɣ and agonists and inverse agonists. In this review, we focus on the role of dynamic properties and plasticity of the RORɣ orthosteric and allosteric binding sites by examining structural information from crystal structures and simulated models. We discuss the possible influences of allosteric ligands on the orthosteric binding site. We find that high structural plasticity favors the druggability of RORɣ, especially for allosteric ligands.

Multi-state recognition pathway of the intrinsically disordered protein kinase inhibitor by protein kinase A

In the nucleus, the spatiotemporal regulation of the catalytic subunit of cAMP-dependent protein kinase A (PKA-C) is orchestrated by an intrinsically disordered protein kinase inhibitor, PKI, which recruits the CRM1/RanGTP nuclear exporting complex. How the PKA-C/PKI complex assembles and recognizes CRM1/RanGTP is not well understood. Using NMR, SAXS, fluorescence, metadynamics, and Markov model analysis, we determined the multi-state recognition pathway for PKI. After a fast binding step in which PKA-C selects PKI's most competent conformations, PKI folds upon binding through a slow conformational rearrangement within the enzyme's binding pocket. The high-affinity and pseudo-substrate regions of PKI become more structured and the transient interactions with the kinase augment the helical content of the nuclear export sequence, which is then poised to recruit the CRM1/RanGTP complex for nuclear translocation. The multistate binding mechanism featured by PKA-C/PKI complex represents a paradigm on how disordered, ancillary proteins (or protein domains) are able to operate multiple functions such as inhibiting the kinase while recruiting other regulatory proteins for nuclear export.

Strategies for identifying dynamic regions in protein complexes: Flexibility changes accompany methylation in chemotaxis receptor signaling states

Bacterial chemoreceptors are organized in arrays composed of helical receptors arranged as trimers of dimers, coupled to a histidine kinase CheA and a coupling protein CheW. Ligand binding to the external domain inhibits the kinase activity, leading to a change in the swimming behavior. Adaptation to an ongoing stimulus involves reversible methylation and demethylation of specific glutamate residues. However, the exact mechanism of signal propagation through the helical receptor to the histidine kinase remains elusive. Dynamics of the receptor cytoplasmic domain is thought to play an important role in the signal transduction, and current models propose inverse dynamic changes in different regions of the receptor. We hypothesize that the adaptational modification (methylation) controls the dynamics by stabilizing a partially ordered domain, which in turn modulates the binding of the kinase, CheA. We investigated the difference in dynamics between the methylated and unmethylated states of the chemoreceptor using solid-state NMR. The unmethylated receptor (CF4E) shows increased flexibility relative to the methylated mimic (CF4Q). Methylation helix 1 (MH1) has been shown to be flexible in the methylated mimic receptor. Our analysis indicates that in addition to MH1, methylation helix 2 also becomes flexible in the unmethylated receptor. In addition, we have demonstrated that both states of the receptor have a rigid region and segments with intermediate timescale dynamics. The strategies used in this study for identifying dynamic regions are applicable to a broad class of proteins and protein complexes with intrinsic disorder and dynamics spanning multiple timescales.

Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain

Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their "fuzzy" nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy "hotspots" as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.

HK97 gp74 Possesses an α-Helical Insertion in the ββα Fold That Affects Its Metal Binding, cos Site Digestion, and In Vivo Activities

The last gene in the genome of the bacteriophage HK97 encodes gp74, an HNH endonuclease. HNH motifs contain two conserved His residues and an invariant Asn residue, and they adopt a ββα structure. gp74 is essential for phage head morphogenesis, likely because gp74 enhances the specific endonuclease activity of the HK97 terminase complex. Notably, the ability of gp74 to enhance the terminase-mediated cleavage of the phage cos site requires an intact HNH motif in gp74. Mutation of H82, the conserved metal-binding His residue in the HNH motif, to Ala abrogates gp74-mediated stimulation of terminase activity. Here, we present nuclear magnetic resonance (NMR) studies demonstrating that gp74 contains an α-helical insertion in the Ω-loop, which connects the two β-strands of the ββα fold, and a disordered C-terminal tail. NMR data indicate that the Ω-loop insert makes contacts to the ββα fold and influences the ability of gp74 to bind divalent metal ions. Further, the Ω-loop insert and C-terminal tail contribute to gp74-mediated DNA digestion and to gp74 activity in phage morphogenesis. The data presented here enrich our molecular-level understanding of how HNH endonucleases enhance terminase-mediated digestion of the cos site and contribute to the phage replication cycle.IMPORTANCE This study demonstrates that residues outside the canonical ββα fold, namely, the Ω-loop α-helical insert and a disordered C-terminal tail, regulate the activity of the HNH endonuclease gp74. The increased divalent metal ion binding when the Ω-loop insert is removed compared to reduced cos site digestion and phage formation indicates that the Ω-loop insert plays multiple regulatory roles. The data presented here provide insights into the molecular basis of the involvement of HNH proteins in phage DNA packing.

Untangling the Conformational Polymorphism of Disordered Proteins Associated With Neurodegeneration at the Single-Molecule Level

A large fraction of the human genome encodes intrinsically disordered proteins/regions (IDPs/IDRs) that are involved in diverse cellular functions/regulation and dysfunctions. Moreover, several neurodegenerative disorders are associated with the pathological self-assembly of neuronal IDPs, including tau [Alzheimer's disease (AD)], α-synuclein [Parkinson's disease (PD)], and huntingtin exon 1 [Huntington's disease (HD)]. Therefore, there is an urgent and emerging clinical interest in understanding the physical and structural features of their functional and disease states. However, their biophysical characterization is inherently challenging by traditional ensemble techniques. First, unlike globular proteins, IDPs lack stable secondary/tertiary structures under physiological conditions and may interact with multiple and distinct biological partners, subsequently folding differentially, thus contributing to the conformational polymorphism. Second, amyloidogenic IDPs display a high aggregation propensity, undergoing complex heterogeneous self-assembly mechanisms. In this review article, we discuss the advantages of employing cutting-edge single-molecule fluorescence (SMF) techniques to characterize the conformational ensemble of three selected neuronal IDPs (huntingtin exon 1, tau, and α-synuclein). Specifically, we survey the versatility of these powerful approaches to describe their monomeric conformational ensemble under functional and aggregation-prone conditions, and binding to biological partners. Together, the information gained from these studies provides unique insights into the role of gain or loss of function of these disordered proteins in neurodegeneration, which may assist the development of new therapeutic molecules to prevent and treat these devastating human disorders.

Computational Disorder Analysis in Ethylene Response Factors Uncovers Binding Motifs Critical to Their Diverse Functions

APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein–protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.

Supramolecular Fuzziness of Intracellular Liquid Droplets: Liquid-Liquid Phase Transitions, Membrane-Less Organelles, and Intrinsic Disorder

Cells are inhomogeneously crowded, possessing a wide range of intracellular liquid droplets abundantly present in the cytoplasm of eukaryotic and bacterial cells, in the mitochondrial matrix and nucleoplasm of eukaryotes, and in the chloroplast’s stroma of plant cells. These proteinaceous membrane-less organelles (PMLOs) not only represent a natural method of intracellular compartmentalization, which is crucial for successful execution of various biological functions, but also serve as important means for the processing of local information and rapid response to the fluctuations in environmental conditions. Since PMLOs, being complex macromolecular assemblages, possess many characteristic features of liquids, they represent highly dynamic (or fuzzy) protein–protein and/or protein–nucleic acid complexes. The biogenesis of PMLOs is controlled by specific intrinsically disordered proteins (IDPs) and hybrid proteins with ordered domains and intrinsically disordered protein regions (IDPRs), which, due to their highly dynamic structures and ability to facilitate multivalent interactions, serve as indispensable drivers of the biological liquid–liquid phase transitions (LLPTs) giving rise to PMLOs. In this article, the importance of the disorder-based supramolecular fuzziness for LLPTs and PMLO biogenesis is discussed.

q-Canonical Monte Carlo Sampling for Modeling the Linkage between Charge Regulation and Conformational Equilibria of Peptides

The overall charge content and the patterning of charged residues have a profound impact on the conformational ensembles adopted by intrinsically disordered proteins. These parameters can be altered by charge regulation, which refers to the effects of post-translational modifications, pH-dependent changes to charge, and conformational fluctuations that modify the pKa values of ionizable residues. Although atomistic simulations have played a prominent role in uncovering the major sequence-ensemble relationships of IDPs, most simulations assume fixed charge states for ionizable residues. This may lead to erroneous estimates for conformational equilibria if they are linked to charge regulation. Here, we report the development of a new method we term q-canonical Monte Carlo sampling for modeling the linkage between charge regulation and conformational equilibria. The method, which is designed to be interoperable with the ABSINTH implicit solvation model, operates as follows: For a protein sequence with n ionizable residues, we start with all 2n charge microstates and use a criterion based on model compound pKa values to prune down to a subset of thermodynamically relevant charge microstates. This subset is then grouped into mesostates, where all microstates that belong to a mesostate have the same net charge. Conformational distributions, drawn from a canonical ensemble, are generated for each of the charge microstates that make up a mesostate using a method we designate as proton walk sampling. This method combines Metropolis Monte Carlo sampling in conformational space with an auxiliary Markov process that enables interconversions between charge microstates along a mesostate. Proton walk sampling helps identify the most likely charge microstate per mesostate. We then use thermodynamic integration aided by the multistate Bennett acceptance ratio method to estimate the free energies for converting between mesostates. These free energies are then combined with the per-microstate weights along each mesostate to estimate standard state free energies and pH-dependent free energies for all thermodynamically relevant charge microstates. The results provide quantitative estimates of the probabilities and preferred conformations associated with every thermodynamically accessible charge microstate. We showcase the application of q-canonical sampling using two model systems. The results establish the soundness of the method and the importance of charge regulation in systems characterized by conformational heterogeneity.

Entropy, Fluctuations, and Disordered Proteins

Entropy should directly reflect the extent of disorder in proteins. By clustering structurally related proteins and studying the multiple-sequence-alignment of the sequences of these clusters, we were able to link between sequence, structure, and disorder information. We introduced several parameters as measures of fluctuations at a given MSA site and used these as representative of the sequence and structure entropy at that site. In general, we found a tendency for negative correlations between disorder and structure, and significant positive correlations between disorder and the fluctuations in the system. We also found evidence for residue-type conservation for those residues proximate to potentially disordered sites. Mutation at the disorder site itself appear to be allowed. In addition, we found positive correlation for disorder and accessible surface area, validating that disordered residues occur in exposed regions of proteins. Finally, we also found that fluctuations in the dihedral angles at the original mutated residue and disorder are positively correlated while dihedral angle fluctuations in spatially proximal residues are negatively correlated with disorder. Our results seem to indicate permissible variability in the disordered site, but greater rigidity in the parts of the protein with which the disordered site interacts. This is another indication that disordered residues are involved in protein function.

Entropy and Information within Intrinsically Disordered Protein Regions

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Hybrid Refinement of Heterogeneous Conformational Ensembles Using Spectroscopic Data

Multistructured biomolecular systems play crucial roles in a wide variety of cellular processes but have resisted traditional methods of structure determination, which often resolve only a few low-energy states. High-resolution structure determination using experimental methods that yield distributional data remains extremely difficult, especially when the underlying conformational ensembles are quite heterogeneous. We have therefore developed a method to integrate sparse, multimultimodal spectroscopic data to obtain high-resolution estimates of conformational ensembles. We have tested our method by incorporating double electron-electron resonance data on the soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) protein syntaxin-1a into biased molecular dynamics simulations. We find that our method substantially outperforms existing state-of-the-art methods in capturing syntaxin's open-closed conformational equilibrium and further yields new conformational states that are consistent with experimental data and may help in understanding syntaxin's function. Our improved methods for refining heterogeneous conformational ensembles from spectroscopic data will greatly accelerate the structural understanding of such systems.

Intrinsic Disorder-Based Emergence in Cellular Biology: Physiological and Pathological Liquid-Liquid Phase Transitions in Cells

The visible outcome of liquid-liquid phase transitions (LLPTs) in cells is the formation and disintegration of various proteinaceous membrane-less organelles (PMLOs). Although LLPTs and related PMLOs have been observed in living cells for over 200 years, the physiological functions of these transitions (also known as liquid-liquid phase separation, LLPS) are just starting to be understood. While unveiling the functionality of these transitions is important, they have come into light more recently due to the association of abnormal LLPTs with various pathological conditions. In fact, several maladies, such as various cancers, different neurodegenerative diseases, and cardiovascular diseases, are known to be associated with either aberrant LLPTs or some pathological transformations within the resultant PMLOs. Here, we will highlight both the physiological functions of cellular liquid-liquid phase transitions as well as the pathological consequences produced through both dysregulated biogenesis of PMLOs and the loss of their dynamics. We will also discuss the potential downstream toxic effects of proteins that are involved in pathological formations.

Stochastic modeling of macromolecules in solution. I. Relaxation processes

A framework for the stochastic description of relaxation processes in flexible macromolecules, including dissipative effects, is introduced from an atomistic point of view. Projection-operator techniques are employed to obtain multidimensional Fokker-Planck operators governing the relaxation of internal coordinates and global degrees of freedom and depending upon parameters fully recoverable from classic force fields (energetics) and continuum models (friction tensors). A hierarchy of approaches of different complexity is proposed in this unified context, aimed primarily at the interpretation of magnetic resonance relaxation experiments. In particular, a model based on a harmonic internal Hamiltonian is discussed as a test case.