Various strategies of using residual dipolar couplings in NMR-driven protein docking: Application to Lys48-linked di-ubiquitin and validation against 15N-relaxation data

Mutational and Environmental Effects on the Dynamic Conformational Distributions of Lys48-Linked Ubiquitin Chains

In multidomain proteins, individual domains connected by flexible linkers are dynamically rearranged upon ligand binding and sensing changes in environmental factors, such as pH and temperature. Here, we characterize dynamic domain rearrangements of Lys48-linked ubiquitin (Ub) chains as models of multidomain proteins in which molecular surfaces mediating intermolecular interactions are involved in intramolecular domain-domain interactions. Using NMR and other biophysical techniques, we characterized dynamic conformational interconversions of diUb between open and closed states regarding solvent exposure of the hydrophobic surfaces of each Ub unit, which serve as binding sites for various Ub-interacting proteins. We found that the hydrophobic Ub-Ub interaction in diUb was reinforced by cysteine substitution of Lys48 of the distal Ub unit because of interaction between the cysteinyl thiol group and the C-terminal segment of the proximal Ub unit. In contrast, the replacement of the isopeptide linker with an artificial ethylenamine linker minimally affected the conformational distributions. Furthermore, we demonstrated that the mutational modification allosterically impacted the exposure of the most distal Ub unit in triUb. Thus, the conformational interconversion of Ub chains offers a unique design framework in Ub-based protein engineering not only for developing biosensing probes but also for allowing new opportunities for the allosteric regulation of multidomain proteins.

Modeling and Structure Determination of Homo-Oligomeric Proteins: An Overview of Challenges and Current Approaches

Protein homo-oligomerization is a very common phenomenon, and approximately half of proteins form homo-oligomeric assemblies composed of identical subunits. The vast majority of such assemblies possess internal symmetry which can be either exploited to help or poses challenges during structure determination. Moreover, aspects of symmetry are critical in the modeling of protein homo-oligomers either by docking or by homology-based approaches. Here, we first provide a brief overview of the nature of protein homo-oligomerization. Next, we describe how the symmetry of homo-oligomers is addressed by crystallographic and non-crystallographic symmetry operations, and how biologically relevant intermolecular interactions can be deciphered from the ordered array of molecules within protein crystals. Additionally, we describe the most important aspects of protein homo-oligomerization in structure determination by NMR. Finally, we give an overview of approaches aimed at modeling homo-oligomers using computational methods that specifically address their internal symmetry and allow the incorporation of other experimental data as spatial restraints to achieve higher model reliability.

An Integrative Approach to Determine 3D Protein Structures Using Sparse Paramagnetic NMR Data and Physical Modeling

Paramagnetic nuclear magnetic resonance (NMR) methods have emerged as powerful tools for structure determination of large, sparsely protonated proteins. However traditional applications face several challenges, including a need for large datasets to offset the sparsity of restraints, the difficulty in accounting for the conformational heterogeneity of the spin-label, and noisy experimental data. Here we propose an integrative approach to structure determination combining sparse paramagnetic NMR with physical modelling to infer approximate protein structural ensembles. We use calmodulin in complex with the smooth muscle myosin light chain kinase peptide as a model system. Despite acquiring data from samples labeled only at the backbone amide positions, we are able to produce an ensemble with an average RMSD of ∼2.8 Å from a reference X-ray crystal structure. Our approach requires only backbone chemical shifts and measurements of the paramagnetic relaxation enhancement and residual dipolar couplings that can be obtained from sparsely labeled samples.

A novel recognition site for polyubiquitin and ubiquitin-like signals in an unexpected region of proteasomal subunit Rpn1

The ubiquitin (Ub)–proteasome system is the primary mechanism for maintaining protein homeostasis in eukaryotes, yet the underlying signaling events and specificities of its components are poorly understood. Proteins destined for degradation are tagged with covalently linked polymeric Ub chains and subsequently delivered to the proteasome, often with the assistance of shuttle proteins that contain Ub-like domains. This degradation pathway is riddled with apparent redundancy—in the form of numerous polyubiquitin chains of various lengths and distinct architectures, multiple shuttle proteins, and at least three proteasomal receptors. Moreover, the largest proteasomal receptor, Rpn1, contains one known binding site for polyubiquitin and shuttle proteins, although several studies have recently proposed the existence of an additional uncharacterized site. Here, using a combination of NMR spectroscopy, photocrosslinking, mass spectrometry, and mutagenesis, we show that Rpn1 does indeed contain another recognition site that exhibits affinities and binding preferences for polyubiquitin and Ub-like signals comparable to those of the known binding site in Rpn1. Surprisingly, this novel site is situated in the N-terminal section of Rpn1, a region previously surmised to be devoid of functionality. We identified a stretch of adjacent helices as the location of this previously uncharacterized binding site, whose spatial proximity and similar properties to the known binding site in Rpn1 suggest the possibility of multivalent signal recognition across the solvent-exposed surface of Rpn1. These findings offer new mechanistic insights into signal recognition processes that are at the core of the Ub–proteasome system.

Aromatic Interactions Drive the Coupled Folding and Binding of the Intrinsically Disordered Sesbania mosaic Virus VPg Protein

The plant Sesbania mosaic virus [a (+)-ssRNA sobemovirus] VPg protein is intrinsically disordered in solution. For the virus life cycle, the VPg protein is essential for replication and for polyprotein processing that is carried out by a virus-encoded protease. The nuclear magnetic resonance (NMR)-derived tertiary structure of the protease-bound VPg shows it to have a novel tertiary structure with an α-β-β-β topology. The quaternary structure of the high-affinity protease-VPg complex (≈27 kDa) has been determined using HADDOCK protocols with NMR (residual dipolar coupling, dihedral angle, and nuclear Overhauser enhancement) restraints and mutagenesis data as inputs. The geometry of the complex is in excellent agreement with long-range orientational restraints such as residual dipolar couplings and ring-current shifts. A "vein" of aromatic residues on the protease surface is pivotal for the folding of VPg via intermolecular edge-to-face π···π stacking between Trp²⁷¹ and Trp³⁶⁸ of the protease and VPg, respectively, and for the CH···π interactions between Leu³⁶¹ of VPg and Trp²⁷¹ of the protease. The structure of the protease-VPg complex provides a molecular framework for predicting sites of important posttranslational modifications such as RNA linkage and phosphorylation and a better understanding of the coupled folding upon binding of intrinsically disordered proteins. The structural data presented here augment the limited structural data available on viral proteins, given their propensity for structural disorder.

