Post-translational modification with one of the isoforms of the small ubiquitin-like modifier (SUMO) affects thousands of proteins in the human proteome. The binding of SUMO to SUMO interacting motifs (SIMs) can translate the SUMOylation event into functional consequences. The E3 ubiquitin ligase RNF4 contains multiple SIMs and connects SUMOylation to the ubiquitin pathway. SIM2 and SIM3 of RNF4 were shown to be the most important motifs to recognize SUMO chains. However, the study of SIM-SUMO complexes is complicated by their typically low affinity and variable binding of the SIMs in parallel and antiparallel orientations. We investigated properties of complexes formed by SUMO3 with peptides containing either SIM2 or SIM3 using molecular dynamics simulations. The affinities of the complexes were determined using a state-of-the-art free energy protocol and were found to be in good agreement with experimental data, thus corroborating our method. Long unrestrained simulations allowed a new interpretation of experimental results regarding the structure of the SIM-SUMO interface. We show that both SIM2 and SIM3 bind SUMO3 in parallel and antiparallel orientations and identified main interaction sites for acidic residues flanking the SIM. We noticed unusual SIM-SUMO interfaces in a previously reported NMR structure (PDB: 2mp2) of a complex formed by a SUMO3 dimer with the bivalent SIM2-SIM3 peptide. Computational determination of the individual SIM-SUMO affinities based on these structural arrangements yielded significantly higher dissociation constants. To our knowledge, our approach adds new opportunities to characterize individual SIM-SUMO complexes and suggests that further studies will be necessary to understand these interactions when occurring in multivalent form.

Interaction of the small ubiquitin-related modifier (SUMO) and peptides containing a SUMO-interacting motif (SIM) has attracted a lot of interest in recent years, yet their structural properties and relationships between the composition of the peptide and binding free energy are not completely understood. We perform molecular dynamics simulations of the complex formed by SUMO and a peptide containing the tight-binding SIM of the protein inhibitor of activated STAT. The calculated standard binding free energy of −5.06 kcal/mol is in reasonable agreement with the experimental value of −6.54 kcal/mol. Experimental results for complexes formed by SUMO and SIM dimers indicate the existence of a parallel and an antiparallel binding mode for similar SIM peptides. We find that the parallel binding mode is highly favored in the present case. Furthermore, the simulations show that residues neighboring the SIM core motif contribute strongly to the binding energy. Structurally, the complex shows differences from the picture in which the SIM core motif lies deep in the SUMO binding groove. This also supports the idea that neighboring residues play an important role in binding.

SUMO targeted ubiqutin ligases (STUbLs) like RNF4 or Arkadia/RNF111 recognize SUMO chains through multiple SUMO interacting motifs (SIMs). Typically, these are contained in disordered regions of these enzymes and also the individual SUMO domains of SUMO chains move relatively freely. It is assumed that binding the SIM region significantly restricts the conformational freedom of SUMO chains. Here, we present the results of extensive molecular dynamics simulations on the complex formed by the SIM2−SIM3 region of RNF4 and diSUMO3. Though our simulations highlight the importance of typical SIM−SUMO interfaces also in the multivalent situation, we observe that frequently other regions of the peptide than the canonical SIMs establish this interface. This variability regarding the individual interfaces leads to a conformationally highly flexible complex. Comparison with previous experimental measurements clearly supports our findings and indicates that our observations can be extended to other multivalent SIM−SUMO complexes.

SUMO targeted ubiqutin ligases (STUbLs) like RNF4 or Arkadia/RNF111 recognize SUMO chains through multiple SUMO interacting motifs (SIMs). Typically, these are contained in disordered regions of these enzymes and also the individual SUMO domains of SUMO chains move relatively freely. It is assumed that binding the SIM region significantly restricts the conformational freedom of SUMO chains. Here, we present the results of extensive molecular dynamics simulations on the complex formed by the SIM2-SIM3 region of RNF4 and diSUMO3. Though our simulations highlight the importance of typical SIM-SUMO interfaces also in the multivalent situation, we observe that frequently other regions of the peptide than the canonical SIMs establish this interface. This variability regarding the individual interfaces leads to a conformationally highly flexible complex. Even though this is in contrast to previous models of the RNF4 SUMO chain interaction, we demonstrate that our simulations are clearly consistent with previous experimental results.