Intrinsically disordered proteins: regulation and disease

Deciphering RNA-Recognition Patterns of Intrinsically Disordered Proteins

Intrinsically disordered regions (IDRs) and protein (IDPs) are highly flexible owing to their lack of well-defined structures. A subset of such proteins interacts with various substrates; including RNA; frequently adopting regular structures in the final complex. In this work; we have analysed a dataset of protein–RNA complexes undergoing disorder-to-order transition (DOT) upon binding. We found that DOT regions are generally small in size (less than 3 residues) for RNA binding proteins. Like structured proteins; positively charged residues are found to interact with RNA molecules; indicating the dominance of electrostatic and cation-π interactions. However, a comparison of binding frequency shows that interface hydrophobic and aromatic residues have more interactions in only DOT regions than in a protein. Further; DOT regions have significantly higher exposure to water than their structured counterparts. Interactions of DOT regions with RNA increase the sheet formation with minor changes in helix forming residues. We have computed the interaction energy for amino acids–nucleotide pairs; which showed the preference of His–G; Asn–U and Ser–U at for the interface of DOT regions. This study provides insights to understand protein–RNA interactions and the results could also be used for developing a tool for identifying DOT regions in RNA binding proteins.

Dynamics, Conformational Entropy, and Frustration in Protein-Protein Interactions Involving an Intrinsically Disordered Protein Domain

Intrinsically disordered proteins (IDPs) are abundant in the eukaryotic proteome. However, little is known about the role of subnanosecond dynamics and the conformational entropy that it represents in protein-protein interactions involving IDPs. Using nuclear magnetic resonance side chain and backbone relaxation, stopped-flow kinetics, isothermal titration calorimetry, and computational studies, we have characterized the interaction between the globular TAZ1 domain of the CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2). We show that the TAZ1/TAD-STAT2 complex retains considerable subnanosecond motions, with TAD-STAT2 undergoing only a partial disorder-to-order transition. We report here the first experimental determination of the conformational entropy change for both binding partners in an IDP binding interaction and find that the total change even exceeds in magnitude the binding enthalpy and is comparable to the contribution from the hydrophobic effect, demonstrating its importance in the binding energetics. Furthermore, we show that the conformational entropy change for TAZ1 is also instrumental in maintaining a biologically meaningful binding affinity. Strikingly, a spatial clustering of very high amplitude motions and a cluster of more rigid sites in the complex exist, which through computational studies we found to overlap with regions that experience energetic frustration and are less frustrated, respectively. Thus, the residual dynamics in the bound state could be necessary for faster dissociation, which is important for proteins that interact with multiple binding partners.

Function and Regulation of Phase-Separated Biological Condensates

Achieving functional specificity while minimizing cost to fitness is a key constraint during evolution. Formation of biological condensates by liquid-liquid phase separation (LLPS) appears to serve as an important regulatory mechanism to generate moderate specificity in molecular recognition while maintaining a reasonable cost for fitness in terms of design complexity. Formation of biological condensates serves as a unique mechanism of molecular recognition achieving some level of specificity without a huge cost to fitness. Rapid formation of biological condensates in vivo induced by specific cellular or environmental triggers has been shown to be an important mechanism for increasing cellular fitness. Here we discuss the functions and regulation of biological condensates, especially those formed by LLPS, involving interactions between proteins and nucleic acids. These condensates are spatially isolated within the cytosol or nucleus and can facilitate specific biochemical functions under conditions such as stress. The misregulation of biological condensates resulting in nondynamic aggregates has been implicated in a number of diseases. Understanding the functional importance of biological condensates and their regulation opens doors for development of therapies targeting dysfunctional biological condensates, as well as spatiotemporal engineering of functions in cells.

Convergence of Artificial Protein Polymers and Intrinsically Disordered Proteins

A flurry of research in recent years has revealed the molecular origins of many membraneless organelles to be the liquid phase separation of intrinsically disordered proteins (IDPs). Consequently, protein disorder has emerged as an important driver of intracellular compartmentalization by providing specialized microenvironments chemically distinct from the surrounding medium. Though the importance of protein disorder and its relationship to intracellular phase behavior are clear, a detailed understanding of how such phase behavior can be predicted and controlled remains elusive. While research in IDPs has largely focused on the implications of structural disorder on cellular function and disease, another field, that of artificial protein polymers, has focused on the de novo design of protein polymers with controllable material properties. A subset of these polymers, specifically those derived from structural proteins such as elastin and resilin, are also disordered sequences that undergo liquid-liquid phase separation. This phase separation has been used in a variety of biomedical applications, and researchers studying these polymers have developed methods to precisely characterize and tune their phase behavior. Despite their disparate origins, both fields are complementary as they study the phase behavior of intrinsically disordered polypeptides. This Perspective hopes to stimulate collaborative efforts by highlighting the similarities between these two fields and by providing examples of how such collaboration could be mutually beneficial.

The impact of O-glycan chemistry on the stability of intrinsically disordered proteins

Protein glycosylation is a diverse post-translational modification that serves myriad biological functions. O-linked glycans in particular vary widely in extent and chemistry in eukaryotes, with secreted proteins from fungi and yeast commonly exhibiting O-mannosylation in intrinsically disordered regions of proteins, likely for proteolysis protection, among other functions. However, it is not well understood why mannose is often the preferred glycan, and more generally, if the neighboring protein sequence and glycan have coevolved to protect against proteolysis in glycosylated intrinsically disordered proteins (IDPs). Here, we synthesized variants of a model IDP, specifically a natively O-mannosylated linker from a fungal enzyme, with α-O-linked mannose, glucose, and galactose moieties, along with a non-glycosylated linker. Upon exposure to thermolysin, O-mannosylation, by far, provides the highest extent of proteolysis protection. To explain this observation, extensive molecular dynamics simulations were conducted, revealing that the axial configuration of the C2-hydroxyl group (2-OH) of α-mannose adjacent to the glycan-peptide bond strongly influences the conformational features of the linker. Specifically, α-mannose restricts the torsions of the IDP main chain more than other glycans whose equatorial 2-OH groups exhibit interactions that favor perpendicular glycan-protein backbone orientation. We suggest that IDP stiffening due to O-mannosylation impairs protease action, with contributions from protein-glycan interactions, protein flexibility, and protein stability. Our results further imply that resistance to proteolysis is an important driving force for evolutionary selection of α-mannose in eukaryotic IDPs, and more broadly, that glycan motifs for proteolysis protection likely coevolve with the protein sequence to which they attach.