Computational structure modeling for diverse categories of macromolecular interactions

Computational protein-protein docking is one of the most intensively studied topics in structural bioinformatics. The field has made substantial progress through over three decades of development. The development began with methods for rigid-body docking of two proteins, which have now been extended in different directions to cover the various macromolecular interactions observed in a cell. Here, we overview the recent developments of the variations of docking methods, including multiple protein docking, peptide-protein docking, and disordered protein docking methods.

Machine Learning Driven Analysis of Large Scale Simulations Reveals Conformational Characteristics of Ubiquitin Chains

Understanding the conformational characteristics of protein complexes in solution is crucial for a deeper insight in their biological function. Molecular dynamics simulations performed on high performance computing plants and with modern simulation techniques can be used to obtain large data sets that contain conformational and thermodynamic information about biomolecular systems. While this can in principle give a detailed picture of protein-protein interactions in solution and therefore complement experimental data, it also raises the challenge of processing exceedingly large high-dimensional data sets with several million samples. Here we present a novel method for the characterization of protein-protein interactions, which combines a neural network based dimensionality reduction technique to obtain a two-dimensional representation of the conformational space with a density based clustering algorithm for state detection and a metric which assesses the (dis)similarity between different conformational spaces. This method is highly scalable and therefore makes the analysis of massive data sets computationally tractable. We demonstrate the power of this approach to large scale data analysis by characterizing highly dynamic conformational phase spaces of differently linked ubiquitin (Ub) oligomers from coarse-grained simulations. We are able to extract a protein-protein interaction model for two unlinked Ub proteins which is then used to determine how the Ub-Ub interaction pattern is altered in Ub oligomers by the introduction of a covalent linkage. We find that the Ub chain conformational ensemble depends highly on the linkage type and for some cases also on the Ub chain length. By this, we obtain insight into the conformational characteristics of different Ub chains and how this may contribute to linkage type and chain length specific recognition.

Diubiquitin-Based NMR Analysis: Interactions Between Lys6-Linked diUb and UBA Domain of UBXN1

Ubiquitination is a process in which a protein is modified by the covalent attachment of the C-terminal carboxylic acid of ubiquitin (Ub) to the ε-amine of lysine or N-terminal methionine residue of a substrate protein or another Ub molecule. Each of the seven internal lysine residues and the N-terminal methionine residue of Ub can be linked to the C-terminus of another Ub moiety to form 8 distinct Ub linkages and the resulting differences in linkage types elicit different Ub signaling pathways. Cellular responses are triggered when proteins containing ubiquitin-binding domains (UBDs) recognize and bind to specific polyUb linkage types. To get more insight into the differences between polyUb chains, all of the seven lysine-linked di-ubiquitin molecules (diUbs) were prepared and used as a model to study their structural conformations in solution using NMR spectroscopy. We report the synthesis of diUb molecules, fully ¹⁵N-labeled on the distal (N-terminal) Ub moiety and revealed their structural orientation with respect to the proximal Ub. As expected, the diUb molecules exist in different conformations in solution, with multiple conformations known to exist for K6-, K48-, and K63-linked diUb molecules. These multiple conformations allow structural flexibility in binding with UBDs thereby inducing unique responses. One of the well-known but poorly understood UBD-Ub interaction is the recognition of K6 polyubiquitin by the ubiquitin-associated (UBA) domain of UBXN1 in the BRCA-mediated DNA repair pathway. Using our synthetic ¹⁵N-labeled diUbs, we establish here how a C-terminally extended UBA domain of UBXN1 confers specificity to K6 diUb while the non-extended version of the domain does not show any linkage preference. We show that the two distinct conformations of K6 diUb that exist in solution converge into a single conformation upon binding to this extended form of the UBA domain of the UBXN1 protein. It is likely that more of such extended UBA domains exist in nature and can contribute to linkage-specificity in Ub signaling. The isotopically labeled diUb compounds described here and the use of NMR to study their interactions with relevant partner molecules will help accelerate our understanding of Ub signaling pathways.

Branching via K11 and K48 Bestows Ubiquitin Chains with a Unique Interdomain Interface and Enhanced Affinity for Proteasomal Subunit Rpn1

Post-translational substrate modification with ubiquitin is essential for eukaryotic cellular signaling. Polymeric ubiquitin chains are assembled with specific architectures, which convey distinct signaling outcomes depending on the linkages involved. Recently, branched K11/K48-linked polyubiquitins were shown to enhance proteasomal degradation during mitosis. To better understand the underlying structural mechanisms, we determined the crystal and NMR structures of branched K11/K48-linked tri-ubiquitin and discovered a previously unobserved interdomain interface between the distal ubiquitins. Small-angle neutron scattering and site-directed mutagenesis corroborated the presence of this interface, which we hypothesized to be influential in the physiological role of branched K11/K48-linked chains. Yet, experiments probing polyubiquitin interactions-deubiquitination assays, binding to proteasomal shuttle hHR23A-showed negligible differences between branched K11/K48-linked tri-ubiquitin and related di-ubiquitins. However, significantly stronger binding affinity for branched K11/K48-linked tri-ubiquitin was observed with proteasomal subunit Rpn1, thereby suggesting a functional impact of this interdomain interface and pinpointing the mechanistic site of enhanced degradation.

Enhancing NMR derived ensembles with kinetics on multiple timescales

