Conformational dynamics and structural plasticity play critical roles in the ubiquitin recognition of a UIM domain

A novel RNF125 variant associated with Tenorio syndrome alters ubiquitin chain binding

A key signalling pathway required for clearance of viruses from host cells relies on the receptor protein, retinoic acid-inducible gene I (RIG-I). The activity of RIG-I is tightly controlled, and once bound to viral dsRNA, addition of lysine 63-linked ubiquitin chains activates signalling. Meanwhile, the addition of lysine 48-linked ubiquitin chains to RIG-I is required to terminate signalling when the infection has been resolved. Really interesting new gene (RING) finger protein 125 (RNF125) is the E3 ligase responsible for addition of the ubiquitin chains that terminate signalling, with disruption of its function associated with Tenorio syndrome. Here we describe a novel RNF125 gene variant in an individual with clinical symptoms including intellectual disability, macrocephaly and congenital heart disease, consistent with Tenorio syndrome. The newly identified Tenorio syndrome-associated variant [(NM_017831.4):c.670G>C p.Glu224Gln] is the first to be found in the ubiquitin interaction motif (UIM) of RNF125. While the E3 ligase activity of this RNF125 variant is retained, it has an impaired ability to interact with lysine 63-linked ubiquitin chains. The function of the UIM in RNF125 is uncertain; however, this study suggests that the UIM binds lysine 63-linked ubiquitin chains, and that this interaction is required for the normal function of RNF125.

Transient exposure of a buried phosphorylation site in an autoinhibited protein

Autoinhibition is a mechanism used to regulate protein function, often by making functional sites inaccessible through the interaction with a cis-acting inhibitory domain. Such autoinhibitory domains often display a substantial degree of structural disorder when unbound, and only become structured in the inhibited state. These conformational dynamics make it difficult to study the structural origin of regulation, including effects of regulatory post-translational modifications. Here, we study the autoinhibition of the Dbl Homology domain in the protein Vav1 by the so-called acidic inhibitory domain. We use molecular simulations to study the process by which a mostly unstructured inhibitory domain folds upon binding and how transient exposure of a key buried tyrosine residue makes it accessible for phosphorylation. We show that the inhibitory domain, which forms a helix in the bound and inhibited stated, samples helical structures already before binding and that binding occurs via a molten-globule-like intermediate state. Together, our results shed light on key interactions that enable the inhibitory domain to sample a finely tuned equilibrium between an inhibited and a kinase-accessible state.

Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs

Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome.

Differential proteomics study of postharvest Volvariella volvacea during storage at 4 °C

The postharvest storage of Volvariella volvacea is an important factor limiting the industry development. Low-temperature storage is the traditional storage method used for most edible fungi, but V. volvacea undergoes autolysis at low temperature. To understand the molecular mechanism underlying the low-temperature autolysis of V. volvacea after harvesting, fruiting bodies of V. volvacea strain V23 were stored at 4 °C. Based on our previous study, in which the changes of morphological and physiological indexes during storage for 0, 6, 12, 24, 30, 36, 48 and 60 h were measured; four time points, namely, 0, 12, 24 and 60 h, were selected for this differential proteomics study. The proteomic changes in the postharvest storage samples were studied by isobaric tags for relative and absolute quantification-coupled two-dimensional liquid chromatography-tandem mass spectrometry (2D LC–MS/MS). A total of 2,063 proteins were identified, and 192 differentially expressed proteins (DEPs), including 24 up-regulated proteins and 168 down-regulated proteins, were detected after 12 h of storage. After 24 h of storage, 234 DEPs, including 48 up-regulated and 186 down-regulated proteins, were observed, and after 60 h, 415 DEPs, including 65 up-regulated proteins and 350 down-regulated proteins, were observed. An in-depth data analysis showed that the DEPs participated in various cellular processes, particularly metabolic processes. In this study, we combined Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses, and the results focused on oxidative phosphorylation and ubiquitin mediated proteolysis pathways. In addition, sdh2, uba1 and ubc1 was confirmed by quantitative real-time polymerase chain reaction, and the results showed that the expression of these genes were consistent with their protein level. Based on the literature and our results, it is speculated that the identified DEPs, such as ATP1, SDH2, COR1, UBA1, COX4, UBC1 and SKP1 play a key role in the low-temperature autolysis of V. volvacea.

Recognition of Poly-Ubiquitins by the Proteasome through Protein Refolding Guided by Electrostatic and Hydrophobic Interactions