Life and Work of Stress Granules and Processing Bodies: New Insights into Their Formation and Function

The dynamic formation of stress granules (SGs), processing bodies (PBs), and related RNA organelles regulates diverse cellular processes, including the coordination of functionally connected messengers, the translational regulation at the synapse, and the control of viruses and retrotransposons. Recent studies have shown that pyruvate kinase and other enzymes localize in SGs and PBs, where they become protected from stress insults. These observations may have implications for enzyme regulation and metabolic control exerted by RNA-based organelles. The formation of these cellular bodies is governed by liquid-liquid phase separation (LLPS) processes, and it needs to be strictly controlled to prevent pathogenic aggregation. The intracellular concentration of key metabolites, such as ATP and sterol derivatives, may influence protein solubility, thus affecting the dynamics of liquid organelles. LLPS in vitro depends on the thermal diffusion of macromolecules, which is limited inside cells, where the condensation and dissolution of membrane-less organelles are helped by energy-driven processes. The active transport by the retrograde motor dynein helps SG assembly, whereas the anterograde motor kinesin mediates SG dissolution; a tug of war between these two molecular motors allows transient SG formation. There is evidence that the efficiency of dynein-mediated transport increases with the number of motor molecules associated with the cargo. The dynein-dependent transport may be influenced by cargo size as larger cargos can load a larger number of motors. We propose a model based on this emergent property of dynein motors, which would be collectively stronger during SG condensation and weaker during SG breakdown, thus allowing kinesin-mediated dispersion.

A co-expression strategy to achieve labeling of individual subunits within a dimeric protein for single molecule analysis

A generic co-expression strategy for site-specific incorporation of a single donor-acceptor dye pair into any position in a dimeric protein, allowing single molecule FRET study of proteins previously inaccessible to this technique, such as the intrinsically disordered prion N-domain of Ure2 in the context of its globular dimeric C-domain.

Structural characterization of the intrinsically disordered domain of Mycobacterium tuberculosis protein tyrosine kinase A

Although intrinsically disordered proteins or protein domains (IDPs or IDD) are less abundant in bacteria than in eukaryotes, their presence in pathogenic bacterial proteins is important for protein-protein interactions. The protein tyrosine kinase A (PtkA) from Mycobacterium tuberculosis possesses an 80-residue disordered region (IDD_PtkA ) of unknown function, located N-terminally to the well-folded kinase core domain. Here, we characterize the conformation of IDD_PtkA under varying biophysical conditions and phosphorylation using NMR-spectroscopy. Our results confirm that the N-terminal domain of PtkA exists as an IDD at physiological pH. Furthermore, phosphorylation of IDD_PtkA increases the activity of PtkA. Our findings will complement future approaches in understanding molecular mechanisms of key proteins in pathogenic virulence.

On the Calculation of SAXS Profiles of Folded and Intrinsically Disordered Proteins from Computer Simulations

Solution techniques such as small-angle X-ray scattering (SAXS) play a central role in structural studies of intrinsically disordered proteins (IDPs); yet, due to low resolution, it is generally necessary to combine SAXS with additional experimental sources of data and to use molecular simulations. Computational methods for the calculation of theoretical SAXS intensity profiles can be separated into two groups, depending on whether the solvent is modeled implicitly as continuous electron density or considered explicitly. The former offers reduced computational cost but requires the definition of a number of free parameters to account for, for example, the excess density of the solvation layer. Overfitting can thus be an issue, particularly when the structural ensemble is unknown. Here, we investigate and show how small variations of the contrast of the hydration shell, δρ, severely affect the outcome, analysis and interpretation of computed SAXS profiles for folded and disordered proteins. For both the folded and disordered proteins studied here, using a default δρ may, in some cases, result in the calculation of non-representative SAXS profiles, leading to an overestimation of their size and a misinterpretation of their structural nature. The solvation layer of the different IDP simulations also impacts their size estimates differently, depending on the protein force field used. The same is not true for the folded protein simulations, suggesting differences in the solvation of the two classes of proteins, and indicating that different force fields optimized for IDPs may cause expansion of the polypeptide chain through different physical mechanisms.

Glycation induces conformational changes in the amyloid-β peptide and enhances its aggregation propensity: molecular insights

The cytotoxicity of the amyloid beta (Aβ) peptide, implicated in the pathogenesis of Alzheimer's disease (AD), can be enhanced by its post-translational glycation, a series of non-enzymatic reactions with reducing sugars and reactive dicarbonyls. However, little is known about the underlying mechanisms that potentially enhance the cytotoxicity of the advanced glycation modified Aβ. In this work, fully atomistic molecular dynamics (MD) simulations are exploited to obtain direct molecular insights into the process of early Aβ self-assembly in the presence and absence of glycated lysine residues. Analyses of data exceeding cumulative timescales of 1 microsecond for each system reveal that glycation results in a stronger enthalpy of association between Aβ monomers and lower conformational entropy, in addition to a sharp overall increase in the beta-sheet content. Further analyses reveal that the enhanced interactions originate, in large part, due to markedly stronger, as well as new, inter-monomer salt bridging propensities in the glycated variety. Interestingly, these conformational and energetic effects are broadly reflected in preformed protofibrillar forms of Aβ small oligomers modified with glycation. Our combined results imply that glycation consolidates Aβ self-assembly regardless of its point of occurrence in the pathway. They provide a basis for further mechanistic studies and therapeutic endeavors that could potentially result in novel ways of combating AGE related AD progression.

Computational modeling suggests impaired interactions between NKX2.5 and GATA4 in individuals carrying a novel pathogenic D16N NKX2.5 mutation

NKX2.5, a homeobox containing gene, plays an important role in embryonic heart development and associated mutations are linked with various cardiac abnormalities. We sequenced the NKX2.5 gene in 100 congenital heart disease (CHD) patients and 200 controls. Our analysis revealed a total of 7 mutations, 3 in intronic region, 3 in coding region and 1 in 3’ UTR. Of the above mutations, one mutation was found to be associated with tetralogy of fallot (TOF) and two (rs2277923 and a novel mutation, D16N) were strongly associated with VSD. A novel missense mutation, D16N (p-value =0.009744), located in the tinman (TN) region and associated with ventricular septal defect (VSD), is the most significant findings of this study. Computational analysis revealed that D16N mutation is pathogenic in nature. Through the molecular modeling, docking and molecular dynamics simulation studies, we have identified the location of mutant D16N in NKX2.5 and its interaction map with other partners at the atomic level. We found NKX2.5-GATA4 complex is stable, however, in case of mutant we observed significant conformational changes and loss of key polar interactions, which might be a cause of the pathogenic behavior. This study underscores the structural basis of D16N pathogenic mutation in the regulation of NKX2.5 and how this mutation renders the structural-functional divergence that possibly leading towards the diseased state.

