Kinetic recognition of the retinoblastoma tumor suppressor by a specific protein target

The retinoblastoma tumor suppressor (Rb) plays a key role in cell cycle control and is linked to various types of human cancer. Rb binds to the LxCxE motif, present in a number of cellular and viral proteins such as AdE1A, SV40 large T-antigen and human papillomavirus (HPV) E7, all instrumental in revealing fundamental mechanisms of tumor suppression, cell cycle control and gene expression. A detailed kinetic study of RbAB binding to the HPV E7 oncoprotein shows that an LxCxE-containing E7 fragment binds through a fast two-state reaction strongly favored by electrostatic interactions. Conversely, full-length E7 binds through a multistep process involving a pre-equilibrium between E7 conformers, a fast electrostatically driven association step guided by the LxCxE motif and a slow conformational rearrangement. This kinetic complexity arises from the conformational plasticity and intrinsically disordered nature of E7 and from multiple interaction surfaces present in both proteins. Affinity differences between E7N domains from high- and low-risk types are explained by their dissociation rates. In fact, since Rb is at the center of a large protein interaction network, fast and tight recognition provides an advantage for disruption by the viral proteins, where the balance of physiological and pathological interactions is dictated by kinetic ligand competition. The localization of the LxCxE motif within an intrinsically disordered domain provides the fast, diffusion-controlled interaction that allows viral proteins to outcompete physiological targets. We describe the interaction mechanism of Rb with a protein ligand, at the same time an LxCxE-containing model target, and a paradigmatic intrinsically disordered viral oncoprotein.     

Expanding the Proteome of an RNA Virus by Phosphorylation of an Intrinsically Disordered Viral Protein

The human proteome contains myriad intrinsically disordered proteins. Within intrinsically disordered proteins, polyproline-II motifs are often located near sites of phosphorylation. We have used an unconventional experimental paradigm to discover that phosphorylation by protein kinase A (PKA) occurs in the intrinsically disordered domain of hepatitis C virus non-structural protein 5A (NS5A) on Thr-2332 near one of its polyproline-II motifs. Phosphorylation shifts the conformational ensemble of the NS5A intrinsically disordered domain to a state that permits detection of the polyproline motif by using ¹⁵N-, ¹³C-based multidimensional NMR spectroscopy. PKA-dependent proline resonances were lost in the presence of the Src homology 3 domain of c-Src, consistent with formation of a complex. Changing Thr-2332 to alanine in hepatitis C virus genotype 1b reduced the steady-state level of RNA by 10-fold; this change was lethal for genotype 2a. The lethal phenotype could be rescued by changing Thr-2332 to glutamic acid, a phosphomimetic substitution. Immunofluorescence and transmission electron microscopy showed that the inability to produce Thr(P)-2332-NS5A caused loss of integrity of the virus-induced membranous web/replication organelle. An even more extreme phenotype was observed in the presence of small molecule inhibitors of PKA. We conclude that the PKA-phosphorylated form of NS5A exhibits unique structure and function relative to the unphosphorylated protein. We suggest that post-translational modification of viral proteins containing intrinsic disorder may be a general mechanism to expand the viral proteome without a corresponding expansion of the genome.     

SLiMPrints: conservation-based discovery of functional motif fingerprints in intrinsically disordered protein regions

Large portions of higher eukaryotic proteomes are intrinsically disordered, and abundant evidence suggests that these unstructured regions of proteins are rich in regulatory interaction interfaces. A major class of disordered interaction interfaces are the compact and degenerate modules known as short linear motifs (SLiMs). As a result of the difficulties associated with the experimental identification and validation of SLiMs, our understanding of these modules is limited, advocating the use of computational methods to focus experimental discovery. This article evaluates the use of evolutionary conservation as a discriminatory technique for motif discovery. A statistical framework is introduced to assess the significance of relatively conserved residues, quantifying the likelihood a residue will have a particular level of conservation given the conservation of the surrounding residues. The framework is expanded to assess the significance of groupings of conserved residues, a metric that forms the basis of SLiMPrints (short linear motif fingerprints), a de novo motif discovery tool. SLiMPrints identifies relatively overconstrained proximal groupings of residues within intrinsically disordered regions, indicative of putatively functional motifs. Finally, the human proteome is analysed to create a set of highly conserved putative motif instances, including a novel site on translation initiation factor eIF2A that may regulate translation through binding of eIF4E.     