Specificity of protein degradation by cellular proteasomes comes from tetra-ubiquitin recognition. We carry out molecular dynamics simulations to characterize how the ubiquitin receptor Rpn10 recognizes in the 26S proteasome K48-linked tetra-ubiquitin. In the binding pose, ubiquitin and Rpn10 interact primarily through hydrophobic patches. However, K48-linked tetra-ubiquitin mostly assumes a closed form in solution prior to binding, in which its hydrophobic patches are not exposed to solvent. Likewise, the hydrophobic ubiquitin interacting motifs (UIMs) of Rpn10 are mostly protected prior to binding. As a result, ubiquitin recognition in the proteasome requires refolding of both K48-linked tetra-ubiquitin and Rpn10. Simulations suggest that conserved complementary electrostatic patterns of Rpn10 and ubiquitins guide protein association (stage 1 in the recognition process), which induces refolding (stage 2), and then facilitates formation of hydrophobic contacts (stage 3). The simulations also explain why Rpn10 has a higher affinity for K48-linked tetra-ubiquitin than for mono-ubiquitin and K48-linked di- and tri-ubiquitins. Simulation results expand on the current view that the flexible arm of Rpn10 acts as an extended fragment of α-helices and flexible coils in the recognition process.

Dimeric interactions and complex formation using direct coevolutionary couplings

We develop a procedure to characterize the association of protein structures into homodimers using coevolutionary couplings extracted from Direct Coupling Analysis (DCA) in combination with Structure Based Models (SBM). Identification of dimerization contacts using DCA is more challenging than intradomain contacts since direct couplings are mixed with monomeric contacts. Therefore a systematic way to extract dimerization signals has been elusive. We provide evidence that the prediction of homodimeric complexes is possible with high accuracy for all the cases we studied which have rich sequence information. For the most accurate conformations of the structurally diverse dimeric complexes studied the mean and interfacial RMSDs are 1.95Å and 1.44Å, respectively. This methodology is also able to identify distinct dimerization conformations as for the case of the family of response regulators, which dimerize upon activation. The identification of dimeric complexes can provide interesting molecular insights in the construction of large oligomeric complexes and be useful in the study of aggregation related diseases like Alzheimer's or Parkinson's.

Evaluation of reduced point charge models of proteins through Molecular Dynamics simulations: application to the Vps27 UIM-1-Ubiquitin complex

Reduced point charge models of amino acids are designed, (i) from local extrema positions in charge density distribution functions built from the Poisson equation applied to smoothed molecular electrostatic potential (MEP) functions, and (ii) from local maxima positions in promolecular electron density distribution functions. Corresponding charge values are fitted versus all-atom Amber99 MEPs. To easily generate reduced point charge models for protein structures, libraries of amino acid templates are built. The program GROMACS is used to generate stable Molecular Dynamics trajectories of an Ubiquitin-ligand complex (PDB: 1Q0W), under various implementation schemes, solvation, and temperature conditions. Point charges that are not located on atoms are considered as virtual sites with a nul mass and radius. The results illustrate how the intra- and inter-molecular H-bond interactions are affected by the degree of reduction of the point charge models and give directions for their implementation; a special attention to the atoms selected to locate the virtual sites and to the Coulomb-14 interactions is needed. Results obtained at various temperatures suggest that the use of reduced point charge models allows to probe local potential hyper-surface minima that are similar to the all-atom ones, but are characterized by lower energy barriers. It enables to generate various conformations of the protein complex more rapidly than the all-atom point charge representation.

Development of inhibitors in the ubiquitination cascade

The ubiquitin proteasome system (UPS) is essential in regulating myriad aspects of protein functions. It is therefore a fundamentally important regulatory mechanism that impacts most if not all aspects of cellular processes. Indeed, malfunction of UPS components is implicated in human diseases such as neurodegenerative and immunological disorders and many cancers. The success of proteasome inhibitors in cancer therapy suggests that modulating enzymes in the ubiquitination cascade would be clinically important for therapeutic benefits. In this review, we summarize advances in developing inhibitors of a variety of UPS components. In particular, we highlight recent work done on the protein engineering of ubiquitin as modulators of the UPS, a novel approach that may shed light on innovative drug discovery in the future.

Bacterial expression and purification of the amyloidogenic peptide PAPf39 for multidimensional NMR spectroscopy

PAPf39 is a 39 residue peptide fragment from human prostatic acidic phosphatase that forms amyloid fibrils in semen. These fibrils have been implicated in facilitating HIV transmission. To enable structural studies of PAPf39 by NMR spectroscopy, efficient methods allowing the production of milligram quantities of isotopically labeled peptide are essential. Here, we report the high-yield expression and purification of uniformly (13)C- and (15)N-labeled PAPf39 peptide, through expression as a fusion to ubiquitin at the N-terminus and an intein at the C-terminus. This allows the study of the PAPf39 monomer conformational ensemble by NMR spectroscopy. To this end, we performed the NMR chemical shift assignment of the PAPf39 peptide in the monomeric state at low pH.

Conformational and dynamic changes at the interface contribute to ligand binding by ubiquitin