The intrinsically disordered Tarp protein from chlamydia binds actin with a partially preformed helix

Tarp (translocated actin recruiting phosphoprotein) is an effector protein common to all chlamydial species that functions to remodel the host-actin cytoskeleton during the initial stage of infection. In C. trachomatis, direct binding to actin monomers has been broadly mapped to a 100-residue region (726-825) which is predicted to be predominantly disordered, with the exception of a ~10-residue α-helical patch homologous to other WH2 actin-binding motifs. Biophysical investigations demonstrate that a Tarp726-825 construct behaves as a typical intrinsically disordered protein; within it, NMR relaxation measurements and chemical shift analysis identify the ten residue WH2-homologous region to exhibit partial α-helix formation. Isothermal titration calorimetry experiments on the same construct in the presence of monomeric G-actin show a well defined binding event with a 1:1 stoichiometry and Kd of 102 nM, whilst synchrotron radiation circular dichroism spectroscopy suggests the binding is concomitant with an increase in helical secondary structure. Furthermore, NMR experiments in the presence of G-actin indicate this interaction affects the proposed WH2-like α-helical region, supporting results from in silico docking calculations which suggest that, when folded, this α-helix binds within the actin hydrophobic cleft as seen for other actin-associated proteins.

Common molecular pathogenesis of disease-related intrinsically disordered proteins revealed by NMR analysis

Intrinsically disordered proteins (IDPs) are either completely unstructured or contain large disordered regions in their native state; they have drawn much attention in the field of molecular pathology. Some of them substantially tend to form protein self-assemblies, such as toxic or non-toxic aggregates and fibrils, and have been postulated to relate to diseases. These disease-related IDPs include Aβ(1-42) [Alzheimer's disease (AD)], Tau (AD and tauopathy), α-synuclein (Parkinson's disease) and p53 (cancer). Several studies suggest that these aggregation and/or fibril formation processes are often initiated by transient conformational changes of the IDPs prior to protein self-assembly. Interestingly, the pathological molecular processes of these IDPs share multiple common features with those of protein misfolding diseases, such as transmissible spongiform encephalopathy (PrPsc) and AL-amyloidosis (VL-domain of γ-immunoglobulin). This review provides an overview of solution NMR techniques that can help analyse the early and transient events of conformational equilibrium of IDPs and folded proteins.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

IDPs in macromolecular complexes: the roles of multivalent interactions in diverse assemblies

Intrinsically disordered proteins (IDPs) have critical roles in a diverse array of cellular functions. Of relevance here is that they are components of macromolecular complexes, where their conformational flexibility helps mediate interactions with binding partners. IDPs often interact with their binding partners through short sequence motifs, commonly repeated within the disordered regions. As such, multivalent interactions are common for IDPs and their binding partners within macromolecular complexes. Here we discuss the importance of IDP multivalency in three very different macromolecular assemblies: biomolecular condensates, the nuclear pore, and the cytoskeleton.

Biotinylated Cell-penetrating Peptides to Study Intracellular Protein-protein Interactions

Here we present a protocol to study intracellular protein-protein interactions that is based on the widely used biotin-avidin pull-down system. The modification presented includes the combination of this technique with cell-penetrating sequences. We propose to design cell-penetrating baits that can be incubated with living cells instead of cell lysates and therefore the interactions found will reflect those that occur within the intracellular context. Connexin43 (Cx43), a protein that forms gap junction channels and hemichannels is down-regulated in high-grade gliomas. The Cx43 region comprising amino acids 266-283 is responsible for the inhibition of the oncogenic activity of c-Src in glioma cells. Here we use TAT as the cell-penetrating sequence, biotin as the pull-down tag and the region of Cx43 comprised between amino acids 266-283 as the target to find intracellular interactions in the hard-to-transfect human glioma stem cells. One of the limitations of the proposed method is that the molecule used as bait could fail to fold properly and, consequently, the interactions found could not be associated with the effect. However, this method can be especially interesting for the interactions involved in signal transduction pathways because they are usually carried out by intrinsically disordered regions and, therefore, they do not require an ordered folding. In addition, one of the advantages of the proposed method is that the relevance of each residue on the interaction can be easily studied. This is a modular system; therefore, other cell-penetrating sequences, other tags, and other intracellular targets can be employed. Finally, the scope of this protocol is far beyond protein-protein interaction because this system can be applied to other bioactive cargoes such as RNA sequences, nanoparticles, viruses or any molecule that can be transduced with cell-penetrating sequences and fused to pull-down tags to study their intracellular mechanism of action.

HyRes: a coarse-grained model for multi-scale enhanced sampling of disordered protein conformations

Efficient coarse-grained (CG) models can be coupled with atomistic force fields to accelerate the sampling of atomistic energy landscapes in the multi-scale enhanced sampling (MSES) framework. This approach may be particularly suitable for generating atomistic conformational ensembles of intrinsically disordered proteins (IDPs). While MSES is relatively robust to inherent CG artifacts, achieving optimal sampling efficiency requires CG modeling to generate the local and long-range fluctuations that are largely consistent with those at the atomistic level. Here, we describe a new hybrid resolution CG model (HyRes) for MSES simulations of disordered protein states, which is specifically designed to provide semi-quantitative secondary structure propensities together with a qualitative description of long-range nonspecific interactions. The HyRes model contains an atomistic description of the backbone with intermediate resolution side chains. The secondary structure propensities are tuned by adjusting the backbone hydrogen-bonding strength and the ϕ/ψ torsion profile. The sizes and covalent geometries of the side chains are parameterized to reproduce distributions derived from atomistic simulations. Lennard-Jones parameters for sidechain beads are assigned to reproduce statistical potentials derived from the protein structural database, and then globally parameterized with nonspecific electrostatic interactions to reproduce the free energy profiles of pair wise interactions and the key conformational properties of model peptides. Application of HyRes to MSES simulations of small IDPs suggests that it is capable of driving faster structural transitions at the atomistic level and increasing the convergence rate compared to the Cα-only Gō-like models previously utilized. With further optimization, we believe that the new CG model could greatly improve the efficiency of MSES simulations of the larger and more complex IDPs frequently involved in cellular signalling and regulation.