Nuclear magnetic resonance (NMR) has the unique advantage of elucidating the structure and dynamics of biomolecules in solution at physiological temperatures, where they are in constant movement on timescales from picoseconds to milliseconds. Such motions have been shown to be critical for enzyme catalysis, allosteric regulation, and molecular recognition. With NMR being particularly sensitive to these timescales, detailed information about the kinetics can be acquired. However, nearly all methods of NMR-based biomolecular structure determination neglect kinetics, which introduces a large approximation to the underlying physics, limiting both structural resolution and the ability to accurately determine molecular flexibility. Here we present the Kinetic Ensemble approach that uses a hierarchy of interconversion rates between a set of ensemble members to rigorously calculate Nuclear Overhauser Effect (NOE) intensities. It can be used to simultaneously refine both temporal and structural coordinates. By generalizing ideas from the extended model free approach, the method can analyze the amplitudes and kinetics of motions anywhere along the backbone or side chains. Furthermore, analysis of a large set of crystal structures suggests that NOE data contains a surprising amount of high-resolution information that is better modeled using our approach. The Kinetic Ensemble approach provides the means to unify numerous types of experiments under a single quantitative framework and more fully characterize and exploit kinetically distinct protein states. While we apply the approach here to the protein ubiquitin and cross validate it with previously derived datasets, the approach can be applied to any protein for which NOE data is available.

Towards a molecular basis of ubiquitin signaling: A dual-scale simulation study of ubiquitin dimers

Covalent modification of proteins by ubiquitin or ubiquitin chains is one of the most prevalent post-translational modifications in eukaryotes. Different types of ubiquitin chains are assumed to selectively signal respectively modified proteins for different fates. In support of this hypothesis, structural studies have shown that the eight possible ubiquitin dimers adopt different conformations. However, at least in some cases, these structures cannot sufficiently explain the molecular basis of the selective signaling mechanisms. This indicates that the available structures represent only a few distinct conformations within the entire conformational space adopted by a ubiquitin dimer. Here, molecular simulations on different levels of resolution can complement the structural information. We have combined exhaustive coarse grained and atomistic simulations of all eight possible ubiquitin dimers with a suitable dimensionality reduction technique and a new method to characterize protein-protein interfaces and the conformational landscape of protein conjugates. We found that ubiquitin dimers exhibit characteristic linkage type-dependent properties in solution, such as interface stability and the character of contacts between the subunits, which can be directly correlated with experimentally observed linkage-specific properties.

Post-translational modification of proteins by covalent attachment of ubiquitin is a key cellular process, regulating for example the fate and recycling of proteins. We present a new method to combine multiscale simulation with advanced analysis methods to characterize the states of ubiquitin-ubiquitin conjugates. We found that the linkage position affects the conformational space of ubiquitin dimers, determining the number and stability of relevant states, the character of subunit contacts and the nature of the surface exposed to possible binding partners.

The Ball and Chain of Polyubiquitin Structures

Ubiquitylation is a post-translational modification implicated in several different cellular pathways. The possibility of forming chains through covalent crosslinking between any of the seven lysines, or the initial methionine, and the C terminus of another moiety provides ubiquitin (Ub) with special flexibility in its function in signalling. Here, we review the knowledge accumulated over the past several years about the functions and structural features of polyUb chains. This analysis reveals the need to understand further the functional role of some of the linkages and the structural code that determines recognition of polyUbs by protein partners.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

Cryo-EM Data Are Superior to Contact and Interface Information in Integrative Modeling

Protein-protein interactions carry out a large variety of essential cellular processes. Cryo-electron microscopy (cryo-EM) is a powerful technique for the modeling of protein-protein interactions at a wide range of resolutions, and recent developments have caused a revolution in the field. At low resolution, cryo-EM maps can drive integrative modeling of the interaction, assembling existing structures into the map. Other experimental techniques can provide information on the interface or on the contacts between the monomers in the complex. This inevitably raises the question regarding which type of data is best suited to drive integrative modeling approaches. Systematic comparison of the prediction accuracy and specificity of the different integrative modeling paradigms is unavailable to date. Here, we compare EM-driven, interface-driven, and contact-driven integrative modeling paradigms. Models were generated for the protein docking benchmark using the ATTRACT docking engine and evaluated using the CAPRI two-star criterion. At 20 Å resolution, EM-driven modeling achieved a success rate of 100%, outperforming the other paradigms even with perfect interface and contact information. Therefore, even very low resolution cryo-EM data is superior in predicting heterodimeric and heterotrimeric protein assemblies. Our study demonstrates that a force field is not necessary, cryo-EM data alone is sufficient to accurately guide the monomers into place. The resulting rigid models successfully identify regions of conformational change, opening up perspectives for targeted flexible remodeling.

Comparison of native and non-native ubiquitin oligomers reveals analogous structures and reactivities

Ubiquitin (Ub) chains regulate a wide range of biological processes, and Ub chain connectivity is a critical determinant of the many regulatory roles that this post-translational modification plays in cells. To understand how distinct Ub chains orchestrate different biochemical events, we and other investigators have developed enzymatic and non-enzymatic methods to synthesize Ub chains of well-defined length and connectivity. A number of chemical approaches have been used to generate Ub oligomers connected by non-native linkages; however, few studies have examined the extent to which non-native linkages recapitulate the structural and functional properties associated with native isopeptide bonds. Here, we compare the structure and function of Ub dimers bearing native and non-native linkages. Using small-angle X-ray scattering (SAXS) analysis, we show that scattering profiles for the two types of dimers are similar. Moreover, using an experimental structural library and atomistic simulations to fit the experimental SAXS profiles, we find that the two types of Ub dimers can be matched to analogous structures. An important application of non-native Ub oligomers is to probe the activity and selectivity of deubiquitinases. Through steady-state kinetic analyses, we demonstrate that different families of deubiquitinases hydrolyze native and non-native isopeptide linkages with comparable efficiency and selectivity. Considering the significant challenges associated with building topologically diverse native Ub chains, our results illustrate that chains harboring non-native linkages can serve as surrogate substrates for explorations of Ub function.