The High-Risk HPV16 E7 Oncoprotein Mediates Interaction between the Transcriptional Coactivator CBP and the Retinoblastoma Protein pRb

The oncoprotein E7 from human papillomavirus (HPV) strains that confer high cancer risk mediates cell transformation by deregulating host cellular processes and activating viral gene expression through recruitment of cellular proteins such as the retinoblastoma protein (pRb) and the cyclic-AMP response element binding binding protein (CBP) and its paralog p300. Here we show that the intrinsically disordered N-terminal region of E7 from high-risk HPV16 binds the TAZ2 domain of CBP with greater affinity than E7 from low-risk HPV6b. HPV E7 and the tumor suppressor p53 compete for binding to TAZ2. The TAZ2 binding site in E7 overlaps the LxCxE motif that is crucial for interaction with pRb. While TAZ2 and pRb compete for binding to a monomeric E7 polypeptide, the full-length E7 dimer mediates an interaction between TAZ2 and pRb by promoting formation of a ternary complex. Cell-based assays show that expression of full-length HPV16 E7 promotes increased pRb acetylation and that this response depends both on the presence of CBP/p300 and on the ability of E7 to form a dimer. These observations suggest a model for the oncogenic effect of high-risk HPV16 E7. The disordered region of one E7 molecule in the homodimer interacts with the pocket domain of pRb, while the same region of the other E7 molecule binds the TAZ2 domain of CBP/p300. Through its ability to dimerize, E7 recruits CBP/p300 and pRb into a ternary complex, bringing the histone acetyltransferase domain of CBP/p300 into proximity to pRb and promoting acetylation, leading to disruption of cell cycle control.     

GlobPlot: Exploring protein sequences for globularity and disorder