Long disordered regions of the C-terminal domain of Abelson tyrosine kinase have specific and additive functions in regulation and axon localization

Abelson tyrosine kinase (Abl) is a key regulator of actin-related morphogenetic processes including axon guidance, where it functions downstream of several guidance receptors. While the long C-terminal domain (CTD) of Abl is required for function, its role is poorly understood. Here, a battery of mutants of Drosophila Abl was created that systematically deleted large segments of the CTD from Abl or added them back to the N-terminus alone. The functionality of these Abl transgenes was assessed through rescue of axon guidance defects and adult lethality in Abl loss-of-function, as well as through gain-of-function effects in sensitized slit or frazzled backgrounds that perturb midline guidance in the Drosophila embryonic nerve cord. Two regions of the CTD play important and distinct roles, but additive effects for other regions were also detected. The first quarter of the CTD, including a conserved PxxP motif and its surrounding sequence, regulates Abl function while the third quarter localizes Abl to axons. These regions feature long stretches of intrinsically disordered sequence typically found in hub proteins and are associated with diverse protein-protein interactions. Thus, the CTD of Abl appears to use these disordered regions to establish a variety of different signaling complexes required during formation of axon tracts.

Intrinsic protein disorder reduces small-scale gene duplicability

Whereas the rate of gene duplication is relatively high, only certain duplications survive the filter of natural selection and can contribute to genome evolution. However, the reasons why certain genes can be retained after duplication whereas others cannot remain largely unknown. Many proteins contain intrinsically disordered regions (IDRs), whose structures fluctuate between alternative conformational states. Due to their high flexibility, IDRs often enable protein–protein interactions and are the target of post-translational modifications. Intrinsically disordered proteins (IDPs) have characteristics that might either stimulate or hamper the retention of their encoding genes after duplication. On the one hand, IDRs may enable functional diversification, thus promoting duplicate retention. On the other hand, increased IDP availability is expected to result in deleterious unspecific interactions. Here, we interrogate the proteomes of human, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Arabidopsis thaliana and Escherichia coli, in order to ascertain the impact of protein intrinsic disorder on gene duplicability. We show that, in general, proteins encoded by duplicated genes tend to be less disordered than those encoded by singletons. The only exception is proteins encoded by ohnologs, which tend to be more disordered than those encoded by singletons or genes resulting from small-scale duplications. Our results indicate that duplication of genes encoding IDPs outside the context of whole-genome duplication (WGD) is often deleterious, but that IDRs facilitate retention of duplicates in the context of WGD. We discuss the potential evolutionary implications of our results.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

An intrinsically disordered linker controlling the formation and the stability of the bacterial flagellar hook

Background

In a macro-molecular complex, any minor change may prove detrimental. For a supra-molecular nano-machine like the bacterial flagellum, which consists of several distinct parts with specific characteristics, stability is important. During the rotation of the bacterial flagellar motor, which is located in the membrane, the flagella rotate at speeds between 200 and 2000 rpm, depending on the bacterial species. The hook substructure of the bacterial flagellum acts as a universal joint connecting the motor to the flagellar filament. We investigated the formation of the bacterial flagellar hook and its overall stability between the FlgE subunits that make up the hook and attempted to understand how this stability differs between bacteria.

Results

An intrinsically disordered segment plays an important role for overall hook stability and for its structural cohesion during motor rotation. The length of this linker segment depends on the species of bacteria; for Salmonella enterica and Campylobacter jejuni it is approximately 37 and 54 residues, respectively. Few residues of the linker are conserved and mutating the conserved residues of the linker yields non-flagellated cells. In the case of Campylobacter, which rotates its flagella at a speed much higher than that of Salmonella, shortening the linker leads to a rupture of the hook at its base, decreasing cell motility. Our experiments show that this segment is required for polymerization and stability of the hook, demonstrating a surprising role for a disordered region in one of the most finely tuned and closely studied macromolecular machines.

Conclusions

This study reveals a detailed functional characteristic of an intrinsically disordered segment in the hook protein. This segment evolved to fulfill a specific role in the formation of the hook, and it is at the core of the stability and flexibility of the hook. Its length is important in the case of bacteria with high-speed rotating flagella. Finding a way of disrupting this linker in Campylobacter might help in preventing infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0438-7) contains supplementary material, which is available to authorized users.

Enhanced Sampling of Intrinsic Structural Heterogeneity of the BH3-Only Protein Binding Interface of Bcl-xL

Antiapoptotic Bcl-xL plays central roles in regulating programed cell death. Partial unfolding of Bcl-xL has been observed at the interface upon specific binding to the pro-apoptotic BH3-only protein PUMA, which in turn disrupts the interaction of Bcl-xL with tumor suppressor p53 and promotes apoptosis. Previous analysis of existing Bcl-xL structures and atomistic molecular dynamics (MD) simulations have suggested that substantial intrinsic structure heterogeneity exists at the BH3-only protein binding interface of Bcl-xL to facilitate its conformational transitions upon binding. In this study, enhanced sampling is applied to further characterize the interfacial conformations of unbound Bcl-xL in explicit solvent. Extensive replica exchange with solute tempering (REST) simulations, with a total accumulated time of 16 μs, were able to cover much wider conformational spaces for the interfacial region of Bcl-xL. The resulting structural ensembles are much better converged, with local and long-range structural features that are highly consistent with existing NMR data. These simulations further demonstrate that the BH3-only protein binding interface of Bcl-xL is intrinsically disordered and samples many rapidly interconverting conformations. Intriguingly, all previously observed conformers are well represented in the unbound structure ensemble. Such intrinsic structural heterogeneity and flexibility may be critical for Bcl-xL to undergo partial unfolding induced by PUMA binding, and likely provide a robust basis that allows Bcl-xL to respond sensitively to binding of various ligands in cellular signaling and regulation.

Small Molecule Enhancement of 20S Proteasome Activity Targets Intrinsically Disordered Proteins