Ubiquitin interacts with numerous domains and motifs in its lifetime that vary in structure but bind the same hydrophobic patch. To identify the structural features of ubiquitin that make it an exceptional protein-protein interaction partner, we have studied the interaction of ubiquitin with the signal transducing adaptor molecule-1 ubiquitin interacting motif (UIM) using nuclear magnetic resonance. Our studies bring to light the role of the inherent backbone flexibility of ubiquitin in its interactions with a large array of binding partners, revealed from the changes in C(α) chemical shifts, backbone dynamics, and hydrogen bond lengths upon UIM binding. The crystal structures of ubiquitin complexes lend further support to our findings, underscoring the importance of the unique and flexible hydrogen bond network within ubiquitin and simultaneously providing insights into the nature of the slow motions. Taken together, our studies provide an in-depth view of the molecular changes associated with ligand recognition by ubiquitin.

Expanding the proteome: disordered and alternatively folded proteins

Proteins provide much of the scaffolding for life, as well as undertaking a variety of essential catalytic reactions. These characteristic functions have led us to presuppose that proteins are in general functional only when well structured and correctly folded. As we begin to explore the repertoire of possible protein sequences inherent in the human and other genomes, two stark facts that belie this supposition become clear: firstly, the number of apparent open reading frames in the human genome is significantly smaller than appears to be necessary to code for all of the diverse proteins in higher organisms, and secondly that a significant proportion of the protein sequences that would be coded by the genome would not be expected to form stable three-dimensional (3D) structures. Clearly the genome must include coding for a multitude of alternative forms of proteins, some of which may be partly or fully disordered or incompletely structured in their functional states. At the same time as this likelihood was recognized, experimental studies also began to uncover examples of important protein molecules and domains that were incompletely structured or completely disordered in solution, yet remained perfectly functional. In the ensuing years, we have seen an explosion of experimental and genome-annotation studies that have mapped the extent of the intrinsic disorder phenomenon and explored the possible biological rationales for its widespread occurrence. Answers to the question 'why would a particular domain need to be unstructured?' are as varied as the systems where such domains are found. This review provides a survey of recent new directions in this field, and includes an evaluation of the role not only of intrinsically disordered proteins but also of partially structured and highly dynamic members of the disorder-order continuum.

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.

In silico elucidation of the recognition dynamics of ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain.

Molecular recognition plays a central role in many biological processes, ensuring specific and efficient interaction between binding partners. Various models for describing the mechanisms of molecular recognition have been proposed, but the validation of these models has been traditionally difficult due to the transient and complex nature of the dynamic recognition process. In the present study, we aim at visually characterizing the mutual interplay between human ubiquitin and its ligands via microsecond time scale molecular dynamics simulation, which is validated rigorously against experimental NMR data. Taking advantage of Boltzmann sampling of molecular dynamics snapshots, we statistically reweight the populations of ubiquitin in the presence of its ligand molecule at intermediate distance range (1–2 nm) to examine the population redistribution mechanisms. These results offer new atomistic insights into this vital protein-protein recognition event.

The structure of human ubiquitin in 2-methyl-2,4-pentanediol: a new conformational switch

A new crystal structure of human ubiquitin is reported at 1.8 Å resolution. Compared with the other known crystal structure or the solution NMR structure of monomeric human ubiquitin, this new structure is similar in its overall fold but differs with respect to the conformation of the backbone in a surface-exposed region. The conformation reported here resembles conformations previously seen in complex with deubiquinating enzymes, wherein the Asp52/Gly53 main chain and Glu24 side chain move. This movement exposes the backbone carbonyl of Asp52 to the exterior of the molecule, making it possible to engage in hydrogen-bond contacts with neighboring molecules, rather than in an internal hydrogen bond with the backbone of Glu24. This particular crystal form of ubiquitin has been used in a large number of solid state NMR studies. The structure described here elucidates the origin of many of the chemical shift differences comparing solution and solid state studies.

WD40 repeat propellers define a ubiquitin-binding domain that regulates turnover of F box proteins

WD40-repeat β-propellers are found in a wide range of proteins involved in distinct biological activities. We define a large subset of WD40 β-propellers as a class of ubiquitin-binding domains. Using the β-propeller from Doa1/Ufd3 as a paradigm, we find the conserved top surface of the Doa1 β-propeller binds the hydrophobic patch of ubiquitin centered on residues I44, L8, and V70. Mutations that disrupt ubiquitin binding abrogate Doa1 function, demonstrating the importance of this interaction. We further demonstrate that WD40 β-propellers from a functionally diverse set of proteins bind ubiquitin in a similar fashion. This set includes members of the F box family of SCF ubiquitin E3 ligase adaptors. Using mutants defective in binding, we find that ubiquitin interaction by the F box protein Cdc4 promotes its autoubiquitination and turnover. Collectively, our results reveal a molecular mechanism that may account for how ubiquitin controls a broad spectrum of cellular activities.