Predicted binding site information improves model ranking in protein docking using experimental and computer-generated target structures

Background

Protein-protein interactions (PPIs) mediate the vast majority of biological processes, therefore, significant efforts have been directed to investigate PPIs to fully comprehend cellular functions. Predicting complex structures is critical to reveal molecular mechanisms by which proteins operate. Despite recent advances in the development of new methods to model macromolecular assemblies, most current methodologies are designed to work with experimentally determined protein structures. However, because only computer-generated models are available for a large number of proteins in a given genome, computational tools should tolerate structural inaccuracies in order to perform the genome-wide modeling of PPIs.

Results

To address this problem, we developed eRank^PPI, an algorithm for the identification of near-native conformations generated by protein docking using experimental structures as well as protein models. The scoring function implemented in eRank^PPI employs multiple features including interface probability estimates calculated by eFindSite^PPI and a novel contact-based symmetry score. In comparative benchmarks using representative datasets of homo- and hetero-complexes, we show that eRank^PPI consistently outperforms state-of-the-art algorithms improving the success rate by ~10 %.

Conclusions

eRank^PPI was designed to bridge the gap between the volume of sequence data, the evidence of binary interactions, and the atomic details of pharmacologically relevant protein complexes. Tolerating structure imperfections in computer-generated models opens up a possibility to conduct the exhaustive structure-based reconstruction of PPI networks across proteomes. The methods and datasets used in this study are available at www.brylinski.org/erankppi.

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.

Ubiquitination of the Dishevelled DIX domain blocks its head-to-tail polymerization

Dishevelled relays Wnt signals from the plasma membrane to different cytoplasmic effectors. Its signalling activity depends on its DIX domain, which undergoes head-to-tail polymerization to assemble signalosomes. The DIX domain is ubiquitinated in vivo at multiple lysines, which can be antagonized by various deubiquitinases (DUBs) including the CYLD tumour suppressor that attenuates Wnt signalling. Here, we generate milligram quantities of pure human Dvl2 DIX domain mono-ubiquitinated at two lysines (K54 and K58) by genetically encoded orthogonal protection with activated ligation (GOPAL), to investigate their effect on DIX polymerization. We show that the ubiquitination of DIX at K54 blocks its polymerization in solution, whereas DIX58-Ub remains oligomerization-competent. DUB profiling identified 28 DUBs that cleave DIX-ubiquitin conjugates, half of which prefer, or are specific for, DIX54-Ub, including Cezanne and CYLD. These DUBs thus have the potential to promote Dvl polymerization and signalosome formation, rather than antagonize it as previously thought for CYLD.

Information-driven structural modelling of protein-protein interactions

Protein-protein docking aims at predicting the three-dimensional structure of a protein complex starting from the free forms of the individual partners. As assessed in the CAPRI community-wide experiment, the most successful docking algorithms combine pure laws of physics with information derived from various experimental or bioinformatics sources. Of these so-called "information-driven" approaches, HADDOCK stands out as one of the most successful representatives. In this chapter, we briefly summarize which experimental information can be used to drive the docking prediction in HADDOCK, and then focus on the docking protocol itself. We discuss and illustrate with a tutorial example a "classical" protein-protein docking prediction, as well as more recent developments for modelling multi-body systems and large conformational changes.

Recovering a representative conformational ensemble from underdetermined macromolecular structural data

Structural analysis of proteins and nucleic acids is complicated by their inherent flexibility, conferred, for example, by linkers between their contiguous domains. Therefore, the macromolecule needs to be represented by an ensemble of conformations instead of a single conformation. Determining this ensemble is challenging because the experimental data are a convoluted average of contributions from multiple conformations. As the number of the ensemble degrees of freedom generally greatly exceeds the number of independent observables, directly deconvolving experimental data into a representative ensemble is an ill-posed problem. Recent developments in sparse approximations and compressive sensing have demonstrated that useful information can be recovered from underdetermined (ill-posed) systems of linear equations by using sparsity regularization. Inspired by these advances, we designed the Sparse Ensemble Selection (SES) method for recovering multiple conformations from a limited number of observations. SES is more general and accurate than previously published minimum-ensemble methods, and we use it to obtain representative conformational ensembles of Lys48-linked diubiquitin, characterized by the residual dipolar coupling data measured at several pH conditions. These representative ensembles are validated against NMR chemical shift perturbation data and compared to maximum-entropy results. The SES method reproduced and quantified the previously observed pH dependence of the major conformation of Lys48-linked diubiquitin, and revealed lesser-populated conformations that are preorganized for binding known diubiquitin receptors, thus providing insights into possible mechanisms of receptor recognition by polyubiquitin. SES is applicable to any experimental observables that can be expressed as a weighted linear combination of data for individual states.

Structural analysis of poly-SUMO chain recognition by the RNF4-SIMs domain

The E3 ubiquitin ligase RNF4 (RING finger protein 4) contains four tandem SIM [SUMO (small ubiquitin-like modifier)-interaction motif] repeats for selective interaction with poly-SUMO-modified proteins, which it targets for degradation. We employed a multi-faceted approach to characterize the structure of the RNF4-SIMs domain and the tetra-SUMO2 chain to elucidate the interaction between them. In solution, the SIM domain was intrinsically disordered and the linkers of the tetra-SUMO2 were highly flexible. Individual SIMs of the RNF4-SIMs domains bind to SUMO2 in the groove between the β2-strand and the α1-helix parallel to the β2-strand. SIM2 and SIM3 bound to SUMO with a high affinity and together constituted the recognition module necessary for SUMO binding. SIM4 alone bound to SUMO with low affinity; however, its contribution to tetra-SUMO2 binding avidity is comparable with that of SIM3 when in the RNF4-SIMs domain. The SAXS data of the tetra-SUMO2-RNF4-SIMs domain complex indicate that it exists as an ordered structure. The HADDOCK model showed that the tandem RNF4-SIMs domain bound antiparallel to the tetra-SUMO2 chain orientation and wrapped around the SUMO protamers in a superhelical turn without imposing steric hindrance on either molecule.