The 20S proteasome is the main protease for the degradation of oxidatively damaged and intrinsically disordered proteins. When accumulation of disordered or oxidatively damaged proteins exceeds proper clearance in neurons, imbalanced pathway signaling or aggregation occurs, which have been implicated in the pathogenesis of several neurological disorders. Screening of the NIH Clinical Collection and Prestwick libraries identified the neuroleptic agent chlorpromazine as a lead agent capable of enhancing 20S proteasome activity. Chemical manipulation of chlorpromazine abrogated its D2R receptor binding affinity while retaining its ability to enhance 20S mediated proteolysis at low micromolar concentrations. The resulting small molecule enhancers of 20S proteasome activity induced the degradation of intrinsically disordered proteins, α-synuclein, and tau but not structured proteins. These small molecule 20S agonists can serve as leads to explore the therapeutic potential of 20S activation or as new tools to provide insight into the yet unclear mechanics of 20S-gate regulation.

Methods of probing the interactions between small molecules and disordered proteins

It is generally recognized that a large fraction of the human proteome is made up of proteins that remain disordered in their native states. Despite the fact that such proteins play key biological roles and are involved in many major human diseases, they still represent challenging targets for drug discovery. A major bottleneck for the identification of compounds capable of interacting with these proteins and modulating their disease-promoting behaviour is the development of effective techniques to probe such interactions. The difficulties in carrying out binding measurements have resulted in a poor understanding of the mechanisms underlying these interactions. In order to facilitate further methodological advances, here we review the most commonly used techniques to probe three types of interactions involving small molecules: (1) those that disrupt functional interactions between disordered proteins; (2) those that inhibit the aberrant aggregation of disordered proteins, and (3) those that lead to binding disordered proteins in their monomeric states. In discussing these techniques, we also point out directions for future developments.

To be disordered or not to be disordered: is that still a question for proteins in the cell?

There is ample evidence that many proteins or regions of proteins lack a well-defined folded structure under native-like conditions. These are called intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). Whether this intrinsic disorder is also their main structural characteristic in living cells has been a matter of intense debate. The structural analysis of IDPs became an important challenge also because of their involvement in a plethora of human diseases, which made IDPs attractive targets for therapeutic development. Therefore, biophysical approaches are increasingly being employed to probe the structural and dynamical state of proteins, not only in isolation in a test tube, but also in a complex biological environment and even within intact cells. Here, we survey direct and indirect evidence that structural disorder is in fact the physiological state of many proteins in the proteome. The paradigmatic case of α-synuclein is used to illustrate the controversial nature of this topic.

Predictions of Backbone Dynamics in Intrinsically Disordered Proteins Using De Novo Fragment-Based Protein Structure Predictions

Intrinsically disordaered proteins (IDPs) are a prevalent phenomenon with over 30% of human proteins estimated to have long disordered regions. Computational methods are widely used to study IDPs, however, nearly all treat disorder in a binary fashion, not accounting for the structural heterogeneity present in disordered regions. Here, we present a new de novo method, FRAGFOLD-IDP, which addresses this problem. Using 200 protein structural ensembles derived from NMR, we show that FRAGFOLD-IDP achieves superior results compared to methods which can predict related data (NMR order parameter, or crystallographic B-factor). FRAGFOLD-IDP produces very good predictions for 33.5% of cases and helps to get a better insight into the dynamics of the disordered ensembles. The results also show it is not necessary to predict the correct fold of the protein to reliably predict per-residue fluctuations. It implies that disorder is a local property and it does not depend on the fold. Our results are orthogonal to DynaMine, the only other method significantly better than the naïve prediction. We therefore combine these two using a neural network. FRAGFOLD-IDP enables better insight into backbone dynamics in IDPs and opens exciting possibilities for the design of disordered ensembles, disorder-to-order transitions, or design for protein dynamics.

Influence of Flexible "ω" on the Activity of E. coli RNA Polymerase: A Thermodynamic Analysis