A major challenge in the proteomics and structural genomics era is to predict protein structure and function, including identification of those proteins that are partially or wholly unstructured. Non-globular sequence segments often contain short linear peptide motifs (e.g. SH3-binding sites) which are important for protein function. We present here a new tool for discovery of such unstructured, or disordered regions within proteins. GlobPlot (http://globplot.embl.de) is a web service that allows the user to plot the tendency within the query protein for order/globularity and disorder. We show examples with known proteins where it successfully identifies inter-domain segments containing linear motifs, and also apparently ordered regions that do not contain any recognised domain. GlobPlot may be useful in domain hunting efforts. The plots indicate that instances of known domains may often contain additional N- or C-terminal segments that appear ordered. Thus GlobPlot may be of use in the design of constructs corresponding to globular proteins, as needed for many biochemical studies, particularly structural biology. GlobPlot has a pipeline interface--GlobPipe--for the advanced user to do whole proteome analysis. GlobPlot can also be used as a generic infrastructure package for graphical displaying of any possible propensity.     

A structure filter for the Eukaryotic Linear Motif Resource

Background  

Many proteins are highly modular, being assembled from globular domains and segments of natively disordered polypeptides. Linear motifs, short sequence modules functioning independently of protein tertiary structure, are most abundant in natively disordered polypeptides but are also found in accessible parts of globular domains, such as exposed loops. The prediction of novel occurrences of known linear motifs attempts the difficult task of distinguishing functional matches from stochastically occurring non-functional matches. Although functionality can only be confirmed experimentally, confidence in a putative motif is increased if a motif exhibits attributes associated with functional instances such as occurrence in the correct taxonomic range, cellular compartment, conservation in homologues and accessibility to interacting partners. Several tools now use these attributes to classify putative motifs based on confidence of functionality.     

Prediction of short linear protein binding regions

Short linear motifs in proteins (typically 3-12 residues in length) play key roles in protein-protein interactions by frequently binding specifically to peptide binding domains within interacting proteins. Their tendency to be found in disordered segments of proteins has meant that they have often been overlooked. Here we present SLiMPred (short linear motif predictor), the first general de novo method designed to computationally predict such regions in protein primary sequences independent of experimentally defined homologs and interactors. The method applies machine learning techniques to predict new motifs based on annotated instances from the Eukaryotic Linear Motif database, as well as structural, biophysical, and biochemical features derived from the protein primary sequence. We have integrated these data sources and benchmarked the predictive accuracy of the method, and found that it performs equivalently to a predictor of protein binding regions in disordered regions, in addition to having predictive power for other classes of motif sites such as polyproline II helix motifs and short linear motifs lying in ordered regions. It will be useful in predicting peptides involved in potential protein associations and will aid in the functional characterization of proteins, especially of proteins lacking experimental information on structures and interactions. We conclude that, despite the diversity of motif sequences and structures, SLiMPred is a valuable tool for prioritizing potential interaction motifs in proteins.     

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.     

Systematic discovery of new recognition peptides mediating protein interaction networks

Many aspects of cell signalling, trafficking, and targeting are governed by interactions between globular protein domains and short peptide segments. These domains often bind multiple peptides that share a common sequence pattern, or “linear motif” (e.g., SH3 binding to PxxP). Many domains are known, though comparatively few linear motifs have been discovered. Their short length (three to eight residues), and the fact that they often reside in disordered regions in proteins makes them difficult to detect through sequence comparison or experiment. Nevertheless, each new motif provides critical molecular details of how interaction networks are constructed, and can explain how one protein is able to bind to very different partners. Here we show that binding motifs can be detected using data from genome-scale interaction studies, and thus avoid the normally slow discovery process. Our approach based on motif over-representation in non-homologous sequences, rediscovers known motifs and predicts dozens of others. Direct binding experiments reveal that two predicted motifs are indeed protein-binding modules: a DxxDxxxD protein phosphatase 1 binding motif with a K_D of 22 μM and a VxxxRxYS motif that binds Translin with a K_D of 43 μM. We estimate that there are dozens or even hundreds of linear motifs yet to be discovered that will give molecular insight into protein networks and greatly illuminate cellular processes.     

Strategies for bacterial expression of protein-peptide complexes: application to solubilization of papillomavirus E6

E6 is a small oncoprotein involved in tumorigenesis induced by papillomaviruses (PVs). E6 often recognizes its cellular targets by binding to short motifs presenting the consensus LXXLL. E6 proteins have long resisted structural analysis. We found that bovine papillomavirus type 1 (BPV1) E6 binds the N-terminal LXXLL motif of the cellular protein paxillin with significantly higher affinity as compared to other E6/peptide interactions. Although recombinant BPV1 E6 was poorly soluble in the free state, provision of the paxillin LXXLL peptide during BPV1 E6 biosynthesis greatly enhanced the protein's solubility. Expression of BPV1 E6/LXXLL peptide complexes was carried out in bacteria in the form of triple fusion constructs comprising, from N- to C-terminus, the soluble carrier protein maltose binding protein (MBP), the LXXLL motif and the E6 protein. A TEV protease cleavage site was placed either between MBP and LXXLL motif or between LXXLL motif and E6. These constructs allowed us to produce highly concentrated samples of BPV1 E6, either covalently fused to the C-terminus of the LXXLL motif (intra-molecular complex) or non-covalently bound to it (inter-molecular complex). Heteronuclear NMR measurements were performed and showed that the E6 protein was folded with similar conformations in both covalent and non-covalent complexes. These data open the way to novel structural and functional studies of the BPV1 E6 in complex with its preferential target motif.     

The HPV16 E7 viral oncoprotein self-assembles into defined spherical oligomers

Despite the fact that E7 is a major transforming oncoprotein in papillomavirus, its structure and precise molecular mechanism of action remain puzzling to date. E7 proteins share sequence homology and proteasome targeting properties of tumor suppressors with adenovirus E1A and SV40 T antigen, two other paradigmatic oncoproteins from DNA tumor viruses. High-risk HPV16 E7, a nonglobular dimer with some properties of intrinsically disordered proteins, is capable of undergoing pH-dependent conformational transitions that expose hydrophobic surfaces to the solvent. We found that treatment with a chelating agent produced a protein that can readily assemble into homogeneous spherical particles with an average molecular mass of 790 kDa and a diameter of 50 nm, as determined from dynamic light scattering and electron microscopy. The protein undergoes a substantial conformational transition from coil to beta-sheet structure, with concomitant consolidation of tertiary structure as judged by circular dichroism and fluorescence. The assembly process is very slow, in agreement with a substantial energy barrier caused by structural rearrangements. The resulting particles are highly stable, cooperatively folded, and capable of binding both Congo Red and thioflavin T, reporters of repetitive beta-sheet structures similar to those found in amyloids, although no fibrillar or insoluble material was observed under our experimental conditions.     

Profile-based short linear protein motif discovery

ABSTRACT: BACKGROUND: Short linear protein motifs are attracting increasing attention as functionally independent sites, typically 3-10 amino acids in length that are enriched in disordered regions of proteins. Multiple methods have recently been proposed to discover over-represented motifs within a set of proteins based on simple regular expressions. Here, we extend these approaches to profile-based methods, which provide a richer motif representation. RESULTS: The profile motif discovery method MEME performed relatively poorly for motifs in disordered regions of proteins. However, when we applied evolutionary weighting to account for redundancy amongst homologous proteins, and masked out poorly conserved regions of disordered proteins, the performance of MEME is equivalent to that of regular expression methods. However, the two approaches returned different subsets within both a benchmark dataset, and a more realistic discovery dataset. CONCLUSIONS: Profile-based motif discovery methods complement regular expression based methods. Whilst profile-based methods are computationally more intensive, they are likely to discover motifs currently overlooked by regular expression methods.     

Expanding the Paradigm: Intrinsically Disordered Proteins and Allosteric Regulation

Allosteric regulatory processes are implicated at all levels of biological function. Recent advances in our understanding of the diverse and functionally significant class of intrinsically disordered proteins have identified a multitude of ways in which disordered proteins function within the confines of the allosteric paradigm. Allostery within or mediated by intrinsically disordered proteins ensures robust and efficient signal integration through mechanisms that would be extremely unfavorable or even impossible for globular protein interaction partners. Here, we highlight recent examples that indicate the breadth of biological outcomes that can be achieved through allosteric regulation by intrinsically disordered proteins. Ongoing and future work in this rapidly evolving area of research will expand our appreciation of the central role of intrinsically disordered proteins in ensuring the fidelity and efficiency of cellular regulation.     

Limited proteolysis of natively unfolded protein 4E‐BP1 in the presence of trifluoroethanol

Natively unfolded (intrinsically disordered (ID) proteins) have been attracting an increasing attention due to their involvement in many regulatory processes. Natively unfolded proteins can fold upon binding to their metabolic partners. Coupled folding and binding events usually involve only relatively short motifs (binding motifs). These binding motifs which are able to fold should have an increased propensity to form a secondary structure. The aim of the present work was to probe the conformation of the intrinsically disordered protein 4E‐BP1 in the native and partly folded states by limited proteolysis and to reveal regions with a high propensity to form an ordered structure. Trifuoroethanol (TFE) in low concentrations (up to 15 vol%) was applied to increase the helical population of protein regions with a high intrinsic propensity to fold. When forming helical structures, these regions lose mobility and become more protected from proteases than random/unfolded protein regions. Limited proteolysis followed by mass spectrometry analysis allows identification of the regions with decreased mobility in TFE solutions. Trypsin and V8 proteases were used to perform limited proteolysis of the 4E‐BP1 protein in buffer and in solutions with low TFE concentrations at 37°C and at elevated temperatures (42 and 50°C). Comparison of the results obtained with the previously established 4E‐BP1 structure and the binding motif illustrates the ability of limited proteolysis in the presence of a folding assistant (TFE) to map the regions with high and low propensities to form a secondary structure revealing potential binding motifs inside the intrinsically disordered protein. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 591–602, 2014.     

Close encounters of the third kind: disordered domains and the interactions of proteins

Protein–protein interactions are thought to be mediated by domains, which are autonomous folding units of proteins. Recently, a second type of interaction has been suggested, mediated by short segments termed linear motifs, which are related to recognition elements of intrinsically disordered regions. Here, we propose a third kind of protein–protein recognition mechanism, mediated by disordered regions longer than 20–30 residues. Bioinformatics predictions and well‐characterized examples, such as the kinase‐inhibitory domain of Cdk inhibitors and the Wiskott–Aldrich syndrome protein (WASP)‐homology domain 2 of actin‐binding proteins, show that these disordered regions conform to the definition of domains rather than motifs, i.e., they represent functional, evolutionary, and structural units. Their functions are distinct from those of short motifs and ordered domains, and establish a third kind of interaction principle. With these points, we argue that these long disordered regions should be recognized as a distinct class of biologically functional protein domains.     

Impaired transport of intrinsically disordered proteins through the Sec61 and SecY translocon; implications for prion diseases

The prion protein (PrP) is composed of two major domains of similar size. The structured C-terminal domain contains three alpha-helical regions and a short two-stranded beta-sheet, while the N-terminal domain is intrinsically disordered. The analysis of PrP mutants with deletions in the C-terminal globular domain provided the first hint that intrinsically disordered domains are inefficiently transported into the endoplasmic reticulum through the Sec61 translocon. Interestingly, C-terminally truncated PrP mutants have been linked to inherited prion disease in humans and are characterized by inefficient ER import and the formation of neurotoxic PrP conformers. In a recent study we found that the Sec61 translocon in eukaryotic cells as well as the SecY translocon in bacteria is inherently deficient in translocating intrinsically disordered proteins. Moreover, our results suggest that translocon-associated components in eukaryotic cells enable the Sec61 complex to transport secretory proteins with extended unstructured domains such as PrP and shadoo.     

Multi-PDZ domain protein MUPP1 is a cellular target for both adenovirus E4-ORF1 and high-risk papillomavirus type 18 E6 oncoproteins

A general theme that has emerged from studies of DNA tumor viruses is that otherwise unrelated oncoproteins encoded by these viruses often target the same important cellular factors. Major oncogenic determinants for human adenovirus type 9 (Ad9) and high-risk human papillomaviruses (HPV) are the E4-ORF1 and E6 oncoproteins, respectively, and although otherwise unrelated, both of these viral proteins possess a functional PDZ domain-binding motif that is essential for their transforming activity and for binding to the PDZ domain-containing and putative tumor suppressor protein DLG. We report here that the PDZ domain-binding motifs of Ad9 E4-ORF1 and high-risk HPV-18 E6 also mediate binding to the widely expressed cellular factor MUPP1, a large multi-PDZ domain protein predicted to function as an adapter in signal transduction. With regard to the consequences of these interactions in cells, we showed that Ad9 E4-ORF1 aberrantly sequesters MUPP1 within the cytoplasm of cells whereas HPV-18 E6 targets this cellular protein for degradation. These effects were specific because mutant viral proteins unable to bind MUPP1 lack these activities. From these results, we propose that the multi-PDZ domain protein MUPP1 is involved in negatively regulating cellular proliferation and that the transforming activities of two different viral oncoproteins depend, in part, on their ability to inactivate this cellular factor.     

Protein intrinsic disorder and human papillomaviruses: increased amount of disorder in E6 and E7 oncoproteins from high risk HPVs

It is recognized now that many functional proteins or their long segments are devoid of stable secondary and/or tertiary structure and exist instead as very dynamic ensembles of conformations. They are known by different names including natively unfolded, intrinsically disordered, intrinsically unstructured, rheomorphic, pliable, and different combinations thereof. Many important functions and activities have been associated with these intrinsically disordered proteins (IDPs), including molecular recognition, signaling, and regulation. It is also believed that disorder of these proteins allows function to be readily modified through phosphorylation, acetylation, ubiquitination, hydroxylation, and proteolysis. Bioinformatics analysis revealed that IDPs comprise a large fraction of different proteomes. Furthermore, it is established that the intrinsic disorder is relatively abundant among cancer-related and other disease-related proteins and IDPs play a number of key roles in oncogenesis. There are more than 100 different types of human papillomaviruses (HPVs), which are the causative agents of benign papillomas/warts, and cofactors in the development of carcinomas of the genital tract, head and neck, and epidermis. With respect to their association with cancer, HPVs are grouped into two classes, known as low (e.g., HPV-6 and HPV-11) and high-risk (e.g., HPV-16 and HPV-18) types. The entire proteome of HPV includes six nonstructural proteins [E1, E2, E4, E5, E6, and E7 (the latter two are known to function as oncoproteins in the high-risk HPVs)] and two structural proteins (L1 and L2). To understand whether intrinsic disorder plays a role in the oncogenic potential of different HPV types, we have performed a detailed bioinformatics analysis of proteomes of high-risk and low-risk HPVs with the major focus on E6 and E7 oncoproteins. The results of this analysis are consistent with the conclusion that high-risk HPVs are characterized by the increased amount of intrinsic disorder in transforming proteins E6 and E7.