Modeling protein-protein complexes using the HADDOCK webserver "modeling protein complexes with HADDOCK"

Protein-protein interactions lie at the heart of most cellular processes. Determining their high-resolution structures by experimental methods is a nontrivial task, which is why complementary computational approaches have been developed over the years. To gain structural and dynamical insight on an atomic scale in these interactions, computational modeling must often be complemented by low-resolution experimental information. For this purpose, we developed the user-friendly HADDOCK webserver, the interface to our biomolecular docking program, which can make use of a variety of low-resolution data to drive the docking process. In this chapter, we explain the use of the HADDOCK webserver based on the real-life Lys48-linked di-ubiquitin case, which led to the 2BGF PDB model. We demonstrate the use of chemical shift perturbation data in combination with residual dipolar couplings and further highlight a few other cases where our software was successfully used. The HADDOCK webserver is available to the science community for free at haddock.science.uu.nl/services/HADDOCK.

Discovery of a novel small molecule inhibitor targeting the frataxin/ubiquitin interaction via structure-based virtual screening and bioassays

Friedreich's ataxia (FRDA) is an autosomal recessive neuro- and cardiodegenerative disorder for which there are no proven effective treatments. FRDA is caused by decreased expression and/or function of the mitochondrial protein frataxin. Here, we report findings that frataxin is degraded via the ubiquitin-proteasomal pathway and that it is ubiquitinated at residue K(147) in Calu-6 cells. A theoretical model of the frataxin-K(147)/Ub complex, constructed by combining bioinformatics interface predictions with information-driven docking, revealed a hitherto unnoticed, potential ubiquitin-binding domain in frataxin. Through structure-based virtual screening and cell-based assays, we discovered a novel small molecule (compound (+)-11) able to prevent frataxin ubiquitination and degradation. (+)-11 was synthesized and tested for specific binding to frataxin by an UF-LC/MS based ligand-binding assay. Follow-up scaffold-based searches resulted in the identification of a lead series with micromolar activity in disrupting the frataxin/Ub interaction. This study also suggests that frataxin could be a potential target for FRDA drug development.

Enzymes of ubiquitination and deubiquitination

Ubiquitination, the covalent attachment of the small protein modifier ubiquitin to a substrate protein is involved in virtually all cellular processes by mediating the regulated degradation of proteins. Aside from proteasomal degradation, ubiquitination plays important roles in transcriptional regulation, protein trafficking, including endocytosis and lysosomal targeting, and activation of kinases involved in signalling processes. A three-tiered enzymatic cascade consisting of E1 or ubiquitin-activating enzyme, E2 or ubiquitin-conjugating enzyme, and E3, or ubiquitin ligases, is necessary to achieve the many forms of ubiquitination known to date. In this chapter, we summarize the current knowledge on the enzymatic machinery necessary for ubiquitin activation and ligation, as well as its removal, and provide some insight into the complexity of regulatory processes governed by ubiquitination.

Polymerization behavior of a bifunctional ubiquitin monomer as a function of the nucleophile site and folding conditions

Biopolymers with repeating modules composed of either folded peptides or tertiary protein domains are considered some of the basic biomaterials that nature has evolved to optimize for energy efficient synthesis and unique functions. Such biomaterials continue to inspire scientists to mimic their exceptional properties and the ways that nature adopts to prepare them. Ubiquitin chains represent another example of nature's approach to use a protein-repeating module to prepare functionally important biopolymers. In the current work, we utilize a novel synthetic strategy to prepare bifunctional ubiquitin monomers having a C-terminal thioester and a nucleophilic 1,2-aminothiol at a desired position to examine their polymerization products under different conditions. Our study reveals that such analogues, when subjected to polymerization conditions under different folding states, afford distinct patterns of polymerization products where both the dynamic and the tertiary structures of the chains play important roles in such processes. Moreover, we also show that the presence of a specific ubiquitin-binding domain, which binds specifically to some of these chains, could interfere selectively with the polymerization outcome. Our study represents the first example of examining the polymerization of designed and synthetic repeating modules based on tertiary protein domains and affords early lessons in the design and synthesis of biomaterial. In regards to the ubiquitin system, our study may have implications on the ease of synthesis of ubiquitin chains with varying lengths and types for structural and functional analyses. Importantly, such an approach could also assist in understanding the enzymatic machinery and the factors controlling the assembly of these chains with a desired length.

NMR structure of a heterodimeric SAM:SAM complex: characterization and manipulation of EphA2 binding reveal new cellular functions of SHIP2

The sterile alpha motif (SAM) for protein-protein interactions is encountered in over 200 proteins, but the structural basis for its interactions is just becoming clear. Here we solved the structure of the EphA2-SHIP2 SAM:SAM heterodimeric complex by use of NMR restraints from chemical shift perturbations, NOE and RDC experiments. Specific contacts between the protein surfaces differ significantly from a previous model and other SAM:SAM complexes. Molecular dynamics and docking simulations indicate fluctuations in the complex toward alternate, higher energy conformations. The interface suggests that EphA family members bind to SHIP2 SAM, whereas EphB members may not; correspondingly, we demonstrate binding of EphA1, but not of EphB2, to SHIP2. A variant of EphB2 SAM was designed that binds SHIP2. Functional characterization of a mutant EphA2 compromised in SHIP2 binding reveals two previously unrecognized functions of SHIP2 in suppressing ligand-induced activation of EphA2 and in promoting receptor coordinated chemotactic cell migration.

Structural and biochemical studies of the open state of Lys48-linked diubiquitin