The Escherichia coli RNA polymerase (RNAP) is a multisubunit protein complex containing the smallest subunit, ω. Despite the evolutionary conservation of ω and its role in assembly of RNAP, E. coli mutants lacking rpoZ (codes for ω) are viable due to the association of RNAP with the global chaperone protein GroEL. With an aim to get better insight into the structure and functional role of ω, we isolated a dominant negative mutant of ω (ω₆), which is predominantly α-helical, in contrast to largely unstructured native ω, and then studied its assembly with reconstituted core1 (α₂ββ') by a biophysical approach. The mutant showed higher binding affinity compared to native ω. We observed that the interaction between core1 and ω₆ is driven by highly negative enthalpy and a small but unfavorable negative entropy term. Extensive structural alteration in ω₆ makes it more rigid, the plasticity of the interacting domain formed by ω₆ and core1 is compromised, which may be responsible for the entropic cost. Such tight binding of the structured mutant (ω₆) affects initiation of transcription. However, once preinitiated, the complex elongates the RNA chain efficiently. The initiation of transcription requires recognition of appropriate σ-factors by the core enzyme (core2: α₂ββ'ω). We found that the altered core enzyme (α₂ββ'ω₆) with mutant ω showed a decrease in binding affinity to the σ-factors (σ⁷⁰, σ³² and σ³⁸) compared to that of the core enzyme containing native ω. In the absence of unstructured ω, the association of σ-factors to the core is less efficient, suggesting that the flexible native ω plays a direct role in σ-factor recruitment.

The contribution of intrinsically disordered regions to protein function, cellular complexity, and human disease

In the 1960s, Christian Anfinsen postulated that the unique three-dimensional structure of a protein is determined by its amino acid sequence. This work laid the foundation for the sequence–structure–function paradigm, which states that the sequence of a protein determines its structure, and structure determines function. However, a class of polypeptide segments called intrinsically disordered regions does not conform to this postulate. In this review, I will first describe established and emerging ideas about how disordered regions contribute to protein function. I will then discuss molecular principles by which regulatory mechanisms, such as alternative splicing and asymmetric localization of transcripts that encode disordered regions, can increase the functional versatility of proteins. Finally, I will discuss how disordered regions contribute to human disease and the emergence of cellular complexity during organismal evolution.

A pH-dependent switch promotes β-synuclein fibril formation via glutamate residues

α-Synuclein (αS) is the primary protein associated with Parkinson's disease, and it undergoes aggregation from its intrinsically disordered monomeric form to a cross-β fibrillar form. The closely related homolog β-synuclein (βS) is essentially fibril-resistant under cytoplasmic physiological conditions. Toxic gain-of-function by βS has been linked to dysfunction, but the aggregation behavior of βS under altered pH is not well-understood. In this work, we compare fibril formation of αS and βS at pH 7.3 and mildly acidic pH 5.8, and we demonstrate that pH serves as an on/off switch for βS fibrillation. Using αS/βS domain-swapped chimera constructs and single residue substitutions in βS, we localized the switch to acidic residues in the N-terminal and non-amyloid component domains of βS. Computational models of βS fibril structures indicate that key glutamate residues (Glu-31 and Glu-61) in these domains may be sites of pH-sensitive interactions, and variants E31A and E61A show dramatically altered pH sensitivity for fibril formation supporting the importance of these charged side chains in fibril formation of βS. Our results demonstrate that relatively small changes in pH, which occur frequently in the cytoplasm and in secretory pathways, may induce the formation of βS fibrils and suggest a complex role for βS in synuclein cellular homeostasis and Parkinson's disease.

Eukaryotic transcription factors: paradigms of protein intrinsic disorder

Gene-specific transcription factors (TFs) are key regulatory components of signaling pathways, controlling, for example, cell growth, development, and stress responses. Their biological functions are determined by their molecular structures, as exemplified by their structured DNA-binding domains targeting specific cis-acting elements in genes, and by the significant lack of fixed tertiary structure in their extensive intrinsically disordered regions. Recent research in protein intrinsic disorder (ID) has changed our understanding of transcriptional activation domains from ‘negative noodles’ to ID regions with function-related, short sequence motifs and molecular recognition features with structural propensities. This review focuses on molecular aspects of TFs, which represent paradigms of ID-related features. Through specific examples, we review how the ID-associated flexibility of TFs enables them to participate in large interactomes, how they use only a few hydrophobic residues, short sequence motifs, prestructured motifs, and coupled folding and binding for their interactions with co-activators, and how their accessibility to post-translational modification affects their interactions. It is furthermore emphasized how classic biochemical concepts like allostery, conformational selection, induced fit, and feedback regulation are undergoing a revival with the appreciation of ID. The review also describes the most recent advances based on computational simulations of ID-based interaction mechanisms and structural analysis of ID in the context of full-length TFs and suggests future directions for research in TF ID.

CtIP/Ctp1/Sae2, molecular form fit for function

Vertebrate CtIP, and its fission yeast (Ctp1), budding yeast (Sae2) and plant (Com1) orthologs have emerged as key regulatory molecules in cellular responses to DNA double strand breaks (DSBs). By modulating the nucleolytic 5'-3' resection activity of the Mre11/Rad50/Nbs1 (MRN) DSB repair processing and signaling complex, CtIP/Ctp1/Sae2/Com1 is integral to the channeling of DNA double strand breaks through DSB repair by homologous recombination (HR). Nearly two decades since its discovery, emerging new data are defining the molecular underpinnings for CtIP DSB repair regulatory activities. CtIP homologs are largely intrinsically unstructured proteins comprised of expanded regions of low complexity sequence, rather than defined folded domains typical of DNA damage metabolizing enzymes and nucleases. A compact structurally conserved N-terminus forms a functionally critical tetrameric helical dimer of dimers (THDD) region that bridges CtIP oligomers, and is flexibly appended to a conserved C-terminal Sae2-homology DNA binding and DSB repair pathway choice regulatory hub which influences nucleolytic activities of the MRN core nuclease complex. The emerging evidence from structural, biophysical, and biological studies converges on CtIP having functional roles in DSB repair that include: 1) dynamic DNA strand coordination through direct DNA binding and DNA bridging activities, 2) MRN nuclease complex cofactor functions that direct MRN endonucleolytic cleavage of protein-blocked DSB ends and 3) acting as a protein binding hub targeted by the cell cycle regulatory apparatus, which influences CtIP expression and activity via layers of post-translational modifications, protein-protein interactions and DNA binding.

Effect of different γ-subunit isoforms on the regulation of AMPK

AMP-activated protein kinase (AMPK) plays a key role in integrating metabolic pathways in response to energy demand. AMPK activation results in a wide range of downstream responses, many of which are associated with improved metabolic outcome, making AMPK an attractive target for the treatment of metabolic diseases. AMPK is a heterotrimeric complex consisting of a catalytic subunit (α) and two regulatory subunits (β and γ). The γ-subunit harbours the nucleotide-binding sites and plays an important role in AMPK regulation in response to cellular energy levels. In mammals, there are three isoforms of the γ-subunit and these respond differently to regulation by nucleotides, but there is limited information regarding their role in activation by small molecules. Here, we determined the effect of different γ-isoforms on AMPK by a direct activator, 991. In cells, 991 led to a greater activation of γ2-containing AMPK complexes compared with either γ1 or γ3. This effect was dependent on the long N-terminal region of the γ2-isoform. We were able to rule out an effect of Ser108 phosphorylation, since mutation of Ser108 to alanine in the β2-isoform had no effect on activation of AMPK by 991 in either γ1- or γ2-complexes. The rate of dephosphorylation of Thr172 was slower for γ2- compared with γ1-complexes, both in the absence and presence of 991. Our studies show that activation of AMPK by 991 depends on the nature of the γ-isoform. This finding may have implications for the design of isoform-selective AMPK activators.

Mycobacterium tuberculosis Complex Exhibits Lineage-Specific Variations Affecting Protein Ductility and Epitope Recognition

The advent of whole-genome sequencing has provided an unprecedented detail about the evolution and genetic significance of species-specific variations across the whole Mycobacterium tuberculosis Complex. However, little attention has been focused on understanding the functional roles of these variations in the protein coding sequences. In this work, we compare the coding sequences from 74 sequenced mycobacterial species including M. africanum, M. bovis, M. canettii, M. caprae, M. orygis, and M. tuberculosis. Results show that albeit protein variations affect all functional classes, those proteins involved in lipid and intermediary metabolism and respiration have accumulated mutations during evolution. To understand the impact of these mutations on protein functionality, we explored their implications on protein ductility/disorder, a yet unexplored feature of mycobacterial proteomes. In agreement with previous studies, we found that a Gly71Ile substitution in the PhoPR virulence system severely affects the ductility of its nearby region in M. africanum and animal-adapted species. In the same line of evidence, the SmtB transcriptional regulator shows amino acid variations specific to the Beijing lineage, which affects the flexibility of the N-terminal trans-activation domain. Furthermore, despite the fact that MTBC epitopes are evolutionary hyperconserved, we identify strain- and lineage-specific amino acid mutations affecting previously known T-cell epitopes such as EsxH and FbpA (Ag85A). Interestingly, in silico studies reveal that these variations result in differential interaction of epitopes with the main HLA haplogroups.

Structural Characterization of the Intrinsically Disordered Protein p53 Using Raman Spectroscopy

The intrinsically disordered protein p53 has attracted a strong interest for its important role in genome safeguarding and potential therapeutic applications. However, its disordered character makes difficult a full characterization of p53 structural architecture. A deep knowledge of p53 structural motifs could significantly help the understanding of its functional properties, in connection with its complex binding network. We have applied Raman spectroscopy to investigate the structural composition and the conformational heterogeneity of both full-length p53 and its DNA binding domain (DBD), in different solvent environments. In particular, a careful analysis of the Amide I Raman band, which is highly sensitive to protein secondary structure elements such as α-helices, β-sheets and random coils, has revealed the presence of extended random coils in p53 and predominant β-sheet regions in its DBD. In addition, this analysis has allowed us to explore the ensemble of interchanging conformations in both p53 and its DBD, highlighting a higher conformational heterogeneity in p53 than in its DBD. Furthermore, by applying a principal components analysis, we have identified the principal spectral markers in both p53 and DBD samples. The combination of the two approaches could be insightful for the study of intrinsically disordered proteins, by offering increased versatility and wide application as a label-free, real-time and non-invasive detection method.

Human EGF-derived direct and reverse short linear motifs: conformational dynamics insight into the receptor-binding residues

Short linear motifs (SLiMs) have been recognized to perform diverse functions in a variety of regulatory proteins through the involvement in protein-protein interactions, signal transduction, cell cycle regulation, protein secretion, etc. However, detailed molecular mechanisms underlying their functions including roles of definite amino acid residues remain obscure. In our previous studies, we demonstrated that conformational dynamics of amino acid residues in oligopeptides derived from regulatory proteins such as alpha-fetoprotein (AFP), carcino-embryonic antigen (CEA), and pregnancy specific β1-glycoproteins (PSGs) contributes greatly to their biological activities. In the present work, we revealed the 22-member linear modules composed of direct and reverse AFP_14-20-like heptapeptide motifs linked by CxxGY/FxGx consensus motif within epidermal growth factor (EGF), growth factors of EGF family and numerous regulatory proteins containing EGF-like modules. We showed, first, the existence of similarity in amino acid signatures of both direct and reverse motifs in terms of their physicochemical properties. Second, molecular dynamics (MD) simulation study demonstrated that key receptor-binding residues in human EGF in the aligned positions of the direct and reverse motifs may have similar distribution of conformational probability densities and dynamic behavior despite their distinct physicochemical properties. Third, we found that the length of a polypeptide chain (from 7 to 53 residues) has no effect, while disulfide bridging and backbone direction significantly influence the conformational distribution and dynamics of the residues. Our data may contribute to the atomic level structure-function analysis and protein structure decoding; additionally, they may provide a basis for novel protein/peptide engineering and peptide-mimetic drug design.

The antiproliferative ELF2 isoform, ELF2B, induces apoptosis in vitro and perturbs early lymphocytic development in vivo

Background

ELF2 (E74-like factor 2) also known as NERF (new Ets-related factor), a member of the Ets family of transcription factors, regulates genes important in B and T cell development, cell cycle progression, and angiogenesis. Conserved ELF2 isoforms, ELF2A, and ELF2B, arising from alternative promoter usage can exert opposing effects on target gene expression. ELF2A activates, whilst ELF2B represses, gene expression, and the balance of expression between these isoforms may be important in maintaining normal cellular function.

Methods

We compared the function of ELF2 isoforms ELF2A and ELF2B with other ELF subfamily proteins ELF1 and ELF4 in primary and cancer cell lines using proliferation, colony-forming, cell cycle, and apoptosis assays. We further examined the role of ELF2 isoforms in haemopoietic development using a Rag1^-/-murine bone marrow reconstitution model.

Results

ELF2B overexpression significantly reduced cell proliferation and clonogenic capacity, minimally disrupted cell cycle kinetics, and induced apoptosis. In contrast, ELF2A overexpression only marginally reduced clonogenic capacity with little effect on proliferation, cell cycle progression, or apoptosis. Deletion of the N-terminal 19 amino acids unique to ELF2B abrogated the antiproliferative and proapoptotic functions of ELF2B thereby confirming its crucial role. Mice expressing Elf2a or Elf2b in haemopoietic cells variously displayed perturbations in the pre-B cell stage and multiple stages of T cell development. Mature B cells, T cells, and myeloid cells in steady state were unaffected, suggesting that the main role of ELF2 is restricted to the early development of B and T cells and that compensatory mechanisms exist. No differences in B and T cell development were observed between ELF2 isoforms.

Conclusions

We conclude that ELF2 isoforms are important regulators of cellular proliferation, cell cycle progression, and apoptosis. In respect to this, ELF2B acts in a dominant negative fashion compared to ELF2A and as a putative tumour suppressor gene. Given that these cellular processes are critical during haemopoiesis, we propose that the regulatory interplay between ELF2 isoforms contributes substantially to early B and T cell development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13045-017-0446-7) contains supplementary material, which is available to authorized users.

Systemic and network functions of the microtubule-associated protein tau: Implications for tau-based therapies

Tau is a microtubule-associated neuronal protein, whose primary role was long thought to regulate axonal microtubule assembly. Tau is subject to many posttranslational modifications and can aggregate into neurofibrillary tangles, which are considered to be a hallmark of several neurodegenerative diseases collectively called "tauopathies". The most common tauopathy is Alzheimer's disease, where tau pathology correlates with sites of neurodegeneration. Tau belongs to the class of intrinsically disordered proteins, which are known to interact with many partners and are considered to be involved in various signaling, regulation and recognition processes. Thus more recent evidence indicates that tau functionally interacts with many proteins and different cellular structures, which may have an important physiological role and may be involved in neurodegenerative processes. Furthermore, tau can be released from neurons and exert functional effects on other cells. This review article weighs the evidence that tau has subtle but important systemic effects on neuronal network function by maintaining physiological neuronal transmission and synaptic plasticity, which are possibly independent from tau's microtubule modulating activities. Implications for tau-based therapeutic approaches are discussed.