Ubiquitin (Ub) is a small protein highly conserved among eukaryotes and involved in practically all aspects of eukaryotic cell biology. Polymeric chains assembled from covalently-linked Ub monomers function as molecular signals in the regulation of a host of cellular processes. Our previous studies have shown that the predominant state of Lys48-linked di- and tetra-Ub chains at near-physiological conditions is a closed conformation, in which the Ub-Ub interface is formed by the hydrophobic surface residues of the adjacent Ub units. Because these very residues are involved in (poly)Ub interactions with the majority of Ub-binding proteins, their sequestration at the Ub-Ub interface renders the closed conformation of polyUb binding incompetent. Thus the existence of open conformation(s) and the interdomain motions opening and closing the Ub-Ub interface is critical for the recognition of Lys48-linked polyUb by its receptors. Knowledge of the conformational properties of a polyUb signal is essential for our understanding of its specific recognition by various Ub-receptors. Despite their functional importance, open states of Lys48-linked chains are poorly characterized. Here we report a crystal structure of the open state of Lys48-linked di-Ub. Moreover, using NMR, we examined interactions of the open state of this chain (at pH4.5) with a Lys48-linkage-selective receptor, the UBA2 domain of a shuttle protein hHR23a. Our results show that di-Ub binds UBA2 in the same mode and with comparable affinity as the closed state. Our data suggest a mechanism for polyUb signal recognition, whereby Ub-binding proteins select specific conformations out of the available ensemble of polyUb chain conformations. This article is part of a Special Issue entitled: Ubiquitin Drug Discovery and Diagnostics.

Synergistic applications of MD and NMR for the study of biological systems

Modern biological sciences are becoming more and more multidisciplinary. At the same time, theoretical and computational approaches gain in reliability and their field of application widens. In this short paper, we discuss recent advances in the areas of solution nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations that were made possible by the combination of both methods, that is, through their synergistic use. We present the main NMR observables and parameters that can be computed from simulations, and how they are used in a variety of complementary applications, including dynamics studies, model-free analysis, force field validation, and structural studies.

Structural model of the TRPP2/PKD1 C-terminal coiled-coil complex produced by a combined computational and experimental approach

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in TRPP2 and PKD1, which form an ion channel/receptor complex containing three TRPP2 and one PKD1. A TRPP2 C-terminal coiled-coil trimer, critical for the assembly of this complex, associates with a single PKD1 C-terminal coiled-coil. Many ADPKD pathogenic mutations result in the abolishment of the TRPP2/PKD1 coiled-coil complex. To gain molecular and functional insights into this heterotetrameric complex, we computationally constructed a structural model by using a two-step docking strategy, based on a known crystal structure of the TRPP2 coiled-coil trimer. The model shows that this tetrameric complex has a novel di-trimer configuration: An upstream trimer made of three TRPP2 helices and a downstream trimer made of two TRPP2 helices and one PKD1 helix. Mutagenesis and biochemical analysis identified critical TRPP2/PKD1 interface contacts essential for the heteromeric coiled-coil complex. Mutation of these interface positions in the full-length proteins showed that these interactions were critical for the assembly of the full-length complex in cells. Our results provide a means to specifically weaken the TRPP2 and PKD1 association, thus facilitating future in vitro and in vivo studies on the functional importance of this association.

Fast approximations of the rotational diffusion tensor and their application to structural assembly of molecular complexes

We present and evaluate a rigid-body, deterministic, molecular docking method, called ELMDOCK, that relies solely on the three-dimensional structure of the individual components and the overall rotational diffusion tensor of the complex, obtained from nuclear spin-relaxation measurements. We also introduce a docking method, called ELMPATIDOCK, derived from ELMDOCK and based on the new concept of combining the shape-related restraints from rotational diffusion with those from residual dipolar couplings, along with ambiguous contact/interface-related restraints obtained from chemical shift perturbations. ELMDOCK and ELMPATIDOCK use two novel approximations of the molecular rotational diffusion tensor that allow computationally efficient docking. We show that these approximations are accurate enough to properly dock the two components of a complex without the need to recompute the diffusion tensor at each iteration step. We analyze the accuracy, robustness, and efficiency of these methods using synthetic relaxation data for a large variety of protein-protein complexes. We also test our method on three protein systems for which the structure of the complex and experimental relaxation data are available, and analyze the effect of flexible unstructured tails on the outcome of docking. Additionally, we describe a method for integrating the new approximation methods into the existing docking approaches that use the rotational diffusion tensor as a restraint. The results show that the proposed docking method is robust against experimental errors in the relaxation data or structural rearrangements upon complex formation and is computationally more efficient than current methods. The developed approximations are accurate enough to be used in structure refinement protocols.

Structural analysis of a viral ovarian tumor domain protease from the Crimean-Congo hemorrhagic fever virus in complex with covalently bonded ubiquitin

Crimean-Congo hemorrhagic fever (CCHF) virus is a tick-borne, negative-sense, single-stranded RNA [ssRNA(-)] nairovirus that produces fever, prostration, and severe hemorrhages in humans. With fatality rates for CCHF ranging up to 70% based on several factors, CCHF is considered a dangerous emerging disease. Originally identified in the former Soviet Union and the Congo, CCHF has rapidly spread across large sections of Europe, Asia, and Africa. Recent reports have identified a viral homologue of the ovarian tumor protease superfamily (vOTU) within its L protein. This protease has subsequently been implicated in downregulation of the type I interferon immune response through cleavage of posttranslational modifying proteins ubiquitin (Ub) and the Ub-like interferon-simulated gene 15 (ISG15). Additionally, homologues of vOTU have been suggested to perform similar roles in the positive-sense, single-stranded RNA [ssRNA(+)] arteriviruses. By utilizing X-ray crystallographic techniques, the structure of vOTU covalently bound to ubiquitin propylamine, a suicide substrate of the enzyme, was elucidated to 1.7 Å, revealing unique structural elements that define this new subclass of the OTU superfamily. In addition, kinetic studies were carried out with aminomethylcoumarin (AMC) conjugates of monomeric Ub, ISG15, and NEDD8 (neural precursor cell expressed, developmentally downregulated 8) substrates in order to provide quantitative insights into vOTU's preference for Ub and Ub-like substrates.

Nuclear magnetic resonance analysis of protein-DNA interactions