GARP2 accelerates retinal degeneration in rod cGMP-gated cation channel β-subunit knockout mice

The Cngb1 locus-encoded β-subunit of rod cGMP-gated cation channel and associated glutamic acid rich proteins (GARPs) are required for phototransduction, disk morphogenesis, and rod structural integrity. To probe individual protein structure/function of the GARPs, we have characterized several transgenic mouse lines selectively restoring GARPs on a Cngb1 knockout (X1^−/−) mouse background. Optical coherence tomography (OCT), light and transmission electron microscopy (TEM), and electroretinography (ERG) were used to analyze 6 genotypes including WT at three and ten weeks postnatal. Comparison of aligned histology/OCT images demonstrated that GARP2 accelerates the rate of degeneration. ERG results are consistent with the structural analyses showing the greatest attenuation of function when GARP2 is present. Even 100-fold or more overexpression of GARP1 could not accelerate degeneration as rapidly as GARP2, and when co-expressed GARP1 attenuated the structural and functional deficits elicited by GARP2. These results indicate that the GARPs are not fully interchangeable and thus, likely have separate and distinct functions in the photoreceptor. We also present a uniform murine OCT layer naming nomenclature system that is consistent with human retina layer designations to standardize murine OCT, which will facilitate data evaluation across different laboratories.

The thermodynamics of protein aggregation reactions may underpin the enhanced metabolic efficiency associated with heterosis, some balancing selection, and the evolution of ploidy levels

Identifying the physical basis of heterosis (or "hybrid vigor") has remained elusive despite over a hundred years of research on the subject. The three main theories of heterosis are dominance theory, overdominance theory, and epistasis theory. Kacser and Burns (1981) identified the molecular basis of dominance, which has greatly enhanced our understanding of its importance to heterosis. This paper aims to explain how overdominance, and some features of epistasis, can similarly emerge from the molecular dynamics of proteins. Possessing multiple alleles at a gene locus results in the synthesis of different allozymes at reduced concentrations. This in turn reduces the rate at which each allozyme forms soluble oligomers, which are toxic and must be degraded, because allozymes co-aggregate at low efficiencies. The model developed in this paper can explain how heterozygosity impacts the metabolic efficiency of an organism. It can also explain why the viabilities of some inbred lines seem to decline rapidly at high inbreeding coefficients (F > 0.5), which may provide a physical basis for truncation selection for heterozygosity. Finally, the model has implications for the ploidy level of organisms. It can explain why polyploids are frequently found in environments where severe physical stresses promote the formation of soluble oligomers. The model can also explain why complex organisms, which need to synthesize aggregation-prone proteins that contain intrinsically unstructured regions (IURs) and multiple domains because they facilitate complex protein interaction networks (PINs), tend to be diploid while haploidy tends to be restricted to relatively simple organisms.

Computer Simulations of Intrinsically Disordered Proteins

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

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

Binding of Amphipathic Cell Penetrating Peptide p28 to Wild Type and Mutated p53 as studied by Raman, Atomic Force and Surface Plasmon Resonance spectroscopies

BACKGROUND

Mutations within the DNA binding domain (DBD) of the tumor suppressor p53 are found in >50% of human cancers and may significantly modify p53 secondary structure impairing its function. p28, an amphipathic cell-penetrating peptide, binds to the DBD through hydrophobic interaction and induces a posttranslational increase in wildtype and mutant p53 restoring functionality. We use mutation analyses to explore which elements of secondary structure may be critical to p28 binding.

METHODS

Molecular modeling, Raman spectroscopy, Atomic Force Spectroscopy (AFS) and Surface Plasmon Resonance (SPR) were used to identify which secondary structure of site-directed and naturally occurring mutant DBDs are potentially altered by discrete changes in hydrophobicity and the molecular interaction with p28.

RESULTS

We show that specific point mutations that alter hydrophobicity within non-mutable and mutable regions of the p53 DBD alter specific secondary structures. The affinity of p28 was positively correlated with the β-sheet content of a mutant DBD, and reduced by an increase in unstructured or random coil that resulted from a loss in hydrophobicity and redistribution of surface charge.

CONCLUSIONS

These results help refine our knowledge of how mutations within p53-DBD alter secondary structure and provide insight on how potential structural alterations in p28 or similar molecules improve their ability to restore p53 function.

GENERAL SIGNIFICANCE

Raman spectroscopy, AFS, SPR and computational modeling are useful approaches to characterize how mutations within the p53DBD potentially affect secondary structure and identify those structural elements prone to influence the binding affinity of agents designed to increase the functionality of p53.

Detection of correlated conformational fluctuations in intrinsically disordered proteins through paramagnetic relaxation interference

Functionally relevant conformational states of intrinsically disordered proteins (IDPs) are typically concealed in a vast space of fast interconverting structures. Here we present a novel methodology, NMR-based paramagnetic relaxation interference (PRI), that allows for direct observation of concerted motions and cooperatively folded sub-states in IDPs. The proposed NMR technique is based on the exploitation of cross correlated electron-nuclear dipolar relaxation interferences in doubly spin-labeled proteins and probes the transient spatial encounter of electron-nucleus spin pairs.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (R_h ) when measured under normal conditions. To investigate sequence-modulation of IDP R_h , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental R_h from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PP_II ), however, modulated R_h with intricate dependence on the simulated propensities. PP_II propensities at levels expected in IDPs correlated with heightened R_h sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP R_h can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.

A sequence-based computational method for prediction of MoRFs

Molecular recognition features (MoRFs) are relatively short segments (10–70 residues) within intrinsically disordered regions (IDRs) that can undergo disorder-to-order transitions during binding to partner proteins.

The transient manifold structure of the p53 extreme C-terminal domain: insight into disorder, recognition, and binding promiscuity by molecular dynamics simulations

The extreme C-terminus of the p53 tumour suppressor (p53-CTD) is a 30 residue long intrinsically disordered region, responsible for regulating the p53 DNA binding activity. Extensive conformational sampling through MD simulations of a p53-CTD derived peptide in solution highlights its propensity to form short and stable secondary structure motifs, specifically localized within the sequence.