Recent methodological and instrumental advances in solution-state nuclear magnetic resonance have opened up the way to investigating challenging problems in structural biology such as large macromolecular complexes. This review focuses on the experimental strategies currently employed to solve structures of protein-DNA complexes and to analyse their dynamics. It highlights how these approaches can help in understanding detailed molecular mechanisms of target recognition.

NMR reveals a different mode of binding of the Stam2 VHS domain to ubiquitin and diubiquitin

The VHS domain of the Stam2 protein is a ubiquitin binding domain involved in the recognition of ubiquitinated proteins committed to lysosomal degradation. Among all VHS domains, the VHS domain of Stam proteins is the strongest binder to monoubiqiuitin and exhibits preferences for K63-linked chains. In the present paper, we report the solution NMR structure of the Stam2-VHS domain in complex with monoubiquitin by means of chemical shift perturbations, spin relaxation, and paramagnetic relaxation enhancements. We also characterize the interaction of Stam2-VHS with K48- and K63-linked diubiquitin chains and report the first evidence that VHS binds differently to these two chains. Our data reveal that VHS enters the hydrophobic pocket of K48-linked diubiquitin and binds the two ubiquitin subunits with different affinities. In contrast, VHS interacts with K63-linked diubiquitin in a mode similar to its interaction with monoubiquitin. We also suggest possible structural models for both K48- and K63-linked diubiquitin in interaction with VHS. Our results, which demonstrate a different mode of binding of VHS for K48- and K63-linked diubiquitin, may explain the preference of VHS for K63- over K48-linked diubiquitin chains and monoubiquitin.

A new crystal form of Lys48-linked diubiquitin

A new crystal form of Lys48-linked diubiquitin was obtained and its structure was determined by X-ray crystallography to 1.6 Å resolution.

Lys48-linked polyubiquitin chains are recognized by the proteasome as a tag for the degradation of the attached substrates. Here, a new crystal form of Lys48-linked diubiquitin (Ub₂) was obtained and the crystal structure was refined to 1.6 Å resolution. The structure reveals an ordered isopeptide bond in a trans configuration. All three molecules in the asymmetric unit were in the same closed conformation, in which the hydrophobic patches of both the distal and the proximal moieties interact with each other. Despite the different crystallization conditions and different crystal packing, the new crystal structure of Ub₂ is similar to the previously published structure of diubiquitin, but differences are observed in the conformation of the flexible isopeptide linkage.

Structural assembly of molecular complexes based on residual dipolar couplings

We present and evaluate a rigid-body molecular docking method, called PATIDOCK, that relies solely on the three-dimensional structure of the individual components and the experimentally derived residual dipolar couplings (RDCs) for the complex. We show that, given an accurate ab initio predictor of the alignment tensor from a protein structure, it is possible to accurately assemble a protein-protein complex by utilizing the RDCs' sensitivity to molecular shape to guide the docking. The proposed docking method is robust against experimental errors in the RDCs and computationally efficient. We analyze the accuracy and efficiency of this method using experimental or synthetic RDC data for several proteins, as well as synthetic data for a large variety of protein-protein complexes. We also test our method on two protein systems for which the structure of the complex and steric-alignment data are available (Lys48-linked diubiquitin and a complex of ubiquitin and a ubiquitin-associated domain) and analyze the effect of flexible unstructured tails on the outcome of docking. The results demonstrate that it is fundamentally possible to assemble a protein-protein complex solely on the basis of experimental RDC data and the prediction of the alignment tensor from 3D structures. Thus, despite the purely angular nature of RDCs, they can be converted into intermolecular distance/translational constraints. Additionally, we show a method for combining RDCs with other experimental data, such as ambiguous constraints from interface mapping, to further improve structure characterization of protein complexes.

Building macromolecular assemblies by information-driven docking: introducing the HADDOCK multibody docking server

Over the last years, large scale proteomics studies have generated a wealth of information of biomolecular complexes. Adding the structural dimension to the resulting interactomes represents a major challenge that classical structural experimental methods alone will have difficulties to confront. To meet this challenge, complementary modeling techniques such as docking are thus needed. Among the current docking methods, HADDOCK (High Ambiguity-Driven DOCKing) distinguishes itself from others by the use of experimental and/or bioinformatics data to drive the modeling process and has shown a strong performance in the critical assessment of prediction of interactions (CAPRI), a blind experiment for the prediction of interactions. Although most docking programs are limited to binary complexes, HADDOCK can deal with multiple molecules (up to six), a capability that will be required to build large macromolecular assemblies. We present here a novel web interface of HADDOCK that allows the user to dock up to six biomolecules simultaneously. This interface allows the inclusion of a large variety of both experimental and/or bioinformatics data and supports several types of cyclic and dihedral symmetries in the docking of multibody assemblies. The server was tested on a benchmark of six cases, containing five symmetric homo-oligomeric protein complexes and one symmetric protein-DNA complex. Our results reveal that, in the presence of either bioinformatics and/or experimental data, HADDOCK shows an excellent performance: in all cases, HADDOCK was able to generate good to high quality solutions and ranked them at the top, demonstrating its ability to model symmetric multicomponent assemblies. Docking methods can thus play an important role in adding the structural dimension to interactomes. However, although the current docking methodologies were successful for a vast range of cases, considering the variety and complexity of macromolecular assemblies, inclusion of some kind of experimental information (e.g. from mass spectrometry, nuclear magnetic resonance, cryoelectron microscopy, etc.) will remain highly desirable to obtain reliable results.

Exploring the linkage dependence of polyubiquitin conformations using molecular modeling

Posttranslational modification of proteins by covalent attachment of a small protein ubiquitin (Ub) or a polymeric chain of Ub molecules (called polyubiquitin) is involved in controlling a vast variety of processes in eukaryotic cells. The question of how different polyubiquitin signals are recognized is central to understanding the specificity of various types of polyubiquitination. In polyubiquitin, monomers are linked to each other via an isopeptide bond between the C-terminal glycine of one Ub and a lysine of the other. The functional outcome of polyubiquitination depends on the particular lysine involved in chain formation and appears to rely on linkage-dependent conformation of polyubiquitin. Thus, K48-linked chains, a universal signal for proteasomal degradation, under physiological conditions adopt a closed conformation where functionally important residues L8, I44, and V70 are sequestered at the interface between two adjacent Ub monomers. By contrast, K63-linked chains, which act as a nonproteolytic regulatory signal, adopt an extended conformation that lacks hydrophobic interubiquitin contact. Little is known about the functional roles of the so-called "noncanonical" chains (linked via K6, K11, K27, K29, or K33, or linked head-to-tail), and no structural information on these chains is available, except for information on the crystal structure of the head-to-tail-linked diubiquitin (Ub(2)). In this study, we use molecular modeling to examine whether any of the noncanonical chains can adopt a closed conformation similar to that in K48-linked polyubiquitin. Our results show that the eight possible Ub(2) chains can be divided into two groups: chains linked via K6, K11, K27, or K48 are predicted to form a closed conformation, whereas chains linked via K29, K33, or K63, or linked head-to-tail are unable to form such a contact due to steric occlusion. These predictions are validated by the known structures of K48-, K63-, and head-to-tail-linked chains. Our study also predicts structural models for Ub(2) chains linked via K6, K11, or K27. The implications of these findings for linkage-selective recognition of noncanonical polyubiquitin signals by various receptors are discussed.

Insights into the DNA cleavage mechanism of human LINE-1 retrotransposon endonuclease

The human LINE-1 endonuclease (L1-EN) contributes in defining the genomic integration sites of the abundant human L1 and Alu retrotransposons. LINEs have been considered as possible vehicles for gene delivery and understanding the mechanism of L1-EN could help engineering them as genetic tools. We tested the in vitro activity of point mutants in three L1-EN residues--Asp145, Arg155, Ile204--that are key for DNA cleavage, and determined their crystal structures. The L1-EN structure remains overall unaffected by the mutations, which change the enzyme activity but leave DNA cleavage sequence specificity mostly unaffected. To better understand the mechanism of L1-EN, we performed molecular dynamics simulations using as model the structures of wild type EN-L1, of two betaB6-betaB5 loop exchange mutants we have described previously to be important for DNA recognition, of the R155A mutant from this study, and of the homologous TRAS1 endonuclease: all confirm a rigid scaffold. The simulations crucially indicate that the betaB6-betaB5 loop shows an anticorrelated motion with the surface loops betaA6-betaA5 and betaB3-alphaB1. The latter loop harbors N118, a residue that alters DNA cleavage specificity in homologous endonucleases, and implies that the plasticity and correlated motion of these loops has a functional importance in DNA recognition and binding. To further explore how these loops are possibly involved in DNA binding, we docked computationally two DNA substrates to our structure, one involving a flipped-out nucleotide downstream the scissile phosphodiester; and one not. The models for both scenarios are feasible and agree with the hypotheses derived from the dynamic simulations. The reduced cleavage activity we have observed for the I204Y mutant above however, favors the flipped out nucleotide model.

NMR structural analysis of DNA recognition by a novel Myb1 DNA-binding domain in the protozoan parasite Trichomonas vaginalis

The transcription regulator, tvMyb1, is the first Myb family protein identified in Trichomonas vaginalis. Using an electrophoretic mobility shift assay, we defined the amino-acid sequence from Lys(35) to Ser(141) (tvMyb1(35-141)) as the minimal DNA-binding domain, encompassing two Myb-like DNA-binding motifs (designated as R2 and R3 motifs) and an extension of 10 residues at the C-terminus. NMR solution structures of tvMyb1(35-141) show that both the R2 and R3 motifs adopt helix-turn-helix conformations while helix 6 in the R3 motif is longer than its counterpart in vertebrate Myb proteins. The extension of helix 6 was then shown to play an important role in protein stability as well as in DNA-binding activity. The structural basis for the tvMyb1(35-141)/DNA interaction was investigated using chemical shift perturbations, residual dipolar couplings, DNA specificity data and data-driven macromolecular docking by HADDOCK. Our data indicate that the orientation between R2 and R3 motifs dramatically changes upon binding to DNA so as to recognize the DNA major groove through a number of key contacts involving residues in helices 3 and 6. The tvMyb1(35-141)/DNA complex model furthers our understanding of DNA recognition by Myb proteins and this approach could be applied in determining the complex structures involving proteins with multiple domains.

Solution NMR structure of putidaredoxin-cytochrome P450cam complex via a combined residual dipolar coupling-spin labeling approach suggests a role for Trp106 of putidaredoxin in complex formation

The 58-kDa complex formed between the [2Fe-2S] ferredoxin, putidaredoxin (Pdx), and cytochrome P450cam (CYP101) from the bacterium Pseudomonas putida has been investigated by high-resolution solution NMR spectroscopy. Pdx serves as both the physiological reductant and effector for CYP101 in the enzymatic reaction involving conversion of substrate camphor to 5-exo-hydroxycamphor. In order to obtain an experimental structure for the oxidized Pdx-CYP101 complex, a combined approach using orientational data on the two proteins derived from residual dipolar couplings and distance restraints from site-specific spin labeling of Pdx has been applied. Spectral changes for residues in and near the paramagnetic metal cluster region of Pdx in complex with CYP101 have also been mapped for the first time using (15)N and (13)C NMR spectroscopy, leading to direct identification of the residues strongly affected by CYP101 binding. The new NMR structure of the Pdx-CYP101 complex agrees well with results from previous mutagenesis and biophysical studies involving residues at the binding interface such as formation of a salt bridge between Asp38 of Pdx and Arg112 of CYP101, while at the same time identifying key features different from those of earlier modeling studies. Analysis of the binding interface of the complex reveals that the side chain of Trp106, the C-terminal residue of Pdx and critical for binding to CYP101, is located across from the heme-binding loop of CYP101 and forms non-polar contacts with several residues in the vicinity of the heme group on CYP101, pointing to a potentially important role in complex formation.